Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,392,281 T→C 6.7% intergenic (‑253/+177) 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,2810TC 6.7% 263.1 / 6.6 91intergenic (‑253/+177)tldD/yhdPputative peptidase/DUF3971‑AsmA2 domains protein
Reads supporting (aligned to +/- strand):  ref base T (37/46);  new base C (3/3);  total (42/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTA  >  NC_000913/3392135‑3392425
                                                                                                                                                  |                                                                                                                                                
aTAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCTCATCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTAt                                                                                                                                                <  6:213386/149‑1 (MQ=255)
  aaGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATc                                                                                                                                               <  2:234336/148‑1 (MQ=37)
    gACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCGTATCCGACCTACGGATATGTTCCGTGGGTAGGCCTGCTAAGACGCACAGCGTCGCATCATGCACCGGCTGCCGGCTGCGGCGTAAACACCAAATCCGa                                                                                                                                            >  8:366612/1‑149 (MQ=11)
           aGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACgg                                                                                                                                     <  7:41852/149‑1 (MQ=37)
           aGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACgg                                                                                                                                     <  8:208825/149‑1 (MQ=37)
             cGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGtt                                                                                                                                   <  7:444762/149‑1 (MQ=37)
              gTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTa                                                                                                                                  <  5:324310/149‑1 (MQ=37)
                  cATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTCTAGGCCTGATAAGACGCACAGCGTCGCATGAGGCAACGGCTGCCGGATGAGGCGTTAACGCCTTATCCGACCTACGGTTAAGtt                                                                                                                              >  5:191616/1‑149 (MQ=11)
                   aTCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAATACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTAt                                                                                                                                                >  6:314813/1‑130 (MQ=38)
                   aTCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTc                                                                                                                             >  7:21452/1‑149 (MQ=37)
                   aTCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTc                                                                                                                             >  7:14586/1‑149 (MQ=37)
                   aTCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGGTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTAt                                                                                                                                                <  5:314813/130‑1 (MQ=38)
                       ggCAACGGCTGTCGGATGCGGCTTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGt                                                                                                                          <  7:195236/148‑1 (MQ=255)
                        gCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCAGATAACACCCACAACCACGCCTCAGGCAACGCCTGCCGGATGCGGAGAAACCGCCTAATCCGACGAACGGATATCTTCCgagt                                                                                                                        >  6:139485/1‑146 (MQ=11)
                          aaCGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATAAGGCGTTTACGCCGCATCCggca                                                                                                                                          <  5:350615/129‑4 (MQ=11)
                          aaCGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATAAGGCGTTTACGCCGCATCCggca                                                                                                                                          >  6:350615/1‑126 (MQ=11)
                          aaCGGCTGTCGGATGCGGAGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGt                                                                                                                      <  2:10613/149‑1 (MQ=255)
                             ggCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATCAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATAAGGCGTTTACGCCGCATCCggca                                                                                                                                          <  4:170239/126‑4 (MQ=11)
                             ggCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATAAGGCGTTTACGCCGCATCCggca                                                                                                                                          >  3:170239/1‑123 (MQ=11)
                             ggCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCAGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAgg                                                                                                                   >  8:205676/1‑149 (MQ=255)
                               cTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCATAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGcc                                                                                                                 <  2:222554/149‑1 (MQ=255)
                                  tCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGa                                                                                                              >  4:370758/1‑149 (MQ=255)
                                      aTGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTa                                                                                                                                  >  7:199258/1‑125 (MQ=35)
                                      aTGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTa                                                                                                                                  <  8:199258/125‑1 (MQ=35)
                                         cggctTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCTTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGAc                                                                                                        <  8:47267/148‑1 (MQ=255)
                                           gcgTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTCTAGGCCAGATAAGACGCACAGCGACGCATCAGGCAACGGCTGCCGCATGCCGCGTAAACGCCCTATCCGACCTACGGTTATCTACAGTTAGTAGGCCTGATAACACGca                                                                                                     >  5:91503/1‑149 (MQ=14)
                                             gTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACATCGTCGGATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGcaca                                                                                                   <  6:224486/149‑1 (MQ=255)
                                                 aCGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCAGACCTACGGTTATGTTCGGTTTGTAGGCCTGATAAGACGCACAGCGt                                                                                               >  7:222964/1‑149 (MQ=255)
                                                  cGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTc                                                                                              <  8:14586/149‑1 (MQ=255)
                                                  cGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTc                                                                                              <  8:21452/149‑1 (MQ=255)
                                                    ccTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCAGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGc                                                                                            <  7:205676/149‑1 (MQ=255)
                                                     cTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCATATCCGACCTACGGTTATGTTCAGTTTGTATGCCAGATAACACCCACAGCGTctc                                                                                            >  1:203362/1‑146 (MQ=14)
                                                        ttCCGACGTACGCTAATGTTCCGTTTGTTGGCTTTATAAGACGCAGAGCTTCTCTTCAGGCAACGGCTTGCGGATGCGGCGTAAACGCCTTTTCCGTCCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATc                                                                                         <  8:72276/147‑1 (MQ=11)
                                                        aTCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCATGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACACACAGCGTCGCATc                                                                                         <  2:2053/148‑1 (MQ=255)
                                                         tCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATACGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAg                                                                                       <  4:277484/149‑1 (MQ=255)
                                                           cGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGTAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGa                                                                                                                                            >  7:359725/1‑94 (MQ=37)
                                                           cGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGTAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGa                                                                                                                                            <  8:359725/94‑1 (MQ=37)
                                                              ccTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTAc                                                                                                                                       <  8:43268/96‑1 (MQ=37)
                                                              ccTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTAc                                                                                                                                       >  7:43268/1‑96 (MQ=37)
                                                              ccTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGtt                                                                                                                                   >  8:418774/1‑100 (MQ=37)
                                                              ccTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGtt                                                                                                                                   <  7:418774/100‑1 (MQ=37)
                                                                 aCGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGc                                                                               >  4:76156/1‑149 (MQ=35)
                                                                   ggTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGc                                                                                                                  >  4:75054/1‑112 (MQ=38)
                                                                   ggTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGc                                                                                                                  <  3:75054/112‑1 (MQ=38)
                                                                   ggTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTg                                                                             >  8:278238/1‑149 (MQ=35)
                                                                   ggTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTg                                                                             >  5:7349/1‑149 (MQ=35)
                                                                       aTGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCg                                                                          >  8:147340/1‑148 (MQ=35)
                                                                          ttCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGAt                                                                       <  6:454369/148‑1 (MQ=35)
                                                                           tCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCgg                                                                         >  5:409465/1‑145 (MQ=34)
                                                                           tCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCgg                                                                         <  6:409465/145‑1 (MQ=34)
                                                                             cGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgcgg                                                                   <  4:369722/149‑1 (MQ=32)
                                                                                  gTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTg                                                               >  5:100097/1‑148 (MQ=32)
                                                                                  cgttgCCTGATGCGACGCTGTGCGTCTTATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGAt                                                                                                             >  2:336239/5‑102 (MQ=11)
                                                                                  cgttgCCTGATGCGACGCTGTGCGTCTTATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGAt                                                                                                             <  1:336239/98‑1 (MQ=11)
                                                                                     ggCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACg                                                           >  3:115651/1‑149 (MQ=34)
                                                                                      gCCTGATAAGACGCACAGCGTCGCATCAGGAAACGGCTGCCGGATGCGGCGTAAACGCCATAAACGACCTACGGTTATGTTCCGTATGTAGGCCTCAAAACACGCACAGCGTCTCATCAGGAAACGGCTGCCGGAAGAGGCCTGAActc                                                          >  6:385425/1‑147 (MQ=11)
                                                                                       ccTGATAAGACGCACAGAGTCGAATCAGGCAACGGCTTCCGGATGCGTCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGcc                                                         <  6:224239/149‑1 (MQ=255)
                                                                                        cTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGccc                                                        >  7:408368/1‑149 (MQ=35)
                                                                                         tGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGccc                                                        >  3:96952/1‑148 (MQ=35)
                                                                                            tAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATc                                                    >  4:151265/1‑149 (MQ=37)
                                                                                             aaGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGa                                                                        >  3:409946/1‑128 (MQ=32)
                                                                                             aaGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGa                                                                        <  4:409946/128‑1 (MQ=32)
                                                                                               gACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGa                                                 <  6:100097/149‑1 (MQ=37)
                                                                                                  gCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGAcc                                               >  4:417842/1‑148 (MQ=37)
                                                                                                   cacaGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTa                                             <  6:7349/149‑1 (MQ=37)
                                                                                                     caGCGTCGCTTCAGGCATCGGCTGGCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACg                                           <  7:366612/149‑1 (MQ=255)
                                                                                                        cGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGtt                                        <  3:370758/149‑1 (MQ=37)
                                                                                                          tCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTAt                                      <  5:385425/149‑1 (MQ=37)
                                                                                                             cATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATg                                     >  5:192444/1‑147 (MQ=37)
                                                                                                             cATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATg                                     <  6:192444/147‑1 (MQ=37)
                                                                                                              aTCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACAGTTATGTTc                                  <  8:408368/149‑1 (MQ=37)
                                                                                                                  ggCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGc                              <  4:96952/149‑1 (MQ=40)
                                                                                                                     aaCGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGCTGCGGCGTGAACGCCCGATCCGACCTACGGTTAGGTTCTGGCTCg                           >  2:304973/1‑149 (MQ=255)
                                                                                                                       cGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACTGCGTCGCATCATGCAACGTCTGCCGGATGCGGCGTTAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGtt                         <  8:222964/149‑1 (MQ=255)
                                                                                                                        ggCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTg                        <  3:76156/149‑1 (MQ=255)
                                                                                                                          cTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGc                      >  3:368265/1‑149 (MQ=255)
                                                                                                                             acGGATGCGGCGTAAACGCCTTCTCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGcc                     <  6:191616/146‑1 (MQ=255)
                                                                                                                               ggATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTg                 <  3:151265/149‑1 (MQ=255)
                                                                                                                                gataaggcgTTTACGCCGCATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGt                          >  8:79025/6‑139 (MQ=14)
                                                                                                                                gataaggcgTTTACGCCGCATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGt                          <  7:79025/134‑1 (MQ=14)
                                                                                                                                     ggcgTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAAc            >  2:394798/1‑148 (MQ=255)
                                                                                                                                     ggcgTAAACGCCTTATCCGACCTACGGTTATGTTCCGTATGTAGGCCTGATAAGACGCACAGCGTCGCATCATGCAACGGCTGCCGGATGCTCCTAGAACGCCCTATCCGACCTAGGGTGATGTTCTGGCTCGTTGGCCTTTGGCa              >  6:162792/1‑146 (MQ=17)
                                                                                                                                      gcgTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGa          >  4:137515/1‑149 (MQ=255)
                                                                                                                                       cgTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGAt         >  1:145512/1‑149 (MQ=255)
                                                                                                                                           aaCGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGcc                                                         <  2:219513/97‑1 (MQ=33)
                                                                                                                                           aaCGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATc     <  1:394798/149‑1 (MQ=255)
                                                                                                                                           aaCGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATc     <  4:368265/149‑1 (MQ=255)
                                                                                                                                            aCGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGa          <  4:115246/143‑1 (MQ=255)
                                                                                                                                            aCGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGa          >  3:115246/1‑143 (MQ=255)
                                                                                                                                              gCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCTGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTa  <  2:145512/149‑1 (MQ=255)
                                                                                                                                              gCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTa  >  6:359380/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
ATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTA  >  NC_000913/3392135‑3392425

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: