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_pgi3,560,4551.G93.9% 86.1 / 0.3 34intergenic (‑2/+1)glpR/glpRpseudogene, 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/1);  new base G (14/17);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

GCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGT  >  NC_000913_3_pae_pgi/3560307‑3560594
                                                                                                                                                      |                                                                                                                                           
gCCGAAATCGAAAAAAACCGCCGAGCCATTGGGAATTTTCTCCGCCACTTTTCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGTCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑c                                                                                                                                              <  2:118417/149‑1 (MQ=255)
         gATAAACAGCGTCGAGACATTGGGGATATGCTCCGCCACTTTGCGGGCGAAGCGCTCTTTTTCTTCGGTCTGGGGGGCCTTGCGATCGTGCCACGGGGTGTTAACCGAACTGGAAGGCAGCGCCCCACCCCCATGCTGG‑CACGGGATCa                                                                                                                                     >  2:339608/1‑149 (MQ=255)
              aCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTTTTTTGTTTTGATGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAAATGGAAGGCAGCGCGGCACCCCCATGATGG‑CGCAGGATGAGGttt                                                                                                                                >  2:30437/1‑149 (MQ=255)
                        gCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAATCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATGAGGTTTTGCTCCGCCa                                                                                                                      <  1:697235/149‑1 (MQ=255)
                            ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTTTTCTGTTTGTGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGGTCAGGTTTTGCTCCGCCAGCTc                                                                                                                  >  2:649671/1‑149 (MQ=255)
                             tGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCTATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCa                                                                                                                 <  1:226988/149‑1 (MQ=255)
                               gggATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTTCACCGGCCAGTTAACCGAACTGGAAGGCAGCGCCGCAACGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCAtt                                                                                                               >  1:301962/1‑149 (MQ=255)
                                ggATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTTTTCGGTCTGGTTGGCCTTGCGATCGTGCCACGGCGTGTTAACAGAACTGGAAGGCAGCGCCGCACCGCAATGATGGAG‑CAGGATACGGTTTTGCTCCGCCAGCTCATTg                                                                                                              >  2:116124/1‑149 (MQ=255)
                                  aTTTGCTCCGCCACTTTGCGGGCGATGCGCGCTCTTTCGTCGGGCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAg                                                                                                            <  2:522081/149‑1 (MQ=255)
                                                     gggCGATGCGCTCTTTTTCTTCGGTCTGTGTGGCCTAGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGTTTTTGCTCCGCCAGCTCATTAAGGTCTCGGCGAATAGTCTGc                                                                                         >  2:691653/1‑149 (MQ=255)
                                                     gggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGc                                                                                         <  1:542405/149‑1 (MQ=255)
                                                         gATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCAtt                                                                                                               >  2:72495/1‑123 (MQ=255)
                                                         gATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCAtt                                                                                                               <  1:72495/123‑1 (MQ=255)
                                                                tctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAAATGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGga                                                                              >  1:342939/1‑149 (MQ=255)
                                                                            gTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAc                                                                  <  2:646440/149‑1 (MQ=255)
                                                                                gggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGc                                                              <  2:342939/149‑1 (MQ=255)
                                                                                 ggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGAAAGGCAGCGCCGCACCGCCATGATGG‑CGAAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGct                                                             >  1:92824/1‑149 (MQ=255)
                                                                                    ggCCTGGTGATCGGGCTACGGGGTGTTAGCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCCCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCtt                                                          <  2:449871/149‑1 (MQ=255)
                                                                                     gCCGTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTc                                                         <  2:567503/149‑1 (MQ=255)
                                                                                      ccgtGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGTAGCGGCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCGCCGCCAGCTCAATGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCg                                                        <  2:471120/146‑1 (MQ=255)
                                                                                      ccTTGCTATCATGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGCGAACTGCTCTACCAGCCCTTCg                                                        >  2:189925/1‑149 (MQ=255)
                                                                                            gATCGTGCCACGGCTTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTg                                                  <  2:662763/149‑1 (MQ=255)
                                                                                             aTCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGACATGCTCTACCAGCTCTTCGGTACTGa                                                 <  2:244785/149‑1 (MQ=255)
                                                                                                gTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTCCTGACAt                                              >  2:628245/1‑149 (MQ=255)
                                                                                                 tgtctgggCGTGTTAACCGGACTGGAAGTCAGCGGCGCGCCTCCATGAGGGTG‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGGCGCGGCGAATAGTCTGCGGGCTGATGGAGAAATGCTCTTCCCGCTCTGCGGTACTGACATa                                             <  2:758644/143‑1 (MQ=255)
                                                                                                    cACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGGCATAAcc                                          >  2:262380/1‑149 (MQ=255)
                                                                                                                aCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGTCTGACTGTGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAc                              <  1:116124/149‑1 (MQ=255)
                                                                                                                aCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGAGGGTTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAc                              <  1:314162/149‑1 (MQ=255)
                                                                                                                   gAACTGGAAGGCAGCGCCGCACGGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAg                           <  2:4598/149‑1 (MQ=255)
                                                                                                                   gAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTATTACTCCGCCAACTCATAGAGGTCGCGGCGAATAGTCTGCGGCCTGACGCAAAAAGGCTCTACCAGCTCTTCGGTACTGACATAACCCAGCTGTTTAACCAg                           >  1:179044/1‑149 (MQ=255)
                                                                                                                    aaCTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTCCGCTCCTGACATAACCCTGCTGTTTAACCAgg                          >  2:375336/1‑148 (MQ=255)
                                                                                                                           aGGCAGCGCCGCACCGCCATGATGG‑CGTAGGATCAGGTTTTGCTCCCGCAGCTTATTGAGGTCGCGGCGAATAGTCTGCGGTCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACAGAACCCTGCTGTTTAACCAGTTCGataa                   <  1:57064/149‑1 (MQ=255)
                                                                                                                                       atcggcATGCTGG‑CGCAGGGTCAGGTTTTGCCCCGCCAGCTCAGTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACGGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTGGTGAc       <  2:403811/144‑1 (MQ=255)
                                                                                                                                            cATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGAAGGAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCTGTTCGATAATAACATTGTGACTTtgt  >  2:56160/1‑149 (MQ=255)
                                                                                                                                            cATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAAGTCGCGACGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgt  >  1:448464/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                           
GCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGT  >  NC_000913_3_pae_pgi/3560307‑3560594

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