Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,392,273 A→T 7.6% intergenic (‑245/+185) 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,2730AT 7.6% 244.7 / 7.0 87intergenic (‑245/+185)tldD/yhdPputative peptidase/DUF3971‑AsmA2 domains protein
Reads supporting (aligned to +/- strand):  ref base A (37/43);  new base T (4/3);  total (41/46)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGAT  >  NC_000913/3392126‑3392418
                                                                                                                                                   |                                                                                                                                                 
gTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTa                                                                                                                                                   >  3:369722/1‑148 (MQ=32)
 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: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                                                                                                                           <  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:21452/149‑1 (MQ=255)
                                                           cGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTc                                                                                       <  8:14586/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                                                                                                                                     <  8:359725/94‑1 (MQ=37)
                                                                    cGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGTAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGa                                                                                                                                     >  7:359725/1‑94 (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                                                                      >  5:7349/1‑149 (MQ=35)
                                                                            ggTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTg                                                                      >  8:278238/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                                                                 <  4:409946/128‑1 (MQ=32)
                                                                                                      aaGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGa                                                                 >  3:409946/1‑128 (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                              >  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                   <  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)
                                                                                                                                                   |                                                                                                                                                 
GTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGAT  >  NC_000913/3392126‑3392418

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: