Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,392,274 A→T 6.1% intergenic (‑246/+184) 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,2740AT 6.1% 260.5 / 4.1 87intergenic (‑246/+184)tldD/yhdPputative peptidase/DUF3971‑AsmA2 domains protein
Reads supporting (aligned to +/- strand):  ref base A (38/43);  new base T (3/3);  total (41/46)
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

TAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATC  >  NC_000913/3392127‑3392422
                                                                                                                                                   |                                                                                                                                                    
gttgCCTGATGCGACGCTGTGCGTCTTATCGGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGGTATGTTCCGTTTGTAGGCCTGATAAGACGCGCAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGGaac                                                                                                                                                     >  3:277484/4‑148 (MQ=11)
    ccTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGcc                                                                                                                                                 <  1:410796/149‑1 (MQ=34)
     cTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCATAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCt                                                                                                                                                >  1:222554/1‑149 (MQ=25)
        aTAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCTCATCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTAt                                                                                                                                             <  6:213386/149‑1 (MQ=255)
          aaGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATc                                                                                                                                            <  2:234336/148‑1 (MQ=37)
            gACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCGTATCCGACCTACGGATATGTTCCGTGGGTAGGCCTGCTAAGACGCACAGCGTCGCATCATGCACCGGCTGCCGGCTGCGGCGTAAACACCAAATCCGa                                                                                                                                         >  8:366612/1‑149 (MQ=11)
                   aGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACgg                                                                                                                                  <  8:208825/149‑1 (MQ=37)
                   aGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACgg                                                                                                                                  <  7:41852/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:14586/1‑149 (MQ=37)
                           aTCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTc                                                                                                                          >  7:21452/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                                                                                                                                       >  6:350615/1‑126 (MQ=11)
                                  aaCGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATAAGGCGTTTACGCCGCATCCggca                                                                                                                                       <  5:350615/129‑4 (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                                                                                                                               <  8:199258/125‑1 (MQ=35)
                                              aTGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTa                                                                                                                               >  7:199258/1‑125 (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                                                                                                                                    >  7:43268/1‑96 (MQ=37)
                                                                      ccTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTAc                                                                                                                                    <  8:43268/96‑1 (MQ=37)
                                                                      ccTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGtt                                                                                                                                <  7:418774/100‑1 (MQ=37)
                                                                      ccTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGtt                                                                                                                                >  8:418774/1‑100 (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                                                                                                          <  1:336239/98‑1 (MQ=11)
                                                                                          cgttgCCTGATGCGACGCTGTGCGTCTTATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGAt                                                                                                          >  2:336239/5‑102 (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)
                                                                                                            acaGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGc                                                                                                                                                  <  2:102288/44‑1 (MQ=38)
                                                                                                             caGCGTCGCTTCAGGCATCGGCTGGCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACg                                        <  7:366612/149‑1 (MQ=255)
                                                                                                                cGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGtt                                     <  3:370758/149‑1 (MQ=37)
                                                                                                                  tCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTAt                                   <  5:385425/149‑1 (MQ=37)
                                                                                                                     cATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATg                                  <  6:192444/147‑1 (MQ=37)
                                                                                                                     cATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATg                                  >  5:192444/1‑147 (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                       <  7:79025/134‑1 (MQ=14)
                                                                                                                                        gataaggcgTTTACGCCGCATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGt                       >  8:79025/6‑139 (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  <  4:368265/149‑1 (MQ=255)
                                                                                                                                                   aaCGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATc  <  1:394798/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATC  >  NC_000913/3392127‑3392422

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: