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.G98.0% 148.9 / ‑4.0 50pseudogene (151/758 nt)glpRDNA‑binding transcriptional repressor GlpR
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (26/23);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

ATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGA  >  NC_000913/3560317‑3560588
                                                                                                                                            |                                                                                                                                     
gtAATCAGCGTCCAGCTGTTGGGGATGTGCTCCGCCCCTTTGCGGGCGATGCGCTCTTTTCCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGCTAACCGAACTGGAAGGCAGCGCCGCACCGCCATCATGGcggtg                                                                                                                                    <  1:762699/142‑5 (MQ=255)
  aaaCAGCGTCGCGCCATTGGGGCTTTGCTCCGCCTCTTTGCGGGGGATGCGCTCTTTTGCTTGGGTCTGGGTGGCTTTGACATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGa                                                                                                                                  <  2:757266/143‑1 (MQ=255)
    aCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAAGCAGCGCCGCACCGCCAAGATGGCGCAGGATc                                                                                                                                >  1:1003667/1‑143 (MQ=255)
     cggCGTCGAGCCATTGGGGATTTGCTCCGCCACGTTGCGGGCGATGCGCTCTTTTTCTTAGGTGTGGTGTCCCTGGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCa                                                                                                                               <  1:843165/141‑1 (MQ=255)
          tCGGGCCGTTGGGGCTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCCGGGTGACCGTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGttt                                                                                                                          <  2:1077834/143‑1 (MQ=255)
          tCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGCTGGCGCAGGGTCAGGatt                                                                                                                          >  2:799111/1‑143 (MQ=255)
            gAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGGCTGGGTGGCCTAGCGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGAAGCGCAGAACCGCCATGATGGCGCAGGATCAGGtgctg                                                                                                                        >  2:861863/1‑139 (MQ=255)
                  ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGACTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCGGGTGTTGCTCCGc                                                                                                                  >  1:530659/1‑143 (MQ=255)
                   cagggATTTGCTCCGCCCCTGTGCGGGCGCTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGcc                                                                                                                 <  1:519078/141‑1 (MQ=255)
                    ggggGTGTGCCCCGCCACGGTGCGGGCGATGCGCTCGTTTTCTTCGGGCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTGAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCa                                                                                                                <  1:149691/143‑1 (MQ=255)
                        aTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTc                                                                                                            >  2:338714/1‑143 (MQ=255)
                            gCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTTGCCTTGCTAGCGTGCCACGGCGTGTGAACCGAGCTGGCAGGCAGCGCCGCACCGCCATGATGGCGCAGGAGAAGGTCTTGATGCGGCACGTCATTg                                                                                                        >  2:69082/1‑143 (MQ=255)
                            gCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGCAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCTGGTTTTGCTCCGCCAGCTCGTTg                                                                                                        >  2:852958/1‑143 (MQ=255)
                            gCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCGCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                                                                                        >  1:561852/1‑143 (MQ=255)
                              tCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAg                                                                                                      >  1:227243/1‑143 (MQ=255)
                              tCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGCTTTGCTCCGCCAGCGCATTgcg                                                                                                      >  2:1111285/1‑141 (MQ=255)
                              tCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCAGGGAGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCCTGATCGCGCAGGATCCGGTTTTGCTCCGCCCGCTcctcgcc                                                                                                      >  1:217191/1‑137 (MQ=255)
                                     cTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATGGAGGTCgcgg                                                                                               >  1:335129/1‑143 (MQ=255)
                                      tttGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAGCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggc                                                                                              <  1:627550/143‑1 (MQ=255)
                                         gCGGGCGACGCGCGCGGGTGCGTCGGTCTGGGTGGCCGTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGaa                                                                                           <  1:231887/143‑1 (MQ=255)
                                            ggCGTTGCGCTCTTGTGTCTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTGTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAg                                                                                        <  2:348539/143‑1 (MQ=255)
                                                     ctcgTTTTCTGCGGTCTGGGTGGGCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGc                                                                               <  2:464488/143‑1 (MQ=255)
                                                              ttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGgagc                                                                      >  2:887360/1‑142 (MQ=255)
                                                               gtggTCTGGGTGGGCTTGCGTGCGTGCCACGGCGTGTTGACCGACCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaa                                                                     <  2:582522/141‑1 (MQ=255)
                                                                    cTGGGTGGCCTTGCGAGCGTGCCACGGCGTGGTAACCGAACTGGGAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                                <  1:545794/143‑1 (MQ=255)
                                                                       ggTGGGCTGGCGCCCGGGTCACGGCGTGTTAACCGAACGGGAAGGCCGCGCCGCACCGCCATGATGGCGCAGGAGCAGGTTTTGCTCCGCCAGCTCATGGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTa                                                             <  2:1077738/143‑1 (MQ=255)
                                                                           tccTGGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAg                                                         <  2:739808/142‑1 (MQ=255)
                                                                           gCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAg                                                         <  1:991071/143‑1 (MQ=255)
                                                                                     cgggCCGCGGCGTGTTGACCGGACTGGAAGGCAGCGCCGCACCGCCATGGTGGCGCAGGATCAGGTTGTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTa                                               <  2:716986/140‑1 (MQ=255)
                                                                                     cGTGCCACGGCGGGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCTGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTa                                               <  1:579100/143‑1 (MQ=255)
                                                                                          cACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCAAGATGGCGCAGGAACAGGTTTAGCTCCGCCAACTCACTGCGGTCGCGACGAATAAGCTGCAGGCTCACCGAGAAATGCTATACCAACTCTTCAGCACTGAc                                          >  2:821683/1‑143 (MQ=255)
                                                                                           aCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACa                                         >  1:319545/1‑143 (MQ=255)
                                                                                            cGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTGTACCAGCTCTTCGGTACTGACAt                                        >  1:1035519/1‑143 (MQ=255)
                                                                                            cGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACAt                                        <  2:1055316/143‑1 (MQ=255)
                                                                                                  tttAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCACTGAGGTCGCGGCGAATAGTCTGCCGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACCAAACCct                                  >  2:1100591/2‑143 (MQ=255)
                                                                                                     aaCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAACTCTTCGGTAAAGACATAAcccacct                               >  2:393073/1‑139 (MQ=255)
                                                                                                      aCCGAACTGGAAGGCAGCGCAGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCCGCTCATAGAGGTCGCGGCGAATAGTCTGCGGGCTGACGCAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgctg                              >  2:1012340/1‑143 (MQ=255)
                                                                                                        cGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGtt                            >  2:988879/1‑143 (MQ=255)
                                                                                                         gAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGttt                           <  1:835347/143‑1 (MQ=255)
                                                                                                         gAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGGCGCGGCGAATAGTCTGCGGGCTGACGGGGAGATGCTCTACCAGCTCTTCGGTACTGACATAACActgctaaga                           >  1:2200/1‑139 (MQ=255)
                                                                                                          aaCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTCTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCACTTCGGTACTGACATAACCCAGCTGTTTa                          >  2:1076428/1‑143 (MQ=255)
                                                                                                           aCTGGAGGGCATCTCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTaa                         <  2:475685/143‑1 (MQ=255)
                                                                                                            gtGGGAGGCGGCGCCGAACCGCCATGTTGGTGCAGTATCAGGTTTTGCTGCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAc                        <  2:282237/142‑1 (MQ=255)
                                                                                                            gcggACGGCAGCACCTCACCGCCAGGATGGCGCAGGCTCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAc                        <  2:393245/141‑1 (MQ=255)
                                                                                                             tGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGCGCTGACGGAGAAATGCTCTACCAGCTCTTCGCTACTGACATAACACTGCTGATTaaat                       >  1:492576/1‑141 (MQ=255)
                                                                                                                 aGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGtt                   <  2:227243/143‑1 (MQ=255)
                                                                                                                     aGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAGTAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGATCTTCGGTACTGACATAACCATGCTGTTTAACCAGATCGat               >  1:401319/1‑143 (MQ=255)
                                                                                                                       cgccgcCCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataa             <  1:11502/143‑1 (MQ=255)
                                                                                                                         ccgcACGGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCGGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataata           <  2:406199/143‑1 (MQ=255)
                                                                                                                               cGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGAGGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtg     >  2:344705/1‑143 (MQ=255)
                                                                                                                                  cATGATGGCGCAAGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATACTACCGTAGTGa  >  2:436632/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                     
ATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGA  >  NC_000913/3560317‑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: