Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,560,455 +G pseudogene (151/758 nt) glpR ← DNA‑binding transcriptional repressor GlpR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,560,4551.G96.2% 77.0 / ‑3.5 26pseudogene (151/758 nt)glpRDNA‑binding transcriptional repressor GlpR
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (17/8);  total (17/9)
Fisher's exact test for biased strand distribution p-value = 3.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTG  >  NC_000913/3560318‑3560587
                                                                                                                                           |                                                                                                                                    
tAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGGGAGCTTGCGATCGTGCCACGGCGTGTTAACCGAGCTGGAAGGCAGCGCCGCACCGCCATGAT‑GGcgcagg                                                                                                                                  >  1:49558/1‑138 (MQ=255)
      gCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAgg                                                                                                                            <  1:336007/143‑1 (MQ=255)
       cGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTAGGTCTGGGTGGCCTTGCGACCGTGCGACGGCGTGTTACCCGAACTGGAAGGCCGCGCCGCACCGCCATGAG‑GGCGCGGGATCAGGt                                                                                                                           >  2:88776/1‑143 (MQ=255)
         tCGGGCCATTGGGGATTCGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGttt                                                                                                                         <  1:446982/143‑1 (MQ=255)
          cGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGtttt                                                                                                                        >  1:545672/1‑143 (MQ=255)
            aGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAACGCAGCGCCGCACCGCCATGAT‑GGCGCAGGAAAAGGTTTTGc                                                                                                                      >  1:200886/1‑143 (MQ=255)
                aTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCCTGAT‑GGCGCAGGATCAGGTTTTGCTcct                                                                                                                  >  1:136801/1‑142 (MQ=255)
                 ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGc                                                                                                                 >  1:536742/1‑143 (MQ=255)
                           gCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                                                                                       <  1:541059/143‑1 (MQ=255)
                                              gATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCCCATTGAGGTCGCGGCGCATAGTca                                                                                    >  1:347001/1‑142 (MQ=255)
                                                                 gTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATg                                                                 >  2:635873/1‑143 (MQ=255)
                                                                      ggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGACGAATAGTATGCGGGCTGACCGAGAAATGCTCTa                                                            >  2:128286/1‑143 (MQ=255)
                                                                           ccTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCGGc                                                       >  2:244999/1‑143 (MQ=255)
                                                                            cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCCTTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGGAATGCTCTACCAGct                                                      >  2:159521/1‑143 (MQ=255)
                                                                            cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGGTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGct                                                      >  1:54162/1‑143 (MQ=255)
                                                                            cTTGCGATCGTGCCACGGCGTGTTAACAGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGATGTGAACGACCGTGTCATTCAGGTCGCGGCGCATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGct                                                      >  1:104259/1‑143 (MQ=255)
                                                                                cGATGGGGCCACGGTGTGTTGAGCGGACTGGAAGGCGGGGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTc                                                  <  2:200886/143‑1 (MQ=255)
                                                                                 gATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCg                                                 <  1:765422/143‑1 (MQ=255)
                                                                                 gATAGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGAGGGTTGACGGAGAAATGCTCTACCAGCTCTTCg                                                 <  1:564451/143‑1 (MQ=255)
                                                                                   tCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCGTGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGATCTACCAGCTCTTCGGt                                               >  1:8774/1‑143 (MQ=255)
                                                                                      tGCCACGGCGTGTGAACCGAGCTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACt                                            <  2:318230/143‑1 (MQ=255)
                                                                                           cGGCGTGTTAACCGAACTGGAAGGCACCGCCGCACCGCCATGCT‑GGCGCAGCATCAGGTTTTGCTCTACCAACTCATTAAGGTCGCAGCGAATAGTCTGCGGGCTGACGAAGAAATGCTCTACCAACTCTCCCGTACtcccat                                       >  1:406417/1‑138 (MQ=255)
                                                                                            ggCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATa                                      <  1:750948/143‑1 (MQ=255)
                                                                                                 gtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCAATGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGCTACCGACATAACCct                                 >  2:70821/1‑143 (MQ=255)
                                                                                                       cGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAACAGCTCTTCGGTACAGAACTAACCctgctgcg                           >  2:444977/1‑141 (MQ=255)
                                                                                                                             ccGCCATGATGGGG‑CGGGAGCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGtt     <  2:558012/143‑1 (MQ=255)
                                                                                                                                ccATGAT‑GGCGCAGGATCAGGTTTTGCTCCGCCAGCCCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtgtg  >  1:744495/1‑143 (MQ=255)
                                                                                                                                           |                                                                                                                                    
TAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAT‑GGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTG  >  NC_000913/3560318‑3560587

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