Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,704,138 A→T 5.4% F208L (TTT→TTA dppC ← dipeptide/heme ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,704,1380AT 5.4% 681.5 / 20.0 185F208L (TTT→TTAdppCdipeptide/heme ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (92/83);  new base T (5/5);  total (97/88)
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.78e-01

Only 100 of 203 total aligned reads displayed.

CGGAGAGCATGGTGCCCCACTCAGGCGTTGGCGGCTGTGCCCCCATGCCGAGGAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACG  >  NC_000913/3703993‑3704286
                                                                                                                                                 |                                                                                                                                                    
cGGAGAGCATGGTGCCCCACTCAGGCGTTGGCGGCTGTGCCCCCATGCCGAGGAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAAc                                                                                                                                                   <  4:1307800/149‑1 (MQ=255)
    gagCATGGTGCCCCACTCAGGCGTTGGCGGCTGTGCCCCCATGCCGAGGAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCt                                                                                                                                               >  2:49165/1‑149 (MQ=255)
      gCATGGTGCCCCACTCAGGCGTTGGCGGCTGTGCCCCCATGCCGAGGAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGa                                                                                                                                             <  3:593475/149‑1 (MQ=255)
       cATGGTGCCCCACTCAGGCGTTGGCGGCTGTGCCCCCATGCCGAGGAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGa                                                                                                                                             <  4:711607/148‑1 (MQ=255)
            tGCCCCACTCAGGCGTTGGCGGCTGTGCCCCCATGCCGAGGAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATc                                                                                                                                       <  1:875646/149‑1 (MQ=255)
               cccACTCAGGCGTTGGCGGCTGTGCCCCCATGCCGAGGAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACAt                                                                                                                                                 >  3:849386/1‑136 (MQ=255)
                ccACTCAGGCGTTGGCGGCTGTGCCCCCATGCCGAGGAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGccc                                                                                                                                   >  4:845863/1‑149 (MQ=255)
                  aCTCAGGCGTTGGCGGCTGTGCCCCCATGCCGAGGAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCg                                                                                                                                 <  1:590179/149‑1 (MQ=255)
                  aCTCAGGCGTTGGCGGCTGTGCCCCCATGCCGAGGAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCg                                                                                                                                 <  1:1347711/149‑1 (MQ=255)
                    tCAGGCGTTGGCGGCTGTGCCCCCATGCCGAGGAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACa                                                                                                                                                  >  2:353866/1‑130 (MQ=255)
                    tCAGGCGTTGGCGGCTGTGCCCCCATGCCGAGGAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGc                                                                                                                               >  4:176694/1‑149 (MQ=255)
                       ggCGTTGGCGGCTGTGCCCCCATGCCGAGGAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACAt                                                                                                                                                 >  3:687951/1‑128 (MQ=255)
                        gCGTTGGCGGCTGTGCCCCCATGCCGAGGAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCAccc                                                                                                                           <  4:902756/149‑1 (MQ=255)
                          gTTGGCGGCTGTGCCCCCATGCCGAGGAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGc                                                                                                                         <  3:1510451/149‑1 (MQ=255)
                            tGGCGGCTGTGCCCCCATGCCGAGGAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCa                                                                                                                       <  4:265603/149‑1 (MQ=255)
                                    gtgCCCCCATGCCGAGGAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGAc                                                                                                               >  3:918831/1‑149 (MQ=255)
                                        ccccATGCCGAGGAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCa                                                                                                                       >  1:1711453/1‑137 (MQ=255)
                                        ccccATGCCGAGGAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCa                                                                                                                       <  2:1711453/137‑1 (MQ=255)
                                        ccccATGCCGAGGAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCgg                                                                                                           >  2:1462982/1‑149 (MQ=255)
                                         cccATGCCGAGGAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCAt                                                                                                                                        <  3:1088858/119‑1 (MQ=255)
                                         cccATGCCGAGGAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGt                                                                                                          >  3:739803/1‑149 (MQ=255)
                                                 gAGGAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAAt                                                                                                  <  3:1204499/149‑1 (MQ=255)
                                                 gAGGAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAAt                                                                                                  >  2:5533/1‑149 (MQ=255)
                                                 gAGGAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAAt                                                                                                  >  2:1811755/1‑149 (MQ=255)
                                                  aGGAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATc                                                                                                 >  2:455917/1‑149 (MQ=255)
                                                   ggAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCgaga                                                                                                                <  1:1681816/133‑1 (MQ=255)
                                                   ggAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCgaga                                                                                                                >  2:1681816/1‑133 (MQ=255)
                                                   ggAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCga                                                                                                                  >  1:1080554/1‑131 (MQ=255)
                                                    gAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGAc                                                                                                       >  4:525342/1‑141 (MQ=255)
                                                    gAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATcgc                                                                                               >  2:1270746/1‑149 (MQ=255)
                                                     aaaCCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATcgcg                                                                                              <  1:1270746/149‑1 (MQ=255)
                                                      aaCCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATcgcgc                                                                                             >  2:257233/1‑148 (MQ=255)
                                                      aaCCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCgg                                                                                             >  4:1668380/1‑149 (MQ=255)
                                                        ccAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTTATTAAc                                                                                                                                                   <  1:833258/93‑1 (MQ=255)
                                                         cAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGAc                                                                                                       <  4:1506746/136‑1 (MQ=255)
                                                         cAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGt                                                                                                     >  2:1406620/1‑138 (MQ=255)
                                                           agagCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCg                                                                                                                                      >  1:477405/1‑103 (MQ=255)
                                                            gagCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACt                                                                                       <  4:592795/149‑1 (MQ=255)
                                                                 gcCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTTATTAACATCttc                                                                                                                                             <  4:1095165/90‑3 (MQ=255)
                                                                  ccATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTTATTAACATCtt                                                                                                                                              <  1:405449/88‑2 (MQ=255)
                                                                   cATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGa                                                                                                        <  2:611734/125‑1 (MQ=255)
                                                                    atatCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCgg                                                                                             >  4:674381/1‑135 (MQ=255)
                                                                    atatCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAg                                                                               >  3:1263097/1‑149 (MQ=255)
                                                                       tCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCa                                                                             <  2:631519/148‑1 (MQ=255)
                                                                               ggCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGgcgc                                                                    >  3:1722690/1‑149 (MQ=255)
                                                                                      gagaAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCgg                                                                  <  4:409205/144‑1 (MQ=255)
                                                                                        gaAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTcc                                                                                    >  3:1420009/1‑124 (MQ=255)
                                                                                        gaAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAgc                                                           >  3:400274/1‑149 (MQ=255)
                                                                                          aaCCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATaaa                                                                                                                                                    >  1:852601/1‑58 (MQ=255)
                                                                                           aCCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGcaca                                                        >  4:1221795/1‑149 (MQ=255)
                                                                                           aCCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGca                                                        <  1:1811755/149‑1 (MQ=255)
                                                                                                  gACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTg                                                 >  1:1523217/1‑149 (MQ=255)
                                                                                                   aCGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGc                                                <  1:1440177/149‑1 (MQ=255)
                                                                                                      ccTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCggc                                             <  3:1142090/149‑1 (MQ=255)
                                                                                                        tGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGAc                                                                                                       >  2:856864/1‑89 (MQ=255)
                                                                                                        tGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCggg                                                                 <  1:511870/127‑1 (MQ=255)
                                                                                                        tGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCggcaa                                           >  1:1212271/1‑149 (MQ=255)
                                                                                                        tGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCggcaa                                           >  1:677517/1‑149 (MQ=255)
                                                                                                        tGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCggcaa                                           >  3:529046/1‑149 (MQ=255)
                                                                                                          aaCAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCggcaagg                                         <  3:577095/149‑1 (MQ=255)
                                                                                                             aaTCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAg                                                   <  4:355045/136‑1 (MQ=255)
                                                                                                             aaTCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAg                                                   >  3:355045/1‑136 (MQ=255)
                                                                                                               tCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACg                                    >  3:299133/1‑149 (MQ=255)
                                                                                                                  gcggcgCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATcgcg                                                                                              >  2:1066008/1‑88 (MQ=255)
                                                                                                                  gcggcgCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCggcaag                                          <  2:841143/140‑1 (MQ=255)
                                                                                                                  gcggcgCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAg                                 <  2:694711/149‑1 (MQ=255)
                                                                                                                      cgcAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGc                                                                              >  2:1037015/1‑100 (MQ=255)
                                                                                                                      cgcAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGc                                                                              <  1:1037015/100‑1 (MQ=255)
                                                                                                                         aaGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAg                                 >  4:534226/1‑142 (MQ=255)
                                                                                                                           ggCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTccac                                                                                  >  2:673521/1‑91 (MQ=255)
                                                                                                                             cAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAg                      >  4:221799/1‑149 (MQ=255)
                                                                                                                             cAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAg                      <  3:885950/149‑1 (MQ=255)
                                                                                                                               gTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAg                                 <  2:1602981/136‑1 (MQ=255)
                                                                                                                               gTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAg                                 >  1:1602981/1‑136 (MQ=255)
                                                                                                                                ttCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAgcgc                   <  2:1212271/149‑1 (MQ=255)
                                                                                                                                   gggAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCgcggc                >  1:1051587/1‑149 (MQ=255)
                                                                                                                                    gtcaGATGTTTATTAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGaa                                  <  3:621794/127‑1 (MQ=255)
                                                                                                                                    gtcaGATGTTTATTAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGaa                                  >  4:621794/4‑130 (MQ=255)
                                                                                                                                     gAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGaa                                  <  4:1043009/129‑1 (MQ=255)
                                                                                                                                     gAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGt              >  3:970769/1‑149 (MQ=255)
                                                                                                                                      aaGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGtt             <  1:613680/149‑1 (MQ=255)
                                                                                                                                        gATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTAc           <  2:1523217/149‑1 (MQ=255)
                                                                                                                                        gATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTAc           <  3:1013023/149‑1 (MQ=255)
                                                                                                                                        gATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTAc           >  3:500940/1‑149 (MQ=255)
                                                                                                                                           gTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGTGGCAAGGCAACGAAGGTCAGTTCCAGCGCGGCGTTACCaa        <  1:547314/149‑1 (MQ=255)
                                                                                                                                            ttAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCggcaag                                          <  3:246354/114‑1 (MQ=255)
                                                                                                                                            ttAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCggcaag                                          <  2:522499/114‑1 (MQ=255)
                                                                                                                                            ttAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAAt       <  4:855733/149‑1 (MQ=255)
                                                                                                                                            ttAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAAt       <  3:1807047/149‑1 (MQ=255)
                                                                                                                                              aataaACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCg     >  1:1475166/1‑149 (MQ=255)
                                                                                                                                               ataaACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCggg                                                                 <  2:1840488/88‑1 (MQ=255)
                                                                                                                                               ataaACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCggg                                                                 >  1:1840488/1‑88 (MQ=255)
                                                                                                                                               ataaACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGTGCGGCGTTAc           >  4:1597036/1‑142 (MQ=255)
                                                                                                                                               ataaACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTAc           <  3:1597036/142‑1 (MQ=255)
                                                                                                                                                taaACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGt                           >  3:1870227/1‑125 (MQ=255)
                                                                                                                                                taaACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGt                           <  4:1870227/125‑1 (MQ=255)
                                                                                                                                                 aaaCATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGaa                                  <  1:649859/117‑1 (MQ=255)
                                                                                                                                                 aaaCATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCaa        <  1:1748248/143‑1 (MQ=255)
                                                                                                                                                 aaaCATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCaa        >  2:1748248/1‑143 (MQ=255)
                                                                                                                                                 aaaCATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACg  >  1:301692/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CGGAGAGCATGGTGCCCCACTCAGGCGTTGGCGGCTGTGCCCCCATGCCGAGGAAACCAAGAGCAGCCATATCGAGAATGGCGTTAGAGAAACCGAGCGACGCCTGAACAATCAGCGGCGCAAGGCAGTTCGGGAAGATGTTAATAAACATCTGACGCATCGCCCCGGCACCCGCCACGCGAGACGCGGTGACGTAATCGCGGTTAACTTCCACCAGCACGGCGGCGCGGGTTAAGCGCACATAGTGCGGCAAGGCAACGAAGGTCAGTGCCAGCGCGGCGTTACCAATCGACG  >  NC_000913/3703993‑3704286

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: