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.G96.2% 75.1 / ‑2.8 26pseudogene (151/758 nt)glpRDNA‑binding transcriptional repressor GlpR
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (15/10);  total (15/11)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATA  >  NC_000913/3560315‑3560579
                                                                                                                                              |                                                                                                                            
cGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGCTGGcgc                                                                                                                             >  1:3243/1‑141 (MQ=255)
           gTCGGGCCAGTGGGGAGTTGCTCCGCCACTTAGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGtt                                                                                                                  <  1:65661/143‑1 (MQ=255)
            tCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGCCTAGCTTGCCTTGCGAACGTGCCACGGCGTGTTAACAGGAAATGAAGGCAGAGACACTCCTCCAACATGGCGCAGGAACAGGttt                                                                                                                 >  2:63496/1‑143 (MQ=255)
                   aTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCg                                                                                                          >  1:532293/1‑143 (MQ=255)
                    ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGcgccgc                                                                                                         >  1:67464/1‑143 (MQ=255)
                       gggATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGCTCGGGTTTTGCTCCGCCAg                                                                                                      >  1:319466/1‑143 (MQ=255)
                             tGCTCCGCCACTTGGCGGGCGCGGCGCTCTGTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCAtt                                                                                                <  2:223668/143‑1 (MQ=255)
                                       cTTTGCGGGCGATGCGCTCCTTTTCTTCGCTCTGCGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcgg                                                                                      <  1:370032/143‑1 (MQ=255)
                                             gggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGCAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCCTCGAGGTCGCGGCGAGTa                                                                                >  1:123245/1‑143 (MQ=255)
                                              ggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGACCGTGCCACGGCGTGTTAACCGAACTGGAGGGCAGCGCCGCACCGCCATGCTGGCGCAGGATCAGGTTTTGCTCCGCGAGCTCATTGAGGGCGCGGCGAATAg                                                                               >  2:347418/1‑143 (MQ=255)
                                                      gcTCTTTTTCTTCGGTCTGGGTTGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCggg                                                                       <  2:2457/143‑1 (MQ=255)
                                                        tctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCt                                                                     <  1:128940/143‑1 (MQ=255)
                                                           ttttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGAc                                                                  >  1:433062/1‑143 (MQ=255)
                                                                 tcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaa                                                            >  1:240524/1‑143 (MQ=255)
                                                                 tcGGGCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaaa                                                           <  2:450937/143‑1 (MQ=255)
                                                                    ggcTGGGTGGCCTTGCGATCGTGCGACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATg                                                         <  1:229605/141‑1 (MQ=255)
                                                                    gTCTGGGTGGCCGTGCGATCGTGCCACGGCGTGGTAGCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATg                                                         <  2:151043/143‑1 (MQ=255)
                                                                           gggCCTTGCGGTCGTGCTACGGCGTGGTAGCCGACCTGGAAGGCGGCGCGGCACCGCCCTGATGGCGCACGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGTCTGACGGAGAAATGCTCTAcc                                                  <  2:462593/142‑1 (MQ=255)
                                                                                       cGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCGCATTGAGGTCGCAGCGAATAGTCTGGGGGGTGACGGAGAAATGCTCTACCAGCTCTTAGGTa                                      >  2:140506/1‑143 (MQ=255)
                                                                                        gtaccACGGCGTGTTAACCGAACTGGAAGACAGCGCCGCACCGCCATGATAGCGCAGGATCAGGTTGTGCTCCGCCAGCCCATTGAGGTCGCGGCGAATAGTCTGCGGACTGACGGACAAATGCTCTACCAGCTCTTCTGTAc                                     >  2:43450/4‑143 (MQ=255)
                                                                                            cACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGAc                                 >  1:450864/1‑143 (MQ=255)
                                                                                             aCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAAGGCTCTACCAGCTCTTCGGTACTGACa                                >  2:263449/1‑143 (MQ=255)
                                                                                                       aaCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgct                      <  1:149391/143‑1 (MQ=255)
                                                                                                              cTGGAAGGCGGCGCCGCGCGGCCATGGTGGGGCAGGATCAGGTTTTGCTCCGCCAGGTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAc               <  2:123245/143‑1 (MQ=255)
                                                                                                                 gAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAg            >  1:250310/1‑143 (MQ=255)
                                                                                                                 gAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATGGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGCAATGCTCTACGAGCTCTTCGGTACGGACATGACCCCGCTGGTTaaggcg            >  2:284123/1‑139 (MQ=255)
                                                                                                                           ccgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTGTTCGGTACTGACATAACCCTGCTGGTTAACCACTTCGataata  >  2:440058/1‑143 (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: