Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,392,280 T→G 7.2% intergenic (‑252/+178) 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,2800TG 7.2% 242.5 / 8.4 91intergenic (‑252/+178)tldD/yhdPputative peptidase/DUF3971‑AsmA2 domains protein
Reads supporting (aligned to +/- strand):  ref base T (35/46);  new base G (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 = 9.66e-01

CTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTA  >  NC_000913/3392132‑3392425
                                                                                                                                                    |                                                                                                                                                 
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                                                                                                                                     <  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)
                                                                                                                                                    |                                                                                                                                                 
CTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTA  >  NC_000913/3392132‑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: