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.G94.9% 112.1 / ‑0.2 39pseudogene (151/758 nt)glpRDNA‑binding transcriptional repressor GlpR
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (19/18);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTG  >  NC_000913/3560323‑3560587
                                                                                                                                       |                                                                                                                                    
cgcggcgAGCCTTGGGGGGTTTGCGCCGACTCGCTGCGGGCGAGGCGCTCTCTGTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAg                                                                                                                             <  2:614577/142‑1 (MQ=255)
    tCGAGCCATTGGCGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGttt                                                                                                                         <  1:129401/143‑1 (MQ=255)
         ccATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTc                                                                                                                    >  1:548677/1‑143 (MQ=255)
                 ggtttGCCCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCCGGCTGGCCTTGCGATCGTGCCACCGCGTCTTAACCGAACTGGAACGCTGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCt                                                                                                            <  1:53948/141‑1 (MQ=255)
                  aTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTGGGGATCGTGCCTGGGCGAGTTGACCGCGCGGGACTGCAGCGCCGCACCGCCATGATGGCGCGGGATCGTGCTTTGCTCCGCCAGccc                                                                                                           >  2:151659/1‑141 (MQ=255)
                    ttGCTCCGCCACGTGGCGGGCGATGCGCTCGTTTTCTTCGGTCTGGGTGGGCTGGCGAGCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCAt                                                                                                         <  1:1355/143‑1 (MQ=255)
                       cTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGa                                                                                                      >  1:325481/1‑143 (MQ=255)
                        tCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGGGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTGATGGaa                                                                                                     >  1:362501/1‑142 (MQ=255)
                          cGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGt                                                                                                   >  1:206827/1‑143 (MQ=255)
                           gCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATGAGGTTTTGCTCCGCCAGCTCATTGAGGTc                                                                                                  >  1:166666/1‑143 (MQ=255)
                              aCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCGGGTTTTGCTCCGCGCGCTCACGGAGGTcgcg                                                                                               >  1:503633/1‑143 (MQ=255)
                               cTTTGCCGGCCATGCGCTCCTTTTCTTCCGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCACGTTTTGCCCCGCCAGCTCATTGAGGTCgcgg                                                                                              <  2:408277/143‑1 (MQ=255)
                                 ttGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggcg                                                                                            >  1:533556/1‑143 (MQ=255)
                                  tGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGa                                                                                           <  1:463365/143‑1 (MQ=255)
                                     gggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATa                                                                                        <  2:183428/143‑1 (MQ=255)
                                      ggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGAGCAGGTTGTGCTCCGCCAGCTCATTGAGGTCGCGGCGAGTAg                                                                                       >  2:523352/1‑143 (MQ=255)
                                             cgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAGGGCGCAGGATCAGGGTTTGCTCCGCCAGCGCATTGAGGGCGCGGCGAAGAGGCTGCgg                                                                                >  1:178039/1‑143 (MQ=255)
                                              gcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGGTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTATGAggg                                                                               >  1:131819/1‑143 (MQ=255)
                                                 ctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGGCGCGGCGAATAGTCTGCGGGCTg                                                                            >  1:471796/1‑143 (MQ=255)
                                                   tgggctgcGGGCTGGGTGGCCGTGCGGTCGTGCCGCGGCGTGGTAGCGGAAGTGGAAGGCAGCGGCGCACGGCCATGGGGGGGCAGGAGCAGGTTTTGCTCCGCCCGCTCATTGGGGTCGCGGCGAATAGTCTGCGGGCTGAc                                                                          <  1:596659/136‑1 (MQ=255)
                                                    tttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACg                                                                         <  1:476268/143‑1 (MQ=255)
                                                      gcttcGGTCTGGGTGGCCTTGCGCTCGTGCCTCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGga                                                                       <  2:417180/142‑1 (MQ=255)
                                                         tcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaa                                                                    <  1:325258/143‑1 (MQ=255)
                                                          cGGTCTGGGTGGCCTTGCGCTCGTGGCACAACGCGTTAACCGAACTGGAACGCGGCGCCGCACCGCCGTGGTGGCGGAGGAGCAGGGTTTGCGCCGCCAGCCCAGTCGGATCGGGGGGAATCGTCTGCCGGCCGACAgacgag                                                                   >  2:249844/1‑139 (MQ=255)
                                                           ggTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAAt                                                                  <  1:219735/143‑1 (MQ=255)
                                                              cTGGGTGGCTTTTCGATGGTGCCACGGCGTGTTAACTGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCATGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGct                                                               <  2:169075/143‑1 (MQ=255)
                                                                  ggggCCTTGCGGCCGTGTCACGGCGTGTTAACCGCGCTGGAAGGCAGCGCCGCACCGTCCTGATGGGGCTGGATCAGGTTTTGCTCCGCGAGCTCTTTGCGGACGCCGCTGTTGCTCTGCGGGCTGACTGAGAAATGCTCTAc                                                           <  2:445080/141‑1 (MQ=255)
                                                                           cGATGGGGCCCCGGCGTGGTGACCGGACTGGGACGCAGCGCGGCACCGCCATGATGGCGCACAATCAGGTTTTACTCCGTCGGTGTGTTGAGTTCGCGCCCAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTc                                                  <  2:267752/143‑1 (MQ=255)
                                                                             aTCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCgg                                                <  2:206827/143‑1 (MQ=255)
                                                                                   ccACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGGGAGACGCTCTGCCAGCTCTTCGGTACTGa                                          >  2:65638/1‑143 (MQ=255)
                                                                                    cACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGAc                                         <  1:398966/143‑1 (MQ=255)
                                                                                     aCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACa                                        <  1:194820/143‑1 (MQ=255)
                                                                                     aCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGAGCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACa                                        >  2:128434/1‑143 (MQ=255)
                                                                                                         gAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCGGGATTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAg                    >  1:345486/1‑143 (MQ=255)
                                                                                                         gAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCTGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAg                    <  2:325481/143‑1 (MQ=255)
                                                                                                            tgCAGCGCCGCGCCGCCATGCTTGCGCCGTATCCGGTTTTGCGCCGCCAGCTCATTTAGGTCGCGGCGCACAGTCGGCTGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAGCCCTGCTGTTTAACCAGTTc                 <  2:662583/142‑1 (MQ=255)
                                                                                                                 cgccgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGagaa            >  2:186009/1‑143 (MQ=255)
                                                                                                                     gcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTAGTGACCTAACCCTGCTGTTTCACCCGTTCGATAATAcc        >  2:425003/1‑143 (MQ=255)
                                                                                                                           ccATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtgtg  >  1:543380/1‑143 (MQ=255)
                                                                                                                                       |                                                                                                                                    
AGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTG  >  NC_000913/3560323‑3560587

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: