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

CGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGC  >  NC_000913_3_pae_tpiA/263952‑264228                                                                                                                                                   |                                                                                                                                   cGGGCGCCCAGGCGAAGTGTAGTGTGCCGCGGCAGGCCCCGGGTAGGGTTGCTCATTATGGTGCTCCACTATTGTCAGTTCAGGGTGGTGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAg                                                                                                                                  <  2:43085/149‑1 (MQ=255)          aGGCGCAGTGTAACCTCCCGCGGCAGGCCCCGGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGAGGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGc                                                                                                                         <  2:360895/149‑1 (MQ=255)                             gCTGCAGGCCCGTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATTCCCGGAACGTCCCTg                                                                                                      >  1:450775/1‑149 (MQ=255)                                    gCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgtgt                                                                                               <  2:484619/149‑1 (MQ=255)                                      cccTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGCACGTCCCTGCCGGTGTaa                                                                                             >  2:536149/1‑149 (MQ=255)                                      cccTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTaa                                                                                             >  1:592098/1‑149 (MQ=255)                                       ccTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAAt                                                                                            >  2:170977/1‑149 (MQ=255)                                         gtgTAGGGTTGCTTATTATAGTTCTCCACTATTGTCAGTTCAGGGTGAGGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGa                                                                                          <  1:159432/147‑1 (MQ=255)                                           gTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATAAGGCAGTGATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTCATGAcc                                                                                        >  1:139717/1‑148 (MQ=37)                                            tAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGACAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAg                                                                                       >  2:623523/1‑149 (MQ=255)                                                 ttGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGGCGGTGTACAGACAGACGCt                                                                                  >  2:457880/1‑149 (MQ=255)                                                 ttGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGCGTAATGCCCGACGCt                                                                                  >  2:157658/1‑149 (MQ=255)                                                 ttGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCCGTGTAATGACAGACGCt                                                                                  >  1:65587/1‑149 (MQ=255)                                                         tACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAGTGACAGCCTCTAGCGTggc                                                                          >  2:469663/1‑147 (MQ=255)                                                              ttCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATTAGGAACGTGCCTTCCGGTGTAATGACAGACGCTACCGTGCCACTgg                                                                     >  2:608490/1‑149 (MQ=25)                                                                 tcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTc                                                                                                          >  2:554981/1‑109 (MQ=255)                                                                 tcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTc                                                                                                          <  1:554981/109‑1 (MQ=255)                                                                 tcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTAt                                                                  >  2:258869/1‑149 (MQ=255)                                                                   cACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGGGTAATGACAGACGCTACCCTGCCACTGGTAttt                                                                >  2:368500/1‑149 (MQ=255)                                                                           tCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAAACGCTACCGTGCCACTGGTAATTGCCGGtca                                                        >  1:639025/1‑147 (MQ=255)                                                                            cAGTTCAGGGTGAGGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCa                                                       <  1:608490/149‑1 (MQ=255)                                                                                 cAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGCCAGAAACCACCCTGCCACTGAGAGTTGCCGGTGTAGTAGc                                                  >  2:149588/1‑149 (MQ=255)                                                                                      tGAGGCTCAGCAGGCAGGTATAGGGCCCATTGCGGTTCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAg                                             <  2:27729/149‑1 (MQ=255)                                                                                       gATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGa                                            <  1:582955/149‑1 (MQ=255)                                                                                        aTGCTCATCAGGCAGGTATAGGGCCCATTGCGGTACTGGCGGCGTTAGGCGTATACCGCTAGGGCTCCGGTGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCAACTGGTATTTCCCGGTGCAGTCGCGAAAGAt                                           >  2:288104/1‑149 (MQ=21)                                                                                          gCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTc                                         >  2:316989/1‑149 (MQ=255)                                                                                           cTCATCAGGCATGTATATGGCCTATTGCGGTTCTGGCGGTGTTCGTGGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTATCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCg                                        <  1:11672/149‑1 (MQ=18)                                                                                              atcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGCCAGACGCTACCGTGCCACTGGTGTTTGCCGGTGCCGTAGCGAAAGATTCGGGa                                     >  1:189333/1‑149 (MQ=255)                                                                                                  ggcaggGATAGGTCCCATTGCTGTGCTCGCGGCGTACGGCATATACCGCCAGGACTCCGGAGATACCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCa                                 <  1:149588/149‑1 (MQ=17)                                                                                                   gcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGGGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGCTGCAGTAGCGAAAGATCCAGGACTcgg                                >  2:144699/1‑147 (MQ=18)                                                                                                   gcaggTATAGGGCCACTATCGGTCCTGACGGCGTTCGGCATATGCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCGGTAGCGAAAGATTCGGGACTCAg                                >  2:246814/1‑149 (MQ=17)                                                                                                     aggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTACTAACAGACGTGACCGTGCCACTGGTCTTTGCCGGTGCAGTAGCGACAGATTCGGGATCCAGGa                              >  1:166135/1‑149 (MQ=255)                                                                                                        tataGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATCCCGCTAGGATTCGCGCGATATCCGGCACGGCCCTGCCGGCCTGATCCCAGACGCTACCGTGCCACACGTATTTTCCAGTGCAGTATCGAAAGATTCGGGGCTCAGGCTGc                           >  2:527514/1‑149 (MQ=255)                                                                                                        tataGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACGGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGc                           >  2:502438/1‑149 (MQ=255)                                                                                                                ccATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGta                   >  1:12676/1‑149 (MQ=255)                                                                                                                    tGCGGTCCTGGCGGGTTGCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATc               <  1:362088/149‑1 (MQ=255)                                                                                                                       ggTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCt                          >  2:330424/1‑135 (MQ=255)                                                                                                                       ggTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCt                          <  1:330424/135‑1 (MQ=255)                                                                                                                        gTCCTGTCGGCGTTCGGCTGAGACCGCTAGTGTTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCa           <  2:450775/149‑1 (MQ=17)                                                                                                                            tGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGc       >  2:406070/1‑149 (MQ=255)                                                                                                                                 gcgTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGc  >  1:230977/1‑149 (MQ=255)                                                                                                                                 gcgTTCGGCGTATACCGCTAGGACTCCGGCGATATCAGGAGCGGCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGGATATCGTCATGTCCCTGc  >  1:353515/1‑149 (MQ=17)                                                                                                                                                   |                                                                                                                                   CGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGC  >  NC_000913_3_pae_tpiA/263952‑264228

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