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	96.5%	81.0 / NA	29	L50L (CTC→CTA)	ykfH	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (15/13);  minor base C (0/1);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTA  >  NC_000913_3_pae_tpiA/263954‑264233                                                                                                                                                 |                                                                                                                                        ggCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCGGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGa                                                                                                                                     <  2:388008/149‑1 (MQ=255)    tCCCAGGCGCAGTGTAATCACCCGCTGCAGGCCCCTGGTAGAGCTCCTCTGGACAGAACTAGCATATTGTCAGTACAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGACGTTCGGCGTATACCGCTAGGACTc                                                                                                                                  >  1:689449/1‑149 (MQ=255)      ccAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCg                                                                                                                                >  1:390895/1‑149 (MQ=255)       cAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCgg                                                                                                                               <  2:500935/149‑1 (MQ=255)              aGTGTAATCTCCCGCTCAAGGCCCCTGGTAGGGTTGCCCAAAAGAGTTCTCCACAATTGTCATGGCGAGGGGATGCTTATCAGGCAGGTATGGGGCCCATTTCGGTCCTGGCGGCGTTCGGCGTTTACCGCTAGGCCCCCGGCGATATc                                                                                                                        >  2:148646/1‑149 (MQ=255)                tgtAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTGCGCGTATACCGCTAGGACTCCGGCGATATCCg                                                                                                                      >  1:507955/1‑149 (MQ=255)                     tctcCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAg                                                                                                                                       <  2:599942/127‑1 (MQ=255)                     tctcCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAg                                                                                                                                       >  1:599942/1‑127 (MQ=255)                                    cccTGGTAGGGTTGCTCATTACAGATCTCCACTATTGTCAGGTCAGGGTGATGCTCATCCGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTaa                                                                                                  >  1:605323/1‑149 (MQ=255)                                     ccTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAAt                                                                                                 <  1:191189/149‑1 (MQ=255)                                               ttGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGTTATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCt                                                                                       >  2:311388/1‑149 (MQ=255)                                               ttGATCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGCTCGGCGTATACCACCAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCt                                                                                       <  1:551455/149‑1 (MQ=255)                                                      ttACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTgg                                                                                >  1:60069/1‑148 (MQ=255)                                                          aGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGGCAGACGCTACCGTGCCACt                                                                            >  1:454022/1‑149 (MQ=255)                                                            ttCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTgg                                                                          <  1:83498/149‑1 (MQ=255)                                                                 cACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAg                                                                                                                                       <  2:550359/83‑1 (MQ=255)                                                                  aCTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGCTGTAATGACAGACGCTCCCGTGCCCATGGTGTTTg                                                                    >  2:693804/1‑149 (MQ=255)                                                                      ttGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCgg                                                                <  1:483/149‑1 (MQ=255)                                                                         tCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTAt                                                                       >  1:548203/1‑139 (MQ=255)                                                                         tCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTAt                                                                       <  2:548203/139‑1 (MQ=255)                                                                                     gATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTg                                                              <  1:208855/136‑1 (MQ=255)                                                                                     gATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTg                                                              >  2:208855/1‑136 (MQ=255)                                                                                         cTCATCAGGCAGGTCTAGGGCCCTTTGCGGTCCTGGCGGCGTACGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGGAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCg                                             <  2:689449/149‑1 (MQ=21)                                                                                          tcatcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGCCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGCCAGACGCTACCGTGCCGCTGGTATTTGCCGGTGCAGTAGCGAAAGCTTCgg                                            >  2:81888/1‑149 (MQ=255)                                                                                                ggcaggTATAGGGCCCGTTGCGGTCCCGGCGGCGTTCGGCGTATACCGGTAGGAATTCTGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCa                                      <  2:353524/149‑1 (MQ=17)                                                                                                                         cTGGCGGCGTTCGTCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTCCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCAttt             >  2:618698/1‑148 (MQ=255)                                                                                                                          acgcggcgTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTCATGACCGCCGCTCCCGTGCCCCTGCTTTTTGCCGATGCAGTAGCGACAGATCCTGGACTCCGGCTCTTCGCGGCCTATCCTCATTGc            >  2:487596/3‑149 (MQ=255)                                                                                                                               gcgTGCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGc       <  2:655633/149‑1 (MQ=255)                                                                                                                                    cGGCGTATCCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTa  <  1:673200/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                        GGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTA  >  NC_000913_3_pae_tpiA/263954‑264233

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