Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,714,358 A→G 6.0% intergenic (+128/+81) srmB → / ← yfiE ATP‑dependent RNA helicase SrmB/putative LysR‑type DNA‑binding transcriptional regulator YfiE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,714,3580AG 6.0% 254.4 / 5.1 83intergenic (+128/+81)srmB/yfiEATP‑dependent RNA helicase SrmB/putative LysR‑type DNA‑binding transcriptional regulator YfiE
Reads supporting (aligned to +/- strand):  ref base A (43/35);  new base G (2/3);  total (45/38)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

ACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAACCGCTTTTCC  >  NC_000913/2714209‑2714507
                                                                                                                                                     |                                                                                                                                                     
aCCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATa                                                                                                                                                       <  2:1870836/150‑1 (MQ=255)
     aaaCGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCg                                                                                                                                                  <  2:699082/150‑1 (MQ=255)
      aaCGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGt                                                                                                                                                 <  1:1951255/150‑1 (MQ=255)
        cGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGGTGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTc                                                                                                                                               <  1:131156/150‑1 (MQ=255)
        cGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTc                                                                                                                                               <  1:76662/150‑1 (MQ=255)
        cGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTc                                                                                                                                               <  1:1935237/150‑1 (MQ=255)
        cGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTc                                                                                                                                               <  1:1384310/150‑1 (MQ=255)
           cAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcg                                                                                                                                            >  1:601024/1‑150 (MQ=255)
           cAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcg                                                                                                                                            >  2:958341/1‑150 (MQ=255)
             gaagaGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAg                                                                                                                                                     <  2:1179670/139‑1 (MQ=255)
             gaagaGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAg                                                                                                                                                     >  1:1179670/1‑139 (MQ=255)
              aatagTTATCTCAATTCCAGTTTTAAGCCTGTTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcgccg                                                                                                                                         <  1:846874/147‑1 (MQ=21)
              aagaGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcgccg                                                                                                                                         <  1:733713/150‑1 (MQ=255)
                   tAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                    <  2:928797/150‑1 (MQ=25)
                     aTCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                  <  1:162911/150‑1 (MQ=25)
                       ctcAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGTAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                                                                >  2:990521/1‑150 (MQ=255)
                        tcAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTATGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcct                                                                                                                               <  1:936101/150‑3 (MQ=21)
                         cAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATg                                                                                                                              >  2:703700/1‑150 (MQ=255)
                         cAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATg                                                                                                                              <  2:375199/150‑1 (MQ=255)
                          aaTGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGa                                                                                                                             >  1:411591/1‑150 (MQ=255)
                          aaTGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGACGCATCCGGCATGa                                                                                                                             <  2:1462161/150‑1 (MQ=255)
                               cAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATaaa                                                                                                                        >  1:1507733/1‑150 (MQ=255)
                                 ggTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAgc                                                                                                                      >  1:1283419/1‑150 (MQ=255)
                                            ggTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCg                                                                                                           <  1:1301353/150‑1 (MQ=255)
                                              tATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                 <  2:1240092/126‑1 (MQ=255)
                                              tATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                 >  1:1240092/1‑126 (MQ=255)
                                               aTTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACa                                                                                                        >  1:455262/1‑150 (MQ=255)
                                                 tAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAg                                                                                                      >  2:1344056/1‑150 (MQ=255)
                                                   aaaGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGtt                                                                                                    <  2:236753/150‑1 (MQ=255)
                                                       tGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAAt                                                                                                <  1:112276/150‑1 (MQ=255)
                                                         cTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGc                                                                                              >  1:1163709/1‑150 (MQ=255)
                                                         cTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGc                                                                                              >  1:1954641/1‑150 (MQ=255)
                                                           ggCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCg                                                                                            >  2:1356256/1‑150 (MQ=255)
                                                           ggCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCg                                                                                            >  2:1448487/1‑150 (MQ=255)
                                                               aaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGtt                                                                                        >  2:356031/1‑150 (MQ=255)
                                                                  gCAAAACTGCCTGATGAGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCAGCATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTaa                                                                                     >  2:5300/1‑150 (MQ=255)
                                                                   cAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAAt                                                                                    >  2:815502/1‑150 (MQ=255)
                                                                      aaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATccc                                                                                 <  2:1302874/150‑1 (MQ=255)
                                                                          gCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACa                                                                             >  2:380805/1‑150 (MQ=255)
                                                                          gCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACa                                                                             >  2:1800102/1‑150 (MQ=255)
                                                                             tGATGCGATACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACatta                                                                          >  2:1890170/1‑150 (MQ=255)
                                                                                 gcgcTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTAtttt                                                                      >  1:736676/1‑150 (MQ=255)
                                                                                  cgcTACGCTTATCAGGCCTACATGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTAttttt                                                                     >  1:235089/1‑150 (MQ=255)
                                                                                    cTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTAt                                                                   <  2:1163709/150‑1 (MQ=255)
                                                                                    cTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTAt                                                                   >  2:1832928/1‑150 (MQ=255)
                                                                                     tACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATc                                                                  >  1:1982154/1‑150 (MQ=255)
                                                                                     tACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATc                                                                  <  2:411591/150‑1 (MQ=255)
                                                                                      aCGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATcc                                                                 >  1:609327/1‑150 (MQ=255)
                                                                                        gCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAg                                                               <  2:1379648/150‑1 (MQ=255)
                                                                                        cgctATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAg                                                               >  2:536374/4‑150 (MQ=255)
                                                                                           tATCAGGCCTACGTGGCTCATGTAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCca                                                            <  1:990521/150‑1 (MQ=255)
                                                                                             tCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCcaca                                                          <  2:609327/150‑1 (MQ=255)
                                                                                               aGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACaaa                                                        >  1:354132/1‑150 (MQ=255)
                                                                                                 gCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAAc                                                      <  2:736676/150‑1 (MQ=255)
                                                                                                   cTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAActc                                                    >  2:1652515/1‑150 (MQ=255)
                                                                                                      cGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAActctct                                                 >  1:1053399/1‑150 (MQ=255)
                                                                                                       gTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCtt                                                <  1:1946628/150‑1 (MQ=255)
                                                                                                          gCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCaa                                             <  2:584748/150‑1 (MQ=255)
                                                                                                           cTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAac                                            >  1:878263/1‑150 (MQ=255)
                                                                                                            tCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAACGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAaca                                           <  2:1507733/150‑1 (MQ=255)
                                                                                                               tGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAacacac                                        >  2:1511157/1‑150 (MQ=255)
                                                                                                                 cAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTg                                      <  1:1517798/150‑1 (MQ=255)
                                                                                                                         tGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAg                              <  2:354132/150‑1 (MQ=255)
                                                                                                                            aTTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAAt                                   >  2:401860/1‑142 (MQ=255)
                                                                                                                            aTTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAAt                                   <  1:401860/142‑1 (MQ=255)
                                                                                                                             tttGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGtgtg                          <  1:1890170/150‑1 (MQ=255)
                                                                                                                                      tCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCgg                 >  2:1842739/1‑150 (MQ=255)
                                                                                                                                       cTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCggg                <  1:280736/150‑1 (MQ=255)
                                                                                                                                       cTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCggg                >  2:62108/1‑150 (MQ=255)
                                                                                                                                         tGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCt              >  2:547263/1‑150 (MQ=255)
                                                                                                                                         tGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCt              >  1:952769/1‑150 (MQ=255)
                                                                                                                                          tcatgCCGGATGCGGCGTGAACGCCTTATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTa             >  1:1390241/5‑150 (MQ=38)
                                                                                                                                           tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTaa            >  2:511983/1‑150 (MQ=255)
                                                                                                                                           catgCCGGATGCGGCGTGAACGCCTTATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTaa            >  2:1112985/4‑150 (MQ=38)
                                                                                                                                            aGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAAc           <  1:703700/150‑1 (MQ=255)
                                                                                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATTCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAAcc          <  2:1954641/150‑1 (MQ=255)
                                                                                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAAcc          <  2:455262/150‑1 (MQ=255)
                                                                                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAAcc          <  1:958341/150‑1 (MQ=255)
                                                                                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAAcc          >  2:1432878/1‑150 (MQ=255)
                                                                                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAACCg         <  1:1511157/150‑1 (MQ=255)
                                                                                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCAGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAACCg         >  1:1008095/1‑150 (MQ=255)
                                                                                                                                                  gATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAACCGCttt     >  2:993145/1‑150 (MQ=255)
                                                                                                                                                   aTAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAACCGCtttt    >  2:429382/1‑150 (MQ=255)
                                                                                                                                                    tAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAACCGCTTTTc   >  1:415177/1‑150 (MQ=255)
                                                                                                                                                     aaGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAACCGCTTTTcc  <  1:543534/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAACCGCTTTTCC  >  NC_000913/2714209‑2714507

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: