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
                                                                                                                                        |                                                                                                                                         
gATCTGGCAGAAGAAAAGCGACGCGCGGTGGATAAAGGGGGTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGGCTTTAc                                                                                                                                    <  2:238406/143‑1 (MQ=255)
        agaagaTATGCGCAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATGGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGGCTTTACGCACtgtc                                                                                                                            >  2:482147/1‑142 (MQ=255)
                 gcgcAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCGCTGCTCTCTCAGGCGCGGGAAAAGGGCCGTGAGAACCATAGATTAACCTTTTGTGGCGAACCGCTGCAACAGACGCAAGGCATTACGCACTGTGTCTTGgagc                                                                                                                   >  1:1729/1‑143 (MQ=255)
                  cgcAGCGCGATGGATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGGCTTTACGCACTGTGGCTTGGcgcg                                                                                                                  >  1:124327/1‑143 (MQ=255)
                              gATAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGcc                                                                                                      >  2:58481/1‑143 (MQ=255)
                                tAAAGGTGCTCGCGGCGTGATCGCCCTGCTCTCTCAGGCGCTGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCaa                                                                                                    >  1:430680/1‑143 (MQ=255)
                                           gcggcgTGATCGCCCTGTTCTCTCCGGCGCTGGAAAATGGCCGGGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCa                                                                                         <  1:281632/143‑1 (MQ=255)
                                                            ctctctCAGGCGCTGGAAAATGGCCGTGAGAACCATTGCTCCACCTATCGTAGCGAACCTCTGCAACAGACACAAAGCTCTCCCCACTGTCGCTTGACGTGAATCTGCAGTCCCAAATTTCACGCAGATTctacaaccgtgaa                                                                        >  1:367466/1‑131 (MQ=255)
                                                                     gcgcTGGAACATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTg                                                               <  2:429397/143‑1 (MQ=255)
                                                                         tGGAAAATGGCCGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGccc                                                           >  1:489643/1‑143 (MQ=255)
                                                                             aaaTGGCCGTGAGAACCATTGTTGAACCTTTTGTGGCGAACCGCGGCAACAGGCACAAGGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGt                                                       <  2:315989/143‑1 (MQ=255)
                                                                                   ccGTGAGAACCATTGTTTAACCTTTTGTGGCGAACCGCTGCAACAGGCACAAGGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCCCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTaaaaaa                                                 >  2:306530/1‑143 (MQ=255)
                                                                                   cagTGAGAACCATGGGTAAACCATTTGGGGCGAACCGATGCAACAGGCAGAAGGCTTTACGGACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGAAGTTTCGAGCCACTGGAAAACGCGCTGGAACATGTaaaaaa                                                 <  1:131413/141‑1 (MQ=255)
                                                                                             cATTGTTTGGCGTATAGTGTCGGACGGCTCCAACAGGCACAAGGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCg                                       <  2:137566/143‑1 (MQ=255)
                                                                                                     aaCCTTTTGTGGCGAACCGCTGCAACAGGCACAAGGCTTTACGCACTGTGGCTTGGCACGAATCTGCAGGCCCAAATTTCGCGCAGTTTCGAGCCACTGGAAACCGCGCTGGCCCATGTAACAAACATTATTCCGACGCCCGc                               >  1:498400/1‑143 (MQ=255)
                                                                                                           tgggggCGAACCGCTGCAACAGGCACAAGGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTa                         <  1:17508/139‑1 (MQ=255)
                                                                                                            tgtgGCGAACCGCTGCAACAGGCACAAGGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCCCAGTTTCGAGCCACTGCAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAg                        <  2:517017/143‑1 (MQ=255)
                                                                                                               ggCGAACCGCGGCAACAGGCACAAGGCTTTACGCACTGTGGCTGGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAg                     <  1:120674/143‑1 (MQ=255)
                                                                                                                      cGCTGCAACAGCCACAAGGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCAtttt              <  2:375238/143‑1 (MQ=255)
                                                                                                                       tcgacAACATGCACACGGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCAttttt             <  1:408896/139‑1 (MQ=255)
                                                                                                                         tGCAACAGGCACAAGGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTa           <  2:457274/143‑1 (MQ=255)
                                                                                                                             aCAGGCACAAGGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTCGAGCCACGGGAAAACGCGCTGGCCCATGTAAAAGACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTGTCTCAc       >  2:52795/1‑143 (MQ=255)
                                                                                                                             aCAGGCACAAGGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGATTCGAGCCACTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTATCAc       <  2:261727/143‑1 (MQ=255)
                                                                                                                                ggCACAAGGCTTTACGCACTGTGGCTTGGCGCGAATCTGCAGGCCAAAATTTCGCGCAGTTTACACGCGCTGGAAAACGCGCTGGCCCATGTAAAAAACATTATTGCGACGCCTGCCGTTTAGCAGGCATTTTTTGTCACCAg    >  1:86153/1‑143 (MQ=255)
                                                                                                                                  cacaAGGCTTTACGCACTGTGGCTTGGCGCGACTCTACAGGCCAAAATTTCACGCACTTTCAAGCCACTACAAAATGCACTGGAGCATGTATAAAACATAATTACAACGCCAGGCGTCTAGCAATCATTAATTGTAACCAgcg  >  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: