Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,392,328 A→C 6.7% intergenic (‑300/+130) tldD ← / ← yhdP putative peptidase/DUF3971‑AsmA2 domains protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,392,3280AC 6.7% 1188.5 / 11.1 434intergenic (‑300/+130)tldD/yhdPputative peptidase/DUF3971‑AsmA2 domains protein
Reads supporting (aligned to +/- strand):  ref base A (263/121);  new base C (2/39);  total (266/168)
Fisher's exact test for biased strand distribution p-value = 4.32e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 447 total aligned reads displayed.

GTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTAC  >  NC_000913/3392180‑3392477
                                                                                                                                                     |                                                                                                                                                      
gTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAg                                                                                                                                                      >  2:3494062/1‑150 (MQ=255)
   aaCGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATACGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCGGATCCGACCTACGGGTATGTTCCGGGTGTAGGCCTGATCAGACGCACAGCgg                                                                                                                                                   >  2:4616266/1‑149 (MQ=14)
    aCGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTc                                                                                                                                                  >  2:4692218/1‑150 (MQ=255)
    aCGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTACGCCTGATAAGACGCACAGCGTc                                                                                                                                                  >  2:2571009/1‑150 (MQ=255)
    aCGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCATATCCGACCTACGGTTATGTTCAGTATGTAGGCCTGATAAGACGcacatagtc                                                                                                                                                   >  2:2448328/1‑145 (MQ=255)
     ccccTTATCCGACCTACGGTTCTGTTCCGTTTGGTGGCCTTATAAGACGCACTGGGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCg                                                                                                                                                 <  2:1124037/148‑1 (MQ=14)
     cGCCTTAACCGACCTACGTTTATGTTCCGTTTGTCGGCCTTATAAGACGCCCAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCg                                                                                                                                                 <  1:2039003/150‑1 (MQ=14)
      gCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCATATCCGACCTACGGTTATGTTCCGTATCTATGCCTTATAAGAAGCACAgacaaga                                                                                                                                                >  1:1451597/1‑144 (MQ=14)
         ttATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATc                                                                                                                                             >  2:6228661/1‑150 (MQ=255)
         ttATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATc                                                                                                                                             >  1:4272226/1‑150 (MQ=255)
         ttATCCGACCTACGGTTATGTTCCGTATGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTCCGGTTATGTTCCGTTTGTAGGCCTGTTAAGACGCACAGCGTCGCATc                                                                                                                                             >  2:5032569/1‑150 (MQ=255)
         ttATCCGACCTACGGTTATGTGCCGTTTGTAGGCCTGATAAGACGCACAGCGGCGCATCAGGCCACGGCTGCCGGAGGCGGCGTAAACGCCGTAGCCGACCTACGGTTAGGTTCCGTGTGTAGGCCTGATCAGACGCACCGCGGCGCATc                                                                                                                                             >  2:769826/1‑150 (MQ=14)
          tttCCGACATTGGGTTGTGTTCCGTTTTTTGCCCTTAAAAGACGCCCAGCGTCGCATATGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCa                                                                                                                                            <  1:49103/148‑1 (MQ=11)
          tATCCGACCTCCGTTTAGGTTCCGTTTGTAGCCCTGATAAGACGCACAGCGTCGCATCCGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCa                                                                                                                                            <  1:925272/150‑1 (MQ=14)
          tATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGCCGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCa                                                                                                                                            <  2:3261439/150‑1 (MQ=255)
           aTCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAg                                                                                                                                           >  1:2355934/1‑150 (MQ=40)
            tCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAgg                                                                                                                                          >  1:5550243/1‑150 (MQ=40)
                       gTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGcc                                                                                                                               >  1:6195756/1‑150 (MQ=35)
                       gTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGcc                                                                                                                               >  2:6053867/1‑150 (MQ=35)
                       gTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCAGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGcc                                                                                                                               >  1:551433/1‑150 (MQ=25)
                             ttCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgc                                                                                                                         >  2:1784631/1‑150 (MQ=35)
                             ttCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAAAGGCTGCCGGAGGCGGCGTAAACGCCGTATCCGACCTACGGTTATGTTCCGTGTGGAGGCCGGATACGACGCACTGCGTCGCATCAGTCCACGGCTGCCGGATgc                                                                                                                         >  2:4126841/1‑150 (MQ=11)
                             ttCCGTTTGTAGGCCGGATACGACGCACAGCGTCGCATCAGGAAACGGCTGCCGTATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCTGGCAACGGCTGCCGGATgc                                                                                                                         <  2:430761/150‑1 (MQ=11)
                                        ggCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGc                                                                                                              <  2:4832063/150‑1 (MQ=34)
                                         gCCCGATAAGACGCACAGCGTCGCATCAGGCACCGGCTGCCGGATGCGGCGTAAACGCCTTACCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGcc                                                                                                             <  2:3809711/150‑1 (MQ=255)
                                              ttAAGACGCTCAGCGTCGCCTCTGTCAACGGGTGCCGGATGCGGCGTCACCGCCTTAGCCGCCCTACGTTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATc                                                                                                        <  2:3159110/149‑1 (MQ=11)
                                               tAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTCTGTTCCGGGTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATcc                                                                                                       >  2:26444/1‑150 (MQ=255)
                                               tAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTCTGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATcc                                                                                                       >  1:1909439/1‑150 (MQ=37)
                                                 acacccacaGCCTCGCAACAGCAAACGGCTGCCGGATGCGGCGTAAACGCATTATCCGACATACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATAAGGAAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGa                                                                                                     <  2:1932686/145‑1 (MQ=11)
                                                   aCGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGAcc                                                                                                   >  2:5435449/1‑150 (MQ=37)
                                                    cGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGCAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCt                                                                                                  >  2:219746/1‑150 (MQ=37)
                                                        caGCTTCGCATCTGGAAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACgg                                                                                              <  2:6260018/150‑1 (MQ=255)
                                                            gTCGAATCTGGAAACGGCTGCCGGATGCGGCTTAAACGCCTTATCCTACCTAGGGTTTTTTTCCTTTTGTCGGCCTGATAAGACGCACTGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTAt                                                                                          <  1:4312108/150‑1 (MQ=11)
                                                              cGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGt                                                                                        >  1:5204141/1‑150 (MQ=37)
                                                               gCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGtt                                                                                       >  1:6013261/1‑150 (MQ=37)
                                                               gCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGtt                                                                                       >  1:3380887/1‑150 (MQ=37)
                                                                 aTCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCt                                                                                     >  1:4986463/1‑150 (MQ=38)
                                                                      gCATCGGCTGCCGTATGCGGCTTAAACGCGTTATCCTCCCTAGGGTTTTGTTCCGTTTGTAGGCCTGATAAGACTCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTc                                                                                <  1:3986033/150‑1 (MQ=14)
                                                                      gCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCTACCTACGGTTAGGTTCCGTTTGTAGGCCTGATAAGACACACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTc                                                                                <  1:857654/150‑1 (MQ=255)
                                                                       cAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCg                                                                               >  1:6503373/1‑150 (MQ=255)
                                                                       cAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCg                                                                               >  2:1292234/1‑150 (MQ=255)
                                                                            gcggCCGTATGCGGCGTAAACGCCTTATCCGACCTACGTTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGc                                                                          <  2:3405476/147‑1 (MQ=255)
                                                                                ccGGTTGCGGCGTAAACGCCTTTTCCTACCTACGGTTTTTTTCCGTTTGTAGGCCTGATAAGACGCACTGCGTCGCATCAGGCAAGGGCTGCCGTATGCGGCGTGAACGCCCTATCCGACCTACGTTTTTGTTCTGGCTCGTTGGCCttt                                                                      <  2:4384216/150‑1 (MQ=17)
                                                                                ccGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTTGGCCTGATAAGACGCACCGCGTCGCATCTGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCttt                                                                      <  2:5391846/150‑1 (MQ=255)
                                                                                ccGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCttt                                                                      >  2:1519284/1‑150 (MQ=255)
                                                                                ccGGACGCGGCGTAAACGCCTTAACCGCCCTACGGTTTTGTTCCGTTTGTAGGCCTGATAAGTCGCCCAGCGGCGCACCAGGCATCTGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCttt                                                                      <  2:6465657/150‑1 (MQ=17)
                                                                                 cggtggcggcgTAACCGCCTTACCCGACCTCCGTTTATGTTCCGTTTGTCGGCCTGATAAGACGCACCGGGTCGCATCTGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTg                                                                     <  2:2826143/145‑1 (MQ=17)
                                                                                 cGGAGGCGGCTTAAACGCCTTATCCGACCTAGGGTTATGTTCCGTTTGTAGGCCTGATAAGACCCACAGCGTCGAATCAGCCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTg                                                                     <  1:5936746/150‑1 (MQ=17)
                                                                                 cGGACGCGGCGAAAACGCCTTTCCCGACCGCCGGTTCTGTTCGGTTTGTAGGCATGATAAGACGCACAGCGTCGCATCCGGCAACGGGTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTg                                                                     <  2:235474/150‑1 (MQ=17)
                                                                                  ggATGCGGCGTAAACGCCTTATCCGACCTAGGGTTTTGTTCGGTTTGTCGGCCTGATAAGACGCACAGCGTCGCATCAGGACACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTgg                                                                    <  1:5645988/150‑1 (MQ=17)
                                                                                  ggATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTgg                                                                    <  2:3567477/150‑1 (MQ=255)
                                                                                  ggAGGCGGCTTAAACGCATTAACCGACATCCGTTTATGGGCGTTTTGTTGGCCTGAAATGACGACATGCGTCTCATCTGCCTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTgg                                                                    <  1:170926/150‑1 (MQ=12)
                                                                                     tGCGGCGTAAACGCCTTATCCGACATCGGTTTTTTTTCGTTTTGTCGGCCTTATAAGACGCACCGCGTCGCAACAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCaa                                                                 <  1:1219087/150‑1 (MQ=17)
                                                                                         gcgTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGAt                                                             >  1:5677321/1‑150 (MQ=255)
                                                                                             aaaCGCCTTATCCTACCTACGGTTTTGTTCGTTTTGTTGGCCTGATAAGACGCACCGCGTCGCAACTGGAAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATc                                                         <  1:6389543/150‑1 (MQ=18)
                                                                                             aaaCGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACAGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATc                                                         <  1:2740489/150‑1 (MQ=255)
                                                                                                  ccTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATg                                                    >  1:1023286/1‑150 (MQ=255)
                                                                                                   cTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCACCGATTATCCTATga                                                   >  2:4804692/1‑150 (MQ=255)
                                                                                                   cTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATga                                                   >  2:3654160/1‑150 (MQ=255)
                                                                                                   cTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCATTGGCAACGATTATCCTATga                                                   >  2:5710254/1‑150 (MQ=255)
                                                                                                    ttATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGTGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATgag                                                  >  1:785588/1‑150 (MQ=255)
                                                                                                    ttATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATgag                                                  >  2:4494807/1‑150 (MQ=255)
                                                                                                    ttATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCATATgag                                                  >  1:626539/1‑150 (MQ=255)
                                                                                                    ttATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCATATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCAATgag                                                  >  1:2490012/1‑150 (MQ=255)
                                                                                                        ccGACCTCGGGTTTTGTCCGGTTTTTAGTCCTGTTAAGACGCCATGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTg                                              <  1:3149789/150‑1 (MQ=21)
                                                                                                        ccGAACTACGGTTATGTTCCGTTTGTAGGCCTGTTAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTg                                              <  1:2196322/150‑1 (MQ=255)
                                                                                                         cgcccTACGGTTATTTTCCGTTTGTAGGCCTGATAAGACGCCCCGCGTCGCTTCTGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTgg                                             <  1:5159613/147‑1 (MQ=21)
                                                                                                         cGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTgg                                             >  2:1540962/1‑150 (MQ=255)
                                                                                                          gACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTggg                                            >  1:4565596/1‑150 (MQ=255)
                                                                                                          gACATACGGTTTTTTTCCGTTTGTAGGCCTGATAAGCCGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTggg                                            <  2:2958154/150‑1 (MQ=255)
                                                                                                           aCCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTgggg                                           >  1:2632178/1‑150 (MQ=255)
                                                                                                            ccaaCGGTTATTTACCGTTTGTCGGCCTGATAAGACGCACACCGTCGAAACAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGc                                          <  1:3498645/147‑1 (MQ=21)
                                                                                                                cGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAAtt                                      >  1:3786857/1‑150 (MQ=255)
                                                                                                                cGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAAtt                                      >  1:6055818/1‑150 (MQ=255)
                                                                                                                  gTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTAc                                    >  1:1245008/1‑150 (MQ=255)
                                                                                                                   ttATGTCCCGTTTGTTCTCCTGTGAAGCCGCACAGCGTCGGATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTAcg                                   <  2:3995488/150‑1 (MQ=25)
                                                                                                                    tttGTTCCGTTTGTCGGCCTTATAAGACGCACAGGGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTAcgc                                  <  2:557854/148‑1 (MQ=255)
                                                                                                                    tATGTTCCGTTTGTTGGCCTGATAAGCCGCACCGCGTCGCATCTGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTAcgc                                  <  1:1855945/150‑1 (MQ=255)
                                                                                                                    tATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTAcgc                                  <  2:717043/150‑1 (MQ=255)
                                                                                                                     aTGTTCGGTTTGTAGGCCTGATAAGACGCACCGCGTGGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTAcgcg                                 <  1:4132042/150‑1 (MQ=255)
                                                                                                                        ttCCGTTTGTAGGCCTGCTAAGACGCACAGCGGCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGccc                              <  1:5999575/150‑1 (MQ=255)
                                                                                                                        ttCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGccc                              >  1:2144523/1‑150 (MQ=255)
                                                                                                                        tccggTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTTAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGccc                              <  1:6486539/146‑1 (MQ=255)
                                                                                                                         tCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCt                             >  1:5820218/1‑150 (MQ=255)
                                                                                                                            gTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGt                          <  1:4757905/150‑1 (MQ=255)
                                                                                                                                 gTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCaaa                      <  1:5394180/150‑1 (MQ=255)
                                                                                                                                     gCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCAt                  >  2:5235891/1‑150 (MQ=255)
                                                                                                                                      ccTTATAAGCCGCACCGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCAtt                 <  1:1750619/150‑1 (MQ=255)
                                                                                                                                         gATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGGTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTgcg              >  1:2891647/1‑150 (MQ=255)
                                                                                                                                          ctAAGCCGCACCGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTgcgc             <  2:4165235/149‑1 (MQ=255)
                                                                                                                                              gACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCtttt         >  2:5776155/1‑150 (MQ=255)
                                                                                                                                              gACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCtttt         >  1:4637273/1‑150 (MQ=255)
                                                                                                                                               aCGCACCGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCttttt        <  1:3040585/150‑1 (MQ=255)
                                                                                                                                                 gCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCt      >  2:2056041/1‑150 (MQ=255)
                                                                                                                                                 gCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCt      <  1:1642116/150‑1 (MQ=255)
                                                                                                                                                 gCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCGTTGGCAACGATTCTCCTATGAGAGTGGGGCACTTACGCGCCCGCGGCACATCATGGCGCTTTTTCt      >  2:1258382/1‑150 (MQ=25)
                                                                                                                                                  cacaGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCACATCATTGCGCTTTTTCtt     >  1:947884/1‑150 (MQ=255)
                                                                                                                                                   acaGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCttt    <  2:4210462/150‑1 (MQ=255)
                                                                                                                                                   acaGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCttt    <  2:5635802/150‑1 (MQ=255)
                                                                                                                                                     aGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTAc  >  1:4506088/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTAC  >  NC_000913/3392180‑3392477

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