Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	421,027	C→T	R14R (CGC→CGT)	proY →	putative transporter ProY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	421,027	0	C	T	98.5%	247.2 / NA	68	R14R (CGC→CGT)	proY	putative transporter ProY
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (40/27);  minor base G (0/1);  total (40/28)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CGTCCTGGGCGCAGCGTTTACCGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTACCGGCTGGACCTACTG  >  NC_000913/420780‑421272                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      cGTCCTGGGCGCAGCGTTTACCGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTcac                                                                                                                                                                                                                                                    >  1:406745/1‑251 (MQ=255)  gTCCTGGGCGCAGCGTTTACCGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTcaca                                                                                                                                                                                                                                                   <  1:464940/251‑1 (MQ=255)    ccTGGGCGCAGCGTTTACCGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACAt                                                                                                                                                                                                                                                  <  2:356988/250‑1 (MQ=255)    ccTGGGCGCAGCGTTTACCGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACAt                                                                                                                                                                                                                                                  >  2:500085/1‑250 (MQ=255)       gggCGCAGCGTTTACCGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGc                                                                                                                                                                                                                                              >  2:505258/1‑251 (MQ=255)                gTTTACCGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCa                                                                                                                                                                                                                                     <  2:116551/251‑1 (MQ=255)                       gCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATTGCACTGGGtt                                                                                                                                                                                                                              <  2:377266/251‑1 (MQ=255)                              gAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAAt                                                                                                                                                                                                                       <  1:238089/251‑1 (MQ=255)                               aaCAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAAtt                                                                                                                                                                                                                      >  1:45064/1‑251 (MQ=255)                                  aaGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGc                                                                                                                                                                                                                   >  1:242314/1‑251 (MQ=255)                                                         aCAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGc                                                                                                                                                                                            >  1:65425/1‑251 (MQ=255)                                                          cAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTACGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCa                                                                                                                                                                                           <  2:164416/251‑1 (MQ=255)                                                             tggtgATGGTGGTTCTGGCAATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTCCCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGAc                                                                                                                                                                                        >  2:234966/1‑251 (MQ=255)                                                              ggtgATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACg                                                                                                                                                                                       >  2:317241/1‑251 (MQ=255)                                                                 gATGGTGGTTATGGGGATTATCTGGGATCCCGCGGCAGCTCGTCAGGTGACATCCCGCGCTCACTCAATCACTGAACATTTGTTTTAACCACGGGGCTGCGAGGCACCGCGGTTTTTTCTTGGGTTGATGGGTTAGGAATTGTTGGTAAGTAAGAACAAGCTAAAGCGTGAGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGccc                                                                                                                                                                                    >  2:230096/1‑250 (MQ=255)                                                                    ggtggtTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATc                                                                                                                                                                                  >  1:362343/1‑250 (MQ=255)                                                                    ggtggtTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATc                                                                                                                                                                                  >  2:65506/1‑250 (MQ=255)                                                                             ggCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGcc                                                                                                                                                                         <  1:96288/250‑1 (MQ=255)                                                                             ggCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGc                                                                                                                                                                          <  1:560685/249‑1 (MQ=255)                                                                                 aTTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTcc                                                                                                                                                                    >  1:175112/1‑251 (MQ=255)                                                                                          gATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGTGCTAAGTACCCGTCACATACGCTGTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGATCGtgtt                                                                                                                                                           <  2:520856/251‑1 (MQ=255)                                                                                             cgcgcgGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGGAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGGACCGGGCTGTTTTACGGTTCGGCATACGCCACCAAAATGGCCGGGCCGAGCGtgttgt                                                                                                                                                         <  2:238356/250‑1 (MQ=255)                                                                                             agcgcgGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTGATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCGAAAGCGTGGGCTAAGTACCCGGCACATACGCTGTATGGCACTGGGTTCAGCAATTGGTACCGGGCTGTGTTACGGTGCGGCAGACGCCATCAAAATGGCCGGTCCGTGAGtgttgt                                                                                                                                                         <  2:57863/249‑1 (MQ=255)                                                                                                   gCAGGTCGTCGGGTGGCCTCCAGCGCTCACTAAATCACTGAACATTTGTATTAACCACGGGGCTTCGAGGGCCCGTGGTTTGTTATTGTGTTGATGGGTTAGGAATGTATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGTCACTGGGTTCAGCATTTGTCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAACGGCCGGTCCGAGctggctgtttgccta                                                                                                                                                  <  2:72423/251‑15 (MQ=255)                                                                                                         cgtcAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTAtattat                                                                                                                                            >  1:270545/1‑251 (MQ=255)                                                                                                         cgtcAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTAtattat                                                                                                                                            >  2:270475/1‑251 (MQ=255)                                                                                                          gtcAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTCACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATc                                                                                                                                           >  1:526851/1‑251 (MQ=255)                                                                                                                        aGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCgcggc                                                                                                                             <  1:575459/251‑1 (MQ=255)                                                                                                                           gcTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGt                                                                                                                           >  2:432305/1‑250 (MQ=255)                                                                                                                                  aaaTCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACACTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCAtt                                                                                                                   <  2:406745/251‑1 (MQ=255)                                                                                                                                      cACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTAtgc                                                                                                               >  2:130816/1‑251 (MQ=255)                                                                                                                                            aCATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTgcgc                                                                                                         <  1:14614/251‑1 (MQ=255)                                                                                                                                             cATTTGTTTTAACCCAGGGTCTGCGAGGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGTTGGAAAGTAAGAACAAGCTAAAGCGTGGTCTAAGTACCCGTCACATACGCTTTATGGGACTGGGTTCATCAAGTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTAGCGTTGGTATCGCGGCGTATATCATTATGCGTgcgcg                                                                                                        <  2:563277/251‑2 (MQ=255)                                                                                                                                             cATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGTATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTCTATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCt                                                                                                        >  2:414434/1‑251 (MQ=255)                                                                                                                                              aTTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTAGTGTGTTGATGGGTTAGGAATTGATGAAAAGTAAGAACAAGCTAAAGCGTGTGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCCATTGGAACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGTCGTATATCATTATGCGTGCGCTg                                                                                                       <  1:409566/251‑1 (MQ=255)                                                                                                                                                  gTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTgg                                                                                                                                                                                                                                 >  2:499693/1‑125 (MQ=255)                                                                                                                                                  gTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTgg                                                                                                                                                                                                                                 <  1:499693/125‑1 (MQ=255)                                                                                                                                                  gTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTggggg                                                                                                   >  1:85855/1‑251 (MQ=255)                                                                                                                                                  gTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTggggg                                                                                                   <  1:25132/251‑1 (MQ=255)                                                                                                                                                    tttAACCGCGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGa                                                                                                  <  1:446018/250‑1 (MQ=255)                                                                                                                                                       aaCCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAAt                                                                                               <  1:317241/250‑1 (MQ=255)                                                                                                                                                       aaCCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATg                                                                                              >  2:120645/1‑251 (MQ=255)                                                                                                                                                              gggCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTAc                                                                                       >  1:66211/1‑251 (MQ=255)                                                                                                                                                                          cccGTGGTTTTTTATTGTGTTGCTGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGgccg                                                                           >  1:500436/1‑251 (MQ=255)                                                                                                                                                                               ggTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGc                                                                      >  1:50264/1‑251 (MQ=255)                                                                                                                                                                                 ttttttATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGctc                                                                    >  1:526379/1‑251 (MQ=255)                                                                                                                                                                                        ttGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGGTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTATTTATcg                                                             >  1:450099/1‑251 (MQ=255)                                                                                                                                                                                            gtTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGtt                                                         <  1:259623/251‑1 (MQ=255)                                                                                                                                                                                                 tGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATc                                                                                                                                  >  2:186253/1‑173 (MQ=255)                                                                                                                                                                                                       aGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGaa                                              <  1:503159/251‑1 (MQ=255)                                                                                                                                                                                                       aGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGa                                               <  2:526851/250‑1 (MQ=255)                                                                                                                                                                                                        ggAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGaa                                              >  2:520932/1‑250 (MQ=255)                                                                                                                                                                                                         gAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGaaaa                                            <  1:500085/251‑1 (MQ=255)                                                                                                                                                                                                          aaTTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAAc                                           >  1:14565/1‑251 (MQ=255)                                                                                                                                                                                                                    aaGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCAGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGc                                 >  1:22785/1‑251 (MQ=255)                                                                                                                                                                                                                     aGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCt                                >  2:274752/1‑251 (MQ=255)                                                                                                                                                                                                                     aGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGCCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCt                                >  1:62222/1‑251 (MQ=255)                                                                                                                                                                                                                      gTAAGAACAAGCTAAAGCGTGGGCTCAGTACCCGTCACATACGCTTTATGGCACTGGGCTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCTGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTg                               >  2:387949/1‑251 (MQ=255)                                                                                                                                                                                                                       tAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTgg                              >  1:461210/1‑251 (MQ=255)                                                                                                                                                                                                                       tAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTgg                              >  1:165924/1‑251 (MQ=255)                                                                                                                                                                                                                           aaCAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAgg                          >  2:176943/1‑251 (MQ=255)                                                                                                                                                                                                                               aGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTAc                      >  1:554994/1‑251 (MQ=255)                                                                                                                                                                                                                               aGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGAAGGTTAc                      >  2:481652/1‑251 (MQ=255)                                                                                                                                                                                                                                  tAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACAtt                   >  2:19122/1‑251 (MQ=255)                                                                                                                                                                                                                                          gggCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGGGTTGTGGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTACCGGCTg           <  2:74745/251‑1 (MQ=255)                                                                                                                                                                                                                                              tAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTACCGGCTGGAcc       >  1:96769/1‑251 (MQ=255)                                                                                                                                                                                                                                               aaGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTACCGGCTGGAcc       <  1:58240/250‑1 (MQ=255)                                                                                                                                                                                                                                                   ccccGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTACCGGCTGGACCTACTg  <  2:121797/250‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      CGTCCTGGGCGCAGCGTTTACCGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTACCGGCTGGACCTACTG  >  NC_000913/420780‑421272

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