Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,392,212 C→G 5.0% intergenic (‑184/+246) 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,2120CG 5.0% 1047.0 / 2.5 361intergenic (‑184/+246)tldD/yhdPputative peptidase/DUF3971‑AsmA2 domains protein
Reads supporting (aligned to +/- strand):  ref base C (175/154);  new base G (4/21);  total (186/175)
Fisher's exact test for biased strand distribution p-value = 3.00e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 365 total aligned reads displayed.

TTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGC  >  NC_000913/3392064‑3392361
                                                                                                                                                    |                                                                                                                                                     
ttttttACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCg                                                                                                                                                      >  1:1030506/1‑150 (MQ=255)
  ttttACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGGTTGTAGGCCTGAGAAGACGCACAGCGTCGCAGCAGGCACCGGCTGACCGCTGGGGCGTAAACGCAGAGTCCGACCTACGGGGGTGTTCCggg                                                                                                                                                    >  2:3207213/1‑148 (MQ=17)
   tttACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTAGTAGGCCTGAAAATACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCAACGGTTATGTTCCGgtt                                                                                                                                                   >  1:4911523/1‑150 (MQ=17)
      aCGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTa                                                                                                                                                <  2:3614048/150‑1 (MQ=255)
        ggCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAgg                                                                                                                                              <  2:1806396/150‑1 (MQ=255)
         gcggCCGGAGGGGGCGTGAACGCCCTATCCAGCCTAGGGTTTTGTTCCGTTTTTTAGGCCTTATAAGACCCACCGCGTCGCAACAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGc                                                                                                                                             <  1:5061333/147‑1 (MQ=14)
         gCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTAGGTTCCGTGTGTAGGc                                                                                                                                             >  1:1329777/1‑150 (MQ=255)
          cTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGcc                                                                                                                                            >  1:428120/1‑150 (MQ=255)
           tGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTTATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCt                                                                                                                                           <  1:3355174/150‑1 (MQ=255)
           tGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCt                                                                                                                                           <  2:4002263/150‑1 (MQ=255)
           tGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGGTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGACGGATGAGGCGAAAACGCCTTAACCGACCTAAGCTTAGGTTCCGATTGTAGGCCt                                                                                                                                           >  1:267485/1‑150 (MQ=14)
             ccGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCAGCAGGCAACGGCTGTCGGATGCGGCGTACACGCCGTAGCCGACCTACGGTTATGGTCGGTGTGGAGGCCTgc                                                                                                                                         >  1:4937191/1‑149 (MQ=14)
              cGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGAt                                                                                                                                        <  2:4443512/150‑1 (MQ=255)
                gATGCGGCGTGAATGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAATGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATaa                                                                                                                                      <  1:5489937/150‑1 (MQ=255)
                 aTGCGGCGTGACCGCCCTATCCAGCCTCCGGTTTTGTTCCGTTTTGTAGGCCTGATAAGACGCACCGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAg                                                                                                                                     <  1:2294590/150‑1 (MQ=14)
                     ggcgTGACCGCCCTATCCAGCCTACGGTTATGTTCCGTTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGc                                                                                                                                 <  1:2178440/150‑1 (MQ=255)
                      gcgTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGca                                                                                                                                <  2:5481280/150‑1 (MQ=255)
                      gcgTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGACGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGca                                                                                                                                >  2:3470292/1‑150 (MQ=255)
                      gcgTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGCCAACGGCTGTCGGATGCGGCGTAAACGCCATATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGca                                                                                                                                >  1:4182213/1‑150 (MQ=255)
                        gTGAACGCCCTATCCAGCCTACGTTTTTGTTCCGGTTTGTAGGCCTTATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGcaca                                                                                                                              <  1:892469/150‑1 (MQ=255)
                        gTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTCCGGTTATGTTCAGTTTGTAGGCCTGATAAGACGcaca                                                                                                                              >  2:5244706/1‑150 (MQ=255)
                        gTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTGTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGCCGTAAACGCCATATCCGACCTACGGTTATGATCCGTTTGTAGGCCTGATAAGACGcaca                                                                                                                              >  1:1900371/1‑150 (MQ=255)
                         ggAACGCCCTATCCAGCCTACGGTTCTGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAg                                                                                                                             >  1:4771829/2‑150 (MQ=255)
                           aaCGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCg                                                                                                                           >  2:1373444/1‑150 (MQ=255)
                            aCGCCCTAGCCTGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGGTCCGTTTGTAGGCCTGATAAGACGCACAGCGt                                                                                                                          <  2:1547868/150‑1 (MQ=255)
                             cGCCCTATCCAGCCTACGGATATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGGCGGATGCGGCGTAAACGCCGTATCCGACCTACGGTTATGTTCAGTTTCTAGGCCTGTTAAGACCCACAGCGTc                                                                                                                         >  2:3892122/1‑150 (MQ=14)
                             cGCCCTATCCAGCCACCGGTTTTGTTCCGGTTTGTCGGCCTGATAAGACGCCCCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTc                                                                                                                         <  2:4605512/150‑1 (MQ=14)
                               cccTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCAGTTTGTAGGCCTGATAAGACGCACAGCGTCGc                                                                                                                       >  2:4646127/1‑150 (MQ=255)
                               cccTAGCCAGCCTAGGGTTATGTTCGGGTTTGTCGGCCTGATAAGACGCACAGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGc                                                                                                                       <  2:2092156/150‑1 (MQ=14)
                                ccTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCa                                                                                                                      >  1:3221211/1‑150 (MQ=255)
                                  tacccAGCCTACGGTTAGGTTCCGGTTTGTCGGCCTGATATGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATc                                                                                                                    <  1:1447340/147‑1 (MQ=255)
                                  tATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATc                                                                                                                    >  1:3957944/1‑150 (MQ=255)
                                  tATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATc                                                                                                                    <  2:3766349/150‑1 (MQ=255)
                                  tATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATc                                                                                                                    <  1:1707346/150‑1 (MQ=255)
                                  tATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATc                                                                                                                    >  2:1190908/1‑150 (MQ=255)
                                   aTCCAGCCTAGGGTTTTGTTCCGGTTTGTTGGCCTTATAAGACGCCCAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCa                                                                                                                   <  1:2416109/150‑1 (MQ=14)
                                         ccTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAAc                                                                                                             <  1:308963/150‑1 (MQ=37)
                                          cTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACg                                                                                                            >  1:2776247/1‑150 (MQ=37)
                                           tACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACgg                                                                                                           >  2:479333/1‑150 (MQ=37)
                                            aCGGTTATGTTCCGGTTTGTAGGCCTGATTAGTCGCTCAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTTAACGCCTTATACGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGc                                                                                                          <  2:2234171/150‑1 (MQ=11)
                                            aCGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGc                                                                                                          >  1:4485137/1‑150 (MQ=37)
                                             cgtttATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCt                                                                                                         <  2:588913/147‑1 (MQ=37)
                                               gTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGc                                                                                                       <  2:4945607/150‑1 (MQ=37)
                                                ttATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGcc                                                                                                      >  1:536261/1‑150 (MQ=37)
                                                 tATGTTCCGGTTTGTAGGCCTTATATGACGCACAGCGTCGCTTCAGGCAACGGCTGTGGGATGCGGCTTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCg                                                                                                     <  2:1978534/150‑1 (MQ=11)
                                                  aTGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCgg                                                                                                    <  2:450485/150‑1 (MQ=37)
                                                     ttCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTGCGTTTGTTGGCGTGTTAATACTCACAGCGTCGAATCAGGCAGCGGCTGCCGGATg                                                                                                 >  1:3491156/1‑150 (MQ=11)
                                                       ccGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTTGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgcg                                                                                               >  2:4136654/1‑150 (MQ=21)
                                                       ccGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgcg                                                                                               >  1:3820014/1‑150 (MQ=34)
                                                        cGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCaa                                                                                                              <  2:2895406/134‑1 (MQ=33)
                                                        cGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgcgg                                                                                              >  1:2603390/1‑150 (MQ=33)
                                                        cGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgcgg                                                                                              >  1:2889459/1‑150 (MQ=33)
                                                         ggTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgcggc                                                                                             >  1:1053498/1‑150 (MQ=32)
                                                         ggTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGAGAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgcggc                                                                                             >  1:1289825/1‑150 (MQ=17)
                                                           tttGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGt                                                                                           >  2:2167485/1‑150 (MQ=32)
                                                               tAGGCCTGATAAGACGCACAGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTAGGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAAc                                                                                       <  2:5350667/150‑1 (MQ=21)
                                                                 ggCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACGGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGc                                                                                     >  2:3094368/1‑150 (MQ=21)
                                                                 ggCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGc                                                                                     >  1:3311481/1‑150 (MQ=34)
                                                                   ccTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCt                                                                                   >  2:4701405/1‑150 (MQ=35)
                                                                       aTAAGACGCACCGCGTCGCATCTGGCATCGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTAGGGTTTTGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATc                                                                               <  2:5372408/150‑1 (MQ=11)
                                                                          cgccGCACTGCGTCGCTTATGTCAACGGCGGTCGGAGGCGGGGTAAACGCCTCATCCGACCTACGTTTATGTTCCGTTTGTCGGCCTGATAAGACGCACTGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGa                                                                            <  2:2508499/147‑1 (MQ=11)
                                                                          aGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGa                                                                            >  1:525162/1‑150 (MQ=37)
                                                                              gCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTa                                                                        >  1:608172/1‑150 (MQ=37)
                                                                                acaGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACg                                                                      >  2:4977619/1‑150 (MQ=37)
                                                                                 caGCGTCGCACCTGGCTTCGGGTGTCGGTTGCGGGGTAAAGGCGTTAAGCGCCCTACGGTTATGTTCCGTTTGTTGTCCTGTTTAGACGCACAGCGTCGCTTAATGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACgg                                                                     <  2:2768336/150‑1 (MQ=11)
                                                                                    cGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGGTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTa                                                                  >  1:807763/1‑150 (MQ=37)
                                                                                     gTCGCATCTGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTAt                                                                 <  2:3315570/150‑1 (MQ=255)
                                                                                         cATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTATGTTc                                                             >  1:3876930/1‑150 (MQ=17)
                                                                                          aTCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCAACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCATATCCGACCTACGGTTATGttac                                                            >  1:1639691/1‑148 (MQ=19)
                                                                                             aGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGtt                                                         >  2:2450018/1‑150 (MQ=40)
                                                                                             aGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGtt                                                         >  1:2700732/1‑150 (MQ=40)
                                                                                              gcaaaCGGGTGTCTGTTGCGGCGTAAACGCCTTATCCGACCTACGGTTTTGTTCCGTTTGTTGGCCTGATAAGACGCACTGCGTCGCATCTGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGttt                                                        <  1:3472937/147‑1 (MQ=14)
                                                                                                  aCGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAg                                                    <  2:3221211/150‑1 (MQ=255)
                                                                                                     gCTGTCGTATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGcc                                                 <  2:1224563/150‑1 (MQ=255)
                                                                                                      cTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGTTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCt                                                <  2:1326615/150‑1 (MQ=255)
                                                                                                             ttGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTCGGCCTGATAAGACGAACCGCGTCGCATCTGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGa                                         <  1:3588005/149‑1 (MQ=255)
                                                                                                               gcggcgTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACg                                       >  1:1548245/1‑150 (MQ=255)
                                                                                                                cggcgTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCAACCTACGGTTATGTTCAGTATCTATGCCAGATAAGAccc                                      >  2:125627/1‑148 (MQ=14)
                                                                                                                 ggcgTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGca                                     >  2:4363729/1‑150 (MQ=255)
                                                                                                                  gcgTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGcac                                    >  1:4583616/1‑150 (MQ=255)
                                                                                                                  gcgTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTGATGTTCCGTGTGTAGGCCTGATAAGACGcac                                    >  2:2849993/1‑150 (MQ=255)
                                                                                                                   cgTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGcaca                                   >  1:3685608/1‑150 (MQ=255)
                                                                                                                   cgTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGcaca                                   >  2:1404428/1‑150 (MQ=255)
                                                                                                                    ttAAACGCATTACCCGCCCTAGGTTTTTTTTCGTTTTGTCGGCCTGATAAGACGCACCGCGTCGCATCAGGCCACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAg                                  <  1:517686/149‑1 (MQ=11)
                                                                                                                        cggCCTTAGCCGACCTACGGTTTGGTGCGGTTTGTCGGCCTGTTTAGCCGCACGGCGCCGCTTCCGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTc                              <  1:5106081/148‑1 (MQ=11)
                                                                                                                        aCGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTc                              >  1:314960/1‑150 (MQ=255)
                                                                                                                         cGCCTTCTCCGCCCTCGGGTTCGGTTCGGTTGGTCGGCCTGATAAGACGCCCAGCGTCGCCTCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCg                             <  2:5437045/150‑1 (MQ=11)
                                                                                                                           ccTTATCCGACATACGTTTATGTTCGTTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCa                           <  2:419461/150‑1 (MQ=255)
                                                                                                                            cTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCAt                          >  2:840405/1‑150 (MQ=255)
                                                                                                                             ttATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCTCATc                         >  2:2319054/1‑150 (MQ=255)
                                                                                                                               aTCCGACCTCCGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAg                       <  2:357921/150‑1 (MQ=255)
                                                                                                                               aTCCGACCTACTTTTATGTTCCGTTTGTCGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAg                       <  1:2075172/150‑1 (MQ=255)
                                                                                                                                tCCGACCTCCGTTTATGTTCCGTTTGTAGGCCTGATATGACGCCCAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAgg                      <  1:2875631/150‑1 (MQ=255)
                                                                                                                                      ctccgtttATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACgg                <  2:3371559/144‑1 (MQ=25)
                                                                                                                                           gTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGcc           >  2:2129050/1‑150 (MQ=35)
                                                                                                                                             tATGTTCCGTTTGTAGACCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGGAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCgg         >  2:4286088/1‑150 (MQ=25)
                                                                                                                                               tGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGAt       >  1:5351316/1‑150 (MQ=35)
                                                                                                                                               tGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCATATCCGACCTACGGTTATGTTCAGTGTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCACCGGCTGCCGGAt       >  1:711821/1‑150 (MQ=255)
                                                                                                                                                 ttCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgc     >  2:1017296/1‑150 (MQ=35)
                                                                                                                                                 ttCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgc     >  1:2927995/1‑150 (MQ=35)
                                                                                                                                                    |                                                                                                                                                     
TTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGC  >  NC_000913/3392064‑3392361

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: