Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,460,477 A→C 8.1% intergenic (‑113/+45) nrdG ← / ← nrdD anaerobic ribonucleoside‑triphosphate reductase activating protein/anaerobic ribonucleoside‑triphosphate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,460,4770AC 8.1% 497.8 / 30.6 171intergenic (‑113/+45)nrdG/nrdDanaerobic ribonucleoside‑triphosphate reductase activating protein/anaerobic ribonucleoside‑triphosphate reductase
Reads supporting (aligned to +/- strand):  ref base A (87/70);  new base C (2/12);  total (89/82)
Fisher's exact test for biased strand distribution p-value = 4.03e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 173 total aligned reads displayed.

GACGATGTCGACAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGCACAC  >  NC_000913/4460329‑4460626
                                                                                                                                                    |                                                                                                                                                     
 aCGATGTCGAGATGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCGGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAATTGTAGGCCGGATAAgg                                                                                                                                                     <  2:369951/150‑1 (MQ=255)
        cGACAGGATAGTACTGATGATAATTCATTATTTCACCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTAc                                                                                                                                              >  1:1345727/1‑150 (MQ=255)
         gACAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTAcg                                                                                                                                             >  1:296338/1‑150 (MQ=255)
         gACAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTAcg                                                                                                                                             <  2:1809793/150‑1 (MQ=255)
          aCAGGATAGTACTGATGATAATTCATTATTTCACCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                            <  2:1345727/150‑1 (MQ=255)
             ggATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTAcgccgc                                                                                                                                         >  2:39951/1‑150 (MQ=255)
                 aGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                     <  1:1176602/150‑1 (MQ=38)
                 aGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                     <  2:61084/150‑1 (MQ=38)
                 aGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATcc                                                                                                                                     <  2:92862/150‑1 (MQ=255)
                    aCTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                  <  1:1920148/150‑1 (MQ=38)
                    aCTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATCTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                  >  2:2039348/1‑150 (MQ=255)
                     cTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                 <  2:1483474/150‑2 (MQ=38)
                     cTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                 <  1:1149733/150‑2 (MQ=38)
                     cTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                 <  1:1459655/150‑1 (MQ=255)
                     cTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                 >  1:982381/1‑150 (MQ=255)
                       gatgatAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAtt                                                                                                                               >  2:1948757/1‑150 (MQ=255)
                         tgatAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTg                                                                                                                             >  2:1915382/1‑150 (MQ=255)
                         tgatAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTg                                                                                                                             >  2:1737119/1‑150 (MQ=255)
                          gatAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGc                                                                                                                            <  1:2064445/150‑1 (MQ=255)
                           atAATTGATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCTGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGGAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCt                                                                                                                           <  2:1360451/150‑1 (MQ=255)
                           atAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCt                                                                                                                           >  1:390793/1‑150 (MQ=255)
                             aaTTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTa                                                                                                                         >  1:788560/1‑150 (MQ=255)
                              aTTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAg                                                                                                                        >  1:202793/1‑150 (MQ=255)
                               ttCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGc                                                                                                                       >  2:2455723/1‑150 (MQ=255)
                               ttCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGc                                                                                                                       >  1:1707708/1‑150 (MQ=255)
                               ttCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGc                                                                                                                       >  1:1012017/1‑150 (MQ=255)
                               ttCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGGTTACGCCGCATCCGGCATTTGCTTAGc                                                                                                                       <  1:362054/150‑1 (MQ=255)
                                 cATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTc                                                                                                                     >  2:2111857/1‑150 (MQ=255)
                                 cATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTc                                                                                                                     >  1:962104/1‑150 (MQ=255)
                                 cATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTc                                                                                                                     >  2:2004809/1‑150 (MQ=255)
                                 cATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTc                                                                                                                     >  2:1773050/1‑150 (MQ=255)
                                    tattTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCaga                                                                                                                  >  1:838548/1‑150 (MQ=255)
                                     attTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCagag                                                                                                                 <  1:67176/150‑1 (MQ=255)
                                              ggCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTa                                                                                                        >  2:2240243/1‑150 (MQ=255)
                                                cTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAAc                                                                                                      >  1:2052160/1‑150 (MQ=255)
                                                cTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAAc                                                                                                      <  2:2238580/150‑1 (MQ=255)
                                                 ttATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAAcc                                                                                                     <  1:335271/150‑1 (MQ=255)
                                                 ttATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGAATTTGCTTAGCTCAGAGTGAAGATTAAcc                                                                                                     <  2:449597/150‑1 (MQ=255)
                                                  tATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCt                                                                                                    >  2:1579893/1‑150 (MQ=255)
                                                       gTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTg                                                                                               >  1:1441340/1‑150 (MQ=255)
                                                         ccGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGcc                                                                                             <  2:1521389/150‑1 (MQ=255)
                                                          cgcgTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGccc                                                                                            >  2:123760/1‑150 (MQ=255)
                                                               gtAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTcc                                                                                       >  2:2004171/1‑150 (MQ=255)
                                                                  ggCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCa                                                                                    >  1:1326800/1‑150 (MQ=255)
                                                                   gCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCaa                                                                                   >  2:88463/1‑150 (MQ=255)
                                                                   gCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGAAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCaa                                                                                   >  2:1079242/1‑150 (MQ=255)
                                                                      gCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATg                                                                                >  1:847886/1‑150 (MQ=255)
                                                                       cATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGt                                                                               >  2:114180/1‑150 (MQ=255)
                                                                       cATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGt                                                                               >  2:2308115/1‑150 (MQ=255)
                                                                          ccGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTa                                                                            <  1:1096599/150‑1 (MQ=255)
                                                                           cGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTaa                                                                           >  1:1671006/1‑150 (MQ=255)
                                                                            ggTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAAc                                                                          >  2:2199349/1‑150 (MQ=255)
                                                                              tAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAAcgc                                                                        <  1:149536/150‑1 (MQ=255)
                                                                                 tCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGAACTGGAGCAATTGACAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACgcggc                                                                     >  2:852358/1‑150 (MQ=255)
                                                                                       ccTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTATCCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTaa                                                               <  2:962104/150‑1 (MQ=255)
                                                                                          gATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTa                                                                                                        >  1:2209847/1‑106 (MQ=255)
                                                                                          gATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTa                                                                                                        >  1:2218431/1‑106 (MQ=255)
                                                                                          gATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTa                                                                                                        <  2:2218431/106‑1 (MQ=255)
                                                                                          gATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAActt                                                            >  1:530979/1‑150 (MQ=255)
                                                                                            tGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAActtct                                                          >  2:173443/1‑150 (MQ=255)
                                                                                              cGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACttcttc                                                        >  1:1943512/1‑150 (MQ=255)
                                                                                                       cgcgTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTAcc                                               >  2:1190158/1‑150 (MQ=255)
                                                                                                         cgTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAg                                             >  1:592584/1‑150 (MQ=255)
                                                                                                         cgTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTACTGCTTACCAg                                             <  2:1997127/150‑1 (MQ=255)
                                                                                                           tCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGTATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGTAGTTTAAACTATCTGCCAATTCCCCAATTTTTTAAAGAGGAGCTTATCTTCTTACTGGTTACCAgag                                           <  2:982381/150‑3 (MQ=37)
                                                                                                               aTCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTaa                                       >  1:2381112/1‑150 (MQ=255)
                                                                                                                tCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTaaa                                      <  1:994406/150‑1 (MQ=255)
                                                                                                                tCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTaaa                                      >  1:1998623/1‑150 (MQ=255)
                                                                                                                tCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTaaa                                      >  1:2330092/1‑150 (MQ=255)
                                                                                                                 cAGACCTGCAGCAATTGCCAAATGTAGGCCGGATAAGGCGTGTACGCCGCATCCTGCATTTGCTTAGGTCAGAGTTAATATTACCCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAAc                                     <  2:1511355/150‑1 (MQ=255)
                                                                                                                  aGACCTGGAGCAATTGCCAAATGTATGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTACCTGCTTACCAGCGTTAAACg                                    <  1:2039348/150‑1 (MQ=255)
                                                                                                                  aGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACg                                    >  1:2299507/1‑150 (MQ=255)
                                                                                                                   gACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACgg                                   >  1:1808589/1‑150 (MQ=255)
                                                                                                                   gACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACgg                                   >  1:2010463/1‑150 (MQ=255)
                                                                                                                    aCCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGa                                  <  2:672335/150‑1 (MQ=255)
                                                                                                                        ggAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTg                              <  1:1021483/150‑1 (MQ=255)
                                                                                                                          aGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCa                            <  1:2240243/150‑1 (MQ=255)
                                                                                                                           gCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCAt                           >  1:1606235/1‑150 (MQ=255)
                                                                                                                               ttGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCgg                       <  2:838548/150‑1 (MQ=255)
                                                                                                                               ttGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCgg                       <  1:2455723/150‑1 (MQ=255)
                                                                                                                                tGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCggg                      <  2:650859/150‑1 (MQ=255)
                                                                                                                                 gCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGc                     <  1:1755174/150‑1 (MQ=255)
                                                                                                                                 gCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGc                     >  2:1315745/1‑150 (MQ=255)
                                                                                                                                  ccAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCt                    >  1:466645/1‑150 (MQ=255)
                                                                                                                                   cAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTa                   <  2:276101/150‑1 (MQ=255)
                                                                                                                                   cAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTa                   >  2:1912814/1‑150 (MQ=255)
                                                                                                                                    aaaTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTCAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTAc                  >  1:446019/1‑150 (MQ=255)
                                                                                                                                         tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTaaa             <  1:2186156/150‑1 (MQ=255)
                                                                                                                                          aGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGATCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAat            <  1:1403957/150‑1 (MQ=255)
                                                                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAata           >  2:2469995/1‑150 (MQ=255)
                                                                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAata           >  2:1069534/1‑150 (MQ=255)
                                                                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTg                              <  2:240814/130‑1 (MQ=255)
                                                                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTg                              >  1:240814/1‑130 (MQ=255)
                                                                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAatat          >  1:2066157/1‑150 (MQ=255)
                                                                                                                                             ccGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATc         >  1:5526/1‑150 (MQ=255)
                                                                                                                                               ggATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCg       >  1:493070/1‑150 (MQ=255)
                                                                                                                                               ggATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCg       <  2:2124401/150‑1 (MQ=255)
                                                                                                                                               ggATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCg       >  2:1641159/1‑150 (MQ=255)
                                                                                                                                                  tAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGcac    <  1:114180/150‑1 (MQ=255)
                                                                                                                                                    aGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGCacac  >  1:1230435/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GACGATGTCGACAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGCACAC  >  NC_000913/4460329‑4460626

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