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.2%	158.5 / NA	55	L50L (CTC→CTA)	ykfH	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (29/24);  total (30/25)
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

TCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGA  >  NC_000913_3_vch_tpiA/263957‑264237                                                                                                                                              |                                                                                                                                             tCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTc                                                                                                                                      <  2:97466/149‑1 (MQ=255)     aGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCTCTATTTTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATCCCGCTAGGA‑CTCCGGc                                                                                                                                  >  1:237268/1‑149 (MQ=255)     aGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGc                                                                                                                                  >  1:496563/1‑149 (MQ=255)                        gCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATACTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATTTCCGGAACGTCCCTg                                                                                                               >  1:235077/1‑149 (MQ=255)                           gCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCg                                                                                                            <  1:84381/149‑1 (MQ=255)                             aGGCCCCTGGTAGGTGTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGgt                                                                                                          <  1:465576/149‑1 (MQ=255)                             aGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGCCAGTTCAGGGGAATGCACATAATGCAGGGATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGACGATATCCGGAACGTCCCTGCCGgt                                                                                                          >  2:48411/1‑149 (MQ=255)                               gCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGgtgt                                                                                                        >  1:292265/1‑149 (MQ=255)                                ccgcTGGTAGGGTTGCTCATTACAGTGTTCCATGAGTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGCCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTa                                                                                                       <  1:434888/146‑1 (MQ=255)                                        aGGGTTGCTCATTACTGTTCTCCACTATTGTCAGGTTAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGa                                                                                               <  2:619723/149‑1 (MQ=255)                                           gTTGCTTATTACAGTTCTCGCCTATGGTCAGTTCAGGGTGATGCTCATCAGGCCGGTATAGGGCCCATTGCGGTCCTGGCGGCGTGCGGCGGATACCGCTAGGA‑CTGGGGAGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGc                                                                                            <  2:116396/149‑1 (MQ=21)                                              gCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTAc                                                                                         >  1:833365/1‑149 (MQ=255)                                               cTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACAGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGGGAAATGACCGACGCTTcc                                                                                        >  2:435733/1‑149 (MQ=255)                                                   ttACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGc                                                                                    >  1:776511/1‑149 (MQ=255)                                                   ttACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGc                                                                                    >  1:548352/1‑149 (MQ=255)                                                   ttACAGCTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGc                                                                                    <  1:538047/149‑1 (MQ=255)                                                     aCAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCa                                                                                  <  2:754119/149‑1 (MQ=255)                                                       aGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACt                                                                                <  1:420785/149‑1 (MQ=255)                                                         ttCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTgg                                                                              <  2:292265/149‑1 (MQ=255)                                                            tcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTAt                                                                           >  2:556410/1‑149 (MQ=255)                                                                cTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGc                                                                       <  1:137395/149‑1 (MQ=255)                                                                  aTTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCg                                                                     <  2:496563/149‑1 (MQ=255)                                                                   ttGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGCCGCTACCGTGCCACTGGTATTTGCCgg                                                                    >  2:315523/1‑149 (MQ=255)                                                                       cAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCa                                                                <  1:686967/149‑1 (MQ=255)                                                                        aGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCAATTGCGGTCCTGGCGGCGTTCTGCGTCTACCGCTAGGC‑GTCCGCCGATATCTGGCACGTCCCTGCCGGTGGCATGACGGACGCCAGCGCGCCGCGGCTGTTCGCCgggcagg                                                               >  2:717989/1‑144 (MQ=255)                                                                          ttCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATATGCCGGTGCAGTa                                                             >  1:809953/1‑149 (MQ=255)                                                                                gTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGaaa                                                       >  1:305633/1‑149 (MQ=255)                                                                                gTGATGCTCATCAGGCAGGTATAGGGCCCATTGCAGTACGGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGaaa                                                       >  1:617985/1‑149 (MQ=255)                                                                                 tGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAg                                                      <  1:435475/149‑1 (MQ=255)                                                                                   aTGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGata                                                                                                                              >  2:431942/1‑75 (MQ=255)                                                                                   aTGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGata                                                                                                                              <  1:431942/75‑1 (MQ=255)                                                                                     gCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGATAGGG‑CGCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTc                                                  >  1:197182/1‑149 (MQ=255)                                                                                       tcatcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGGACGTCCCTGCCGGGGTAATGACAGACGCTACCGTGCCACCGGTATGTGCCGGTGCAGTAGCGAAAGATTCgg                                                >  2:18923/1‑149 (MQ=21)                                                                                       tcatcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCgg                                                >  2:30944/1‑149 (MQ=255)                                                                                         atcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTGGGA‑CTCCGGCGATGTCGGGAACGTCCCTGCCGTTGTAATGACAGGCACTACCGTGCCACACGTATTTGCCGGTGCAGTAGCGCACGCTTCgggc                                              >  2:308541/1‑148 (MQ=255)                                                                                               caggTATAGGGCCCAGTGCGTGCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTGCGTCGATAGCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAgg                                        <  1:587837/149‑1 (MQ=18)                                                                                                aggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGa                                       >  2:454860/1‑149 (MQ=255)                                                                                                aggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGa                                       <  1:556410/149‑1 (MQ=255)                                                                                                        ggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCg                               >  2:591009/1‑149 (MQ=255)                                                                                                        ggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCg                               <  1:591009/149‑1 (MQ=255)                                                                                                        ggCCCATTGCGGTCCTGGCAGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTGATGACAGACGCTACCGTACCACTGGTCTTTGCCGGTGCAGTAGCCAACGATTCGGGCATCAGGCTCCTGGCg                               >  2:703156/1‑149 (MQ=255)                                                                                                         gCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGAGTCAGGGTGCTGGCgg                              >  2:150251/1‑149 (MQ=255)                                                                                                           ccATTGCGGTCCTGGCGGCGTTCGGCGTAAAAAG‑GTGGATCTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGAATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGta                            >  2:179825/1‑149 (MQ=17)                                                                                                            cATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGtat                           <  1:454860/149‑1 (MQ=255)                                                                                                                  ggTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACAGTACGGTTTTGAGGTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTc                     >  1:817726/1‑149 (MQ=17)                                                                                                                   gTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGaa                                                        >  1:758351/1‑113 (MQ=255)                                                                                                                   gTCCTGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGaa                                                        <  2:758351/113‑1 (MQ=255)                                                                                                                    tCCTGGCGGCGTTCGGCTTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTAAGGATGCTGGCGGTATATCGTCAt                   <  1:450469/149‑1 (MQ=255)                                                                                                                       tGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAGCGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGa                                              <  2:412392/119‑1 (MQ=255)                                                                                                                       tGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAGCGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGa                                              >  1:412392/1‑119 (MQ=255)                                                                                                                       tGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGc                <  2:238817/149‑1 (MQ=255)                                                                                                                       tGGCGGCGTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCACGATGCTGCCGGTATATCGTCATTGc                >  2:564442/1‑149 (MQ=255)                                                                                                                              gTTCGGCGTATACCGCTAGGA‑CTCCGGCGATATCCGGCACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGt         >  1:408518/1‑149 (MQ=255)                                                                                                                                  cgCGTATACCGCCAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTAc     <  2:237268/148‑1 (MQ=255)                                                                                                                                     gTATCCCGCTAGGA‑CTCCGGCGATATCCGGAACGTCGCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGAATTTGACGGCCCAGTAGTGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGa  <  1:452883/149‑1 (MQ=14)                                                                                                                                              |                                                                                                                                             TCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGA‑CTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGA  >  NC_000913_3_vch_tpiA/263957‑264237

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