Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,560,455:1 +G 100% 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.G91.7% 63.8 / ‑1.1 24pseudogene (151/758 nt)glpRDNA‑binding transcriptional repressor GlpR
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (9/12);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATA  >  NC_000913/3560315‑3560579
                                                                                                                                               |                                                                                                                            
cGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGcgc                                                                                                                             >  1:123396/1‑141 (MQ=255)
                    ttGGGGGTTTGCTGCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGGCTGGGTGGCCTTGCGATCGTGCCACGGCGTGCTAACCGAGTTCGAAGCCGGCGCCGCGCCGCCCGGATGGCGCAGGACGAGGTTTAGCTCCGc                                                                                                         >  2:191027/1‑143 (MQ=255)
                    ttGGGGATTTGTGCCGCCGCTTTGCGGGCGATGGGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGCTCGTGCCACGGCGTGTTAACCGAACTGGAAGGCCGCGTCGCACCGCCATGATGGCGCAGGATCAGGTTTTGTTCCGc                                                                                                         <  2:147538/143‑1 (MQ=255)
                    ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTGACCGCACTGAACGGCCGCGGCGCACACTCGTGATGGCACAGGCTCAGGGTGTGcgccgc                                                                                                         >  2:232013/1‑143 (MQ=255)
                    ggggggAGGTGCTCCGCGACTTGGCGGGCGGTGCGCTCGTTTTCTTCGGTCTGGGTGGCGTTGCGATCGCGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGc                                                                                                         <  2:10431/141‑1 (MQ=255)
                                  cGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGt                                                                                           >  1:472459/1‑143 (MQ=255)
                                       cTTTGCGTGCTATGCGCGCTTTTTCTCCGGTCTGGGTGGCTTTGCGATCGTGCCACGGCGTGTTAACCGAGCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcgg                                                                                      <  2:441544/143‑1 (MQ=255)
                                          tGCGGGCGGTGCGCGCGTGTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGAGCAGGTTTTGCTCCGCCAGCTCGTTGAGGTCGCGGCGa                                                                                   <  2:297913/143‑1 (MQ=255)
                                           gcgggcggGGTGCGCTTTTTGTTCGGGCTTGGTGGCGTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGaa                                                                                  <  2:249379/143‑1 (MQ=255)
                                           gCGGGCGATGCGCTCTTTTTCTTCGGTCTGGCTGGCCGTGCTATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGaa                                                                                  <  1:244887/143‑1 (MQ=255)
                                                    gcgcTCTTTTTCTTCGGTCCGGGTGGCCTTGCGGTCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCg                                                                         <  1:460920/143‑1 (MQ=255)
                                                          tttttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTGATTGAGGTCGCGGCGAATAGTCTGCGGGCTGa                                                                   >  2:289482/1‑143 (MQ=255)
                                                                      cTGGGGGGCCTTGCGAGCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                       <  2:308754/143‑1 (MQ=255)
                                                                           tGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                  >  1:472330/1‑143 (MQ=255)
                                                                             gCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGCTGGACCAAATTCATGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAg                                                <  1:111778/143‑1 (MQ=255)
                                                                                        gTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTAc                                     >  1:201859/1‑143 (MQ=255)
                                                                                         tGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAGAGGCTCTACCAGCTCTTCGGTAcg                                    >  1:313369/1‑142 (MQ=255)
                                                                                           ccACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTTGTTGAGGAGGCGCCAAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGa                                  >  1:355280/1‑143 (MQ=255)
                                                                                               tgCGTGTTAACCGAACTGGACGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATGGAGCTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATa                              <  1:261116/142‑1 (MQ=255)
                                                                                                gCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGGGCAGGATCAGGTCCTGCTCAGATCGATCAGAAATGGCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATaa                             >  1:487637/1‑143 (MQ=255)
                                                                                                              cgggAAGGCAGCGCCGCGCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAc               <  2:316849/141‑1 (MQ=255)
                                                                                                                  aaGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTCAACCAGt           >  2:169095/1‑143 (MQ=255)
                                                                                                                    ggCAGCGCCGCACGGCCATGATGGCGCAGGGTCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc         <  1:416163/143‑1 (MQ=255)
                                                                                                                    ggCAGCGCCGCACCGCCAGGGTGGCGCAGGATCAGGTTCTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc         <  1:143660/143‑1 (MQ=255)
                                                                                                                           ccgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCCAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGCCATAACCCTGCTGTCTAACCAGTTCGatagtc  >  2:309087/1‑140 (MQ=255)
                                                                                                                                               |                                                                                                                            
CGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATA  >  NC_000913/3560315‑3560579

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: