Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_vch_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_vch_tpiA	264,098	0	G	T	98.0%	146.7 / NA	52	L50L (CTC→CTA)	ykfH	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (28/22);  total (29/23)
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

GGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGT‑GCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAATTT  >  NC_000913_3_vch_tpiA/263954‑264243                                                                                                                                                 |                                                                                                                                                   ggCTCCCAGGCGCAGTGTCATCTCCCTCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGa                                                                                                                                                <  1:141777/149‑1 (MQ=255)   cTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTACTCCACTTTTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACt                                                                                                                                              >  1:111907/1‑149 (MQ=255)               gtgtAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATcc                                                                                                                                  <  2:136287/149‑1 (MQ=255)                tgtAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCg                                                                                                                                 <  2:502369/149‑1 (MQ=255)                tgtAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACAATTGTAAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGTCTCCGGCGATATCCg                                                                                                                                 >  2:736788/1‑149 (MQ=255)                  tAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGGCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTACTGTCGGAGTTCGGCGTATACCGCTGGGACTCAGGCGGTGTCCGGa                                                                                                                               >  1:153160/1‑149 (MQ=255)                  tAATCTCCCGCTGCAGGCCACTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCGGGGTGATGCTCATCACGCAGGTATAGGGCCCATTGCGGTCGTGGCGGCGTTCGGCCTACACCGCACAGCCTACGGCCATATCCgca                                                                                                                               >  2:662829/1‑147 (MQ=255)                          cGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGGGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAAAGTCCCt                                                                                                                       >  1:10936/1‑149 (MQ=255)                                  gCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTGCCTGCCGgtgt                                                                                                               >  1:226569/1‑149 (MQ=255)                                   ccccTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTa                                                                                                              <  1:638499/149‑1 (MQ=255)                                    cccTGGTAGGGTTGCTCATTACAGTTCTCCACTCTTGTTAGTGCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTaa                                                                                                             >  1:775738/1‑149 (MQ=255)                                         gTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTGCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACa                                                                                                        >  1:614953/1‑149 (MQ=255)                                            gggTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGAc                                                                                                     <  2:406617/149‑1 (MQ=255)                                             ggTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTCTAGGGCCCATTGCGGTCCTGGCGGCGTACGGCGTATACCGCTAGGACTCCGGCGATATCCGGGACGTCCCCGCCGGGGTCATGACAGCCg                                                                                                    >  1:261546/1‑149 (MQ=37)                                             ggTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCg                                                                                                                   <  2:144723/134‑1 (MQ=255)                                             ggTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCg                                                                                                                   >  1:144723/1‑134 (MQ=255)                                               ttGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCt                                                                                                  >  1:558045/1‑149 (MQ=255)                                               ttGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGAGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCt                                                                                                  >  1:598154/1‑149 (MQ=255)                                                  cTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAAGGACAGACGCTAcc                                                                                               >  1:598140/1‑149 (MQ=255)                                                     aTTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGACCGTCCCTGCCGGTGTAATGACAGACGCTGCCGT‑g                                                                                           >  1:442657/1‑149 (MQ=255)                                                      ttACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGAATAGGGCCCATTATGATCATGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCAGGGGTAATTACAGACGCTACCGT‑Gc                                                                                          >  1:772189/1‑149 (MQ=255)                                                             tctcCCCTATTGTAAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGT‑GCCACTGGt                                                                                   <  1:462545/149‑1 (MQ=255)                                                                ccACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGT‑GCCCCTGGTAtt                                                                                >  2:54913/1‑149 (MQ=255)                                                                 cACTATTGTCAGTTCAGGGTGATGGTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGCCTCCGGCGATATCCGGAACGTCCCTGCCGGGGTAATGACAGACGCTACCGT‑GCCACTGGTAttt                                                                               <  1:714617/149‑1 (MQ=255)                                                                 cACTATTGTCAGTTCAGGGTGATGCTCATCAGACAGGTATACGGCCCATTGCGGTCCTGGCGGCGATCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGT‑GCCACTGGTAttt                                                                               >  2:744547/1‑149 (MQ=255)                                                                  aCTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGT‑GCCACTGGTATTTg                                                                              >  2:511498/1‑149 (MQ=255)                                                                   cTATTGTCAGGTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGT‑GCCACTGGTATTTGc                                                                             <  2:636233/149‑1 (MQ=255)                                                                     tttGTCAGTTCAGGGTGATGGTCCTCAGGCAGGTATGGGGCCCATGGTGGTCCTGGTGGCGGGCGGCGTCTACAGCTAGGTGGCGGGTGATCTGCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGT‑GCCACTGGTATTTGCCg                                                                           <  2:715393/148‑1 (MQ=18)                                                                     aTTGTCAGTTCAGGGTGATGCTCATCAGGAAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCTGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGT‑GCCACTGGTATTTGCCg                                                                           <  1:139322/149‑1 (MQ=255)                                                                      ttGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCCTTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGGTATCCGGAACGTCCCTGCCGGTGTAATGCCAGACCCTACCGT‑GCCACTGGTATTTTCCgg                                                                          >  2:524071/1‑149 (MQ=255)                                                                        gTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGT‑GCCACTGGTATTTGCCGGTg                                                                        >  2:424920/1‑149 (MQ=255)                                                                         tCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGT‑GCCACTGGTATTTGCCGGTGc                                                                       >  1:4771/1‑149 (MQ=255)                                                                           aGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCGTGGCGGCGTTCGGCGTATACCGCTCGGGCTCGGGCGAAGTCCGGGAAGTCCCTGCGGGTGTGATGAAAGACGCTTGAGT‑GCAACTGGTATTTGCCGGTGAAg                                                                     >  2:274420/1‑149 (MQ=255)                                                                             ttCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGT‑GCCACTGGTAt                                                                                 >  2:238655/1‑135 (MQ=255)                                                                             ttCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGT‑GCCACTGGTAt                                                                                 <  1:238655/135‑1 (MQ=255)                                                                             ttCAGGGGGATGCGCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGT‑GCCACTGGTATTTGCCGGTGCAGTa                                                                   <  1:142456/149‑1 (MQ=255)                                                                             ttCAGGGAGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTGCCTGCCGGTGGAATGACAGACGCTACCGTGGCCAC‑GGTTTTTGACGGTGCGGtt                                                                   >  1:210090/1‑148 (MQ=255)                                                                                        gccccgcaggcaggTATAGGGCCCATTGCGGCACTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCTGGTGTAATGACAGACGCTACCGT‑GCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTc                                                        <  2:111907/143‑1 (MQ=21)                                                                                        gcccatcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGT‑GCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTc                                                        <  1:567515/146‑1 (MQ=255)                                                                                          tcatcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGT‑GCCACTGGTATTTGc                                                                             <  1:3163/126‑1 (MQ=255)                                                                                          tcatcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGT‑GCCACTGGTATTTGc                                                                             >  2:3163/1‑126 (MQ=255)                                                                                          gcaggaggctggTATAGGGCCCTTGTGGGTGCTGGGGGCGTTCGGCGTATATCGCGAGGACTCCGGCGATATCCGGCACGTCCCTGCCGGTGTAATGATAGACGCTACCGT‑GCCACTGGTATTTGGCGGTGCAGTAGCGAAAGATTCgg                                                      <  2:639519/144‑1 (MQ=17)                                                                                             gcaggcaggTATAGGCCCCATTGCGGTCCGAGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGT‑GCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGAc                                                   <  2:321871/148‑1 (MQ=255)                                                                                                      tataGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGCAACTTTTCTGCCGGTGTAATCACAGACGCTACCGT‑GCCACTGGTATTTGCCGGTGCAGTAGTGGAAGATTCGGGACGCAGGAGGc                                          >  1:381598/1‑149 (MQ=255)                                                                                                          gggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAAGGTCCCTGCCGGTGTAATGACAGACGCTACCGT‑GCCACTGGTATTTACCGGTGCAGTAGCGAAAGATTCGGGACGCAGGATGCTGGc                                      >  1:20696/1‑149 (MQ=255)                                                                                                            gCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGT‑GCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCgg                                    <  1:511498/149‑1 (MQ=255)                                                                                                            gCCCATTGCGGTCCTGGCGGCGTAAGACGTATACAGCTCACCCTCCGGCAGTATCCGGAACGTCCCTGCCGGTGTGAGGACAGACGCTACCGT‑GCCCCTGGTATTTGCAGGGGCAGTCGCGACAGCTTCGGCACTCCGGATCCTGGCgg                                    >  2:722238/1‑149 (MQ=255)                                                                                                             cccATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGT‑GCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGCCTCAGGATGCTGGCGGt                                   >  1:93077/1‑149 (MQ=255)                                                                                                                   gCGGTCCGGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCAGGTGTAATGACAGACGCTACCGT‑GCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCg                             <  2:696186/149‑1 (MQ=255)                                                                                                                        cgtGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGT‑GCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCAtt                        <  1:362036/147‑1 (MQ=255)                                                                                                                                        gTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGATACCGT‑GCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGAGTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGa        <  2:614953/149‑1 (MQ=255)                                                                                                                                          ataCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGT‑GCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATa      <  1:424920/149‑1 (MQ=255)                                                                                                                                              cgctAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGT‑GCCACTGGTATTGGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTCCGTTTACTGATAAttt  >  2:611441/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                   GGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGT‑GCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAATTT  >  NC_000913_3_vch_tpiA/263954‑264243

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