Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	46,661	T→C	N285N (AAT→AAC)	yaaU →	putative transporter YaaU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	46,661	0	T	C	95.1%	112.9 / ‑2.0	42	N285N (AAT→AAC)	yaaU	putative transporter YaaU
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (19/20);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.22e-01

AAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGGGAATGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTAATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAAT  >  NC_000913/46414‑46886                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                  acgggCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGATTTCGATGAAGAGCAGCCGAAGCAAACCCGTTTTCGCGATATGTTTAATCGCAGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGCTCCCAATGTTCGCCATTTACCCCTTTGGCCCGCAAAACGTTGGATTGTTGGGTTTGGGTGTTGGCCAAAACGCGGCGCTGGGGAAcgtg                                                                                                                                                                                                                                >  2:458736/3‑251 (MQ=255)         gTAAAAGAGTGCGAAGAGATGATGATAAACCTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAg                                                                                                                                                                                                                        <  1:156947/251‑1 (MQ=255)           aaaaGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGCGTTGGCACACACGCGGCGCTGGGGAACGTGGTGATTAGcc                                                                                                                                                                                                                      >  2:132070/1‑251 (MQ=255)                agTGCGCACAGATGATGATAAAACTGATTGGCAAACCGGAAGCGTTCGATGAAGAGCAGCCGCAGCAAACCAGATTTCCAGCTATGTTTAATCGCCACAATTTTCCTCTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATACCAATGTTCGCCATTTACACCGTTGGCCCGTAAATCGTTGGTTTTTTGGGATTGGGTGGTGGCGAAAACGCGGCGCTGGTGAGCGTGGCGATTAGCCTGttt                                                                                                                                                                                                                 >  2:42263/1‑250 (MQ=255)                   gCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCtt                                                                                                                                                                                                               <  1:458258/250‑1 (MQ=255)                                 ataaaaCTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTAt                                                                                                                                                                                                >  2:157666/4‑251 (MQ=255)                                 aaaaaaCGGTTGGGCGAACCGGGGGCTTGCGATGAAGACCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTTCACCTTTGGCCCGCAAATCGTTGGTTTTTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTAt                                                                                                                                                                                                <  2:54306/248‑1 (MQ=255)                                         gTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGCACGTGGTGCTTAGCCTGTTCTTTATGCTCGGCTGTATTccgccgg                                                                                                                                                                                        >  2:384537/1‑250 (MQ=255)                                                      gTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTAt                                                                                                                                                                                                            >  1:467654/1‑218 (MQ=255)                                                                     gagCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCg                                                                                                                                                            <  1:384537/251‑1 (MQ=255)                                                                               agcaATCCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGGGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCAttgttg                                                                                                                                                  >  1:390676/1‑251 (MQ=255)                                                                                ccaAATCCGTTTTCGGGTTTTGTCTAATCGCCGCCATTTTCCTTTGGTGCTTGTTGTTGCCGCCATCTGGGCCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTTCATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGa                                                                                                                                                 <  2:457891/250‑1 (MQ=255)                                                                                   aaCCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTg                                                                                                                                              <  1:114213/251‑1 (MQ=255)                                                                                        gTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGc                                                                                                                                         >  1:182811/1‑251 (MQ=255)                                                                                                      tttAATCGCCGCCATTCTCCTTTTGTTCTGTTTGTTGCCGCCATCTGTCCCTGCCAGGTGATCCCAATGTTCTCCATTTCCACCTTTGGCCCGCACATCGTTGGTTTGTTGGGCTTGGGTGTTGGCCCAAACGCGGCGCTGGGGAACGTGGTGATTCGCCTGTTCTTTATGCTCGGCTGTCTTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCTTCCATTGTTGATTGGCTGCTTTGCCCTGATGAc                                                                                                                           >  2:27094/1‑251 (MQ=255)                                                                                                              cctccATTTTCCTTTTGTTCTGTTTGTTGCCACCATCTGGACCTGCCAGGTGCTCCCAATGTTCGCCATTTACACCTTTTGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTAGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGCTGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTTGCAGCTTGTCCATGATGAcgctggcgc                                                                                                                   >  1:5975/4‑251 (MQ=255)                                                                                                               cgccATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGt                                                                                                                                                                                                                   >  2:121519/1‑154 (MQ=255)                                                                                                               cgccATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGt                                                                                                                                                                                                                   <  1:121519/154‑1 (MQ=255)                                                                                                                        ccTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGtt                                                                                                                                                                              <  1:175816/182‑1 (MQ=255)                                                                                                                        ccTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGtt                                                                                                                                                                              >  2:175816/1‑182 (MQ=255)                                                                                                                                   gtttgttGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTAt                                                                                                                                                                                                            >  1:311505/1‑141 (MQ=255)                                                                                                                                   gtttgttGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTAt                                                                                                                                                                                                            <  2:311505/141‑1 (MQ=255)                                                                                                                                   gtttgtgGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTGGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCTATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCtgat                                                                                                <  1:2297/249‑4 (MQ=255)                                                                                                                                           ccgccATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTaa                                                                                                                                                                            <  1:197827/165‑1 (MQ=255)                                                                                                                                           ccgccATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTaa                                                                                                                                                                            >  2:197827/1‑165 (MQ=255)                                                                                                                                                aTCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGGTTTgggggtggttccggatattggg                                                                                   >  2:196207/1‑231 (MQ=255)                                                                                                                                                       ccTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGGGAATGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGACTGGCAGCTTTGcc                                                                                                                                   >  1:413892/1‑194 (MQ=255)                                                                                                                                                       ccTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGGGAATGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGACTGGCAGCTTTGcc                                                                                                                                   <  2:413892/194‑1 (MQ=255)                                                                                                                                                       catggcaggggttCCCAAGGTTCGTCATTTACACCTTTGGCCCGCAAATCGTGGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCGCTGGGGAACGGGGTGATTAGCCTGTTCTTTAGGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATctggctg                                                                          <  2:350564/239‑1 (MQ=255)                                                                                                                                                                       aaTGTTCGCCATTTACACCTTTGGCCCGCAAATCGTGGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTAGGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATGTGGGGATCTGGCTGGTAGTGATGGCCTTTg                                                          <  2:390676/251‑1 (MQ=255)                                                                                                                                                                         tGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGCTTTGTTGGGATTGGGTGTTGGCCAAAACGCGGCGCTGGGGACCGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGCGGCTGGGGTTAAACACTGCCGGACGGAGTCCATTGTTGATTGGCAGCTTTGCCATGATGCCGCTGCCGCGGGCGGTTTTGGGGGTGGTCCCGGGTGTGGGGAttggggtggggtgggtgggctttgtg                                                        >  2:37219/1‑226 (MQ=255)                                                                                                                                                                             cGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGgtgt                                                    <  1:157666/251‑1 (MQ=255)                                                                                                                                                                                   ttacacctgtgcggcGCAAATCGTTGGTTTGTTGGGATTGGGTGTTAGCAATAACGGGGCGTTGGGGAGCGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGGATGCCGCCGATGCTCTGGTTAACCACTGGCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGATGGCGGTTTTGTGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTGGCGGTGTATGcc                                               <  1:475279/236‑1 (MQ=255)                                                                                                                                                                                     acacCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTGATGCTCGGCTGTATTCCTCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCttt                                            <  1:390018/251‑1 (MQ=255)                                                                                                                                                                                               cccGCAAATTGTCGGATTGGTGGGATTGGGGGTTGGCAAAAACGCGGCGCTGGGGGACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTATGCCATGATGACGCTGGCGCTGGCGGATTTGGGGGTGAGCCCGGGTATGGGGACCTGGTGGGGGGGGAGTGCCTTTGTGGTGTATGCCtttttgttctgcg                                  >  2:83036/1‑243 (MQ=255)                                                                                                                                                                                                        cGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGAt                                                                                                                                                >  2:140602/1‑132 (MQ=255)                                                                                                                                                                                                        cGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGAt                                                                                                                                                <  1:140602/132‑1 (MQ=255)                                                                                                                                                                                                         gTTGGTTTGTTGGGATTGTGTGTTGGCAAAAACGCGTGGCTGGGGAAGGTGGTGAGTAGCCTGTTCTTTATTCTCGGCTGTATTCCGCCGATGCGTTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCTGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTaa                        <  2:388897/251‑1 (MQ=255)                                                                                                                                                                                                           tgggttgttgGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATgg                                                                <  1:368741/205‑1 (MQ=255)                                                                                                                                                                                                           tGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTTGGGCTGATCCCGGGTATGGGGGTCTGGCTGGTGGTGATGGGCTTTTCGGTGTAAGCCTTTTTTTCTGGCGGGGCGGGGAAtt                      >  1:136036/1‑251 (MQ=255)                                                                                                                                                                                                             gTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTg                    <  1:403237/251‑1 (MQ=255)                                                                                                                                                                                                                               ttGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGGTATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGGGTATGCCTTTTTTTCTGGCGGGCCGGGGAATTTGCAGTGGCGCTATCCTAAt  >  1:354038/1‑251 (MQ=255)                                                                                                                                                                                                                                tGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAAt  <  1:333785/250‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                  AAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGGGAATGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTAATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAAT  >  NC_000913/46414‑46886

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