Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,159,068 G→A 30.8% M193I (ATG→ATA oxyR → DNA‑binding transcriptional dual regulator OxyR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,159,0680GA30.8% 27.9 / 15.1 29M193I (ATG→ATAoxyRDNA‑binding transcriptional dual regulator OxyR
Reads supporting (aligned to +/- strand):  ref base G (7/13);  new base A (7/2);  total (14/15)
Fisher's exact test for biased strand distribution p-value = 5.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

GCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGC  >  NC_000913/4158928‑4159210
                                                                                                                                            |                                                                                                                                              
gcgcTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCt                                                                                                                                              <  2:289940/143‑1 (MQ=255)
              tttctaaGAATTCATTGAAGAGCAGGGGAATGAAGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACt                                                                                                                                <  1:42498/138‑1 (MQ=255)
                      cagtcattGGAGTGCCGGTGTTTGATGAGCCAATGTTGCTGGCTATCGATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATACTGGAAGATGGTCACTGTTTgcgc                                                                                                                        <  1:425880/140‑1 (MQ=255)
                               aaGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGACCGATCTGGCAGGGGAAAAACTGCTGATACTGGAAGATGGTCACTATTTGCGCGATCAGGCa                                                                                                               >  1:14752/1‑143 (MQ=255)
                                  tGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGTTCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGCTTGCGCGATCAGGCAATg                                                                                                            >  2:353821/1‑143 (MQ=255)
                                     cggggttCGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGTAAGATGGTTATTGTTTGCGCGATCAGGCAATGGGt                                                                                                         <  1:90311/139‑1 (MQ=255)
                                     cGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAGACTGCTGATACTGGAAGATGGTCACGATGAGCGCGGTCAGGCAATggga                                                                                                         >  2:68663/1‑142 (MQ=255)
                                          tttGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTg                                                                                                    >  2:45605/1‑143 (MQ=255)
                                             catgcGCAAAGGTGGCCGGGTATCTAGGAAGATCACCCGCGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATACTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGttt                                                                                                 <  1:271955/142‑1 (MQ=255)
                                                  cgctcggTTGCCGGCTAGCTACGCAGATCACCCGTGGGCGCACCGCGAATGCGTACCGATGGCCGATCCCGCAGGGGAAACACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAg                                                                                            <  2:178513/137‑1 (MQ=255)
                                                   ccAATGTGGCGGGGTATCGATGAAGATCAGCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGc                                                                                           <  1:394092/143‑1 (MQ=255)
                                                            cTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGCTGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCggg                                                                                  >  1:115442/1‑142 (MQ=255)
                                                                  gtcgatGAGGACCACCCGTGGCCGAACCGCGAAGGCGTACCGATGGCCGGTCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGa                                                                            <  2:247796/139‑1 (MQ=255)
                                                                  atctagGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGa                                                                            <  1:317102/143‑1 (MQ=255)
                                                                      atGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATACTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATaca                                                                        >  1:42719/1‑143 (MQ=255)
                                                                        gAATATCACGCGTGGGCGAACCGCGAATGCGTACCGATGGCCTTTCTGGCAGTGTAATTTTTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATacaca                                                                      <  1:17761/143‑1 (MQ=255)
                                                                         aaGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGCAAAACTGCTGATACTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATacacac                                                                     >  2:107818/1‑143 (MQ=255)
                                                                            aTCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATACTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTc                                                                  >  1:394658/1‑143 (MQ=255)
                                                                             tCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTcc                                                                 >  2:260581/1‑143 (MQ=255)
                                                                                   gTGGGCGAAACGCGAATGCGTACCGATGGCCGATCTGGCAGCGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGa                                                           <  1:40302/143‑1 (MQ=255)
                                                                                          aaCCGCGAATGCTTACCGATGGCCGATCTGGCAGGGCAAAAACTGCTGATGCTGGAGGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCt                                                    <  2:212559/143‑1 (MQ=255)
                                                                                           aCCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGCAAAACTGCTGACGCTGGCAGATGGTCACTGTTTGCGCGCTCAGGCAACGGGCGTCTGTCTTCAAGCCGGCGCGGATGAAGACACCCACCTCCGCCCGACCCGCCTg                                                   >  2:77006/1‑143 (MQ=255)
                                                                                               cgAATGCGTACCGATGGGCGATCTGGCAGGGGAAACAATGCTGATGCGATAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGaaa                                               <  2:457554/143‑1 (MQ=255)
                                                                                               cgAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGaaa                                               <  1:45605/143‑1 (MQ=255)
                                                                                               cgAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGGCACTGTTTGCGCGATCAGGCAATGGGTATCTGTTTTGAAGCCGGTGCGGATGAAGATACACACTTCCGCACGACCAGCCGGGaaa                                               >  2:29686/1‑143 (MQ=255)
                                                                                                         ccGACGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGACGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCaa                                     <  1:397406/143‑1 (MQ=255)
                                                                                                                  gATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCCGCGCTACAtgttgcc                            >  1:183313/1‑141 (MQ=255)
                                                                                                                          aGGGGAAAAACTGCTGATACTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTa                    >  1:331106/1‑143 (MQ=255)
                                                                                                                               aaaaaCTGCTGATACTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACGAGCCTGGAAACTATGCGCAACATGGTGGCGGCAAGTAGAggg               >  2:125614/1‑143 (MQ=255)
                                                                                                                                            gCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGAGCCCTTTACTGc  >  1:242824/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                              
GCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCTGGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCGGGGCGGATGAAGATACACACTTCCGCGCGACCAGCCTGGAAACTCTGCGCAACATGGTGGCGGCAGGTAGCGGGATCACTTTACTGC  >  NC_000913/4158928‑4159210

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: