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.G97.7% 142.6 / ‑3.0 44pseudogene (151/758 nt)glpRDNA‑binding transcriptional repressor GlpR
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (22/21);  total (22/22)
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.18e-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.

CGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGT  >  NC_000913/3560315‑3560586
                                                                                                                                              |                                                                                                                                   
cGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCAGCACCGCCATGATGGcgc                                                                                                                                    >  1:529072/1‑141 (MQ=255)
    aaaCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGa                                                                                                                                <  1:707371/143‑1 (MQ=255)
      aCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATc                                                                                                                              >  1:874933/1‑143 (MQ=255)
         gCGTCGCGCCATTGGGGATTTGCCCGCCCACTGTGCGTGCGATGCGCTCTTTTTCTTCGTTCTGGGTGGCCTTGCGATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAgg                                                                                                                           <  2:914865/143‑1 (MQ=255)
         gCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAgg                                                                                                                           <  1:345452/143‑1 (MQ=255)
              gAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCGTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTg                                                                                                                      <  1:119023/143‑1 (MQ=255)
                    ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTGCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGc                                                                                                                <  2:711108/143‑1 (MQ=255)
                         ggtttGCTCCGCCACTTGGCGGGCGATGCGCTCGTGTTCTGCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCt                                                                                                           <  1:319875/141‑1 (MQ=255)
                         gATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCt                                                                                                           >  1:541754/1‑143 (MQ=255)
                               cTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTGAACCGAACTGGAAGGCAGCGCCGCACCGACATGATGGCGCAGGATCAGGGGTCGCTCCGCCAGCTCATTgc                                                                                                     >  1:207357/1‑142 (MQ=255)
                                tCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCCTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCCCCGCACCGCCCTGATGGCGCAGCCTCAGGGTTTGCCCCGCCCGCTCAtgcag                                                                                                    >  2:426054/1‑139 (MQ=255)
                                  cGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCGGCCAGCGCATTGAGGt                                                                                                  >  1:243595/1‑143 (MQ=255)
                                             gggCGATGCGCTCTTTTTCTACGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGACAGCTCATTGAGGTCTCGGAGAATa                                                                                       >  2:206583/1‑143 (MQ=255)
                                                cGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGTATCAGGTTTTGCTCAGCCAACTCATTGAGGTCGCGGCGAATAGTc                                                                                    >  1:552729/1‑143 (MQ=255)
                                                   tGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAGTAGTCTGc                                                                                 >  1:829238/1‑143 (MQ=255)
                                                      gcTCTTTTTCTTCGGTCTGGGTGGNCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCNCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTNATTGAGGTCGCGGCGGATCGTCGCCggg                                                                              >  1:76987/1‑143 (MQ=255)
                                                         ctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTg                                                                           >  1:207100/1‑143 (MQ=255)
                                                         ctTTTTCCTCGGCCTGGGTGGCCTTGCGATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTg                                                                           <  2:442788/143‑1 (MQ=255)
                                                          tttttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCAGTGAGGTCGCGGCGAAGAGTCTGCGGGCTGa                                                                          >  1:265501/1‑143 (MQ=255)
                                                           ttttCTTCGGTCTGGGGGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGAc                                                                         <  1:240016/143‑1 (MQ=255)
                                                            tttCTTCGGGCGTTGTTGCCTTTCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCAGGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACg                                                                        <  1:3235/143‑1 (MQ=255)
                                                                ttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGgaga                                                                    >  2:337828/1‑143 (MQ=255)
                                                                ttcGGTCTGGGTAGCCTTGCGATCGTGCCACGGCGTGTTCACCGAACTGGAAGGCAGCGCCGCACCGCCATGCTGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCACTGAGGACGCGGCCAATAGTCCGCGGGCCGACGgaga                                                                    >  2:446327/1‑143 (MQ=255)
                                                                 tctgTCTGAGTGGCCTTGCGATCGTGCCACGGCTTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCAGGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaa                                                                   <  1:305072/140‑1 (MQ=255)
                                                                    gTCTGGGTGGCCGTGCGATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATg                                                                <  1:64918/143‑1 (MQ=255)
                                                                       tGGGTGGCCTTGCGATCGTGCCACGGCGTGGTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGGGCAGGATCAGGTTTTGCTCCGGCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGctc                                                             <  1:525278/143‑1 (MQ=255)
                                                                       cgggTGGCCTTGCGATCGCGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGCTGGCGCCGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGctc                                                             <  2:799729/142‑1 (MQ=255)
                                                                         ggTGGCCTTGGGGTCGTGCCACGTCGGGTTAACCGAACTGTCAAGCAGCACCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATGGAGGTCGCGGCGAATCGTCTGCGGGCTGACGGAGAAATGCTCTa                                                           <  2:363465/143‑1 (MQ=255)
                                                                               cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGCCGGAGAAATGCTCTCCCCGct                                                     >  1:29776/1‑143 (MQ=255)
                                                                                tggcgcTCGTGCCACGGCGTGCTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGctc                                                    <  1:827119/141‑1 (MQ=255)
                                                                                   cGATCGTGCCACGGCGTGTCAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTc                                                 <  1:785759/143‑1 (MQ=255)
                                                                                    gATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCg                                                >  1:220726/1‑143 (MQ=255)
                                                                                    gATCGGGCCAGGGGGTGGTAACCGAACTGGAAGGCAGCGCCGCGCCGCCATGAGGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCg                                                <  2:817700/143‑1 (MQ=255)
                                                                                        gTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGCAATGCTCTACCAGCTCTTCGGTAc                                            >  1:915913/1‑143 (MQ=255)
                                                                                            cACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGAc                                        >  1:785188/1‑143 (MQ=255)
                                                                                              cGGCGTGTTATCGGATGTGTAAGGCGGCGCCGCACCGCCATGATGGCGCAGGGTCAGGTTTTGCGCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACAt                                      <  2:91614/143‑1 (MQ=255)
                                                                                                           gAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCCGCTGgtt                         >  2:562432/1‑143 (MQ=255)
                                                                                                            aaCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCGGTTTa                        >  1:111222/1‑143 (MQ=255)
                                                                                                              cTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAc                      <  2:207100/143‑1 (MQ=255)
                                                                                                                ggAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCGCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGGAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCa                    <  2:212626/143‑1 (MQ=255)
                                                                                                                 gAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAg                   >  2:455503/1‑143 (MQ=255)
                                                                                                                   aGGCAGCGCCGCACCGCCATGGTGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGtt                 <  2:621100/143‑1 (MQ=255)
                                                                                                                         cgccgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTATACAAGCTCTTCGGTACTGCCAGAACCTCGCTGTTTAACCAGTTCGataa           >  1:571717/1‑143 (MQ=255)
                                                                                                                                ccGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGtt    <  1:881944/143‑1 (MQ=255)
                                                                                                                                  gCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTACGAGAGTACCGTagt  >  2:583373/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                   
CGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGT  >  NC_000913/3560315‑3560586

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: