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.2% 55.4 / ‑3.6 24pseudogene (151/758 nt)glpRDNA‑binding transcriptional repressor GlpR
Reads supporting (aligned to +/- strand):  ref base . (1/0);  major base G (13/8);  minor base A (1/1);  total (15/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGA  >  NC_000913/3560319‑3560588
                                                                                                                                         |                                                                                                                                      
aaaCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCGTGCGATCGTGCCACGGCGTCTTACCCGAACTGGAAGGCCGGGCAGCACCGCCCTGGTGGcacaggc                                                                                                                                   >  2:264088/1‑137 (MQ=255)
   cgggtgcGAGCCATTGGGGATTTGCTCCGCCACGTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCa                                                                                                                               <  2:240146/137‑1 (MQ=255)
                ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGc                                                                                                                  <  1:118486/143‑1 (MQ=255)
                     gATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCGTGCGATCGTGCCACGGCGTGTTAACCGGACTGGAAGGCAGCGCCGCAACGCGATGATGGCGCAGGATCAGCATAGGCTCCGCCAGcg                                                                                                             >  2:457472/1‑142 (MQ=255)
                             ccGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAgg                                                                                                     >  1:233966/1‑143 (MQ=255)
                              cggctcgttGCGGGCGGTGCGCTCTTGTGCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTATCCGTACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGt                                                                                                    <  2:540683/136‑1 (MQ=255)
                              cGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGt                                                                                                    <  1:641293/143‑1 (MQ=255)
                                    tttGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCCGGCATTGCACCTCCAGCGCATTGAGGTCgcggc                                                                                              >  2:321208/1‑143 (MQ=255)
                                         gggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATCGAGGTCGCGGCGAATa                                                                                         >  2:137885/1‑143 (MQ=255)
                                         gggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCCCATCGAGGTCGCGGCGactc                                                                                         >  2:667270/1‑140 (MQ=255)
                                                 cgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCgg                                                                                 >  2:542330/1‑143 (MQ=255)
                                                      tttttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAGTCGTCTGCGGGCTGa                                                                            >  1:5814/1‑143 (MQ=255)
                                                              cGGTCCGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaaa                                                                    <  1:368022/143‑1 (MQ=255)
                                                                gTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCCCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGCCGGAGAGATg                                                                  >  2:381870/1‑143 (MQ=255)
                                                                           cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGGAATGCTCTACCAGct                                                       >  1:500884/1‑143 (MQ=255)
                                                                                   cGTGCAACGGCGTGTTGACCGAGCGGGAAGGCGGTGCCGCACCGCCATGATGGCGCAGGCTCAGGTTTTGCGCCGCCAGCTCATTGAGGTCGCGGCGAATCGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTa                                               <  2:400055/143‑1 (MQ=255)
                                                                                    gTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCAACGCCAAGACAATGTAAATGCTTGTTATGGACTAGTTACTCATTGAGGTCGCGGCGAATAGACTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTAc                                              >  1:676176/1‑143 (MQ=255)
                                                                                            gCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCACAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATaa                                      >  1:679610/1‑143 (MQ=255)
                                                                                             cGTGTTACCCGACCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAc                                     <  1:34525/143‑1 (MQ=255)
                                                                                                                 gCAGCGCCGCACCGCCATGATGGCGCAGGATAAGGTGTTGCTCGGCCAGCTCATTGAGGTAGCGGAGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACTAGATGTTCGGTACTGAGCTAACCATGATGTGTAAACAGCTCg                 >  2:372937/1‑143 (MQ=255)
                                                                                                                     cgccggACCGCCATGATGGCGCAGGATCAGGTCTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCCGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACCTAACCCTGCTGTTTAACCAGTTCGataa             <  1:264088/143‑1 (MQ=255)
                                                                                                                     cgccgcACCGCCATGCTGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAGAGTCTGCGGGCTGACGGCGGAATGCTCTACCCGCTCTTCGGTACTGACATAACCCTGCTGTCTAGCCAGTTCGacaa             >  1:593701/1‑143 (MQ=255)
                                                                                                                      gccgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGGTTAACCAGTCCGataat            >  2:175464/1‑143 (MQ=255)
                                                                                                                        cgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGCTAATAc          >  1:132364/1‑143 (MQ=255)
                                                                                                                                cATGGTGGCACAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGa  <  1:417688/143‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                      
AAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGA  >  NC_000913/3560319‑3560588

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: