Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,714,359 A→C 6.0% intergenic (+129/+80) 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,3590AC 6.0% 261.7 / 7.8 83intergenic (+129/+80)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 C (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 = 7.70e-01

AAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAACCGCTTTTCC  >  NC_000913/2714214‑2714507
                                                                                                                                                 |                                                                                                                                                    
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:1384310/150‑1 (MQ=255)
   cGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTc                                                                                                                                               <  1:1935237/150‑1 (MQ=255)
      cAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcg                                                                                                                                            >  1:601024/1‑150 (MQ=255)
      cAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcg                                                                                                                                            >  2:958341/1‑150 (MQ=255)
        gaagaGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAg                                                                                                                                                     >  1:1179670/1‑139 (MQ=255)
        gaagaGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAg                                                                                                                                                     <  2:1179670/139‑1 (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:375199/150‑1 (MQ=255)
                    cAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATg                                                                                                                              >  2:703700/1‑150 (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                                                                                                                                 >  1:1240092/1‑126 (MQ=255)
                                         tATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                 <  2:1240092/126‑1 (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:1800102/1‑150 (MQ=255)
                                                                     gCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACa                                                                             >  2:380805/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)
                                                                                                                                                 |                                                                                                                                                    
AAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTAAATGCCGGGTTTAATCCCGACATTATTTTTATCCAGCCACAAAACTCTCTTCAACACACTGAATAAAAGTGTGCATCGCCGGGCTAACCGCTTTTCC  >  NC_000913/2714214‑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: