Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	246,747	G→T	S12S (TCG→TCT)	yafL →	NlpC/P60 family protein YafL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	246,747	0	G	T	94.4%	50.9 / NA	18	S12S (TCG→TCT)	yafL	NlpC/P60 family protein YafL
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (10/7);  total (10/8)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGGGGTATCTCCTGTGCTGTATATTTATTCAGCTTGAATTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGCACAACCGGGAGCAGATAAAACAGCATCCAGAATGGTTTCCTGCTCCGCTCAAG  >  NC_000913/246502‑246993                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       ttGGGGTATCTCCTGTGCTGTATATTTATTCAGCTTGAATTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGAtt                                                                                                                                                                                                                                                   >  2:139525/1‑251 (MQ=255)                 gCTGTATATTTATTCAGCTTGAATTGTAGCTCAACTGAGCAACCGCGACACTCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGACATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGAtttgtttg                                                                                                                                                                                                                                   >  1:96897/1‑251 (MQ=255)                                                                     ccATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGATATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTAcg                                                                                                                                                                               <  2:28201/251‑1 (MQ=255)                                                                          aGTAATTTAATAAACGGATGGTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGATTGTTCAGAAACGATCCGTGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACgacgcc                                                                                                                                                                          >  2:223116/1‑251 (MQ=255)                                                                           gTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGccgccc                                                                                                                                                                         <  1:30256/251‑1 (MQ=255)                                                                                                                                 ttATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTg                                                                                                                    >  2:57435/1‑250 (MQ=255)                                                                                                                                           ggTTTATTAAGTTTGTTCAGAAACGACCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGc                                                                                                         <  1:4070/251‑1 (MQ=255)                                                                                                                                                            cAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGa                                                                                          >  1:179849/1‑249 (MQ=255)                                                                                                                                                                       tgggATAAACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTTCTTTCCGTCGATTCCATCCTTTGTATTTACTGGTTTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTTCGCCGCTCGCCTGCGGATGCAAAACCGTGCGCTTTTATTAAAACAGTACAAAACTCCTCTGAAACAGCATGCAAGCTTTATATTGGAAGGCactaccga                                                                              >  2:77859/2‑243 (MQ=255)                                                                                                                                                                        gggATACACTGCAACATTAAAGCATGACCAGCCTTTATTGAGTATGTCCTTGCCGTCGTTTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCTTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGaaa                                                                            <  1:33954/251‑1 (MQ=255)                                                                                                                                                                               acTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCaaa                                                                       >  2:82179/1‑249 (MQ=255)                                                                                                                                                                                                tGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGcac                                                     >  1:197650/1‑250 (MQ=255)                                                                                                                                                                                                   ccAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAACACAGAACCAACCTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGACGAAAGCAAAAGGGCGCTACGCCAGCACAAcc                                                 >  1:236791/1‑251 (MQ=255)                                                                                                                                                                                                        cATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGCACAACCGGGAg                                            <  1:211018/251‑1 (MQ=255)                                                                                                                                                                                                            attGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGCACAACCGGGAGCAg                                         <  1:116011/250‑1 (MQ=255)                                                                                                                                                                                                                             cGTCGATTCCATCCTTTGTATTTTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGCACAACCGGGAGCAGATAAAACAGCATCCAGaa                       <  1:194684/251‑1 (MQ=255)                                                                                                                                                                                                                              gTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCCCATCACATATCCCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGCACAACCGGGAGCAGATAAAACAGCATCCAGAAt                      <  2:245157/251‑1 (MQ=255)                                                                                                                                                                                                                                                  tGTCTGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACCTCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGCACAACCGGGAGCAGATAAAACAGCATCCAGACTGGTTTCCTGCTCCGCTCacg  >  2:84134/1‑249 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       TTGGGGTATCTCCTGTGCTGTATATTTATTCAGCTTGAATTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATGGTTTATTAAGTTTGTTCAGAAACGATCCGGGATACACTGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCCTTGCCGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGTTTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGCACAACCGGGAGCAGATAAAACAGCATCCAGAATGGTTTCCTGCTCCGCTCAAG  >  NC_000913/246502‑246993

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