Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	46,652	A→G	A282A (GCA→GCG)	yaaU →	putative transporter YaaU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	46,652	0	A	G	95.5%	126.5 / ‑2.0	44	A282A (GCA→GCG)	yaaU	putative transporter YaaU
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (20/22);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

CGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGGGAATGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTAATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAAT  >  NC_000913/46410‑46886                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           cGCAAAGGGGGAGTAAAAGAGTGCGAATAGATGATGTGAAAACTGTGTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGTTCTGTTTAATCGCCGCCATTTTCCTTTTTTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTTATCCCAATGTTCGCCATTTACACCTTTGGCTCGCAAATCGTTGGTTTGGTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGaa                                                                                                                                                                                                                                    <  2:195933/251‑1 (MQ=255)     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                                                                                                                                                                                                            <  2:311505/141‑1 (MQ=255)                                                                                                                                       gtttgttGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTAt                                                                                                                                                                                                            >  1:311505/1‑141 (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                                                                                                                                   <  2:413892/194‑1 (MQ=255)                                                                                                                                                           ccTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGGGAATGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGACTGGCAGCTTTGcc                                                                                                                                   >  1:413892/1‑194 (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)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           CGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGGGAATGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTAATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAAT  >  NC_000913/46410‑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: ‑