Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,560,455	+G	intergenic (‑2/+1)	glpR ← / ← glpR	pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,560,455	1	.	G	98.0%	150.1 / ‑3.8	49	intergenic (‑2/+1)	glpR/glpR	pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (27/21);  total (28/21)
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

GATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACG  >  NC_000913/3560316‑3560590                                                                                                                                              |                                                                                                                                        gATAAACCGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGTTCTTCTACGGCTTTTGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGcgca                                                                                                                                       <  1:1194596/143‑4 (MQ=255) gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGcgca                                                                                                                                       <  1:1016339/143‑4 (MQ=255) gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCAGCACCGCCATGAT‑GGcgca                                                                                                                                       >  2:1078149/1‑140 (MQ=255) gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTAAGGGTTACCTTGTGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGcgca                                                                                                                                       <  1:169190/143‑4 (MQ=255)         gccgcGAGCCATCGGGGATTTGCGCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAgg                                                                                                                               <  1:987490/139‑1 (MQ=255)         gCGTCGGGCCCTGGGGGTCGTGCTCCGCCACTTGGTGGGCGAGGCGCTCTTGTTCTGCGGTCTGGGTGGCCTGGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAgg                                                                                                                               <  2:1039800/143‑1 (MQ=255)                 ccATTGGGGATGTGCTGCGCCACGTTGCGGGCGATGCGCGCTTTTTCTCCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTc                                                                                                                       <  1:475558/143‑1 (MQ=255)                    tgggggATTTGCTGCGGCCCTTTGCGGGCGATGCGCTCTTTTTGTTCGGTCCGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGc                                                                                                                    <  2:175123/141‑1 (MQ=255)                    gtGGGGATGTGGTCGGCCACTTTGCGGGAGAGGGGCGCTGTGTCGTCGGTCTGGGCGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGCGGGCCTCCCCACACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGc                                                                                                                    <  2:258086/142‑1 (MQ=255)                          aTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCATGCCACGACGTGTTAACCGACCTGGAAGGCAGCGCCGCACCGCCATGAT‑GCCGCACGGTCAGGCTTTGCTCCGCCAGCTc                                                                                                              >  2:458369/1‑143 (MQ=255)                             tGCTCCGCCACTTTCCGGGCCATCCGCTCTCTTTCTTCCGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCACGATCATGTTTTGCTCCGCCAGCTCAtt                                                                                                           <  1:6769/143‑1 (MQ=255)                                tCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGACCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGACATGATGGGC‑GAGGCGCAGGTCTTGCTCCGCCCGCACATTGAg                                                                                                        >  2:743025/1‑143 (MQ=255)                                  cGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGCACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATCGAGGt                                                                                                      >  2:645885/1‑143 (MQ=255)                                    ccACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAAGAAGCACCGCACCACCATGAT‑GGCGCAGAATCAGGTTTTGCTCCGCCAGCTCATTGAGCTcg                                                                                                    >  1:491668/1‑143 (MQ=255)                                        tttGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCGTGCGATCGTGCCACGGCGTGGTAAACGAACTGGAAGGCAGCGCAGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCACATCGAAGGAgcggc                                                                                                >  2:401110/1‑143 (MQ=255)                                          tGCGGGCGATGCGCTCTTTGTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGa                                                                                              <  1:645987/143‑1 (MQ=255)                                             gggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGCT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAtc                                                                                           >  1:42970/1‑142 (MQ=255)                                               gCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGt                                                                                         >  1:525294/1‑143 (MQ=255)                                                    gcgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCg                                                                                    >  1:371856/1‑143 (MQ=255)                                                      gcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCggg                                                                                  >  1:180794/1‑143 (MQ=255)                                                      gcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCggg                                                                                  >  1:134711/1‑143 (MQ=255)                                                      gcTCTTGTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCggg                                                                                  <  1:408904/143‑1 (MQ=255)                                                            tttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGCACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCAGTGAGGTCGCGGCGAAGAGTCTGCGGGCTGACg                                                                            >  1:1096528/1‑143 (MQ=255)                                                              gctcgggTGTGGGTGGCGTGGCGATCGTGCCACGGCGTGTTAGCCGAACTGGAGGGCAGCGCCCCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGga                                                                          <  2:748058/138‑1 (MQ=255)                                                               ctgcGGTCTGGGTGGCCTTGCGATCGTGCCACTTCGTGGTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGgag                                                                         <  1:498590/140‑1 (MQ=255)                                                               cgtcGGTCTGGGTGGCCTTGGGATCGTGCCACGGCGGGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGgag                                                                         <  1:1033831/141‑1 (MQ=255)                                                                ttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGATCATAGTTCTCGCGGCGAATAGTCTGCGGGCTGACGgaga                                                                        >  1:813005/1‑143 (MQ=255)                                                                 tcGGTCTGGGGGGCCTTGCGATCGCGCCACGGCGTGTTAACCGAACTGGAAGGCAGTGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCTGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaa                                                                       <  1:1148479/143‑1 (MQ=255)                                                                  cGGGCTGGGGGGCCTAGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaaa                                                                      <  1:1033419/143‑1 (MQ=255)                                                                      cggggTGGCCTTGCGGTCGTGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCGCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                                  <  2:268571/141‑1 (MQ=255)                                                                       tGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTTGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGctc                                                                 <  1:776192/143‑1 (MQ=255)                                                                       tGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCACTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGctc                                                                 <  1:282808/143‑1 (MQ=255)                                                                          gTGGTCTTGCGATCGCGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAc                                                              <  2:690744/143‑1 (MQ=255)                                                                          gTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGTTCTAc                                                              >  1:1126720/1‑143 (MQ=255)                                                                          gTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAc                                                              >  1:1137266/1‑143 (MQ=255)                                                                           tGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGCTCAGGTTTTGCTCCGCCAGCTCATTGAAGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                             >  1:36426/1‑143 (MQ=255)                                                                               cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGCGAAATGCTCTACCAgac                                                         >  2:1156173/1‑141 (MQ=255)                                                                                ttGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCCGCGCCGCACCGCAATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGctc                                                        >  1:750900/1‑143 (MQ=255)                                                                                           ccACGGCTTGTCAGCCGAACTGGCAGGCAGCGCCGCACCGCCATGAT‑GGCGCTGGATCAGGTTTTGCTCCGCAAGCCCAGTGAAGACACGTCTAATCGCCTGCGGGCTGACGGAGACATGCTCTACCAGCTCTTCGGTACTGa                                             <  2:255015/143‑1 (MQ=255)                                                                                                 cGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGAGGGAGAAATGCTGTACCAGCTCTTCGGTAGTGACATAAc                                       >  1:1427/1‑143 (MQ=255)                                                                                                    gtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCct                                    >  2:716949/1‑143 (MQ=255)                                                                                                       aaCCGAACTGGAAGGCAGCGCCGCACCGCCAGGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTGTTCGGTCCTGAGATAACAcggct                                 >  1:693393/1‑143 (MQ=255)                                                                                                         ccGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGt                               >  2:586862/1‑143 (MQ=255)                                                                                                          cGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGCCATAACCCTGCTGtt                              >  1:415529/1‑143 (MQ=255)                                                                                                               tGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAcc                         >  1:149570/1‑143 (MQ=255)                                                                                                               tGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAcc                         >  1:971887/1‑143 (MQ=255)                                                                                                                 gAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAAAAGTGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAg                       >  2:252373/1‑143 (MQ=255)                                                                                                                    ggCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTCGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTCTAACCAGtac                    >  2:740190/1‑141 (MQ=255)                                                                                                                            cgcACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAc            >  1:739909/1‑143 (MQ=255)                                                                                                                             gcACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAcc           >  1:663052/1‑143 (MQ=255)                                                                                                                               aCCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGt         <  1:1174269/143‑1 (MQ=255)                                                                                                                                ccGCCATGAT‑GGCGCAGGATCAGGGTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGtt        <  1:1202930/143‑1 (MQ=255)                                                                                                                                      tgat‑gGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACg  <  1:404991/143‑1 (MQ=255)                                                                                                                                              |                                                                                                                                        GATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACG  >  NC_000913/3560316‑3560590

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