Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,354,061 C→T 7.5% G218G (GGG→GGA dcuR ← DNA‑binding response regulator in two‑component regulatory system with DcuS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,354,0610CT 7.5% 164.7 / 3.3 53G218G (GGG→GGAdcuRDNA‑binding response regulator in two‑component regulatory system with DcuS
Reads supporting (aligned to +/- strand):  ref base C (26/23);  new base T (2/2);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

CCGTTATTGGCAATATTGTTTCAGTAGTGAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTATT  >  W3110S.gb/4353992‑4354125
                                                                     |                                                                
ccGTTATTGGCAATATTGTTTCAGTAGTGAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGcc                                                                 >  1:1005130/1‑71 (MQ=255)
     aTTGGCAATATTGTTTCAGTAGTGAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGt                                                             <  1:331823/70‑1 (MQ=255)
     aTTGGCAATATTGTTTCAGTAGTGAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGt                                                             >  1:3117771/1‑70 (MQ=255)
     aTTGGCAATATTGTTTCAGTAGTGAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGt                                                             >  1:4172439/1‑70 (MQ=255)
            atatTGTTTCAGTAGTGAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGTCCCGTGACGCCAt                                                     <  1:1887845/71‑1 (MQ=255)
            atatTGTTTCAGTAGTGAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCAt                                                     <  1:4248137/71‑1 (MQ=255)
             tatTGTTTCAGTAGTGAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCAt                                                     >  1:2437324/1‑70 (MQ=255)
             tatTGTTTCAGTAGTGAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCAt                                                     <  1:1860181/70‑1 (MQ=255)
             tatTGTTTCAGTAGTGAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATa                                                    <  1:2825292/71‑1 (MQ=255)
                 gTTTCAGTAGTGAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATg                                                 <  1:538063/70‑1 (MQ=255)
                        tagtgagtagtgTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTg                                          <  1:2623679/70‑1 (MQ=255)
                         agtgagtagtgTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGAt                                              <  1:2971686/65‑1 (MQ=255)
                             agtagtgTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATaa                                                   >  1:2079369/1‑56 (MQ=255)
                                  tgtTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGa                                                           >  1:941544/1‑43 (MQ=255)
                                  tgtTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGa                                                           >  1:164708/1‑43 (MQ=255)
                                  tgtTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATaa                                                   >  1:3883381/1‑51 (MQ=255)
                                  tgtTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATa                                                    <  1:3488771/50‑1 (MQ=255)
                                  tgtTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATg                                                 <  1:215965/53‑1 (MQ=255)
                                  tgtTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACa                                   <  1:154985/67‑1 (MQ=255)
                                  tgtTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGAt                               >  1:1538015/1‑71 (MQ=255)
                                    ttCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGAtgt                             >  1:3129270/1‑71 (MQ=255)
                                      cTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCAt                                                     >  1:2426907/1‑45 (MQ=255)
                                      cTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGa                                >  1:3467321/1‑66 (MQ=255)
                                      cTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTgg                           >  1:2099502/1‑71 (MQ=255)
                                       tGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGa                                      >  1:170535/1‑59 (MQ=255)
                                         ccTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATa                                             >  1:44288/1‑50 (MQ=255)
                                         ccTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGc                          <  1:4506097/69‑1 (MQ=255)
                                          cTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAg                        <  1:445314/70‑1 (MQ=255)
                                          cTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAg                        <  1:318369/70‑1 (MQ=255)
                                           tGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTgg                                         >  1:4471717/1‑52 (MQ=255)
                                            gAATACGGTAACGGTAAACTGGACGTCCCGTGACGCAATAATGGATACTGGTGAACAAGATGTGGCAGtt                      >  1:3498050/1‑70 (MQ=255)
                                            gAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATa                                                    <  1:1367288/40‑1 (MQ=255)
                                             aaTACGGTAACGGTAAACTGGACGTCCCGTGACGCCATAATGGATACTGGTGAACAAGAtg                              >  1:3510339/1‑61 (MQ=255)
                                             aaTACGGTAACGGTAAACTGGACGCCCCGTGACGCCAt                                                     <  1:781996/38‑1 (MQ=255)
                                             aaTACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCa                         <  1:3067829/66‑1 (MQ=255)
                                              aTACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGTACAAGATGTGGCAGTTg                     <  1:2379276/69‑1 (MQ=255)
                                               tACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGa                    >  1:641636/1‑69 (MQ=255)
                                                acggtaacggtaAACTGGACGCCCCGTGACGCCATAATg                                                 <  1:210466/39‑1 (MQ=255)
                                                 cggtaacggtaaACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAg                        >  1:1902121/1‑63 (MQ=255)
                                                   gtaacggtaaACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGAccagc               >  1:159818/1‑70 (MQ=255)
                                                    taacggtaaACTGGACGCCCCGTGACTCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGAccagcca             <  1:3364557/71‑1 (MQ=255)
                                                     aacggtaaACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGAccagc               >  1:3528762/1‑68 (MQ=255)
                                                      acggtaaACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGAccagccag            <  1:4492900/70‑1 (MQ=255)
                                                         gtaaACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGa           >  1:1838582/1‑68 (MQ=255)
                                                          taaACTGGACGCCCCGTGACGCCATAATGGATACTGGTTAACAAGATGTGGCAGTTGACCAGCCAGATGAg       <  1:926536/71‑1 (MQ=255)
                                                          taaACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGAt          >  1:192337/1‑68 (MQ=255)
                                                           aaaCTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAg                                 >  1:780366/1‑44 (MQ=255)
                                                            aaCTGGCCGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAgg      >  1:62086/1‑70 (MQ=255)
                                                            aaCTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAg       >  1:4637666/1‑69 (MQ=255)
                                                             aCTGGACGTCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGc                          <  1:510727/49‑1 (MQ=255)
                                                              cTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGAtg                              >  1:117031/1‑44 (MQ=255)
                                                              cTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGa           >  1:463150/1‑63 (MQ=255)
                                                               tGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTg                     <  1:605260/52‑1 (MQ=255)
                                                                 gACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTAtt  <  1:2675728/69‑1 (MQ=255)
                                                                     |                                                                
CCGTTATTGGCAATATTGTTTCAGTAGTGAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTATT  >  W3110S.gb/4353992‑4354125

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: