Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_pae_tpiA	264,098	G→T	L50L (CTC→CTA)	ykfH ←	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_pae_tpiA	264,098	0	G	T	95.0%	51.3 / NA	20	L50L (CTC→CTA)	ykfH	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (12/7);  total (12/8)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACG‑TCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTT  >  NC_000913_3_pae_tpiA/263963‑264231                                                                                                                                        |                                                                                                                                       gcgcAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGa                                                                                                                           >  2:195341/1‑149 (MQ=255)     aGTGTAATCTCCCGCTGCAGGCCCCTTGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATc                                                                                                                       <  1:111209/149‑1 (MQ=255)                       aGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGACCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCACTAGGACTCCGCAGATAACCGGAACG‑TCCCTGCCgct                                                                                                    >  2:250891/1‑147 (MQ=38)                        gtccccTGGTAGGGTGGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACG‑TCCCTGCCGgtg                                                                                                   <  2:289660/147‑1 (MQ=255)                            ggtGGTAGGGGTGGTCATTGCGGTGGTCCGCTAGTGTCGGTTCAGGGTGCTGCTGGTCAGGCAGGTATAGGGCCCTTTGCGGTCCTGGCGGCGTTCGGCGTTGACCGCTAGGCCTCGGGCGATATCCGGAACG‑TCCCTGCCGGTGTAAt                                                                                               <  2:204885/147‑1 (MQ=255)                            ccgggTAGGGTTGCTCATTACGGTTCTCCACTATTGGCAGTTCAGGGGGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACG‑TCCCTGCCGGTGTAAt                                                                                               <  2:83392/146‑1 (MQ=255)                                          tCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGTAGTCCGGCGATATCCGGAACG‑TCCCTGCCGGTGTAATGACAGACGCTACCg                                                                                 >  1:383365/1‑149 (MQ=255)                                            aTTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATCTCCGGAACG‑TCCCTGCCGGTGTAATGAAAGACCATACCGTg                                                                               >  1:500303/1‑149 (MQ=37)                                                  gTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACG‑TTCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTg                                                                         <  1:439664/149‑1 (MQ=255)                                                  cttcttcACTATTTTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACG‑TCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTg                                                                         <  2:323042/148‑1 (MQ=255)                                                   ttCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATATCGCTGGGACTCCGGCGATATCCGGAACG‑TACCTGCCGGTGTAATGACAGACGCTACCGTGCCACTgg                                                                        >  2:451461/1‑149 (MQ=255)                                                      tcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACG‑TCCCTGCCGGTGGAATGACAGAGGCTACCGTGCCACTGGTAt                                                                     >  1:174874/1‑149 (MQ=255)                                                        cACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGACCGCAACCCGCAAGT‑CAATGACAGACGCTACCGTGCCACTGGTAttt                                                                   >  1:406119/1‑149 (MQ=25)                                                         aCTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACA‑TACCTGCCGGAGTAATGACAGACGCTACCGTGCCCCTGGTATTTg                                                                  >  2:288864/1‑149 (MQ=21)                                                           tATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACG‑TCCCTGCCGGTGTAATGACAGCCGCTACCCTGCCACTGGTCTTTGcc                                                                >  1:226614/1‑149 (MQ=255)                                                                          gTGATGCTCATCAGGCAGGTATAGGTGCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTATGACTCCGGCGATATCCGCAACG‑TCCCTGCCGGTGTAACGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGaaa                                                 <  1:472039/149‑1 (MQ=21)                                                                                           ggcataGGGCCAATTGCGGGGCTGACGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACG‑TCCCTGGCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGAt                                <  2:24918/146‑1 (MQ=18)                                                                                           ggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACG‑TCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGAt                                >  1:148361/1‑149 (MQ=255)                                                                                                                   gcggcgTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACG‑TCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGGATATCGTCATTGccc        >  2:492748/1‑149 (MQ=255)                                                                                                                         ttCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACG‑TCCCTGCCGGTGTAATGACAGACGCTCCCGTGCCACTGGTCTTTGCCGGTGCAGTAGCGAAAGATTCGGGAGTCAGGATGCTGGCGGTTTATCGTCATTGCCCTGCGtt  >  1:310218/1‑149 (MQ=17)                                                                                                                                        |                                                                                                                                       GCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACG‑TCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTT  >  NC_000913_3_pae_tpiA/263963‑264231

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