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

GCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGT  >  NC_000913_3_pae_pgi/3560307‑3560594
                                                                                                                                                      |                                                                                                                                           
gCCGAAATCGAAAAAAACCGCCGAGCCATTGGGAATTTTCTCCGCCACTTTTCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGTCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGc                                                                                                                                              <  2:118417/149‑1 (MQ=255)
         gATAAACAGCGTCGAGACATTGGGGATATGCTCCGCCACTTTGCGGGCGAAGCGCTCTTTTTCTTCGGTCTGGGGGGCCTTGCGATCGTGCCACGGGGTGTTAACCGAACTGGAAGGCAGCGCCCCACCCCCATGCTGGCACGGGATCa                                                                                                                                     >  2:339608/1‑149 (MQ=255)
              aCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTTTTTTGTTTTGATGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAAATGGAAGGCAGCGCGGCACCCCCATGATGGCGCAGGATGAGGttt                                                                                                                                >  2:30437/1‑149 (MQ=255)
                        gCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAATCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATGAGGTTTTGCTCCGCCa                                                                                                                      <  1:697235/149‑1 (MQ=255)
                            ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTTTTCTGTTTGTGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGGTCAGGTTTTGCTCCGCCAGCTc                                                                                                                  >  2:649671/1‑149 (MQ=255)
                             tGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCTATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCa                                                                                                                 <  1:226988/149‑1 (MQ=255)
                               gggATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTTCACCGGCCAGTTAACCGAACTGGAAGGCAGCGCCGCAACGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCAtt                                                                                                               >  1:301962/1‑149 (MQ=255)
                                ggATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTTTTCGGTCTGGTTGGCCTTGCGATCGTGCCACGGCGTGTTAACAGAACTGGAAGGCAGCGCCGCACCGCAATGATGGAGCAGGATACGGTTTTGCTCCGCCAGCTCATTg                                                                                                              >  2:116124/1‑149 (MQ=255)
                                  aTTTGCTCCGCCACTTTGCGGGCGATGCGCGCTCTTTCGTCGGGCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAg                                                                                                            <  2:522081/149‑1 (MQ=255)
                                                     gggCGATGCGCTCTTTTTCTTCGGTCTGTGTGGCCTAGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGTTTTTGCTCCGCCAGCTCATTAAGGTCTCGGCGAATAGTCTGc                                                                                         >  2:691653/1‑149 (MQ=255)
                                                     gggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGc                                                                                         <  1:542405/149‑1 (MQ=255)
                                                         gATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCAtt                                                                                                               >  2:72495/1‑123 (MQ=255)
                                                         gATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCAtt                                                                                                               <  1:72495/123‑1 (MQ=255)
                                                                tctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAAATGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGga                                                                              >  1:342939/1‑149 (MQ=255)
                                                                            gTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAc                                                                  <  2:646440/149‑1 (MQ=255)
                                                                                gggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGc                                                              <  2:342939/149‑1 (MQ=255)
                                                                                 ggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGAAAGGCAGCGCCGCACCGCCATGATGGCGAAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGct                                                             >  1:92824/1‑149 (MQ=255)
                                                                                    ggCCTGGTGATCGGGCTACGGGGTGTTAGCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCCCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCtt                                                          <  2:449871/149‑1 (MQ=255)
                                                                                     gCCGTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTc                                                         <  2:567503/149‑1 (MQ=255)
                                                                                      ccgtGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGTAGCGGCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCGCCGCCAGCTCAATGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCg                                                        <  2:471120/146‑1 (MQ=255)
                                                                                      ccTTGCTATCATGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGCGAACTGCTCTACCAGCCCTTCg                                                        >  2:189925/1‑149 (MQ=255)
                                                                                            gATCGTGCCACGGCTTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTg                                                  <  2:662763/149‑1 (MQ=255)
                                                                                             aTCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGACATGCTCTACCAGCTCTTCGGTACTGa                                                 <  2:244785/149‑1 (MQ=255)
                                                                                                gTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTCCTGACAt                                              >  2:628245/1‑149 (MQ=255)
                                                                                                 tgtctgggCGTGTTAACCGGACTGGAAGTCAGCGGCGCGCCTCCATGAGGGTGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGGCGCGGCGAATAGTCTGCGGGCTGATGGAGAAATGCTCTTCCCGCTCTGCGGTACTGACATa                                             <  2:758644/143‑1 (MQ=255)
                                                                                                    cACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGGCATAAcc                                          >  2:262380/1‑149 (MQ=255)
                                                                                                                aCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGTCTGACTGTGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAc                              <  1:116124/149‑1 (MQ=255)
                                                                                                                aCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGAGGGTTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAc                              <  1:314162/149‑1 (MQ=255)
                                                                                                                   gAACTGGAAGGCAGCGCCGCACGGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAg                           <  2:4598/149‑1 (MQ=255)
                                                                                                                   gAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTATTACTCCGCCAACTCATAGAGGTCGCGGCGAATAGTCTGCGGCCTGACGCAAAAAGGCTCTACCAGCTCTTCGGTACTGACATAACCCAGCTGTTTAACCAg                           >  1:179044/1‑149 (MQ=255)
                                                                                                                    aaCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTCCGCTCCTGACATAACCCTGCTGTTTAACCAgg                          >  2:375336/1‑148 (MQ=255)
                                                                                                                           aGGCAGCGCCGCACCGCCATGATGGCGTAGGATCAGGTTTTGCTCCCGCAGCTTATTGAGGTCGCGGCGAATAGTCTGCGGTCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACAGAACCCTGCTGTTTAACCAGTTCGataa                   <  1:57064/149‑1 (MQ=255)
                                                                                                                                       atcggcATGCTGGCGCAGGGTCAGGTTTTGCCCCGCCAGCTCAGTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACGGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTGGTGAc       <  2:403811/144‑1 (MQ=255)
                                                                                                                                            cATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGAAGGAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCTGTTCGATAATAACATTGTGACTTtgt  >  2:56160/1‑149 (MQ=255)
                                                                                                                                            cATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAAGTCGCGACGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgt  >  1:448464/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                           
GCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGT  >  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: