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

CAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAATT  >  NC_000913_3_vch_tpiA/263960‑264242                                                                                                                                           |                                                                                                                                                 cAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCgg                                                                                                                                        >  2:866707/1‑149 (MQ=255)  aGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGc                                                                                                                                       >  1:676368/1‑149 (MQ=255)      gcAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATGGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGTCAGGGCTCCGGCGata                                                                                                                                   >  2:590429/1‑149 (MQ=255)          tgtAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGGATACCGCTAGGACTCCGGCGATATCCg                                                                                                                               >  1:605731/1‑149 (MQ=255)             aaTCCCCCGACGCAGCCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGaa                                                                                                                            <  1:661253/149‑1 (MQ=38)              aTCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAAc                                                                                                                           >  1:16132/1‑149 (MQ=255)              aTCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAAc                                                                                                                           >  1:643088/1‑149 (MQ=255)                    cGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCt                                                                                                                     >  2:771310/1‑149 (MQ=255)                         cAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCgg                                                                                                                <  1:261104/149‑1 (MQ=255)                                 tGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGa                                                                                                        >  1:646464/1‑149 (MQ=255)                                  ggTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGAc                                                                                                       >  2:93762/1‑149 (MQ=255)                                   gcaGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCCGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACa                                                                                                      <  1:748471/147‑1 (MQ=255)                                           gCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGTGCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCTCTAGGACTCCGGCGATATCCGGAGCGTCCCTGCAGGTGCAAGGGCAGACGCtcc                                                                                              >  1:215413/1‑147 (MQ=25)                                                ttACAGTTCTCCTATATTGTCAGTTTATGGTGATGCGTATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGc                                                                                         <  1:486473/149‑1 (MQ=21)                                                  aCAGTTCTCCACTATTGTCAGTTCAGGGTGATGCACATCAGGCAGGTATAGGGCCTATTGCGGCCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCa                                                                                       <  1:929853/149‑1 (MQ=255)                                                  aCAGTTCTCCACTAGTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCa                                                                                       <  1:71307/149‑1 (MQ=255)                                                     gTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTg                                                                                    <  1:771310/149‑1 (MQ=255)                                                      ttCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCg                                                                                                                                         >  2:896188/1‑95 (MQ=255)                                                      ttCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCg                                                                                                                                         <  1:896188/95‑1 (MQ=255)                                                           cACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATTACAGACGCTACCGTGCCACTGGTAttt                                                                              <  2:646464/149‑1 (MQ=255)                                                             cgaTTGTAAGTTCTGGTGGATGCTAATCAGGCAACTATGGCACAGAGTGCGGTCTGGGCGGTGTGCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGc                                                                            <  2:893654/147‑1 (MQ=17)                                                             cTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGc                                                                            >  2:1037451/1‑149 (MQ=255)                                                                   tCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCTTGCCACTGGTAGTTGCCGGTGc                                                                      >  2:444163/1‑149 (MQ=255)                                                                   tCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTCCCGGTGc                                                                      >  1:656570/1‑149 (MQ=255)                                                                   tCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAAACGCTACCGTGCCACTGGTATTTGCCGGTGc                                                                      >  1:1000358/1‑149 (MQ=255)                                                                       ttCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTa                                                                  >  1:423186/1‑149 (MQ=255)                                                                         cAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGc                                                                <  1:93762/149‑1 (MQ=255)                                                                            ggTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGaa                                                             >  1:349977/1‑149 (MQ=255)                                                                                  gCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTgg                                                                                   <  1:452820/121‑1 (MQ=255)                                                                                  gCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTgg                                                                                   >  2:452820/1‑121 (MQ=255)                                                                                  gCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGAAAGCCGCTACCGTGCCACTGGTATTTGCCGGTGCCGTAAGGaaaacttc                                                       >  1:176759/1‑144 (MQ=255)                                                                                    tcatcaTGCAGGTATAGGGCCCATTACGGTCCTGGCGGCGTTCGGCGTACACCGCTAGGACTCCTGCGATATCCGGAACGTAACTGCCGGTGTAATGACAGATGCTACCGTGCCACTGGTATTCGCCGGTGCAGTAGCGAAAGATTCgg                                                     <  1:944317/149‑1 (MQ=18)                                                                                           gcaggTATGGGGCCCATTGCGGTGGTGGCAGCGTTCGGCGTATACCGCCAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAg                                              <  2:155309/149‑1 (MQ=18)                                                                                               gTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTTCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCCGTAGCGCAAGATTCGGGATTCAGGATg                                          >  2:755044/1‑149 (MQ=18)                                                                                                  taGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTc                                                       >  1:343840/1‑133 (MQ=255)                                                                                                  taGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTc                                                       <  2:343840/133‑1 (MQ=255)                                                                                                    gggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGc                                     >  1:861020/1‑149 (MQ=255)                                                                                                    gggCCCATTGCGGTCCTGGCGGCGCTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGc                                     >  2:891595/1‑149 (MQ=255)                                                                                                     ggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGCCTCAGGATGCTGGCg                                    >  1:696141/1‑149 (MQ=255)                                                                                                        ccATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTATGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGta                                 <  1:590429/149‑1 (MQ=255)                                                                                                        ccATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGta                                 <  2:861020/149‑1 (MQ=255)                                                                                                         ctaggCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGtat                                <  1:60730/145‑1 (MQ=255)                                                                                                         cactGCGGGCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGtat                                >  2:130774/4‑149 (MQ=255)                                                                                                          aTTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTGCCTGCCGGTGTAATGACAGTCGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGtata                               >  2:792591/1‑149 (MQ=255)                                                                                                            tGCGGTCCTGGCGGCGTTGGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATc                             <  1:628759/149‑1 (MQ=255)                                                                                                            tGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATc                             <  2:676368/149‑1 (MQ=255)                                                                                                                gTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCg                                                                          <  1:735293/100‑1 (MQ=255)                                                                                                                gTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCg                                                                          >  2:735293/1‑100 (MQ=255)                                                                                                                 tCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCAGGAGATATGTGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGGCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATAGCGTcgt                        >  1:516735/1‑147 (MQ=17)                                                                                                                      gcggcgTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGccc                   <  1:124368/149‑1 (MQ=255)                                                                                                                              cGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTTGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGtttc           >  2:21286/1‑148 (MQ=255)                                                                                                                                       ccGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTGCCGTGCCACTGGTATTTGCCGGTGCAGTCGCGAAAGATTCGGGACTCAGGATGCCGGCGGTAGATCGTCATGGCGCTGCGTTTGCTGATAAtt  >  2:17270/1‑149 (MQ=14)                                                                                                                                           |                                                                                                                                                 CAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAATT  >  NC_000913_3_vch_tpiA/263960‑264242

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑