Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	995,498	0	A	G	32.2%	95.1 / 11.8	59	F164S (TTC→TCC)	ssuD	FMNH2‑dependent alkanesulfonate monooxygenase
Reads supporting (aligned to +/- strand):  ref base A (11/29);  new base G (19/0);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 5.34e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCGCTCGGCGGCCTGCCACGCTTCGTCGTTAGTTTCACGAACAATCACATGCAGACGAATACCGAAACGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACATGAATATGTTTACCGTT‑GAAATCG‑AC‑‑GGTTTCTCTCT‑GCAATAAACGCCGCCAGACCTGGGTAAATTCCGCCGAGGCTTCGTAGCGCTCGCTATGATCAAGGAACACTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGG  >  NC_000913/995254‑995736                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     ccGCTCGGCGGCCTGCCACGCTTCGTCGTTAGTTTCACGAACAATCACATGCAGACGAATACCGAAACGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGGTG‑AATCGCCGGG‑Aagagg                                                                                                                                                                                                                                                 >  1:187067/1‑249 (MQ=255)      cggcggCCTGCCACGCTTCGTCGTTAGTTTCACGAACAATCACATGCAGACGAATACCGAAACGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGGGGCGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGtttt                                                                                                                                                                                                                                          <  2:260968/251‑1 (MQ=255)         cggcCTGCCACGCTTCGTCGTTAGTTTCACGAACAATCACATGCAGACGAATCCCGAAACGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGCTTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAACGTGGCGGATACGGCTGTTG‑GATCGCCGGG‑GAAAGG‑AGtttttct                                                                                                                                                                                                                                       >  1:416730/1‑248 (MQ=255)              tGCCACGCTTCGTCGTTAGTTTCACGAACAATCACATGCAGACGAATACCGAAACGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGCTACGGGTGTTG‑AATCTCCGGG‑GAGAGC‑AGGTTTGCTCcgcg                                                                                                                                                                                                                                  >  1:386157/1‑251 (MQ=255)                  aagCTTCGTCGTTAGTTTCACGAACAATCCCATGCAGACGAATACCGAACCGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGCCATCTGACGCTCCGCCAAAGTAAAGTGGCGGGTACGGCTGTTG‑ATTCACCGGG‑GAGAGCCAGTTTTgtccatcccc                                                                                                                                                                                                                                >  2:303414/3‑241 (MQ=255)                   cGCTTCGTCGTTAGTTTCACGAACAATCACATGCAGACGAATACCGAAACGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACAt                                                                                                                                                                                                                             <  1:130412/251‑1 (MQ=255)                    gcgtcgtcgTTAGTTTCACGAACAATCACATGCAGACGAATACCGAACCGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACATg                                                                                                                                                                                                                            <  1:221240/248‑1 (MQ=255)                    gCTTCGTCGTTAGTTTCACGAACAATCACATGCAGACGAATACCGAAACGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCCGG‑GAGAGC‑GGGTTTGCTCCGCGCCCATg                                                                                                                                                                                                                            >  1:28432/1‑251 (MQ=255)                        cgtcgtTAGTTTCACGAACAATCACATGCAGACGAATACCGAAACGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACATGAata                                                                                                                                                                                                                        <  1:113820/251‑1 (MQ=255)                         gtcgtTAGTTTCACGAACAATCACATGCAGACGAATACCGAAACGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACATGAatat                                                                                                                                                                                                                       <  1:282060/251‑1 (MQ=255)                         ggcgtTAGTTTCACGAACAATCACATGCAGACGGATACCGAAACGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCTCTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCATCTCCTGTGCTACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACATGAatat                                                                                                                                                                                                                       <  1:261878/249‑1 (MQ=255)                                gTTTCACGAACAATCACATGCAGACGAATACCGAAACGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCCGG‑GAGAGC‑AGGTTTTCTCCGCGCACATGAATATGttttcc                                                                                                                                                                                                                >  1:92512/1‑248 (MQ=255)                                            aTCACATGCAGACGAATACCGAAACGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACATGAATATGTTTACCGTT‑GAAATCG‑Ac                                                                                                                                                                                                  <  2:128022/251‑1 (MQ=255)                                            aTCACATGCAGACGAATACCGAAACGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGACGGG‑AAAAGCCAGTTTTG‑TCCGCGCCCATGAAGATGTTTACCGTT‑TAAATTG‑Ac                                                                                                                                                                                                  >  1:199065/1‑251 (MQ=255)                                                caTGCAGACGAATACCGAAACGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTGCAGTCGC‑GGG‑GAGAGC‑CGGTTTGCGCCGCGCACACGCATCTGTTTTCCGTTTGAAAttaccgtt                                                                                                                                                                                                >  1:96225/1‑244 (MQ=255)                                                 aTGCAGACGAATACCGAAACGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTGGAAT‑GCCGGGGAAGAGA‑GGTTTTTC‑CCGCGCCCCTGAATATGTTTACCCTT‑ggaaacggggggtt                                                                                                                                                                                              >  2:167661/1‑238 (MQ=255)                                                    cAGACGAATACCGAAACGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑GAACGCCGGG‑GAGAGC‑AGTTTTGGTCCGCCCCCAATAGTATGTTTACCGtttagatgggaaaggggcg                                                                                                                                                                                            >  2:29463/1‑234 (MQ=255)                                                    cAGACGAATACCGAAACGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACATGAATATGTTTACCGtt                                                                                                                                                                                                             <  1:29463/234‑1 (MQ=255)                                                     aGACGAATCCCGAAACGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACATGAATATGTTTACCGTT‑GAAATCG‑AC‑‑GGTTtctct                                                                                                                                                                                       <  2:241196/251‑1 (MQ=255)                                                             cccGAAACGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACATGAATATGTTTACCGTT‑GAAATCG‑AC‑‑GGTTTCTCTCT‑GCAATa                                                                                                                                                                              <  2:28432/250‑1 (MQ=255)                                                              ccGAAACGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACATGAATATGTTTACCGTT‑GAAATCG‑AC‑‑GGTTTCTCTCT‑GCAATaa                                                                                                                                                                             <  1:63986/251‑1 (MQ=255)                                                                   aCGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACATGAATATGTTTACCGTT‑GAAATCG‑AC‑‑GGTTTCTCTCT‑GCAATAAAcgcc                                                                                                                                                                        <  2:131189/251‑1 (MQ=255)                                                                     gAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTT‑AATCGCCGGG‑AAAAGC‑AGTTTTGCTCCCCGCCCATGAATATGTTTTCCGTT‑GAAATCGCAG‑‑GGTTT‑TCTCT‑TCAATAAAcgccgc                                                                                                                                                                      >  1:131429/1‑251 (MQ=255)                                                                     gAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGGAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATAGCCGGG‑GAGGGC‑GCATTTGTTCCGCGCACACGCGTATGTTTTCCGTT‑TAAATCG‑TCGAGGTTTCTCTC‑‑‑CAATAAAcgccgc                                                                                                                                                                      >  2:297581/1‑251 (MQ=255)                                                                           tGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCCGTTG‑GAACGCCCGG‑GAGAGC‑AATTTTGCTCCACGGACATTAATATGTTTACCGTT‑GAAAACG‑AC‑‑CGTTTCTCTTT‑GCCATTAACgccgcccgagc                                                                                                                                                                >  1:241663/1‑246 (MQ=255)                                                                             cGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGCAACGGCTGGTG‑AATCGCCGGG‑GAGAGC‑AGTTTTTCTCCCCGCACCTGGAAATGTTTTACGTT‑TAAATCG‑AC‑‑GGGTTCTTTCT‑GCAATAAACtccgcccgagccg                                                                                                                                                              >  2:54563/1‑244 (MQ=255)                                                                              gTCCAGGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACATGAATATGTTTACCGTT‑GAAATCG‑AC‑‑GGTTTCTCTCT‑GCAATAAACGCCGCCAGACCTgg                                                                                                                                                             <  1:261933/251‑1 (MQ=255)                                                                               tCCATGCGCGGCAGCTTTCGCCCGCGCTTGTTCGATTTTCTCTGTAACCAGTTCCGGCGGTTCGCCCCAGGTGTGTTAGAGATCAACCGGTGCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACTGCTGTTG‑AATCGCCCGG‑AAGAGC‑AGTTTTGCTCCGCGCACGTGAATAGCTTTACCGTT‑GAAAGCG‑AC‑‑TGTTTCTCTCT‑GCAATAGAGGCCTCCAGACCTgg                                                                                                                                                             <  2:416730/250‑1 (MQ=255)                                                                                   tGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGGTG‑AAACCCCGGG‑GAGAGA‑CGTTTTGCCCCCCGCCCATGAATATTTTTTCCGTT‑TTAATCG‑AC‑‑GGGTTCTCTCTTGCA‑TAAACGCCCCCAGACCTggggaa                                                                                                                                                         >  2:246838/1‑247 (MQ=255)                                                                                          cAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACATGAATATGTTTACCGTT‑GAAATCG‑AC‑‑GGTTTCTCTCT‑GCAATAAACGCCGCCAGACCTGGGTAAATTccgcc                                                                                                                                                 <  2:92512/251‑1 (MQ=255)                                                                                                  gCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACATGAATATGTTTACCGTT‑GAAATCG‑AC‑‑GGTTTCTCTCT‑GCAATAAACGCCGCCAGACCTGGGTAAATTCCGCCGAGGCTTc                                                                                                                                         <  1:54563/251‑1 (MQ=255)                                                                                                   cccGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑GAGAGC‑AGTTTTTCTCCGCGCCCCTGAATATGTTTTACGTT‑GAAATCG‑AC‑‑CGGTTCTCTCT‑TCCATAAACGCCGCCAGACCTggggaacatcccccccggcttac                                                                                                                                        >  1:284824/1‑231 (MQ=255)                                                                                                    ccGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGGTG‑AATCGCCGGG‑AAGAGG‑AGTTTTTCTCCGCGCACCTGAAAATGTTTACCGTT‑GGAATCG‑AC‑‑GGTTTTTCTCT‑GCAATAAAAGCCCCCCAACCTTGGGAAATTCCGCCGAgggctcgg                                                                                                                                       >  1:105015/1‑245 (MQ=255)                                                                                                     cGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGGTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGGCAAAGTAAAGTGGCGGATACGGCTGTGG‑AAGCGCCGGG‑AAGAGC‑AGTGTTGCGCCGCGCACATGAATATGTTGACCGTT‑GAAATCG‑AC‑‑GGTTGCTCTCT‑GCAATAAACGCCGCCAGACCTGGGTAAATTCCGCCGAGGCTTCGTa                                                                                                                                      <  1:84407/251‑1 (MQ=255)                                                                                                                    ttCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACATGAATATGTTTACCGTT‑GAAATCG‑AC‑‑GGTTTCTCTCT‑GCAATAAACGCCGCCAGACCTGGGTAAATTCCGCCGAGGCTTCGTAGCGCTCGCTATGAt                                                                                                                        <  1:258964/250‑1 (MQ=255)                                                                                                                      ctctTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGCTACGGCTGTTT‑AATCGCCCGG‑AAGAGC‑AGTTTTTCTCCCCGCACATGAATATGTTTTCCCGT‑GAAAACG‑AC‑‑GGTTTCTCTCT‑GCAATAAACGCCGCCAGACCCGGGGAAATTCCGCCCCGGGTTCGTAACGCCCACTATGATCa                                                                                                                      >  1:80747/1‑250 (MQ=255)                                                                                                                                    ccGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACATGAATATGTTTACCGTT‑GAAATCG‑AC‑‑GGTTTCTCTCT‑GCAATAAACGCCGCCAGACCTGGGTAAATTCCGCCGAGGCTTCGTAGCGCTCGCTATGATCAAGGAACACTCCGTCg                                                                                                       <  1:53381/251‑1 (MQ=255)                                                                                                                                      ggcggTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACATGAATATGTTTACCGTT‑GAAATCG‑AC‑‑GGTTTCTCTCT‑GCAATAAACGCCGCCAGACCTGGGTAAATTCCGCCGAGGCTTCGTAGCGCTCGCTATGATCAAGGAACACTCCGTCGcc                                                                                                     <  1:246838/251‑1 (MQ=255)                                                                                                                                       gcggTTCGCCCCAGGGGAGGTAGAGATCAACCTGTTCTGCCGCCAGCGCCTGGGAGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACATGAATATGTTTACCGTT‑GAAATCG‑AC‑‑GGTTTCTCTCT‑GCAATAAACGCCGCCAGACCTGGGTAAATGCCGCCGAGGCTTCGTAGCGCTCGCTATGATCAAGGAACACTCCGTCGCCt                                                                                                    <  2:132873/251‑1 (MQ=255)                                                                                                                                         ggTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑GAGAGC‑AGTTTTGCTCCGCGCACATGAATATGTTTACCGTT‑GAAATCG‑AC‑‑GGTTTCTCTCT‑GCAATAAACGCCGCCAGACCTGGGTAAATTCCGCCGAGGCTTCGTAGCGCCCGCTATGATCAAGGAACACTCCGTCGCCTGc                                                                                                  >  1:188626/1‑251 (MQ=255)                                                                                                                                          gTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑GAGAGC‑AGTTTTGCTCCGCGCACAAGGATATGTTTTCCGTT‑GAAATTG‑AC‑‑GGTTTCTCTCT‑TCAATAAACGCCGCCCGACCTGGGTAAATTCCCCCGGGGCTTCGTAGCGCCCGCTATTGTCAAGGGACACTCCGTCGCCTGcc                                                                                                 >  1:214214/1‑251 (MQ=255)                                                                                                                                          gTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACATGAATATGTTTACCGTT‑GAAATCG‑AC‑‑GGTTTCTCTCT‑GCAATAAACGCCGCCAGACCTGGGTAAATTCCGCCGAGGCTTCGTAGCGCTCGCTATGATCAAGGAACACTCCGTCGCCTGcc                                                                                                 <  1:193799/251‑1 (MQ=255)                                                                                                                                              gCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTGTAAT‑GCCGGG‑GAGAGC‑AGTTTTTGTCCGCGCCCATGAATATGTTTTCCGTT‑GAAAACC‑AC‑‑GGGTTTTCTCT‑GCAATAAACCCCCCCAGACCTGGGGAAAATCCCCCCAGGCCTCGTAGCGCGCCCTATGTTCAAgggaccacacgtcgccagccagct                                                                                             >  2:52705/1‑229 (MQ=255)                                                                                                                                              gCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACATGAATATGTTTACCGTT‑GAAATCG‑AC‑‑GGTTTCTCTCT‑GCAATAAACGCCGCCAGACCTGGGTAAATTCCGCCGAGGCTTCGTAGCGCTCGCTATGATCAAGGAACACTCCGt                                                                                                         <  1:188786/239‑1 (MQ=255)                                                                                                                                              gCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGGTG‑GATCGCCGGG‑GAGAGC‑AGGTTTTCTTCCCGCACATGGATATGTTTTACGTT‑GAAAATG‑CC‑‑GGTTTTTCTCT‑GCAATAAAAGCCGCCCGAACTGGGGAAATTCCGCCGAGGCTTCGTAGCGCTCGCTTTGATCAAgggaaaatccgctgatcgggag                                                                                               >  2:188786/1‑229 (MQ=255)                                                                                                                                                 ccAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑GATCGCCGGG‑GAGAGC‑AGGTTTGCTCCGCGCACATGAATATGTTTACCGTT‑GAAATCG‑AC‑‑GGGTTCTCTCT‑GCAATAAACGCCGCCCGACCTGGGTAAATTCCGCCGCGGCTTCGTAGCGCTCGCTATGATCAAAGAACACTCCCTCGCCTGCCAGCCCtt                                                                                          >  1:298453/1‑251 (MQ=255)                                                                                                                                                  cAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGAAATCTGACGATCCACCAAAGTAAAGTGGCGGCTACGGCTGCTG‑AATCGCCGGG‑GAAAAC‑ACATTTGCTCCGCGCACAAGAATATGTTTACCGTT‑TAAAACG‑GC‑‑GGGTTCTCTCC‑GCAAAAAAAACCGCCAGACCTGGGGAAATTCCGCCCAGGCTTCGTAGCGCTCGCTATGATCAAGGAAAACTCCCTCCCCTGcccgatattt                                                                                         >  1:65642/1‑243 (MQ=255)                                                                                                                                                  cAGGTGAGGTAGAGACCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACATGAATATGTTTACCGTT‑GAAATCG‑AC‑‑GGTTTCTCTCT‑GCAATAAACGCCGCCAGACCTGGGTAAATTCCGCCGAGGCTTCGTAGCGCTCGCTATGATCAAGGAACACTCCGTCGCCTGCCAGCTCTtg                                                                                         <  2:28237/251‑1 (MQ=255)                                                                                                                                                     gTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCAGCCACAGTAAAGTGGCGGATACGGCTGTTT‑AATCGCCGGG‑GAGAGC‑AGTTTTGCTCCGCGCCCATGAATATGATTACCGTT‑GAAATCG‑GC‑‑GGTTTCTCTCT‑GCAATAAACGCCGCCAGACCTGGGTAAAATCCGCCGAGGCCTCGTTGCGCTCGCTATGATTAAGGAACACTCCGTCGCCTTCCAGCCCtttt                                                                                        >  1:311791/1‑247 (MQ=255)                                                                                                                                                                  cAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACATGAATATGTTTACCGTT‑GAAATCG‑AC‑‑GGTTTCTCTCT‑GCAATAAACGCCGCCAGACCTGGGTAAATTCCGCCGAGGCTTCGTAGCGCTCGCTATGATCAAGGAACACTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCtgt                                                                          <  1:134577/250‑1 (MQ=255)                                                                                                                                                                   aaCCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCCAAGTAAAGTGGCGGATACGGCTCTTG‑AATCCCCGCG‑AAGGAC‑CGGTTTGGTCCGCACACAAGAATAAGTTTACCGTT‑GCAATTG‑CC‑‑GGTTTCTCTTT‑GCAATAAACGCCCCCAGACCTGGGGAAATTACGCCGAGGCTTCGGAGAGCTCGCTATGATCAAGGAACCCCCCGTCGCCGACCACCTCTTGTGGGTCCCTGCCtgtgg                                                                        >  1:37474/1‑250 (MQ=255)                                                                                                                                                                                   aGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGGTG‑AATCACCGGG‑AAGAGC‑AGTTTTGCTCCGCGCCCAAGGATATGTTTACCGTT‑GAAACCC‑AC‑‑GGTTTTTCTCT‑GCAATAAAACCCGCCAGAACCGGGGAACATCCCCCGAGGCTTCGTAGAGCTCGCTATGATCAAGGGACCCACCGTCGCCCGCCAGCTCTTGTTGATAGCTACCtgtggcccgggtaaaaacagc                                                        >  1:397508/1‑234 (MQ=255)                                                                                                                                                                                                cATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTTCTCCGCGCACATGAATATGTTTAACGTT‑GAAATCG‑AC‑‑GGTTTCTCTCC‑TCAATAAACGCCGCCCGACCTGGGTAAATTCCCCCGCGGCTTCGTAGCGCTCGCTCTGGTCCAAGAACACTCCGTCGCCTCCCACCTCTTGTGGATCGCTCGCTGTGAcc                                                                      >  1:205866/1‑224 (MQ=255)                                                                                                                                                                                                cATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACATGAATATGTTTACCGTT‑GAAATCG‑AC‑‑GGTTTCTCTCT‑GCAATAAACGCCGCCAGACCTGGGTAAATTCCGCCGAGGCTTCGTAGCGCTCGCTATGATCAAGGAACACTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTTTGAcc                                                                      <  2:205866/224‑1 (MQ=255)                                                                                                                                                                                                                           gCGGATACTGCTGTTG‑AATCGGCCGG‑AATAGC‑AGTTTTGCGCCGTGCAGATTAATAGGTTTACCGTT‑GAAATCG‑AC‑‑GGTTTCTCTCT‑GCAATAAACGCCGCCAGACTTGTGTAAATGCCGCCGAGGCTTCGTAGCGCTCGCTATGATCAAGGAACACTCCTTCGCCTGCCAGCGCTTGTGGATCGCTGCCTGTGACCAGTGTATACAACGCACGTCCAGGTGAGAGACGGTCAAGCGTGGTGGCTTggcgg                <  2:59415/251‑1 (MQ=255)                                                                                                                                                                                                                            cGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACATGAATATGTTTACCGTT‑GAAATCG‑AC‑‑GGTTTCTCTCT‑GCAATAAACGCCGCCAGACCTGGGTAAATTCCGCCGAGGCTTCGTAGCGCTCGCTATGATCAAGGAACACTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTggcggg               <  2:188626/251‑1 (MQ=255)                                                                                                                                                                                                                             ggATACGGCTGTTG‑ACTCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACATGAATATGTTTACCGTT‑GAAATCG‑AC‑‑GGTTTCTCTCT‑GCAATAAACGCCGCCAGACCTGGGTAAATTCCGCCGAGGCTTCGTAGCGCTCGCTATGATCAAGGAACACTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTggcgggc              >  2:27252/1‑251 (MQ=255)                                                                                                                                                                                                                                      tGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACATGAATATGTTTACCGTT‑GAAATCG‑AC‑‑GGTTTCTCTCT‑GCAATAAACGCCGCCAGACCTGGGTAAATTCCGCCGAGGCTTCGTAGCGCTCGCTATGATCAAGGAACACTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGGGGCGGCCTGGCGGGCGGCAACGGt     >  2:139169/1‑251 (MQ=255)                                                                                                                                                                                                                                          g‑AATCGCCGGG‑AAGAGC‑AGTTGTGCTCCGCGCACATGAATATGTTTACCGTT‑GAAATCG‑AC‑‑GGTTTCTCTCT‑GCAATAAACGCCGCCAGACCTGGGTAAATTCCGCCGAGGCTTCGTAGCGCTCGCTATGATCAAGGAACACTCCGGCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATGTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAgg  <  2:257705/250‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     CCGCTCGGCGGCCTGCCACGCTTCGTCGTTAGTTTCACGAACAATCACATGCAGACGAATACCGAAACGAATTTTGCGTCCATGCGCGGCAGCTTTCGCCCGCACTTGTTCGATTTTCTCTTTAACCAGTTCCGGCGGTTCGCCCCAGGTGAGGTAGAGATCAACCTGTTCTGCCGCCAGCTCCTGGGCGACATCTGACGATCCGCCAAAGTAAAGTGGCGGATACGGCTGTTG‑AATCGCCGGG‑AAGAGC‑AGTTTTGCTCCGCGCACATGAATATGTTTACCGTT‑GAAATCG‑AC‑‑GGTTTCTCTCT‑GCAATAAACGCCGCCAGACCTGGGTAAATTCCGCCGAGGCTTCGTAGCGCTCGCTATGATCAAGGAACACTCCGTCGCCTGCCAGCTCTTGTGGATCGCTGCCTGTGACCAGGTTAAACAACGCACGTCCATTTGAGAGACGGTCAAGCGTGGCGGCCTGGCGGGCGGCAACGGTAGG  >  NC_000913/995254‑995736

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