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.1%	154.8 / NA	53	R14R (CGC→CGT)	proY	putative transporter ProY
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (27/24);  minor base G (0/1);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GTCCTGGGCGCAGCGTTTACCGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTACCGGCTGGACCTACTG  >  NC_000913/420781‑421272                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      gTCCTGGGCGCAGCGTTTACCGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTcac                                                                                                                                                                                                                                                    >  1:368894/1‑250 (MQ=255)  tCCTGGGCGCAGCGTTTACCGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACAt                                                                                                                                                                                                                                                  >  1:418741/1‑251 (MQ=255)           cAGCGTTTACCGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTAt                                                                                                                                                                                                                                         >  1:40210/1‑251 (MQ=255)            ggCGTTTTCCGCTGGCCGAACAAGGTCTTGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGGCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATg                                                                                                                                                                                                                                        <  2:368894/250‑1 (MQ=255)                 ttACCGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGACGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACACTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGTGGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGAACCCGTCACAAAAGCTTTATGGcgct                                                                                                                                                                                                                                   >  2:494565/1‑248 (MQ=255)                        tGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGACCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTACGAACAAGCTAAAGCGTGGGCTAAGTACCCGACACATACGCTTTATGGCACTgggatc                                                                                                                                                                                                                             >  1:411623/1‑247 (MQ=255)                               aCAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTg                                                                                                                                                                                                                     >  1:91652/1‑251 (MQ=255)                                    gTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCAcc                                                                                                                                                                                                                >  1:120588/1‑251 (MQ=255)                                                  aTGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACgg                                                                                                                                                                                                  <  1:243574/251‑1 (MQ=255)                                                         cAGTGGTGATGGTTGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCa                                                                                                                                                                                           <  1:85367/251‑1 (MQ=255)                                                                         tcttgCCATTATCTGGGATCGCGCTGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGTGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATg                                                                                                                                                                            <  1:524828/246‑1 (MQ=255)                                                                            ggCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCg                                                                                                                                                                        >  1:493930/1‑251 (MQ=255)                                                                            ggCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCg                                                                                                                                                                        <  1:51523/251‑1 (MQ=255)                                                                             gCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCgg                                                                                                                                                                       <  1:494565/251‑1 (MQ=255)                                                                                      tGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACCTTTGTTTTAACCACGGGGCTGCGATGCCCCGGGGTTTTTTATTGTGTTGATGGGATAGGAATTGATGGCAAGTCAGAACAAGCTAAAGCGTGGGCTAAGTACACGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGAAACGGGCTGTTTTACGGTTCGACCGACACCATCCAAATGGCGGTTCCGAgagt                                                                                                                                                              >  2:3442/1‑248 (MQ=255)                                                                                      gtgtaTCGCGCGGCAGTGCGTCAGGTGTCCTCCAGCGCTGCCTAAATCACTGAACCTTTGTTGTAACCCCGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATTGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCTATCAAAATGGCCGGTCCGAGCgt                                                                                                                                                              <  2:519003/247‑1 (MQ=255)                                                                                      gccgatagcgctgtaGGTCGTCAGGTGATCTCCAGCTCTCACTAAATTTCTGATCAGTTGTTTTAACCACGGGGATGCGATGCCCCGTGGTTTTTTATGGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTACAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGTAGTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGGTCGGCCGTCGCTATAAAAATGGCCGGTCTGAGCg                                                                                                                                                               <  2:448408/236‑1 (MQ=255)                                                                                       gggATCGCGCGGAAGGTCGTCAGGTGACCGCCAGCGCTCACTAAATCAATGAACATTTGTTGTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCg                                                                                                                                                               <  2:418741/249‑1 (MQ=255)                                                                                             gcgcgGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGttgttg                                                                                                                                                       <  1:298002/251‑1 (MQ=255)                                                                                              cgcgGCAGGTCGTCAGGTGACCTCCAGCGCTCACTGAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTCCGGTTCGGCAGACGCCATCAAAATGGCCGGACCGAGCgagagcgga                                                                                                                                                       >  1:448408/1‑242 (MQ=255)                                                                                                            aGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAACGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCCGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCggt                                                                                                                                        >  1:475248/1‑251 (MQ=255)                                                                                                               gtaCCTCAAGCGCTCACTAAAGCACTGAACAGTTGTTTTAACCACGGGGCTGGTATGCCCCGTGGTTTTTTATTGTGTTGAGGGGGTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATGGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTTGCCTATATTATCggtggt                                                                                                                                     <  2:40210/249‑1 (MQ=255)                                                                                                                             cACTAAATCACTGAACATTTGTTTTATCCACGGGGCTGCGATGCCCCTTGGTTTTGTATTGTGTTGATGTGTTAGGAATTGATGGAAAGTAAGAACAAGTTAAAGCGTGGGCTAAGTACCCGTCACATACGTTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTGTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTatat                                                                                                                       <  1:141718/251‑1 (MQ=255)                                                                                                                               cTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATc                                                                                                                      >  1:365494/1‑250 (MQ=255)                                                                                                                                 aaaTCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCa                                                                                                                     <  1:189646/249‑1 (MQ=255)                                                                                                                                        tGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTAtgcgtg                                                                                                            >  1:580539/1‑251 (MQ=255)                                                                                                                                         gaacatgggTTTTAACCACGGGGCTGTGATGCCCCGTTGTTTTTTATTGTGTGGTGGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCCCTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTAGATCATTAtgcgt                                                                                                             <  2:112338/241‑1 (MQ=255)                                                                                                                                         gAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTAtgcgtg                                                                                                            <  1:104495/250‑1 (MQ=255)                                                                                                                                                         cACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTTTATTATCGGTGGTATCGCGGCGTATATCATTATGCGTTCGCTGGGGGGAATggcg                                                                                           >  2:496928/1‑248 (MQ=255)                                                                                                                                                            ggggCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGCGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGt                                                                                         >  1:195020/1‑250 (MQ=255)                                                                                                                                                             gggCTGCGAGGCCCCGTGGTGTTTTATTGTGTGGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTCCCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTAc                                                                                       <  2:21831/251‑1 (MQ=255)                                                                                                                                                                   cGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAAcc                                                                                 >  1:134429/1‑251 (MQ=255)                                                                                                                                                                      ggctccgtgtgttttgATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTCCCCGTCACATACGCTTTATGGCACTGGTTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCAGCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCg                                                                               <  2:381561/239‑1 (MQ=255)                                                                                                                                                                           cggggTTTGTTATTGTGTTGAGGGGTTAGGAATTGATGTAAAGGAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTGTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGgtcgca                                                                         <  2:195020/248‑2 (MQ=255)                                                                                                                                                                               gTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGct                                                                     <  1:411038/251‑1 (MQ=255)                                                                                                                                                                                   tttATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCttt                                                                 >  1:203876/1‑251 (MQ=255)                                                                                                                                                                                           gtTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGTTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTCTATCATTATGCGTGCGCTGGGGGAAATGTCGGTACCGAACCCGGCCGCACGCTCTTTCTCGCGtt                                                         >  2:210063/1‑251 (MQ=255)                                                                                                                                                                                           gtTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGCGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGGATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGtt                                                         >  1:429831/1‑251 (MQ=255)                                                                                                                                                                                               aTGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGACCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTCTgc                                                     >  1:140690/1‑251 (MQ=255)                                                                                                                                                                                                          aatGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATGGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGgcc                                                                            <  2:57073/215‑1 (MQ=255)                                                                                                                                                                                                          aTTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGgcc                                                                            >  1:57073/1‑217 (MQ=255)                                                                                                                                                                                                             gATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCg                                       >  2:86661/1‑251 (MQ=255)                                                                                                                                                                                                             gATGGAAAGTAAGAACAAGATAAGGCTTGGGCCAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAAGTGGCTCCGGGCTGGTTTACGGTTCGGAAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTc                                        <  2:178124/250‑1 (MQ=255)                                                                                                                                                                                                                   aatgaaGAACAAGCTAAAGCGTGGGATAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACAGGGCTGATTTACGGTTCGGCAGACGCCATCAGAATGGCCGGTCCGCGAGTGTTGTATGCCTATCTTATAGGGGGAGACGAGGCGCAAGTCATTGTGCGTGCGCTGGGGGTAAAGTCGGAGAATAACGCGTCCGCCAGCCCTTGCTCGCGTTATGAGCAGGAAAACCTCGGCACGc                                 >  2:403007/5‑251 (MQ=255)                                                                                                                                                                                                                    aGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCt                                >  1:328547/1‑251 (MQ=255)                                                                                                                                                                                                                          aaacatcgaaagagtgtgtaaacaccccGGCACATACGCGTTATGGCACCGGGTTTAGCAATTGGCCCCGGGCTGTTTTACTGTTCGGCAGACGCCATCAAAATGGCCGGTCAGAGTGTGATGTTGGCCTATATTATCGGTGTTATCGCGGCGTATATCGTTATGCTTGCGCTGGGGGAAACGTCGGTACATAACCCGGCCGCCAGCTCGTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAgg                          <  2:181784/226‑1 (MQ=255)                                                                                                                                                                                                                            cAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTg                               >  1:668520/1‑244 (MQ=255)                                                                                                                                                                                                                            cAAGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGACCTGGGTTCAGCAATTGTCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTTGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTg                               <  2:668520/244‑1 (MQ=255)                                                                                                                                                                                                                             aaGCTAAAGCGTGGGCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTATCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGtt                        >  1:632088/1‑250 (MQ=255)                                                                                                                                                                                                                                          gtttAAGTACCCGTCACATACGCTTTAGGGCACTGGGTTCAGCAATTTGCACCGGGCTTTTTTACGGTTCTGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCTCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTACCGGCTgg          <  2:580539/248‑1 (MQ=255)                                                                                                                                                                                                                                          ggCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTTGCAGGTTACATTaagatcggg          >  2:531477/1‑243 (MQ=255)                                                                                                                                                                                                                                          ggCTAAGTACCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTa                  <  1:531477/243‑1 (MQ=255)                                                                                                                                                                                                                                                  aCCCGTCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTACCGGCTGGACCTACTg  >  1:159639/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      GTCCTGGGCGCAGCGTTTACCGCTGGCCGAACAAGGTCTGGCGTGGTTAATGCCAACAGTGGTGATGGTGGTTCTGGCCATTATCTGGGATCGCGCGGCAGGTCGTCAGGTGACCTCCAGCGCTCACTAAATCACTGAACATTTGTTTTAACCACGGGGCTGCGATGCCCCGTGGTTTTTTATTGTGTTGATGGGTTAGGAATTGATGGAAAGTAAGAACAAGCTAAAGCGTGGGCTAAGTACCCGCCACATACGCTTTATGGCACTGGGTTCAGCAATTGGCACCGGGCTGTTTTACGGTTCGGCAGACGCCATCAAAATGGCCGGTCCGAGCGTGTTGTTGGCCTATATTATCGGTGGTATCGCGGCGTATATCATTATGCGTGCGCTGGGGGAAATGTCGGTACATAACCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTACCGGCTGGACCTACTG  >  NC_000913/420781‑421272

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: ‑