Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,392,201 C→G 5.7% intergenic (‑173/+257) 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,2010CG 5.7% 1355.9 / 8.3 450intergenic (‑173/+257)tldD/yhdPputative peptidase/DUF3971‑AsmA2 domains protein
Reads supporting (aligned to +/- strand):  ref base C (229/185);  new base G (2/31);  total (232/218)
Fisher's exact test for biased strand distribution p-value = 1.22e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 456 total aligned reads displayed.

TGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCT  >  NC_000913/3392052‑3392349
                                                                                                                                                     |                                                                                                                                                    
  cagtttcggATTTTTTACGGGTGGCGGTTGGGGGGTGAAAGGCCTTTCCAGACTTCGGTTTTTGTACGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCTGGCAACGGCTGGCGGATGCGTCGTAAACGCCTTATCCGACCTACgg                                                                                                                                                    <  2:1353599/143‑1 (MQ=9)
  cAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACgg                                                                                                                                                    >  1:2641681/1‑150 (MQ=255)
       gagTATTTTTTCCGGCTGCCGGATGCGGCGTGAACGCCCTAACCAGCCTACGTTTATGTTCCGGTTTGTCGGCCTGATACGACGCACAGCGTCGCACCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATg                                                                                                                                               <  1:466437/150‑1 (MQ=17)
       gagGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATg                                                                                                                                               >  2:4338275/1‑150 (MQ=255)
       gagGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTAGCCGACCTAGGGTtgtg                                                                                                                                               >  2:4117565/1‑150 (MQ=255)
       gagGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGGTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATg                                                                                                                                               >  1:1251938/1‑150 (MQ=255)
           cttttttACGGCTGCCGGATGCGGCGTGAACGCCCTAACCAGCCTACGTTTATGTTCCGGTTTGTAGGCCTGATAAGACGCCCAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTcc                                                                                                                                           <  2:2531472/149‑1 (MQ=255)
            ttttttACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCg                                                                                                                                          <  2:908450/150‑1 (MQ=255)
               tttACGGCTGCCGGATGCGGCGTGAACGCCCTCTCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGttt                                                                                                                                       <  2:1430523/150‑1 (MQ=255)
                ttACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTg                                                                                                                                      <  2:5103382/150‑1 (MQ=255)
                  aCGGCTGCCGGCTGGGGCGTGAACGCCCTAACCAGCCTACGGTTATGTTCGGGTTGGTAGGCCTTATAAGCCGCACAGCGTCGCATCCGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTa                                                                                                                                    <  1:22409/150‑1 (MQ=17)
                  aCGGCGGCCGGTTGCGGCGTGACCGCCCTACCCTGCCTACGTTTTTGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTa                                                                                                                                    <  2:504484/150‑1 (MQ=14)
                   cGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCGGGTTTGTTGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAg                                                                                                                                   <  2:2211290/150‑1 (MQ=255)
                    ggCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTGATCCGACCTCCGGTTATGTTCCGTTTGTAgg                                                                                                                                  >  1:3513909/1‑150 (MQ=255)
                    ggCTGCCGGAGGCGGCGTGAACGCCCTATCCAGCCTACGGTTTTGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAgg                                                                                                                                  <  1:3239885/150‑1 (MQ=255)
                     gCTGCCGGATGCGGCGTTAACGCCCTATCCAGCCTACGTTTATGTTCCGGTTTTTTGGCCTTATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGc                                                                                                                                 <  1:1825132/150‑1 (MQ=14)
                     gCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTTTGTTCCGGTTTTTAGGCCTGATAAGACGCACCGCGTCGAATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGc                                                                                                                                 <  1:1735976/150‑1 (MQ=255)
                     gCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGc                                                                                                                                 <  2:3357553/150‑1 (MQ=255)
                       tGCCGGAGGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTCGGCCTGATACGACGCACCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCt                                                                                                                               <  1:250481/150‑1 (MQ=255)
                            gATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTGTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGAGGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTATGTAGGCCTGATaa                                                                                                                          >  2:2889240/1‑150 (MQ=255)
                                cggcgTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACg                                                                                                                      >  1:3474889/1‑150 (MQ=255)
                                 ggcgTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCAGTGTGTAGGCCTGATAAGACGc                                                                                                                     >  2:3416650/1‑150 (MQ=255)
                                  gcgTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCGTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGca                                                                                                                    >  1:789641/1‑150 (MQ=255)
                                       aaCGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCg                                                                                                               >  1:670926/1‑150 (MQ=255)
                                        aCGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGt                                                                                                              <  2:4223766/150‑1 (MQ=255)
                                        aCGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTGATGTTCCGTGAGTAGGCATGATAAGACGCACAGCGt                                                                                                              >  1:2572344/1‑150 (MQ=255)
                                         cGCCCTATCCTGCCTCGGGTTTTGTGCGGGTTTGTCGGCCTTATAAGACGCACAGGGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTc                                                                                                             <  1:4878214/150‑1 (MQ=14)
                                         cGCCCTATCCACCCTACGCTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTc                                                                                                             <  2:4129450/150‑1 (MQ=255)
                                          gCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCg                                                                                                            >  2:1697229/1‑150 (MQ=255)
                                            ccTATCCAGCCAACGTTTATGTTCCGGTTTGTAGGCCTGACAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCa                                                                                                          <  1:2046651/150‑1 (MQ=255)
                                             cTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCAt                                                                                                         <  1:1382623/150‑1 (MQ=255)
                                             cTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGGCGGATGCGGCGGAAACGCCTTATCCGCCCTACGGGGATGTTCCGTTTGTAGGCCTGGTAAGACGCACAGCGTCGCag                                                                                                         >  1:3026270/1‑149 (MQ=14)
                                               aTCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTGATGTTCCGGTTGTAGGCCTGATAAGACGCACAGCGTCGCATCa                                                                                                       >  2:2749259/1‑150 (MQ=255)
                                                tCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCATATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATACGACGCAAAGCGACGCATCAg                                                                                                      >  2:618972/1‑150 (MQ=255)
                                                 ccAGCGTAGGTTTTTGTTCGGTTTTTTTGGCCTGATAAGACGCACAGCTTCGCATCTGGCAACGGCTGTCTGATGCGGCGTAAACGCCTTATCCGACCTACGGTTTTGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAgg                                                                                                     <  1:5058939/150‑1 (MQ=11)
                                                 ccAGCCTACGGTTTTGTTCCGGTTTGTCGGCCTGATAAGACGCCCAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAgg                                                                                                     <  1:5465807/150‑1 (MQ=255)
                                                 ccAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAgg                                                                                                     <  2:3366064/150‑1 (MQ=40)
                                                  cAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGACGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGc                                                                                                    >  1:50044/1‑150 (MQ=37)
                                                   aGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTATGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCa                                                                                                   >  1:5796846/1‑150 (MQ=38)
                                                            ttctGTTCCGGTTTGTAGGCCTGATAAGACGCCCAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGcc                                                                                          <  2:3772897/147‑1 (MQ=255)
                                                               tGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGCCGCATCCGGCAGCCGTTGCCTGa                                                                                       <  2:2276154/150‑1 (MQ=11)
                                                                 ttCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATg                                                                                     >  1:3612264/1‑150 (MQ=35)
                                                                  tCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTGTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGACGGATgc                                                                                    >  2:2001636/1‑150 (MQ=25)
                                                                      gTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCACCGGCTGCCGGAGgcggcg                                                                                >  1:1636963/1‑150 (MQ=17)
                                                                           tAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAAc                                                                           >  1:4240945/1‑150 (MQ=34)
                                                                             ggCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGc                                                                         <  2:4376538/150‑1 (MQ=34)
                                                                              gCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGcc                                                                        >  1:5734793/1‑150 (MQ=34)
                                                                              gCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGcc                                                                        >  2:5609936/1‑150 (MQ=34)
                                                                              gCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGcc                                                                        >  2:5180637/1‑150 (MQ=34)
                                                                              gCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGAAAACGcc                                                                        >  2:4553950/1‑150 (MQ=21)
                                                                                cTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCtt                                                                      >  1:5562694/1‑150 (MQ=35)
                                                                                  taTAAGACGCACTGCGTCGCATCTGGCAACGGCTGTCGTATGCGGCGTAAACGCCTTATCCGACCTACGGTTTTGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTAt                                                                    <  2:1084946/149‑1 (MQ=255)
                                                                                  gATAATACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGTTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCTGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTAt                                                                    <  2:4935022/150‑1 (MQ=255)
                                                                                     taGACGCACAGCGTCGAATCTGGCATCGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                 <  2:4719345/149‑1 (MQ=255)
                                                                                      aGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGa                                                                >  2:4683903/1‑150 (MQ=37)
                                                                                       gACGCACAGCGTCGCATCTGGAAACGGCTGTCGGATGGGGCGTAAACGCCTTATCCGACCTACGTTTTTGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGAc                                                               <  2:2326785/150‑1 (MQ=11)
                                                                                       gACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCGTCAGGCCACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGAc                                                               >  2:5345282/1‑150 (MQ=255)
                                                                                        aCGCACAGCGTCGCATCAGGAAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGAcc                                                              <  2:1893969/150‑1 (MQ=37)
                                                                                          gCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCGTATCCGACCTACGGGGATGTTCCGGTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGAGGCGGCGGAAACGCCTTCTCCGACCTa                                                            >  2:5575537/1‑150 (MQ=11)
                                                                                            acaGCGTCGCATCAGGCCACGGCTGTCGGATGCGGCGTAAACGCCGTATCCGACCTACGGGTATGTTCCGGTTGTAGGCCTGATAAGACGCACAGCGTCGCCTCAGGCACCGGCTGCCGGCGCCGGCGGAAACGCCTTATCCGACCTACg                                                          >  1:1007481/1‑150 (MQ=11)
                                                                                             ccgcgtcgcgtcCGGCAACGGCTGTCGGATGCGGCGTAACCGCCTTATCCGACCTACGGTTAGGTTCCATTTGTCGGCCTGATAAGACGCCCCGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACgg                                                         <  2:5273851/148‑1 (MQ=11)
                                                                                                 gTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCGGTTTGTTGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAAAGCCTTATCCGACCTACGGTTAt                                                     >  1:4397915/1‑150 (MQ=255)
                                                                                                 gTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTAt                                                     <  2:2943026/150‑1 (MQ=37)
                                                                                                 gTCGCATCAGGCAACGGCTGTCGGACGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTGGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTAt                                                     <  2:3511270/150‑1 (MQ=255)
                                                                                                    gCATCTGGAAACGGCTGTCGTTTGCGGCGTAACCGCCTTATCCGACCTACGGTTTTGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGtt                                                  <  2:2663016/150‑1 (MQ=11)
                                                                                                      aTCTGGCCAGGGCTGGGGGTTGCGGGGTAACCGACTTATCCGACCTACGTTTTTGTTCGGTTTGTCGGCCTGATAAGACGCACCGCGTCGCATCTGCCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTcc                                                <  1:4171085/150‑1 (MQ=11)
                                                                                                      aTCTGGCAACGGCTGGGGTTTGCGGCTTAAACGCATTATCCGACCTACGTTTTTTTTCGTTTTGTTGGCCTTATAAGACGCACAGCGTCGCATCTGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTcc                                                <  1:3940779/150‑1 (MQ=11)
                                                                                                         aGGCAACGGCTGTCGGATGCGGCGTTAACGCCGTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGGCGCATCAGGCCACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTGCggg                                             >  1:2131989/1‑148 (MQ=11)
                                                                                                         aGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTACGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATACGACCAACGGTTATGTTCCGtt                                             >  1:5346537/1‑150 (MQ=255)
                                                                                                           gCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTg                                           >  2:4045976/1‑150 (MQ=255)
                                                                                                           gCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCGTATCCGACCTACGGTTATGTTCCGTTTg                                           >  2:1235158/1‑150 (MQ=255)
                                                                                                            cAACGGCTGTCGGATGCGGCGTAAAAGCCTTATCCGACCTACGGTTATGTTCCGTTTTTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACCGCTGCCGGATGCGGCGTAAACGCCTTATACGACCAACGGTTATGTTCCGTTTGt                                          >  1:2883490/1‑150 (MQ=14)
                                                                                                             aaCGGCTGGCGTATGCGGCGTAACCGCCTTATCCGACCTAGGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTa                                         <  2:3216520/150‑1 (MQ=255)
                                                                                                                ggCTGTGGTATGCGGCGTAAACGCCTTATCCGACCTCCGTTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGc                                      <  2:1251938/150‑1 (MQ=255)
                                                                                                                ggCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCGTATCCGACCTACGGTTATGTTCCGTTTGTAGGc                                      >  2:5028832/1‑150 (MQ=255)
                                                                                                                ggCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGGTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTGTGTAGGc                                      >  2:1529325/1‑150 (MQ=255)
                                                                                                                  cTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCt                                    <  2:4583623/150‑1 (MQ=255)
                                                                                                                    gTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGa                                  >  1:3772995/1‑150 (MQ=255)
                                                                                                                    gTCGGATAAGGCGTTTACGCCGCATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCAGa                                  >  2:5504106/1‑150 (MQ=14)
                                                                                                                     tCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGAt                                 >  1:4048863/1‑150 (MQ=255)
                                                                                                                     tCGGATGCGGCGTAAACGCCTTAGCCGCCCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGAt                                 <  2:4834607/150‑1 (MQ=255)
                                                                                                                            cggcgTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGc                          >  2:3046175/1‑150 (MQ=255)
                                                                                                                            cggcgTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCGTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGc                          >  1:4773891/1‑150 (MQ=255)
                                                                                                                               cgTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATGCGACGctgt                       <  1:3348297/150‑4 (MQ=255)
                                                                                                                                 gaaaGGCGTTATCAGACCTTGGTTTATGGTCGGTTTGTCGGCCTGTTAAGACGCCCCGGGTCGCTTATGGCACCGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGc                     <  2:579057/149‑1 (MQ=11)
                                                                                                                                   aaCGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGt                   >  2:3559243/1‑150 (MQ=255)
                                                                                                                                   aaCGCATTATCCGACCTACGTTTTTGTTCCGTTTTTAGGCCTGATATTTCGCACAGCGTCGCATCTGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTTTGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGt                   <  1:2439936/150‑1 (MQ=14)
                                                                                                                                    aCGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCGTATCCGACCTACGGTTATGTTCAGTTTGTAGGCCTGATGAGACGCACAGCGTc                  >  2:2053816/1‑150 (MQ=255)
                                                                                                                                    aCGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACAGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTc                  <  2:4874014/150‑1 (MQ=255)
                                                                                                                                    aCGCCTTATCCCACCTAGGGTTATTTTCCTTTTTTCGGCATTTTAATACGAAAAGCTTCGCATCTGGAAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTc                  <  1:2367942/150‑1 (MQ=14)
                                                                                                                                     cGCCTTATCCGACATACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCACCGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGGTCTGTTCTGTATGGAGGCCTGATTAGACGGACACCGTCg                 >  2:3634257/1‑150 (MQ=14)
                                                                                                                                      gCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGGTTGTAGGCCTGATAAGACGCACAGAGTCGc                >  2:1106580/1‑150 (MQ=255)
                                                                                                                                        ctaagccGACCAACGGTTATGTTCCGTATGTAGGCCTGATAAGACGCACAGCGTCGCATCCGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGCTAAGACGCACAGCGTCGCAt              <  2:778766/145‑1 (MQ=14)
                                                                                                                                        cgtATCAGACCTAGGGTTAGGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCAt              <  2:1820610/148‑1 (MQ=255)
                                                                                                                                        cTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCAt              >  1:5194855/1‑150 (MQ=255)
                                                                                                                                        cTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCAt              >  1:1937463/1‑150 (MQ=255)
                                                                                                                                           caccGCCCTAGGGTTATTTACGGTTTGTCGGCCTGATACGCCGCACAGGGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAg           <  1:1545598/148‑1 (MQ=11)
                                                                                                                                           aTCCGACCTACGTTTATGTTCCGTTTGTAGGCCTCAAAAGACGCACAGCGTCGCATCAGGAAACGGCTGCCGGATGCGGCGTAAACGCCTTAACCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAg           <  2:5662616/150‑1 (MQ=14)
                                                                                                                                                  caaCGTTTTTGTTCGGTTTTTGGGCCTGTTAAGACGCTCAGCGTCGCATCTGGCAACGGCTGCCGTATGCGGCGTAAACGCCTTAACCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACgg    <  2:4514837/148‑1 (MQ=11)
                                                                                                                                                  cTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACgg    >  1:2768793/1‑150 (MQ=35)
                                                                                                                                                     |                                                                                                                                                    
TGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCT  >  NC_000913/3392052‑3392349

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: