Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,395,812 G→C 6.4% L29V (CTG→GTG)  nuoJ ← NADH:ubiquinone oxidoreductase, membrane subunit J

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,395,8120GC 6.4% 160.1 / 2.5 58L29V (CTG→GTG) nuoJNADH:ubiquinone oxidoreductase, membrane subunit J
Reads supporting (aligned to +/- strand):  ref base G (23/31);  new base C (2/2);  total (25/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CCCAGGTTGAGCATCATCACCACGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTCGTTCTCTGCTTC  >  NC_000913/2395667‑2395958
                                                                                                                                                  |                                                                                                                                                  
cccAGGTTGAGCATCATCACCACGAACACGAACAGCACCATTATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAATAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGc                                                                                                                                                 <  1:384607/149‑1 (MQ=255)
        gAGCATCATCACCACGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCg                                                                                                                                         >  2:732634/1‑149 (MQ=255)
        gAGCATCATCACCACGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCg                                                                                                                                         <  2:664863/149‑1 (MQ=255)
         aGCATCATCACCACGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGCACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCgg                                                                                                                                        >  1:133130/1‑149 (MQ=255)
           catcatCACCACGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCCGGTACACGCACTGC‑TGTACCTGAt                                                                                                                                      <  2:596409/149‑1 (MQ=255)
               atcaCCACGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGt                                                                                                                                  <  1:199573/149‑1 (MQ=255)
                tcaCCACGAACACGAACAGCACCATAATGGCCCCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACACCGGAAATCGCCAGCAGCGAAATTATCAGGTACC‑GCTGTGCGTGTACCGGATTGGTa                                                                                                                                 >  2:1076149/1‑149 (MQ=255)
                tcaCCACGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTa                                                                                                                                 >  1:889663/1‑149 (MQ=255)
                 caccacGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCCGGTACACGCACTGC‑TGTACCTGAttatttc                                                                                                                                <  1:260756/149‑6 (MQ=38)
                   ccacGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATgg                                                                                                                              <  2:604642/149‑1 (MQ=255)
                     acgaccacgaacAGCACCATAATCGCACCCGCGTAGACTATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGt                                                                                                                            <  1:149302/149‑1 (MQ=255)
                       gaacacgaccaGCACCATAATCGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGa                                                                                                                          <  2:698773/149‑1 (MQ=255)
                        aacacgtacaGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGAt                                                                                                                         <  1:939796/149‑1 (MQ=255)
                        aacacgaacaGCACCAAAATGGCACCCGCGTAGACGAAAATTTCCATCGCACCGGCGAAGTAAGCGCCCAGTGAGAACAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGAt                                                                                                                         <  1:647612/149‑1 (MQ=255)
                             gaacaGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACTCCCCGGAAATCGCCAGCAGCGAAATATTCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGCTCACTc                                                                                                                    >  1:82371/1‑149 (MQ=255)
                                 aGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCaa                                                                                                                <  2:336280/149‑1 (MQ=255)
                                        tAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGAGATCACTCGCAAGGTCGCa                                                                                                         <  1:232647/149‑1 (MQ=255)
                                        tAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCCGCGAAGTAATCGCCCAGTGAGAAGACCACCCAGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTAGCAAGGTCGCa                                                                                                         >  2:1067892/1‑149 (MQ=255)
                                         aaTGCCACCCGCGTAGACGATAAATTCCAGCGCACCGGCGAAGTAATCGCCCAGTGCGAAGAACACCCCGGAACTCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTCTGGGTGATCACTCGCAAGGTCGCaa                                                                                                        <  2:1167508/149‑1 (MQ=255)
                                                       aGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAgg                                                                                          <  2:565221/149‑1 (MQ=255)
                                                           gATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTCCCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCcac                                                                                      >  1:1021631/1‑149 (MQ=255)
                                                             tAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAg                                                                                    >  2:93202/1‑149 (MQ=255)
                                                             tAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAg                                                                                    >  1:719618/1‑149 (MQ=255)
                                                               aTTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGat                                                                                  <  2:442381/149‑1 (MQ=255)
                                                                   ccAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATaa                                                                              <  1:1086349/149‑1 (MQ=255)
                                                                   ccAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGGGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATaa                                                                              <  1:1069768/149‑1 (MQ=255)
                                                                    cAGCCCACCGGCGAACAAAGCGCCCAGTGAGAAGAACACCCCGGCAATCGACAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGAACCGGATAGGAATCGGTGATCACTAGAAAGGTCGCAAGTATAGATATCAGGCCACATATATaaa                                                                             >  2:785184/1‑149 (MQ=255)
                                                                      gcgcACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATaaaaa                                                                           <  2:374408/149‑1 (MQ=255)
                                                                        gcACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACT‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGc                                                                         <  1:99718/149‑1 (MQ=255)
                                                                          aCCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCTCCAGCAGCGAAATAATCAGGTACA‑GCACTGCGAGTACCGGGTTGGTATGGGTGATCACTCGCAAGGACGCAACTATGGAAATCAGGCCACAGATATAAAAAGcca                                                                       >  1:1043886/1‑147 (MQ=255)
                                                                                    tAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGc                                                             >  2:1169642/1‑149 (MQ=255)
                                                                                    tAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGc                                                             >  2:193428/1‑149 (MQ=255)
                                                                                      aGCGCCCAGTGAGAAGAACCCCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACG‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGccc                                                           <  2:279042/149‑1 (MQ=255)
                                                                                       gcgcCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGcccc                                                          >  2:532744/1‑149 (MQ=255)
                                                                                          cccAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCctct                                                       <  1:382345/149‑1 (MQ=255)
                                                                                             aGTGAGAAGAACACACCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCTCAGATATAAAAAGCGAACTCCATTGCACCTCTACt                                                    >  2:1128302/1‑149 (MQ=255)
                                                                                               tGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTa                                                  >  1:841207/1‑149 (MQ=255)
                                                                                               tGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTa                                                  <  2:224613/149‑1 (MQ=255)
                                                                                                gagaAGAACACCCCGGAAATCGCCAGCAGCGATATAATCTGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCTCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTAc                                                 <  1:601450/149‑1 (MQ=255)
                                                                                                gagaAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCATGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTAc                                                 <  1:606374/149‑1 (MQ=255)
                                                                                                gagaAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTAc                                                 <  1:531436/149‑1 (MQ=255)
                                                                                                      aaCACCCAGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAAc                                           <  1:286980/149‑1 (MQ=255)
                                                                                                           cccGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACTGTAACAGGct                                      >  2:839701/1‑149 (MQ=255)
                                                                                                           cccGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGGGTACCGGATTGGTATGGGTGATCTCTCGCAAGGTCGGATGTATGGCTTTCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGct                                      <  2:1059128/149‑1 (MQ=255)
                                                                                                            ccGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGctc                                     >  2:938955/1‑149 (MQ=255)
                                                                                                                acaTCGCCAGCAGCGAAATAATCAGGTACC‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGGCGCAAGTATGGCTATCAGGCCACAGCTATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGa                                 <  2:388497/147‑1 (MQ=255)
                                                                                                                  aTCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACg                               <  2:1097820/149‑1 (MQ=255)
                                                                                                                   tCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGt                              >  1:116012/1‑149 (MQ=255)
                                                                                                                     gCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTc                                                                   >  2:738419/1‑110 (MQ=255)
                                                                                                                     gCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTc                                                                   <  1:738419/110‑1 (MQ=255)
                                                                                                                      ccAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGa                           <  1:93202/149‑1 (MQ=255)
                                                                                                                      ccAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGATCTTGACGTCGa                           >  2:299028/1‑149 (MQ=255)
                                                                                                                        agcagcGAAATAATAAGGTACA‑GCAGTGCGTGTACCGGATTGGAATGGGTGAACACTCGCAAGGTCCCAAGTATGGCTATCAGGCCACATATATAAAAAGCGAAATCCAATGCCAATCACCTTACGGAAAAAGGCTCTTGACGTCGATa                         >  2:306523/1‑149 (MQ=255)
                                                                                                                           agcGAAATATTGAGGTACA‑GCAGTGCGTGTACCGGATTGGTTTGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCAGAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAggc                      <  1:25855/149‑1 (MQ=255)
                                                                                                                            gcGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAggct                     <  2:14415/149‑1 (MQ=255)
                                                                                                                                  caatCCGGTACACGCACTGC‑TGTACCTGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAggcttggctt               >  1:1080501/2‑149 (MQ=255)
                                                                                                                                   aatCCGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCAATCTGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTATCCTTACGGAAACAGGCTCTTGACGTCTATAGGCTGGGCTTc              >  2:1075020/1‑149 (MQ=255)
                                                                                                                                     tCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTCGt            >  2:354524/1‑149 (MQ=255)
                                                                                                                                            ca‑gcagTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTCGTTCTCTGc     <  2:977806/149‑1 (MQ=255)
                                                                                                                                             a‑gcagTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTCGTTCTCTGCt    >  1:1014404/1‑149 (MQ=255)
                                                                                                                                                cagTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTCGTTCTCTGCtta  >  2:254825/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CCCAGGTTGAGCATCATCACCACGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTCGTTCTCTGCTTC  >  NC_000913/2395667‑2395958

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑