Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,492,834	0	T	G	40.7%	163.1 / 16.6	93	G89G (GGT→GGG)	tsgA	putative transporter
Reads supporting (aligned to +/- strand):  ref base T (12/42);  new base G (37/0);  total (51/42)
Fisher's exact test for biased strand distribution p-value = 2.92e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTA  >  NC_000913/3492687‑3492976                                                                                                                                                     |                                                                                                                                                 gCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGttt                                                                                                                                                <  2:867122/150‑1 (MQ=255)     gTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTTCCGGGTTGATg                                                                                                                                            >  1:170053/1‑150 (MQ=255)     gTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTCGCGGCTTgtt                                                                                                                                             >  4:1220716/1‑147 (MQ=255)       ttCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGGTTGATGtt                                                                                                                                          >  2:349641/1‑150 (MQ=255)        tCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTTCCCGGTTTATGt                                                                                                                                           >  4:444720/1‑148 (MQ=255)         ccAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTTTTGtt                                                                                                                                          >  3:16916/1‑148 (MQ=255)         ccAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGGTTGATGtt                                                                                                                                          >  3:354740/1‑148 (MQ=255)         ccAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTACCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAg                                                                                                                                       <  4:1155976/151‑1 (MQ=255)          cAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGGTTTATGTTTAGc                                                                                                                                      >  4:1662336/1‑151 (MQ=255)            gTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTTTTGTTCAGcc                                                                                                                                     >  4:796027/1‑150 (MQ=255)             tATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGGTTGATGGTTAGCcac                                                                                                                                   >  2:1194745/1‑151 (MQ=255)              aTGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCcaca                                                                                                                                  <  2:284775/151‑1 (MQ=255)              aTGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCcaca                                                                                                                                  <  2:1247989/151‑1 (MQ=255)              aTGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCcaca                                                                                                                                  <  2:905882/151‑1 (MQ=255)               tGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCCGGTTGTTGTTCAGCcac                                                                                                                                   >  1:358942/1‑149 (MQ=255)                gAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAg                                                                                                                                 <  4:1102093/150‑1 (MQ=255)                        ccttcaccttcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCt                                                                                                                        <  2:1174725/151‑1 (MQ=255)                              ccttcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCGCAGCCTGGCGCTGTtctc                                                                                                                  <  2:1496280/151‑1 (MQ=255)                               cttcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTTCCGGTTTGGTGTTTATCCACAACCTGGCGCTGTtttc                                                                                                                  >  2:969763/1‑150 (MQ=255)                               cttcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTTCCGGTTTGATGTTTAGCCACAACCTGGCGCTGtttt                                                                                                                   >  3:44971/1‑147 (MQ=255)                                 tcCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGGTTGATGTTCAGCCCCAGCCTGGCGCTGTtctc                                                                                                                  >  3:177676/1‑148 (MQ=255)                                  ccTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTTCCCGTTTGTTGTTTAGCCCCAAccctgcg                                                                                                                          >  4:1030584/1‑134 (MQ=255)                                  ccTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTc‑ggcg                                                                                                             <  2:711166/151‑1 (MQ=255)                                   cTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTTCCGGGTTTATGTTTAGCCACAGCCTGGCGCTGTTCTt‑ggcg                                                                                                             >  4:108514/1‑150 (MQ=255)                                     cAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑Ggcggcg                                                                                                          <  2:876510/151‑1 (MQ=255)                                     cAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑Ggcggcg                                                                                                          <  2:292886/151‑1 (MQ=255)                                      aaCGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTTCCGGGTTTATGTTCAGCCACAGCCCGGCGCTGTTCTC‑Ggcggcgg                                                                                                         >  4:305327/1‑150 (MQ=255)                                       aCGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGGAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTTGTGTTTAGCCACAGCCTGGCGCTGTTCTc‑ggcgg                                                                                                            >  2:74583/1‑147 (MQ=255)                                       aCGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGAt                                                                                                        <  4:444904/151‑1 (MQ=255)                                       aCGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGGTTGGTGTTTCGCCCCAGCCTGGCGCTGTTCTc‑ggcgg                                                                                                            >  4:76114/1‑147 (MQ=255)                                         gCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGGTGTTCAGCCACACCCCGGCGCTGTTTTT‑Ggcggcgcttt                                                                                                      >  3:69489/1‑147 (MQ=255)                                          ccGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTTCCGGGTTTATGTTTAGCCACAGCCCGGCGCTGTTTTC‑GGCGGCGgtgt                                                                                                      >  2:185996/1‑150 (MQ=255)                                            ggCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGGTTGTTGTTTAGCCCCAGCCTGGCGCTGTTCTC‑GGCGGCGGTGtt                                                                                                     >  4:143943/1‑149 (MQ=255)                                                ttttAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGGTTGATGGTCAGCCCCAACCTGCCGCTGTTCTC‑CGGGGCGATGGTTAtt                                                                                                 >  3:209306/1‑149 (MQ=255)                                                  ttAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATtctc                                                                                              <  3:370386/150‑1 (MQ=255)                                                   tAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGCCTTT‑TCCTGGTGGTGCCGGGGGGTTGCCGGTTTTTTGTTCAGCCCCAACCTGGCGCCGTTTCT‑GGCGGCGGGGTTCttttccgg                                                                                            >  4:1585920/1‑146 (MQ=255)                                                     aTCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGGTTTATGTTTAGCCCCAGCCCGGCGCTGTTCTC‑CGCGGCGATGTTTATTTTTgggg                                                                                          >  3:242853/1‑151 (MQ=255)                                                      tctcTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTTCCGGGTTGATGTTCAGCCACAACCTGGCGCTGTTTTC‑GGCGGCGGTGTTTATTTTCggggg                                                                                         >  4:404242/1‑150 (MQ=255)                                                      tctcTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGGTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGGTGTTCATTTTCggggg                                                                                         >  2:182108/1‑150 (MQ=255)                                                         ctATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTTTCGGGGTGTTGTTTAGGCGCAACCCCGCGCTGTTTTC‑GGAGGGGATGGTGATTTTTGggggggg                                                                                      >  4:1062739/1‑150 (MQ=255)                                                              ttCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTTTC‑GGCGGCGGTGTTCCTTCTCGGGGGGGTCcgcg                                                                                 >  4:406459/1‑151 (MQ=255)                                                                ccTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGGTTTATGTTCAGCCACAGCCTGGCGCTGTTTTC‑GGCGGCGATGTTTATTCTCGGGGGGGGCAGCgg                                                                                >  3:4949/1‑150 (MQ=255)                                                                  tCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTTCCGGGTTGGTGTTTAGCCACACCCTTGCGCCGTTCTT‑TGCCGGGGTGTTTATTCTCGGGGGGGGCCGCGGcgt                                                                             >  3:439137/1‑149 (MQ=255)                                                                   cAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTTCCGGGTTGGTGTTCAGCCCCACCCTCGCGCTGTTTTC‑GGCGGGGGTGGTTTTTTTCGGGGGGGTGCGCggggt                                                                             >  2:1344361/1‑147 (MQ=255)                                                                    aaCGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGGTTTGTGTTTAGCCACAGCCTGGCGCTGTTTTC‑GGCGGCGGTGTTTATTTTCGGGGGGGGGAGCgggctc                                                                            >  4:1612174/1‑146 (MQ=255)                                                                    aaCGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCCGGTTGTTGTTTAGCCACAGCCTGGCGCTGTTCTC‑GGCGGGGGTGTTTTTTCTCGGGGGGGGGAGCGGCCTc                                                                            >  2:1009359/1‑150 (MQ=255)                                                                     aCGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGTCGGGTTTATGTTTCGGCACACCCCCGCGCTGTTCTC‑GGGGGCGATGTTCGCTGTCGGGGGGGGCAACGGCaaa                                                                            >  2:210653/1‑147 (MQ=255)                                                                      tgCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGGTTGTTGTTTAAGCCCAACCTTGCGCTGTTTTC‑GGCGGCGGTGGTCATTCTCGGGGGGGGGAGGGGCATCTcc                                                                         >  2:1153735/2‑151 (MQ=255)                                                                       gCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCa                                                                        <  1:197860/151‑1 (MQ=255)                                                                         cTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGGTTTATGTTTTGCCCCACCCCCGCGCTGTTTTT‑GGCGGCGGTGGTTTTTTTCGGGGGGGGCAGCgggctccacc                                                                        >  2:761819/1‑141 (MQ=255)                                                                          tGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTTCCGGGTTTGTGTTCAGCCCCAACCCGGGGCTGGTCTC‑GGCGGGGGGGGTTTTTTTTGGGGGGGGCAGCGGCATCAccc                                                                        >  2:1070433/1‑147 (MQ=255)                                                                           ggCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTc                                                                    <  4:1395670/151‑1 (MQ=255)                                                                            gCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCg                                                                   <  2:429073/151‑1 (MQ=255)                                                                               gATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGAt                                                                 <  2:1349873/150‑1 (MQ=255)                                                                                  ggAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGGCGGCTTTATTTTTAGACCCAACCTCGCGCTGTTCTCGGGCGGC‑ATGTTTATTCTCGGGGGGGGCAGCGGCGTTACCCTGTCGGTTTGt                                                             >  2:1444186/1‑151 (MQ=255)                                                                                    aaaTCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTAc                                                           <  4:1073036/151‑1 (MQ=255)                                                                                     aaTCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACa                                                          <  4:34753/151‑1 (MQ=255)                                                                                     aaTCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACa                                                          <  4:1298147/151‑1 (MQ=255)                                                                                      aTCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑t                                                        <  3:275097/151‑1 (MQ=255)                                                                                        cGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑tt                                                       <  4:811402/150‑1 (MQ=255)                                                                                             cGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTTCCCGTTTGTTTTTTAGCCCCCACCTCGCGCTGTTTTC‑GGCGGCGATTTTCTTTTTCGGGGGGGTCAGCGGCATCCCCCTTTCGCTTGGTTCA‑TTCCTTGt                                                 >  3:385054/1‑151 (MQ=255)                                                                                             cGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGt                                                 <  4:1490048/151‑1 (MQ=255)                                                                                               ttGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGGTTGATGTTTAGCCCCAGCCCGGCGCTGTTCTC‑GGCGGGGGTGTTTCTTTTCGGGGGGGGCAGCGGCCTCACCCTGTCGGtttggacctttcttg                                                   >  1:308871/1‑135 (MQ=255)                                                                                                tGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGGTTTATGTTCAGCCACAGCCTGGGGCTGGTTTC‑GGGGGGGGTGGTTATTTTCGGGGGGGGGCGGGGGATCACCATGTCGGTTGGGTAA‑TTTCTGGtg                                                >  2:1593243/1‑148 (MQ=255)                                                                                                  aaaaCGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTTCCGGGTTGATGTTCAGCCACAACCCGGCGCTGTTCTC‑GGCGGGGGTGTTTCTTTTCGGGGGGGGCAGCGGCATCACCCTGTCGGTTGGGAAA‑TTTCTGGGaaaa                                             >  4:696087/1‑148 (MQ=255)                                                                                                  aaaaCGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGGTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGGTGGTTCTTCTCGGGGGGGGCAGCGGGATCACCCTGTCGATTGGTACA‑TTCCTGGGaa                                               >  4:602781/1‑148 (MQ=255)                                                                                                      cGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAAt                                        <  4:1010666/151‑1 (MQ=255)                                                                                                      cGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGGTTTATGTTCAGCCACAGCCCGGCGCTGTTTTC‑GGCGGCGGTGTTCATTTTCGGGGGGGGCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACaaa                                         >  4:659363/1‑150 (MQ=255)                                                                                                      cGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCCGGTTGGTGTTTAGCCACAACCTGGCGCTGTTCTC‑GGCGGCGATGTTCCTTTTCGGGGGGGGCAGCGGCATCACCATGTCGATTGGGACA‑TTTCTGGGAAAACaaa                                         >  4:435870/1‑150 (MQ=255)                                                                                                       gCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATg                                       <  1:245680/151‑1 (MQ=255)                                                                                                        cAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGt                                      <  2:334158/151‑1 (MQ=255)                                                                                                         aGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGt                                      <  2:29170/150‑1 (MQ=255)                                                                                                           ttACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGTTTTTTAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGGTGTTCATTTTCGGGGGGGTCAGCGGCCTCACCATGTTGGTTGGTTCA‑ATCCTGGTAACACaaa                                         >  2:117700/1‑145 (MQ=255)                                                                                                              cGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGGTTTCTGTTCCGGCACAACCTCGGGCGGTTTTC‑CGCGGCGGCGTTCGTTTTCGGGGGGGGCCGGGGCCTCACCCTTTCGCGCGGTTCATTTCCTGGttaacaaa                                          >  2:1444213/1‑135 (MQ=255)                                                                                                               gTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCCGTTTTATGTTCAGCCACAGCCTGGCGCTGTTTTC‑GGCGGCGATGTTCATTCTCGGGGGGGGCAGCGGCATCACCCTGTCGGTTGGTACA‑TTCCTGGTAACACAAATGTATGAAgg                               >  4:1166198/1‑151 (MQ=255)                                                                                                                ttttGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGGTTTATGTTTAGCCACAGCCCGGCGCTGTTTTC‑GGCGGGGGTGTTTATTCTCGGGGGGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGGAAAACAAATGTATGAAggg                              >  4:483966/1‑151 (MQ=255)                                                                                                                 tttGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGGTTTATGTTTAGCCACAACCTGGCGCTGTTCTC‑GGCGGGGGTGGTTATTCTCGGGGGGGGGAGCGGCATCACCATGTCGGTTTGGACA‑TTTCTGGTGACACAAAATTATGAAagg                              >  2:1469569/1‑150 (MQ=255)                                                                                                                  ttGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTATGAAGGGc                             <  2:1275484/150‑1 (MQ=255)                                                                                                                   tGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTATGAAGGGCg                            <  2:1661633/150‑1 (MQ=255)                                                                                                                    ggCTTTCTCCTGATGGTG‑CTGGCGGTTTCCGGGTTGTTGGTTAGGCACAACCCCGCGGCGTTCTT‑CGCGGGGGTGTTCATTCTCGGGGGGGGCCGCCGGGTCACCCATACGGTTGGGACA‑ATTCTGGTAACAGAAATTTATGAAgggg                             >  4:1204079/1‑147 (MQ=255)                                                                                                                    ggCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTATGAAGGGCGt                           <  4:1011921/150‑1 (MQ=255)                                                                                                                    ggCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTATGAAGGGCGTc                          <  3:14391/151‑1 (MQ=255)                                                                                                                    ggCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTATGAAGGGCGTc                          <  2:1251349/151‑1 (MQ=255)                                                                                                                        ttCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTATGAAGGGCGTCAgc                       <  2:1610208/150‑1 (MQ=255)                                                                                                                        ttCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTATGAAGGGCGTCAg                        <  2:870440/149‑1 (MQ=255)                                                                                                                         tctcCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTATGAAGGGCGTCAgcg                      <  2:1212360/150‑1 (MQ=255)                                                                                                                         tctcCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTATGAAGGGCGTCAgcg                      <  4:210993/150‑1 (MQ=255)                                                                                                                            ccTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTATGAAGGGCGTCAGcgcg                    <  2:1061383/149‑1 (MQ=255)                                                                                                                               gATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTc                <  4:1374913/150‑1 (MQ=255)                                                                                                                                aTGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTccc              <  4:1595209/151‑1 (MQ=255)                                                                                                                                aTGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTccc              <  4:1006473/151‑1 (MQ=255)                                                                                                                                          gCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTAt     <  1:248494/151‑1 (MQ=255)                                                                                                                                           cGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTAtt    <  4:775590/151‑1 (MQ=255)                                                                                                                                             gTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTa  <  3:308888/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GCCTGTTTCCAGTATGAGTAATACCTTCACCTTCCTCAACGCCGGCATTTTAATCTCTATCTTCCTCAACGCCTGGCTGATGGAAATCGTCCCGTTGAAAACGCAGTTACGTTTTGGCTTTCTCCTGATGGTG‑CTGGCGGTTGCCGGTTTGATGTTCAGCCACAGCCTGGCGCTGTTCTC‑GGCGGCGATGTTCATTCTCGGGGTGGTCAGCGGCATCACCATGTCGATTGGTACA‑TTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTA  >  NC_000913/3492687‑3492976

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑