Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,460,476 A→G 8.2% intergenic (‑112/+46) 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,4760AG 8.2% 494.3 / 26.3 171intergenic (‑112/+46)nrdG/nrdDanaerobic ribonucleoside‑triphosphate reductase activating protein/anaerobic ribonucleoside‑triphosphate reductase
Reads supporting (aligned to +/- strand):  ref base A (87/70);  new base G (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 = 4.94e-01

Only 100 of 172 total aligned reads displayed.

GACGATGTCGACAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGCACA  >  NC_000913/4460329‑4460625
                                                                                                                                                   |                                                                                                                                                     
gACGATGTCGACAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAg                                                                                                                                                     >  1:966084/1‑150 (MQ=255)
 aCGATGTCGAGATGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCGGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAATTGTAGGCCGGATAAgg                                                                                                                                                    <  2:369951/150‑1 (MQ=255)
       tCGACAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCTTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTa                                                                                                                                              <  1:1446851/150‑1 (MQ=255)
        cGACAGGATAGTACTGATGATAATTCATTATTTCACCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTAc                                                                                                                                             >  1:1345727/1‑150 (MQ=255)
         gACAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGtaaacg                                                                                                                                            >  1:410868/1‑145 (MQ=255)
         gACAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTAcg                                                                                                                                            <  2:1809793/150‑1 (MQ=255)
         gACAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTAcg                                                                                                                                            >  1:296338/1‑150 (MQ=255)
          aCAGGATAGTACTGATGATAATTCATTATTTCACCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                           <  2:1345727/150‑1 (MQ=255)
           cAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                                                          <  2:2176555/150‑1 (MQ=255)
           cAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                                                          <  1:520953/150‑1 (MQ=255)
                tAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAAcgccttatc                                                                                                                                     <  2:2247668/150‑6 (MQ=38)
                 aGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                    <  1:1176602/150‑1 (MQ=38)
                 aGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                    <  2:61084/150‑1 (MQ=38)
                 aGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                    <  2:410868/150‑1 (MQ=38)
                 aGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATcc                                                                                                                                    <  2:92862/150‑1 (MQ=255)
                  gTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCg                                                                                                                                   >  2:2010732/1‑150 (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:2473705/150‑2 (MQ=38)
                     cTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                <  1:1149733/150‑2 (MQ=38)
                     cTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                >  1:982381/1‑150 (MQ=255)
                     cTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                <  2:106999/150‑1 (MQ=255)
                     cTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                <  1:1459655/150‑1 (MQ=255)
                      tgatgaTAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcct                                                                                                                               <  1:1564284/150‑3 (MQ=38)
                       gatgatAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAtt                                                                                                                              >  2:1948757/1‑150 (MQ=255)
                       catgatAATTCATTATTTCACCTGTCTTATCAATCCGCGAGTACGCTGCATCCGGTAAACAATAACTAAAGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAtt                                                                                                                              >  2:1021483/2‑150 (MQ=25)
                             aaTTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTa                                                                                                                        >  1:788560/1‑150 (MQ=255)
                              aTTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAg                                                                                                                       >  1:202793/1‑150 (MQ=255)
                               ttCATTATTTCGCCTGGCTTATCAGTCCGCTTGTAGGCTGCATCCGGTTATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGc                                                                                                                      <  1:2215715/150‑1 (MQ=255)
                               ttCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGc                                                                                                                      >  1:1012017/1‑150 (MQ=255)
                               ttCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGc                                                                                                                      >  2:2455723/1‑150 (MQ=255)
                               ttCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGc                                                                                                                      >  1:1707708/1‑150 (MQ=255)
                                 cATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTc                                                                                                                    >  2:1773050/1‑150 (MQ=255)
                                 cATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTc                                                                                                                    >  2:2111857/1‑150 (MQ=255)
                                 cATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTc                                                                                                                    >  1:1608983/1‑150 (MQ=255)
                                  attattTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCa                                                                                                                   <  1:2047278/150‑1 (MQ=255)
                                    tattTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCaga                                                                                                                 >  1:838548/1‑150 (MQ=255)
                                      tttCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGt                                                                                                               >  1:2320034/1‑150 (MQ=255)
                                        tCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGa                                                                                                             <  1:2336533/150‑1 (MQ=255)
                                           ccTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGa                                                                                                          >  2:1755174/1‑150 (MQ=255)
                                             tGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGTATAAGGCGTTTACGCCGCATCCTTCATTTGCTTAGCTCTTAGTGAAGttt                                                                                                        <  2:1225446/150‑1 (MQ=37)
                                                cTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAAc                                                                                                     <  2:2238580/150‑1 (MQ=255)
                                                 ttATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAAcc                                                                                                    <  1:335271/150‑1 (MQ=255)
                                                 ttATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAAcc                                                                                                    <  1:1948757/150‑1 (MQ=255)
                                                 ttATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGAATTTGCTTAGCTCAGAGTGAAGATTAAcc                                                                                                    <  2:449597/150‑1 (MQ=255)
                                                               gtAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTcc                                                                                      >  2:2004171/1‑150 (MQ=255)
                                                                   gCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCaa                                                                                  >  2:88463/1‑150 (MQ=255)
                                                                   gCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCaa                                                                                  >  2:1426787/1‑150 (MQ=255)
                                                                   gCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGAAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCaa                                                                                  >  2:1079242/1‑150 (MQ=255)
                                                                       cATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGt                                                                              >  2:2308115/1‑150 (MQ=255)
                                                                       cATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGt                                                                              >  2:114180/1‑150 (MQ=255)
                                                                          ccGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTa                                                                           <  1:1096599/150‑1 (MQ=255)
                                                                            ggTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAAc                                                                         >  2:2199349/1‑150 (MQ=255)
                                                                            ggTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAAc                                                                         >  1:474099/1‑150 (MQ=255)
                                                                                aTCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACgcgg                                                                     <  1:2118370/150‑1 (MQ=255)
                                                                                 tCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGAACTGGAGCAATTGACAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACgcggc                                                                    >  2:852358/1‑150 (MQ=255)
                                                                                      gCCTGATGCGACGCTGTCTCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATACGGCATTTTCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTa                                                               <  2:1916342/150‑1 (MQ=255)
                                                                                       ccTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTATCCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTaa                                                              <  2:962104/150‑1 (MQ=255)
                                                                                        cTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAAc                                                             <  1:2004171/150‑1 (MQ=255)
                                                                                          gATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTa                                                                                                       >  1:2218431/1‑106 (MQ=255)
                                                                                          gATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTa                                                                                                       <  2:2218431/106‑1 (MQ=255)
                                                                                          gATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAActt                                                           >  1:61149/1‑150 (MQ=255)
                                                                                          gATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAActt                                                           >  1:530979/1‑150 (MQ=255)
                                                                                            tGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAActtct                                                         >  2:173443/1‑150 (MQ=255)
                                                                                                   cTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCt                                                  >  1:1296630/1‑150 (MQ=255)
                                                                                                         cgTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTACTGCTTACCAg                                            <  2:1997127/150‑1 (MQ=255)
                                                                                                           tCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGTATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGTAGTTTAAACTATCTGCCAATTCCCCAATTTTTTAAAGAGGAGCTTATCTTCTTACTGGTTACCAgag                                          <  2:982381/150‑3 (MQ=37)
                                                                                                              tATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTa                                       <  2:1296630/150‑1 (MQ=255)
                                                                                                               aTCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTaa                                      >  1:2381112/1‑150 (MQ=255)
                                                                                                                 cAGACCTGCAGCAATTGCCAAATGTAGGCCGGATAAGGCGTGTACGCCGCATCCTGCATTTGCTTAGGTCAGAGTTAATATTACCCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAAc                                    <  2:1511355/150‑1 (MQ=255)
                                                                                                                    aCCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGa                                 <  2:672335/150‑1 (MQ=255)
                                                                                                                    aCCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGa                                 <  1:2199349/150‑1 (MQ=255)
                                                                                                                      cTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACg                               >  1:1903578/1‑150 (MQ=255)
                                                                                                                        ggAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTg                             <  1:1021483/150‑1 (MQ=255)
                                                                                                                            cAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATc                         >  2:15831/1‑150 (MQ=255)
                                                                                                                               ttGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCgg                      <  2:838548/150‑1 (MQ=255)
                                                                                                                                tGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCggg                     <  2:650859/150‑1 (MQ=255)
                                                                                                                                 gCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGc                    >  1:1395745/1‑150 (MQ=255)
                                                                                                                                 gCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGc                    >  2:1315745/1‑150 (MQ=255)
                                                                                                                                   cAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTa                  >  2:1912814/1‑150 (MQ=255)
                                                                                                                                   cAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTa                  <  2:276101/150‑1 (MQ=255)
                                                                                                                                    aaaTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTCAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTAc                 >  1:446019/1‑150 (MQ=255)
                                                                                                                                      aTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCt               <  2:1565648/150‑1 (MQ=255)
                                                                                                                                        gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTaa             <  1:2010732/150‑1 (MQ=255)
                                                                                                                                         tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTaaa            <  1:1426787/150‑1 (MQ=255)
                                                                                                                                         tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTaaa            <  1:2054769/150‑1 (MQ=255)
                                                                                                                                         tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTaaa            <  1:1782103/150‑1 (MQ=255)
                                                                                                                                         aatgCCGGATGCGGCGTAAACGCCTTATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTaaa            >  1:1732965/4‑150 (MQ=255)
                                                                                                                                          aGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAAAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAat           >  1:951466/1‑150 (MQ=255)
                                                                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAata          >  2:2469995/1‑150 (MQ=255)
                                                                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTg                             <  2:240814/130‑1 (MQ=255)
                                                                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAatat         >  1:2066157/1‑150 (MQ=255)
                                                                                                                                             ccGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATc        >  1:5526/1‑150 (MQ=255)
                                                                                                                                             ccGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATc        >  1:225760/1‑150 (MQ=255)
                                                                                                                                             ccGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATc        >  1:377576/1‑150 (MQ=255)
                                                                                                                                               ggATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCg      >  2:1641159/1‑150 (MQ=255)
                                                                                                                                                 aTAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGca    >  1:68931/1‑150 (MQ=255)
                                                                                                                                                  tAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGcac   <  1:114180/150‑1 (MQ=255)
                                                                                                                                                   aaGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGcaca  <  1:2440403/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GACGATGTCGACAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGCACA  >  NC_000913/4460329‑4460625

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: