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

GGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGAT‑GCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAA  >  NC_000913_3_pae_tpiA/263954‑264240                                                                                                                                                 |                                                                                                                                                ggCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGAGGCTCATCAGGCAGGTATAGGTCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGa                                                                                                                                             <  1:172212/149‑1 (MQ=255)    tCCCCGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCGCTATGGTCAGTTCAGGGTGATGGTCAGCCGGCAGGTATAGGGCCCACTTCGGTCCTCGCGGCGTTCGGCGTATCCCGCTAGGACTc                                                                                                                                          >  2:129781/1‑149 (MQ=255)       cagactcAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCgg                                                                                                                                       >  1:738387/7‑149 (MQ=255)       cAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTCTTGTCAGTTCAGGGTGATGGTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCgg                                                                                                                                       >  1:295689/1‑149 (MQ=255)              aGTGTAATCCCCCGCCGCAGGCCCCTGGTAGGGTTGCTCATTACAGTGCTCCACGATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTAGAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATc                                                                                                                                <  1:465076/149‑1 (MQ=255)                  tAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGa                                                                                                                            >  2:132584/1‑149 (MQ=255)                   aaTCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGaa                                                                                                                           <  2:738387/149‑1 (MQ=255)                    aTCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAAc                                                                                                                          >  2:831012/1‑149 (MQ=255)                        cccGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATACGGCCCATTGCGGTCCTGGCGGCGTCCGGCGTATACCGCTAGGACTCCGGCGATATCCGGGACGTcc                                                                                                                      >  2:732314/1‑149 (MQ=255)                                  gCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCAGGAACGTCCCTGCCGgtgt                                                                                                            >  1:18006/1‑149 (MQ=255)                                    cccTGGTAGGGTTGCTCATTACAGTTGTCCACTATTGTCAGTTAATGACGGTGGTCAGCATGCAGGTATAAGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGCTCTCAGCACCGTCCCTGCCGGTGTaa                                                                                                          >  2:700428/1‑149 (MQ=255)                                        ggTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCGGGGTGATGCTCATCAGGCCGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGAc                                                                                                      >  1:939594/1‑149 (MQ=255)                                               ttGCTCATTACAGTTCTCCACTATTGTCAGTTAAGGGTGATGCACATCAGGCAGGGATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCGTGCCGGAGTAATgacggacgat                                                                                               >  1:658657/1‑147 (MQ=255)                                                 gCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTAc                                                                                             >  1:965221/1‑149 (MQ=255)                                                     aTTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACGGGCGCTACCGTg                                                                                         >  1:712236/1‑149 (MQ=255)                                                     aTTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTg                                                                                         >  1:722884/1‑149 (MQ=255)                                                       tACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGGAATGACACACGCTACCGTAcc                                                                                       >  2:47125/1‑149 (MQ=255)                                                         cAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCg                                                                                                                <  1:757786/122‑1 (MQ=255)                                                         cAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCg                                                                                                                >  2:757786/1‑122 (MQ=255)                                                             tccccACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGt                                                                                 <  1:235491/146‑1 (MQ=255)                                                               tcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTAt                                                                               >  2:634024/1‑149 (MQ=255)                                                                  aCTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTGCCGTGCCACTgcgctttg                                                                            >  2:106667/1‑142 (MQ=255)                                                                     aTTGCCAGATAAGGGCGATGCTCAGCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCCAGGACTCCGGCGAGATCCGGAACGTCCCTGCCGGTGTAATTACAGACGCTACCGTGCCACTGGTATTTGCCg                                                                         <  2:26366/149‑1 (MQ=21)                                                                      ttGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCCCTGGTATTTGCCgg                                                                        >  1:525893/1‑149 (MQ=255)                                                                             ttCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTa                                                                 >  2:456901/1‑149 (MQ=255)                                                                             ttCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTa                                                                 >  1:382529/1‑149 (MQ=255)                                                                              tCAGGGTGATGCTCAGCAGGCAGGTATAGGCCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTACGACTCCGGCGATACCCGGAACGTCCATGCCAGTTTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAg                                                                <  1:354623/149‑1 (MQ=21)                                                                                 gggTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGa                                                             >  1:531351/1‑149 (MQ=255)                                                                                    tGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGaaac                                                          >  2:548233/1‑148 (MQ=255)                                                                                    tGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTAGTATTTGCACGTGCAGTAGCGACAg                                                          >  2:405834/1‑149 (MQ=21)                                                                                        gCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAg                                                                                                                                               <  2:399382/60‑1 (MQ=255)                                                                                                 gcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATGCCGCTAGGACTCCGGCGATATCCGGAACGTGCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTGGCGAAAGATTCGGGACTCAg                                             >  2:387994/1‑149 (MQ=255)                                                                                                 gcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATGTCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAACATTCGGGACTCAg                                             >  1:436200/1‑149 (MQ=255)                                                                                                 gcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAg                                             <  1:166077/149‑1 (MQ=255)                                                                                                     gTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGact                                         >  2:789589/1‑147 (MQ=255)                                                                                                          gggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACCCAGGAT‑GCTGGc                                   >  1:582452/1‑149 (MQ=255)                                                                                                           ggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGAT‑GCTGGCg                                  <  1:156684/149‑1 (MQ=255)                                                                                                            gCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGAT‑GCTGGCgg                                 <  1:831012/149‑1 (MQ=255)                                                                                                             ccctgtGCGGTCCCGGCGGCTCTCGGCGTATACCGCTAGGCCTCCGTCGATAGCCGGAACGTCCCTGCCGGTGAAATGACAGACGCAACCGTGCCACTGGTACTTGCCGGTGCAGTAGCGAAAGATTCGGGACGCAGGAT‑GCTGGCGGt                                <  2:344435/144‑1 (MQ=17)                                                                                                              ccATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGAT‑GCTGGCGGta                               >  1:835475/1‑149 (MQ=255)                                                                                                                      gTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGAT‑GCTGGCGGTATATCGTCa                       <  1:563743/149‑1 (MQ=255)                                                                                                                         cTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGAT‑GCTGGCGGTATATCGTCATTg                    <  1:456901/149‑1 (MQ=255)                                                                                                                           ggcggcGTTCGGCTTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGAT‑GCTGGCGGTATATCGTCATTGcc                  <  2:620632/149‑1 (MQ=255)                                                                                                                               gcgGTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGAT‑GCTGGCGGTATATCGTCATTGCCCTGc              <  2:334713/149‑1 (MQ=255)                                                                                                                                  gtCGGTGCATACCGCTATGCCTCCGGCGATTTCAGGAACGTCCCTGCCGGTGTAAGGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGAT‑GCTGGCGGTATATCGTCATTGCCCTGCGtt           <  2:367272/148‑1 (MQ=17)                                                                                                                                   tCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGAT‑GCTGGCGGTATATCGTCATTGCCCTGCGttt          <  1:634024/149‑1 (MQ=255)                                                                                                                                    cGGCGTATACCGCTAGGACGCCGGAGATATACGGAGCCTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGAT‑GCTGGCGGTATATCGTCATTGCCCTGCGTTTa         <  1:289166/149‑1 (MQ=14)                                                                                                                                      gCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGaaa                                                           <  1:614858/98‑1 (MQ=255)                                                                                                                                      gCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGaaa                                                           >  2:614858/1‑98 (MQ=255)                                                                                                                                        gTATACCGCTAGGACTGCGGGCATAGCCGGAACGTCGCTGCCGGTGTAATGACGGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATCGCGCGGGGT‑GAGGGTCATTGCCCTGCGttaattac     <  2:755855/149‑7 (MQ=14)                                                                                                                                           taCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAg                                             >  2:985187/1‑107 (MQ=255)                                                                                                                                           taCCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAg                                             <  1:985187/107‑1 (MQ=255)                                                                                                                                           caCCGCCAGGACTCCGTCGATATCCGGAACGTCCCTGCCGGTGCAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGAT‑GCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATaa  <  2:815458/148‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                GGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGAT‑GCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAA  >  NC_000913_3_pae_tpiA/263954‑264240

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