Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_pae_pgi	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_pae_pgi	3,560,455	1	.	G	98.0%	162.6 / ‑4.0	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 . (0/1);  new base G (25/23);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTT  >  NC_000913_3_pae_pgi/3560307‑3560592                                                                                                                                                       |                                                                                                                                          gCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGTTGTCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑Gc                                                                                                                                            >  2:117594/1‑149 (MQ=255)         cGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCTCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGAGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGagc                                                                                                                                    >  1:360714/1‑147 (MQ=255)          gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAACGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCa                                                                                                                                   >  1:288916/1‑149 (MQ=255)           aTAAACGGCGTAGGGCCATTGGGGGGTTGCTCCGCCACTTTGGGGGCGATGCGTTGTTTTGCTTCGGTCTGGGGGGCCTTGCGAGCGTGCCACGGCGTGTTAACCGAAATTGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAg                                                                                                                                  <  2:324764/149‑1 (MQ=255)                     tCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTcc                                                                                                                        >  1:202864/1‑149 (MQ=255)                          ccATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAg                                                                                                                   >  2:827028/1‑149 (MQ=255)                             ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTc                                                                                                                >  1:181452/1‑149 (MQ=255)                                gcgATTGGCTCCGCCATTGTGCGGGCGATGGGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCAtt                                                                                                             <  2:8623/147‑1 (MQ=255)                                 ggATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                                                                                            <  1:563221/149‑1 (MQ=255)                                       gCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCTGGTTTTGCTCCGACAGCTCATTGAGGTcg                                                                                                      >  1:981321/1‑149 (MQ=255)                                                      gggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGc                                                                                       >  2:126317/1‑149 (MQ=255)                                                      gggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGc                                                                                       <  1:828015/149‑1 (MQ=255)                                                           aTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCt                                                                                  >  1:974095/1‑149 (MQ=255)                                                              cgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGAc                                                                               <  2:102844/149‑1 (MQ=255)                                                                  ctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGATTTTCTCCGCCAGCTCATTGAGGTCGCAGCGAATAGTCTGCGGGCTGACGgag                                                                           >  1:245281/1‑149 (MQ=255)                                                                       tcttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGCAAGATGCAACGGGGTCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATg                                                                      <  1:246758/149‑1 (MQ=255)                                                                                ggggTGGCCTTGCGATCGTGCCACGGCGTTTTAACCGACCTGGAAGGCAGCGCCGCACCGCCACGCTGCGG‑TAGAATCAGGTTTTCTTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAg                                                             <  1:250661/148‑1 (MQ=255)                                                                                   gTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTATAACTCCGCCACCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGctc                                                          <  1:202090/149‑1 (MQ=255)                                                                                         ttGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGCCGGAGAAATGCTCTACCAGCTCTTCGGt                                                    >  1:656978/1‑149 (MQ=255)                                                                                          tGCGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTa                                                   <  2:667580/149‑1 (MQ=255)                                                                                              aTCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGa                                               <  1:517094/149‑1 (MQ=255)                                                                                                 gTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                                    <  1:502802/125‑1 (MQ=255)                                                                                                 gTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                                    >  2:502802/1‑125 (MQ=255)                                                                                                     cACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGctct                                                                  <  1:436684/123‑1 (MQ=255)                                                                                                     cACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGctct                                                                  >  2:436684/1‑123 (MQ=255)                                                                                                     cACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGCCATAAcc                                        >  2:614572/1‑149 (MQ=255)                                                                                                     cACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAcc                                        >  1:911822/1‑149 (MQ=255)                                                                                                      aCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTg                                                                                        <  1:102605/100‑1 (MQ=255)                                                                                                      aCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTg                                                                                        >  2:102605/1‑100 (MQ=255)                                                                                                      aCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAccc                                       >  2:252265/1‑149 (MQ=255)                                                                                                        ggCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCTCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctg                                     <  1:117594/149‑1 (MQ=255)                                                                                                         gCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgc                                    >  2:553075/1‑149 (MQ=255)                                                                                                             gtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACCTAACCctgcttct                                >  2:109317/1‑146 (MQ=255)                                                                                                                aaCCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGTTGACGGAGAAATGCTGTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTaa                             <  2:754832/149‑1 (MQ=255)                                                                                                                 aCCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAc                            <  2:181452/149‑1 (MQ=255)                                                                                                                      aCTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGtt                       <  2:787590/149‑1 (MQ=255)                                                                                                                        tGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAccc                                       <  2:638708/131‑1 (MQ=255)                                                                                                                        tGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAccc                                       >  1:638708/1‑131 (MQ=255)                                                                                                                         ggAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGCACTGACATAACCCTGCTGTTTAACCAGTTCGa                    <  2:974095/149‑1 (MQ=255)                                                                                                                          gAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAc                            >  2:243436/1‑140 (MQ=255)                                                                                                                          gAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGat                   >  2:159746/1‑149 (MQ=255)                                                                                                                          gAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGat                   <  1:1048962/149‑1 (MQ=255)                                                                                                                          gAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACCGACATAACCCTGCTGTTTAAc                            <  1:243436/140‑1 (MQ=255)                                                                                                                              gCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTCCGataatt               >  2:320223/1‑148 (MQ=255)                                                                                                                                aGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAcc                                        <  2:682080/122‑1 (MQ=255)                                                                                                                                aGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAcc                                        >  1:682080/1‑122 (MQ=255)                                                                                                                                  cgccgcACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCCATAATACCGt           >  2:158502/1‑149 (MQ=255)                                                                                                                                   gccgcACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGtt          <  1:126317/149‑1 (MQ=255)                                                                                                                                           gCCATGATG‑GCGCAGGATCAGGTTTTGGTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGtt  <  2:288916/149‑1 (MQ=255)                                                                                                                                           gCCATGATG‑GCGCACGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTACCCAGTTCCATACTACCGTTGTGACttt  >  2:887759/1‑149 (MQ=255)                                                                                                                                                       |                                                                                                                                          GCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTT  >  NC_000913_3_pae_pgi/3560307‑3560592

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