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

CGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACT‑GGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAATT  >  NC_000913_3_pae_tpiA/263952‑264242                                                                                                                                                   |                                                                                                                                                  ggggCTCCCGGGCGCAGTGTGTTTGGCCGCTGCAGGCTCCTGGTAGGGTTGCTCATTACAGTCCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGGCCTTTGCCATCCTTGCGGCGTTTGGCGTATACCGCTAg                                                                                                                                                 <  2:179022/148‑1 (MQ=255) cGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAg                                                                                                                                                 <  1:233098/149‑1 (MQ=255)  gggCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCTGCGTATACCGCTAgg                                                                                                                                                >  2:10419/1‑149 (MQ=255)    gCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGGGTATACCACTAGGAc                                                                                                                                              >  2:667871/1‑149 (MQ=255)    gCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGAc                                                                                                                                              <  1:456597/149‑1 (MQ=255)        ccAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCCAGGCCTCCg                                                                                                                                          >  2:225009/1‑149 (MQ=255)        ccAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGGTGCTCATTACAGTTCTCCACTATTGTCAGTTCCGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTTCg                                                                                                                                          >  2:225008/1‑149 (MQ=255)         gaGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGAGCACAGTTCGTTCCTGGCGTTCTTCGGCGTATACCGCTAGGACTCCgg                                                                                                                                         <  1:369733/148‑1 (MQ=255)             cgcAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAATTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACTGCTAGGACGCCGGCGat                                                                                                                                     >  1:276996/1‑149 (MQ=255)                attgtAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATc                                                                                                                                  <  1:35266/147‑1 (MQ=255)                     aaTCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGaa                                                                                                                             <  1:360093/149‑1 (MQ=255)                      aTCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGACCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAAc                                                                                                                            >  2:629651/1‑149 (MQ=255)                                     ccccTGGTAGGGTTGCTCATCACAGTTCTCCACTATAGTCAGTTCAGGGGGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGGTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGGGTa                                                                                                             >  2:506167/1‑149 (MQ=255)                                       ccTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAAt                                                                                                           <  2:539122/149‑1 (MQ=255)                                       ccTGGTAGGGTGGCTCATTACAGTTCTCCACTATTGCCCGTACAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAAt                                                                                                           <  1:361155/149‑1 (MQ=255)                                            tAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAg                                                                                                      >  2:128434/1‑149 (MQ=255)                                               ggTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACg                                                                                                   >  2:271400/1‑149 (MQ=255)                                                     tCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCATGTATAGGGCCCATTGCGGTCCTGGGGGCGTTCGGCGTATACCGCTAGGACTCCGGAGATATCCGGAACGTCCCTGCCGGTGTAATGTCCGACCCTACCg                                                                                             >  2:25195/1‑149 (MQ=255)                                                            aGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACt                                                                                      >  2:623039/1‑149 (MQ=255)                                                                 ctcACTATTGTCGGTTCAGGGGGATGCTCATCAGGGAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTTGGCGTATACCGCTAGGACTCCGGCGATATCCGGATCGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACT‑GGTAt                                                                                <  2:337931/147‑1 (MQ=21)                                                                  ccACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACT‑GGTAtt                                                                               >  1:178970/1‑149 (MQ=255)                                                                     cTATTGTCAGTTCAGGGTGATGCTCATCAGCCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACT‑GGTATTTGc                                                                            <  1:626585/149‑1 (MQ=255)                                                                        ttGTCAGTTCAGGGTGATTCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACT‑GGTATTTGCCgg                                                                         <  2:52274/149‑1 (MQ=255)                                                                        ttGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCCCT‑GGTATTTGCCgg                                                                         >  1:128443/1‑149 (MQ=255)                                                                        tggTCAGTTGAGGGTGATGCTCATCAGGCATGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACT‑GGTATTTGCCgg                                                                         <  2:276996/147‑1 (MQ=255)                                                                            cAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGAGATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACT‑GGTATTTGCCGGTGCa                                                                     <  1:163112/149‑1 (MQ=255)                                                                              gTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCACATTGCGGTCCTGGCGGCGGTCGGCGTATACCGCTAGGCCTCCGGCGCTATCCCGAAGATCGCTGCCGCTGGTATGACCCACGATCCCGTG‑CACTCGCTATTTGCCGGtatgac                                                                   >  2:444722/1‑144 (MQ=255)                                                                               ttCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGTTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACT‑GGTATTTGCCGGTGCAGTa                                                                  <  1:670333/149‑1 (MQ=255)                                                                                    ggggATGCTCACCGGGCAGGGGTAGCGCGCATGGCGGTCCTGGCGGCGTTCGTCGTATACCGCTAGGACTCGGGCGATGTCCGCAACGGCCCTGCCGGTGTAATGACGGACGCTACCGTGCCACT‑GGTATTTGCCGGTGCAGTAGCTaa                                                             <  2:146361/146‑1 (MQ=255)                                                                                       gATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCCACCGTGCCACT‑GGTCTTTGCCGGTGCCGTAGCGACAGa                                                          >  1:646007/1‑149 (MQ=25)                                                                                          gCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACT‑GGTATTTGCCGGTGCAGTAGCGAAAGATTc                                                       >  2:530278/1‑149 (MQ=255)                                                                                           gatttcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACT‑GGTATTTGCCGGTGCAGTAGCGAAAGATTCg                                                      <  1:344328/145‑1 (MQ=255)                                                                                              atcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACGGACGCTACCGTGCCACT‑GGTATTTGCCGGTGCGGTAGCGAAAGAttgcgga                                                   >  1:481715/1‑144 (MQ=255)                                                                                                 aggcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACT‑GGTATTTCCCGACGCAGCACCGACAGATTcttctatc                                                >  2:344328/1‑142 (MQ=21)                                                                                                   gcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGAc                                                                                                        <  2:448758/92‑1 (MQ=255)                                                                                                   gcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGAc                                                                                                        >  1:448758/1‑92 (MQ=255)                                                                                                    caggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACG‑GTTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAgg                                             <  1:271400/149‑1 (MQ=255)                                                                                                            gggCCCATTGCGGTCCTGGCGGCGTTCGTCGTATACCGTTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACT‑GGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGc                                     <  1:10419/149‑1 (MQ=255)                                                                                                            gggCACATTGTGGGCCTAGGGGGGTTCGATGGATAGCGATAGGACGCCGGCGATATCAGGAACGGCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACT‑GGTATTGGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGc                                     <  1:315391/149‑1 (MQ=14)                                                                                                              gCCCATTGCGGTCATGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACT‑GGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCgg                                   <  1:662974/149‑1 (MQ=255)                                                                                                                       ggTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACT‑GGTATTTGCCGGTGCCGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTc                          >  1:193315/1‑149 (MQ=255)                                                                                                                        gTCCCGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACT‑GGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCa                         <  1:335288/149‑1 (MQ=255)                                                                                                                          ccTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACT‑GGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGCATATCGTCAtt                       <  1:420264/149‑1 (MQ=255)                                                                                                                                 gcgTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACT‑GGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGc                <  1:477317/149‑1 (MQ=255)                                                                                                                                       ggCGTCTACCGCTAGGACCCCGGCGATGTCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACT‑GGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTAc          <  1:667871/149‑1 (MQ=255)                                                                                                                                       ggCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAAGGACAGACGCTACCGTGCTACT‑GGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTAc          <  1:629651/149‑1 (MQ=255)                                                                                                                                             taCCGCTAGGAGTTCGGCGAGATCCGGAACGTCGCTGCGGGTGTGATGACAGACGCTAGCGTGCCACT‑GGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGAtca    <  2:410828/149‑3 (MQ=14)                                                                                                                                              aCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACT‑GGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAAt   >  1:673206/1‑149 (MQ=255)                                                                                                                                               ccGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACT‑GGTATTTGCCGGTGCAGTAGCGAAAGATTCgg                                                     <  2:234382/98‑1 (MQ=255)                                                                                                                                               ccGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACT‑GGTATTTGCCGGTGCAGTAGCGAAAGATTCgg                                                     >  1:234382/1‑98 (MQ=255)                                                                                                                                               ccGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACT‑GGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCCGGATGCTGGCGGTATATCGTCTTTGCCCGGGCTTTACTGCTAGtt  >  2:243358/1‑149 (MQ=14)                                                                                                                                                   |                                                                                                                                                  CGGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACT‑GGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAATT  >  NC_000913_3_pae_tpiA/263952‑264242

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