Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 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_tpiA3,560,4551.G94.4% 107.5 / ‑0.0 37intergenic (‑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 . (0/2);  new base G (14/20);  total (14/23)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTG  >  NC_000913_3_pae_tpiA/3560316‑3560593
                                                                                                                                              |                                                                                                                                          
gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGCATGGCGTTGCGTTCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGGTCa                                                                                                                                    >  2:797736/1‑149 (MQ=255)
       aGCGTCGCGCCATTGGGGATTTGCTCCGCCCCTGTGCGGGCCATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCTACTGCGTGTTAACCGATCTGGAAGGAAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTg                                                                                                                             <  1:387528/149‑1 (MQ=255)
        gCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGc                                                                                                                            <  1:538932/149‑1 (MQ=255)
                         aTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAg                                                                                                           <  1:463977/149‑1 (MQ=255)
                         aTTTGCTCCGCCACTTGGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAg                                                                                                           <  1:925959/149‑1 (MQ=255)
                          tttGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAgg                                                                                                          >  1:678684/1‑149 (MQ=255)
                                    cACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGaa                                                                                                <  1:699082/149‑1 (MQ=255)
                                      ttttGCTGGTGATGCGCTGTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCTGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATa                                                                                              <  2:19444/148‑1 (MQ=255)
                                                gATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCAtt                                                                                                              >  2:914332/1‑123 (MQ=255)
                                                gATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCAtt                                                                                                              <  1:914332/123‑1 (MQ=255)
                                                          ttttCTTCCGTCTGGGTGGCCTTCGGATCGTGCCACGGCGGGTTAACCGAACTGGAAGGCAGCGCCCCCCCGCCCTGCCAACCCCGGCTCAGGTTGTGCGCCTCCAGCTCATTGAGGTCGCGCCGAATAGCCTGCGGGCTGACGGAGaa                                                                          <  2:433240/149‑1 (MQ=255)
                                                           tttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaaa                                                                         >  1:770027/1‑149 (MQ=255)
                                                             gcttcGGTCTGGGTGGCCGTGCGATCGTGCCACGGCGTGGTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATg                                                                       <  2:858477/148‑1 (MQ=255)
                                                               ttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                                     >  1:88936/1‑149 (MQ=255)
                                                                  ggTCCGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTa                                                                  <  1:205185/149‑1 (MQ=255)
                                                                       gggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGc                                                             <  1:772245/149‑1 (MQ=255)
                                                                         atGGCCTTGCGATCGTGCCACGGCGTGTCAACCGAACTGGAAAGCAGCGCCGCACCGCCCTGATAGCGCACGATCAAGTCTTGCTCCGCCACCAGATTAGGGTCGCGGCGAATAGTCTGCGCGCCGACGGAGACATGCCCTACCAGccc                                                           >  2:703034/2‑147 (MQ=255)
                                                                           ggCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCtt                                                         <  1:189668/149‑1 (MQ=255)
                                                                             tcTTGCGGTCGTGCGGCGGGGTGTTAATCGAACTGGAAGGCAGCGCCGCACCGCTATGGTGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCg                                                       <  2:804719/148‑1 (MQ=255)
                                                                             ccTTGGGATCGTGCCACGGCGTGTTAACCGACCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCg                                                       <  2:132895/149‑1 (MQ=255)
                                                                              cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCCCTACCAGCTCTTCgg                                                      >  1:923656/1‑149 (MQ=255)
                                                                               ttGGGATCGTGCCCCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCAGGATGCCGCAGGAGCAGGGTTTGCGCCGCCAGCTCATTGAGGGCGCCGCGAATTGTCTGCGGCGCGACGGAGATATGATCTACCAGCTCTTCgtt                                                     >  2:442404/1‑147 (MQ=255)
                                                                                     tCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGAc                                               <  2:770027/149‑1 (MQ=255)
                                                                                     tCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGAc                                               >  2:883730/1‑149 (MQ=255)
                                                                                          ccACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAACTCTTCGGTACTGACATAAc                                          >  2:669071/1‑149 (MQ=255)
                                                                                           cccGGCGTTTTAACCGAACTGGAAGGCAGCGCCGCTCCGTCATGACGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAcc                                         <  1:797736/147‑1 (MQ=255)
                                                                                             cGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCgg                                                                                      <  1:846379/102‑1 (MQ=255)
                                                                                             cGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCgg                                                                                      >  2:846379/1‑102 (MQ=255)
                                                                                             cGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCgg                                                                                      >  2:846398/1‑102 (MQ=255)
                                                                                             cGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCgg                                                                                      <  1:846398/102‑1 (MQ=255)
                                                                                               tcgttctAACCAAAGTGGAAGTCTGCGGCGGACAGCCATGATGGGGGAGGTTCAGTTTTTTGTCCGTAAGGTCTTAGAGGTCGAGGCAAATAGTCTTCGCTTTGCCTTAGTAATGCTCTACCAGCACTTCGGTACTGACATAACCctgc                                     <  2:767481/143‑1 (MQ=255)
                                                                                                    gtAACCGGACGGGAAGGCAGGGGCGAACCGCCATGAGGGCGCGGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGttt                                <  2:143366/148‑1 (MQ=255)
                                                                                                       aCCGAGCGGGAAGGCAGCGCCGTGCCGTCATGGTGGCTGAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGGCGCTAATAGTCTGCGGGCTGACGGAGATATGCTCTACCAGCTCTTCGGTACTGCCATAAGCCTGCTGTTTGAc                             <  2:623108/149‑1 (MQ=255)
                                                                                                                         cccgaACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGtt           <  1:44783/148‑1 (MQ=255)
                                                                                                                                  ccATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAAGTCTCTACCAGCTCTTCGGTACTGACATAACCCTGCGGTTTAACAAGTTCGCTAATACCGTTGTGACGTtg  >  2:151019/1‑149 (MQ=255)
                                                                                                                                   cATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATaa                                           <  2:584107/107‑1 (MQ=255)
                                                                                                                                   cATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATaa                                           >  1:584107/1‑107 (MQ=255)
                                                                                                                                              |                                                                                                                                          
GATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTG  >  NC_000913_3_pae_tpiA/3560316‑3560593

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: