Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,726,507 Δ1 bp coding (485/600 nt) nemR → transcriptional repressor for the nemRA‑gloA operon, quinone‑, glyoxal‑, and HOCl‑activated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,726,5070T.96.0% 95.0 / ‑3.1 25coding (485/600 nt)nemRtranscriptional repressor for the nemRA‑gloA operon, quinone‑, glyoxal‑, and HOCl‑activated
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (11/13);  minor base A (1/0);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGAC  >  NC_000913/1726371‑1726644
                                                                                                                                        |                                                                                                                                         
gATCTGGCAGAAGAAAAGCGACGCGCGGTGGATAAAGGGGGTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAG‑GCTTTAc                                                                                                                                    <  2:238406/143‑1 (MQ=255)
        agaagaTATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATGGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAG‑GCTTTACGCACtgtc                                                                                                                            >  2:482147/1‑142 (MQ=255)
                 gcgcAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCGCTGCTCTCTCAGGCGCGGGAAAAGGGCCGTGAGAACCATAGATTAACCTTTTGTGGCGAACCGCTGCAACAGACGCAAG‑GCATTACGCACTGTGTCTTGgagc                                                                                                                   >  1:1729/1‑143 (MQ=255)
                  cgcAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAG‑GCTTTACGCACTGTGGCTTGGcgcg                                                                                                                  >  1:124327/1‑143 (MQ=255)
                              gATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAG‑GCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGcc                                                                                                      >  2:58481/1‑143 (MQ=255)
                                tAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAG‑GCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCaa                                                                                                    >  1:430680/1‑143 (MQ=255)
                                           gcggcgTGATCGCCCTGTTCTCTCCGGCGCTGGAAAATGGCCGGGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAG‑GCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCa                                                                                         <  1:281632/143‑1 (MQ=255)
                                                            ctctctCAGGCGCTGGAAAATGGCCGTGAGAACCATTGCTCCACCTATCGTAGCGAACCTCTGCAACAGACACAA‑AGCTCTCCCCACTGTCGCTTGACGTGAATCTGCAGTCCCAAATTTCACGCAGATTctacaaccgtgaa                                                                        >  1:367466/1‑131 (MQ=255)
                                                                     gcgcTGGAACATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAG‑GCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTg                                                               <  2:429397/143‑1 (MQ=255)
                                                                         tGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAG‑GCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGccc                                                           >  1:489643/1‑143 (MQ=255)
                                                                             aaaTGGCCGTGAGAACCATTGTTGAACCTTTTGTGGCGAACCGCGGCAACAGGCACAAG‑GCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGt                                                       <  2:315989/143‑1 (MQ=255)
                                                                                   ccGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAG‑GCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCCCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTaaaaaa                                                 >  2:306530/1‑143 (MQ=255)
                                                                                   cagTGAGAACCATGGGTAAACCATTTGGGGCGAACCGATGCAACAGGCAGAAG‑GCTTTACGGACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGAAGTTTCGAGCCACTGGAAAACGCGCTGGAACATGTaaaaaa                                                 <  1:131413/141‑1 (MQ=255)
                                                                                             cATTGTTTGGCGTATAGTGTCGGACGGCTCCAACAGGCACAAG‑GCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCg                                       <  2:137566/143‑1 (MQ=255)
                                                                                                     aaCCTTTTGTGGCGAACCGCTGCAACAGGCACAAG‑GCTTTACGCACTGTGGCTTGGCACGAATCTGCAGGCCCAAATTTCGCGCAGTTTCGAGCCACTGGAAACCGCGCTGGCCCATGTAACAAACATTATTCCGACGCCCGc                               >  1:498400/1‑143 (MQ=255)
                                                                                                           tgggggCGAACCGCTGCAACAGGCACAAG‑GCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTa                         <  1:17508/139‑1 (MQ=255)
                                                                                                            tgtgGCGAACCGCTGCAACAGGCACAAG‑GCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCCCAGTTTCGAGCCACTGCAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAg                        <  2:517017/143‑1 (MQ=255)
                                                                                                               ggCGAACCGCGGCAACAGGCACAAG‑GCTTTACGCACTGTGGCTGGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAg                     <  1:120674/143‑1 (MQ=255)
                                                                                                                      cGCTGCAACAGCCACAAG‑GCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCAtttt              <  2:375238/143‑1 (MQ=255)
                                                                                                                       tcgacAACATGCACACG‑GCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCAttttt             <  1:408896/139‑1 (MQ=255)
                                                                                                                         tGCAACAGGCACAAG‑GCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTa           <  2:457274/143‑1 (MQ=255)
                                                                                                                             aCAGGCACAAG‑GCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACGGGAAAACGCGCTGGCCCATGTAAAAGACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTGTCTCAc       >  2:52795/1‑143 (MQ=255)
                                                                                                                             aCAGGCACAAG‑GCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGATTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCAc       <  2:261727/143‑1 (MQ=255)
                                                                                                                                ggCACAAG‑GCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTACACGCGCTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTGTCACCAg    >  1:86153/1‑143 (MQ=255)
                                                                                                                                  cacaAG‑GCTTTACGCACTGTGGCTTGGCGCGACTCTACAGGCCAAAATTTCACGCACTTTCAAGCCACTACAAAATGCACTGGAGCATGTATAAAACATAATTACAACGCCAGGCGTCTAGCAATCATTAATTGTAACCAgcg  >  1:272652/1‑141 (MQ=255)
                                                                                                                                        |                                                                                                                                         
GATCTGTCAGAAGATATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGTGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCACCAGAC  >  NC_000913/1726371‑1726644

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑