Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,392,201 C→G 7.1% 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 7.1% 658.7 / 2.7 217intergenic (‑173/+257)tldD/yhdPputative peptidase/DUF3971‑AsmA2 domains protein
Reads supporting (aligned to +/- strand):  ref base C (110/87);  new base G (2/18);  total (112/105)
Fisher's exact test for biased strand distribution p-value = 7.78e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 219 total aligned reads displayed.

CAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTG  >  NC_000913/3392054‑3392350
                                                                                                                                                   |                                                                                                                                                     
cAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACgg                                                                                                                                                     <  2:1425187/150‑1 (MQ=255)
    agagGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTAt                                                                                                                                                 <  2:508684/150‑1 (MQ=255)
     gagGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATg                                                                                                                                                >  1:276674/1‑150 (MQ=255)
        gATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTAGCCAGCCTACGGTTATGTTCCGGTTTTTAGGCCTGATAAGACGCACAGCGCCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTc                                                                                                                                             <  1:1512092/150‑1 (MQ=255)
           tttttACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTGTGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCACCGGCTGTCGGGTGCGGCGTACTCGCCGTTTCCGACCTACGGGTCTGTTcggg                                                                                                                                          >  2:2264910/1‑147 (MQ=17)
              ttACGGCTGCCGGATGCGGCGTGAACGCCCTACCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTg                                                                                                                                       <  2:291902/150‑1 (MQ=255)
                  ggCTGCGGGATGCGGCGTGAACGGCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAgg                                                                                                                                   <  2:1566755/150‑1 (MQ=255)
                  ggCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAgg                                                                                                                                   <  2:2902947/150‑1 (MQ=255)
                  ggCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAgg                                                                                                                                   <  2:2659468/150‑1 (MQ=255)
                   gCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGTTTTTGTTCCGTTTTGTAGGCCTTATAAGACGCACAGCGTCGCAACAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGc                                                                                                                                  <  1:564429/150‑1 (MQ=14)
                   gCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGc                                                                                                                                  <  2:1641680/150‑1 (MQ=255)
                   gCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGc                                                                                                                                  >  1:2851608/1‑150 (MQ=255)
                   gCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTAGGTTCCGGTTGTAGGc                                                                                                                                  >  2:864253/1‑150 (MQ=255)
                   gCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGACGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGc                                                                                                                                  >  1:1842591/1‑150 (MQ=255)
                      gCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAATGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTg                                                                                                                               <  2:1438514/150‑1 (MQ=255)
                      gCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGGTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTCTCCGACCTACGGTTATGTTCCGTTTGTAGGCCTg                                                                                                                               >  2:2006246/1‑150 (MQ=255)
                       ccGGTTGCGGCGTGACCGCCCTATCCAGCCTACGGTTTTGTGCCGGTTTGTAGGCCTTATAAGCCGCACAGCGTCGCATCAGCCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGa                                                                                                                              <  1:829484/150‑1 (MQ=14)
                       ccGGAGGCGGCGTTACCGCCCTATCCAGCCTCCGGTTTTGTTCCGGTTTGTCGGCCTTATAAGACGCACAGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGa                                                                                                                              <  1:532467/150‑1 (MQ=14)
                            tGCGGCGTGAACGCCCTATCCAGCCTACGTTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGa                                                                                                                         <  2:222646/150‑1 (MQ=255)
                             gcggcgTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGAc                                                                                                                        <  1:1982740/150‑1 (MQ=255)
                             gcggcgTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGAc                                                                                                                        <  2:2660625/150‑1 (MQ=255)
                               ggcgTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGc                                                                                                                      >  1:1887023/1‑150 (MQ=255)
                                 cgTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGTGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGcac                                                                                                                    >  2:2810276/1‑150 (MQ=255)
                                 cgTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGcac                                                                                                                    >  2:1040241/1‑150 (MQ=255)
                                  gTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGcaca                                                                                                                   >  2:1086319/1‑150 (MQ=255)
                                    gAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGc                                                                                                                 >  1:901203/1‑150 (MQ=255)
                                     aaCGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCg                                                                                                                <  2:2388430/150‑1 (MQ=255)
                                      aCGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGt                                                                                                               <  2:351147/150‑1 (MQ=255)
                                      aCGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTGATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGt                                                                                                               >  1:741160/1‑150 (MQ=255)
                                        gCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCg                                                                                                             >  2:1765628/1‑150 (MQ=255)
                                          ccTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGACGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCa                                                                                                           >  1:1855457/1‑150 (MQ=255)
                                            tATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATc                                                                                                         >  1:1120405/1‑150 (MQ=255)
                                            tATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCGGATAAGACGCACAGCGGCGCAGCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATc                                                                                                         <  1:255494/150‑1 (MQ=255)
                                             aTCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCa                                                                                                        <  2:215299/150‑1 (MQ=255)
                                              tCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACTGCGTCGCATCAGGCCACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAg                                                                                                       <  1:1009204/150‑1 (MQ=255)
                                              tCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTCCGGGTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAg                                                                                                       >  2:634849/1‑150 (MQ=255)
                                              tCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAg                                                                                                       >  2:1613025/1‑150 (MQ=40)
                                              tCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAg                                                                                                       >  2:1049147/1‑150 (MQ=40)
                                              tCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTGTGTAGGCCTGATAAGACGCACAGCGACGCATCAg                                                                                                       >  1:1477371/1‑150 (MQ=255)
                                              tCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCAGGTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAg                                                                                                       >  2:1497777/1‑150 (MQ=255)
                                                  gCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCAACAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCaa                                                                                                   <  1:72685/150‑1 (MQ=38)
                                                     tACGGTTATGGTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACgg                                                                                                <  2:2164568/150‑1 (MQ=37)
                                                      aCGGTTCTGGGCCGGTTTGTCGGCCTGATACGACGCCCAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGc                                                                                               <  2:2869464/150‑1 (MQ=11)
                                                      aCGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGc                                                                                               <  2:2070137/150‑1 (MQ=37)
                                                              gTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCACCGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTCGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGAt                                                                                       <  2:185180/150‑1 (MQ=255)
                                                               ttCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGAAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATg                                                                                      >  1:2438301/1‑150 (MQ=25)
                                                                tCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgc                                                                                     >  2:2887475/1‑150 (MQ=35)
                                                                tCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgc                                                                                     >  1:2339455/1‑150 (MQ=35)
                                                                 ccGGTTTGTAGGCCTGATAAGCCGCACAGCGTCGCATCAGGCAACGGCGGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgcg                                                                                    <  2:1833644/150‑1 (MQ=21)
                                                                 ccGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGTCAACGGCTGCCGGATgcg                                                                                    >  2:2011727/1‑150 (MQ=21)
                                                                 ccGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgcg                                                                                    >  2:39576/1‑150 (MQ=34)
                                                                 ccGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgcg                                                                                    >  2:2747203/1‑150 (MQ=34)
                                                                 ccGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgcg                                                                                    >  1:2316092/1‑150 (MQ=34)
                                                                  cGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCACCGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgcgg                                                                                   >  2:1440580/1‑150 (MQ=18)
                                                                  cGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATCAGACGCACAGCGTCGCGTCAGGCGACGGCTGCCGGATgcgg                                                                                   >  2:1785361/1‑150 (MQ=16)
                                                                  cGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgcgg                                                                                   >  2:531007/1‑150 (MQ=33)
                                                                    gtgtGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgcggcg                                                                                 >  1:210010/4‑150 (MQ=17)
                                                                     tttGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGt                                                                                >  1:84330/1‑150 (MQ=32)
                                                                        gttggCCTTATAAGACGCACCGCGGCGCCTCTGGAAAGGGCTGTCGGATGCGGCTTAAACGCCTTATCCGACCTTCGTTTTTGTTCCGTTTGTCGGCCTGATAAGACGCACAGCGTCGCATCTGGCAACGGCTGCCGGATGCGGCGTaaa                                                                             <  2:1059323/147‑1 (MQ=11)
                                                                           ggCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGc                                                                          <  2:730287/150‑1 (MQ=34)
                                                                           ggCCTGATAAGACGCACAGCGTCGCATATGGCAACGGCTGGCGGATGCGGCGTAAACGCCTTATCCGACCTACGTTTTTGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGc                                                                          <  2:2821603/150‑1 (MQ=11)
                                                                              cTGATACGACGCACCGCGTCGCCTCCGGCAACGGCTGTCGGATGCGGCGTAACCGCCTTAACCGACCTACGGTTATGTACCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCtt                                                                       <  2:905755/150‑1 (MQ=11)
                                                                              cTGAGACGACGCACCGCGTCGCATCAGGCAACGGCTGGCGGATGCGGCGTACACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCtt                                                                       <  2:2630788/150‑1 (MQ=11)
                                                                               tGATACGCCGCACCGCGTCGAACCCGGCAACGGCTGTGGTATGCGGCGTAACCGCCTTACCCGACCTAGGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTa                                                                      <  2:2027/150‑1 (MQ=11)
                                                                                     gACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGAc                                                                >  1:1359032/1‑150 (MQ=37)
                                                                                        gCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTa                                                             >  1:164596/1‑150 (MQ=37)
                                                                                         cccaGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTAGGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTAc                                                            <  2:2493175/148‑1 (MQ=37)
                                                                                          acaGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCGTATCCGACCTACg                                                           >  2:2584300/1‑150 (MQ=21)
                                                                                                    aTCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCCTATCCGACCTACGGTTATGTTcc                                                 >  1:2014606/1‑150 (MQ=17)
                                                                                                    aTCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTATGTTTGCCTGATAAGACGCACAGCGTCGCGTCAGGCCCCGGCTGCCGGGTGCGGCGTAAACGCCTGATCCGACCTACGGGTATGTTcg                                                 >  2:1661695/1‑149 (MQ=11)
                                                                                                       aGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGtt                                              >  2:469485/1‑150 (MQ=40)
                                                                                                           aaCGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGGTATGTTCCGGTTGTa                                          >  1:369041/1‑150 (MQ=255)
                                                                                                              ggCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCAGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTATGTAGGc                                       >  2:447771/1‑150 (MQ=255)
                                                                                                                gtGTCGGTGGCGGCGTAAACGCCTTATCCGACCTACTGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATGAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCt                                     <  1:2366540/149‑1 (MQ=255)
                                                                                                                cTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCt                                     >  2:2365622/1‑150 (MQ=255)
                                                                                                                 tGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGGTTGTAGGCCTg                                    >  2:2375394/1‑150 (MQ=255)
                                                                                                                      gaggcggcgTAACCGCCTTCGCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTAGGTTCCGTTTGTAGGCCTGATAAg                               <  2:2625372/147‑1 (MQ=14)
                                                                                                                      gATGCGGCTTAAACGCCTTATCCGCCCTCCGTTTATGTTCCGTTTGTCGGCCTGATAAGACGCACCGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAg                               <  2:1010566/150‑1 (MQ=14)
                                                                                                                       ctGCGGCGTAAACGCCTTACCCGACCTAGGTTTATGTCCCGTTTGTTGGCCTGATAAGACGCACCGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGa                              <  1:103198/149‑1 (MQ=14)
                                                                                                                           ggcgTAAACGCCTTATCCGACCTACGTTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGca                          <  1:2597598/150‑1 (MQ=255)
                                                                                                                             cgTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGcaca                        <  2:2594046/150‑1 (MQ=255)
                                                                                                                             cgTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGcaca                        >  1:1379403/1‑150 (MQ=255)
                                                                                                                               taacccccTTAACCGACCTACGGTTAGGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGc                      <  2:1887023/144‑1 (MQ=255)
                                                                                                                                  aCGCGTTAGCCGACCTACTTTTATGTTCGGTTTGTCGGCCTGATAAGACGCACAGCGTCGCATCTGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTc                   <  2:826226/150‑1 (MQ=14)
                                                                                                                                  aCGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTc                   >  1:1726309/1‑150 (MQ=255)
                                                                                                                                  aCGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTc                   >  1:1205769/1‑150 (MQ=255)
                                                                                                                                  aCGCATTATCCGACCTAGGGTTTTTTTCCGTTTTTTGGCCTTATACGCCGCACAGCGTCGCATCTGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGACCTGATAAGACGCACAGCGTc                   <  1:2689035/150‑1 (MQ=14)
                                                                                                                                     ccTTATCCGACCTAGGGTTTTGTTCCGTTTGTAGGCCTGATAAGCCGCACAGCGTCGCATCTGGCCACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCa                <  2:1314829/150‑1 (MQ=14)
                                                                                                                                     ccTTATCCGACCTACGTTTATGTTCCTTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCTGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCa                <  2:2771676/150‑1 (MQ=255)
                                                                                                                                     ccTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCa                >  1:386153/1‑150 (MQ=255)
                                                                                                                                     ccTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCa                >  2:878317/1‑150 (MQ=255)
                                                                                                                                     ccTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCa                >  1:1594934/1‑150 (MQ=255)
                                                                                                                                       ttATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATc              >  2:125078/1‑150 (MQ=255)
                                                                                                                                         aTCCGACCTAGGTTTTTGTTCGGTTTGTAGGCCTGATAAGACGCCCAGCGTCGCATCTGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAg            <  2:597384/150‑1 (MQ=11)
                                                                                                                                         aTCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAg            >  2:424392/1‑150 (MQ=40)
                                                                                                                                         aTCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGACGCATCAg            >  2:2314960/1‑150 (MQ=255)
                                                                                                                                             gACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCaa        >  1:1554610/1‑150 (MQ=38)
                                                                                                                                               ccTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACg      >  2:1146516/1‑150 (MQ=35)
                                                                                                                                                cTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACgg     >  2:2569942/1‑150 (MQ=35)
                                                                                                                                                   cGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTg  >  1:444073/1‑150 (MQ=35)
                                                                                                                                                   |                                                                                                                                                     
CAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTG  >  NC_000913/3392054‑3392350

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