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% 78.8 / ‑3.7 26pseudogene (151/758 nt)glpRDNA‑binding transcriptional repressor GlpR
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (16/9);  total (16/10)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACC  >  NC_000913/3560313‑3560581
                                                                                                                                                |                                                                                                                              
atCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCCCCGCCCTGCTGGc                                                                                                                                 >  1:612074/1‑143 (MQ=255)
    aTAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGcgcag                                                                                                                             >  1:217464/1‑139 (MQ=255)
       accggCGTGGAGCCATTGTGGAGTTCTTCCGCCACTTGGCGGGCGATGCGCTCTTTTTCTTCGGTCCGGGTGGCCTCGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGAt                                                                                                                          <  2:213940/139‑1 (MQ=255)
              tCGAGCCAGTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGttt                                                                                                                   <  1:527243/143‑1 (MQ=255)
                    cATTGGGGATTGGCTCCGCCACAGTGCGGGCGGTGCGCTCTGTTTCCGCGGGCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTTACCGAACTGGAAGCTAGCGCCGCACCGCCATGTTGGCGCAGGATCAGGTTTTGCTcc                                                                                                             <  1:545314/143‑1 (MQ=255)
                     aTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCg                                                                                                            >  1:735372/1‑143 (MQ=255)
                             tttGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGACGGCGCACGATCAGGTTTTGCTCCGCCAGCTCa                                                                                                    >  2:780409/1‑143 (MQ=255)
                                gCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCCTTg                                                                                                 >  1:520848/1‑143 (MQ=255)
                                 cTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGa                                                                                                >  1:289367/1‑143 (MQ=255)
                                  tgtgCCACTTGGCGGGCGCGTCGGTCTGTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAg                                                                                               <  2:253449/140‑1 (MQ=255)
                                    cGCCACTTTGCCGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTTCCACGGCGGGTTCACCGAACTGGAACGCAGCGCCCCACCGCCCTGACGGCGCAGGCTCAGGTTTTCCTCCACCAGCCCAtcacgca                                                                                             >  2:703877/1‑137 (MQ=255)
                                        aCTTTGCGGGCGATGCGCTCTTTTTCGTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgcg                                                                                         <  2:759189/143‑1 (MQ=255)
                                          tttGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggc                                                                                       >  1:447004/1‑143 (MQ=255)
                                                ggCGATGCGCTCTGTTGCGTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAg                                                                                 <  1:277690/143‑1 (MQ=255)
                                                   gATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCt                                                                              >  1:436760/1‑143 (MQ=255)
                                                          tctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGGCGCGGCGAATGGTCTGCGGGCt                                                                       >  1:584794/1‑143 (MQ=255)
                                                                                   ggCGATCGGGCCGCGGCGTGTGATCCGAACTGGACGGCGGCGCCGCACCGCCATGGTGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGCGTTGTCCCCTGGCTGACGGAGAAATGCTCTACCAGctct                                              <  2:264909/142‑1 (MQ=255)
                                                                                      gATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTCGCGGAGAGATGGTCTACCATCTGttag                                           >  2:497404/1‑141 (MQ=255)
                                                                                        tCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGt                                         >  1:509408/1‑143 (MQ=255)
                                                                                        tCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGt                                         >  1:174286/1‑143 (MQ=255)
                                                                                             ccACGGCGTGTTAACCGAACTGCAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGa                                    <  2:127307/143‑1 (MQ=255)
                                                                                                   cGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAc                              >  1:388408/1‑143 (MQ=255)
                                                                                                     tgtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGCCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAccc                            >  2:747880/1‑143 (MQ=255)
                                                                                                          aCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgccg                       >  1:713141/1‑141 (MQ=255)
                                                                                                          aCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTCTTGCTCCGCCAGCTCATCGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGCGACATGCTCTACCAGCTATTCGCGAACGACATAACGctgatg                       >  2:486203/1‑143 (MQ=255)
                                                                                                                     gggCAGGGCCGCGTCGCCAGGATGGGGCAGGGTCAGGTTTTGCCCCGCCAGCGCATTGAGGTCGCGGCGATTAGTCTGCGGTCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGtt            <  2:748219/142‑1 (MQ=255)
                                                                                                                           cgccgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGACATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataa      >  1:773013/1‑143 (MQ=255)
                                                                                                                               gcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAcc  <  1:248569/143‑1 (MQ=255)
                                                                                                                                                |                                                                                                                              
ATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACC  >  NC_000913/3560313‑3560581

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: