Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,790,027 C→A 13.4% H113N (CAT→AAT)  uhpA → DNA‑binding response regulator in two‑component regulatory system wtih UhpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,790,0270CA13.4% 111.3 / 3.9 37H113N (CAT→AAT) uhpADNA‑binding response regulator in two‑component regulatory system wtih UhpB
Reads supporting (aligned to +/- strand):  ref base C (14/18);  new base A (3/2);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

GCTTAACGCGGGGGCACGCGGCTTTCTTTCCAAACGCTGTAGCCCGGATGAACTCATTGCTGCGGTGCATACGGTTGCCACGGGCGGCTGTTATCTGACGCCGGATATTGCCATTAAACTGGCATCCGG  >  W3110S.gb/3789960‑3790088
                                                                   |                                                             
gcTTAACGCGGGGGCACGCGGCTTTCTTTCCAAACGCTGTAGCCCGGATGAACTCATTGCTGCGGTGCa                                                              <  1:3064988/69‑1 (MQ=255)
 cTTAACGCGGGGGCACGCGGCTTTCTTTCCAAACGCTGTAGCCCGGATGAACTCATTGCTGCGGTGCATAc                                                           <  1:2890754/71‑1 (MQ=255)
          ggggCACGCGGCTTTCTTTCCAAACGCTGTAGCCCGGATGAACTCATTGCTGCGGTGCATACGGTTGCCAc                                                  <  1:3472276/71‑1 (MQ=255)
                     ctttctttCCAAACGCTGTAGCCCGGATGAACTCATTGCTGCGGTGCa                                                              <  1:2191621/48‑1 (MQ=255)
                          tttcCAAACGCTGTAGCCCGGATGAACTCATTGCTGCGGTGCATACGGTTGCCACGGGCGGCTGTTa                                      <  1:556393/67‑1 (MQ=255)
                          tttcCAAACGCTGTAGCCCGGATGAACTCATTGCTGCGGTGCATACGGTTGCCACGGGCGGCTGTTa                                      <  1:3068166/67‑1 (MQ=255)
                          tttcCAAACGCTGTAGCCCGGATGAACTCATTGCTGCGGTGCATACGGTTGCCACGGGCGGCTGTTa                                      <  1:2513949/67‑1 (MQ=255)
                           ttcCAAACGCTGTAGCCCGGATGAACTCATTGCTGCGGTGCATACGGt                                                        <  1:1004987/48‑1 (MQ=255)
                           ttcCAAACGCTGTAGCCCGGATGAACTCATTGCTGCGGTGCATACGGt                                                        <  1:2572852/48‑1 (MQ=255)
                           ttcCAAACGCTGTAGCACGGATGAACTCATTGCTGCGGTGCATACGGt                                                        <  1:4108336/48‑1 (MQ=255)
                                 aCGCTGTAGCCCGGATGAACTCATTGCTGCGGTGCATACGGTTGCCACGGGCGGCTGTTATCTGACGCCgg                           <  1:4159580/71‑1 (MQ=255)
                                  cGCTGTAGCCCGGATGAACTCATTGCTGCGGTGCa                                                              >  1:1542442/1‑35 (MQ=255)
                                  cGCTGTAGCCCGGATGAACTCATTGCTGCGGTGCa                                                              >  1:3525813/1‑35 (MQ=255)
                                    cTGTAGCCCGGATGAACTCATTGCTGCGGTGAATACGGTTGCCACGGGCGGCTGTTATCTGACGCCGGata                        <  1:1576683/71‑1 (MQ=255)
                                       tAGCCCGGATGAACTCATTGCTGCGGTGCATACGGTTGCCACGGGCGGCTGTTATCTGACGCCGGATATTg                     <  1:1275915/71‑1 (MQ=255)
                                            cGGATGAACTCATTGCTGCGGTGCATACGGTTGCCACGGGCGGCTGTTATCTGACGCCGGa                          <  1:242560/61‑1 (MQ=255)
                                             ggATGAACTCATTGCTGCGGTGCATACGGTTGCCACGGGCGGCTGTTATc                                    >  1:1808193/1‑50 (MQ=255)
                                             ggATGAACTCATTGCTGCGGTGCATACGGTTGCCACGGGCGGCTGTTATCTGACGCCgg                           <  1:1668181/59‑1 (MQ=255)
                                             ggATGAACTCATTGCTGCGGTGCATACGGTTGCCACGGGCGGCTGTTATCTGACGCCgg                           <  1:2444478/59‑1 (MQ=255)
                                             ggATGAACTCATTGCTGCGGTGCATACGGTTGCCACGGGCGGCTGTTATCTGACGCCGGATATTGCCAtt                >  1:4528484/1‑70 (MQ=255)
                                             ggATGAACTCATTGCTGCGGTGCATACGGTTGCCACGGGCGGCTGTTATCTGACGCCGGATATTGCCATTa               <  1:51448/71‑1 (MQ=255)
                                              gATGAACTCATTGCTGCGGTGCATACGGTTGCCACGGGCGGCTGTTATCTGACGCCGGATATTGCCATTaa              <  1:4647864/71‑1 (MQ=255)
                                                tGAACTCATTGCTGCGGTGCATACGGTTGCCACGGGCGGCTGTTATCTGACGCCGGATATTGCCATTaaa             >  1:787904/1‑70 (MQ=255)
                                                tGAACTCATTGCTGCGGTGCATACGGTTGCCACGGGCGGCTGTTATCTGACGCCGGATATTGCCATTAAAc            >  1:3772466/1‑71 (MQ=255)
                                                  aaCTCATTGCTGCGGTGCATACGGTTGCCACGGGCGGCTGTTATCTGACGCCGGATATTGCCATTAAACTg          >  1:4358126/1‑71 (MQ=255)
                                                  aaCTCATTGCTGCGGTGAATACGGTTGCCACGGGCGGCTGTTATCTGACGCCGGATATTGCCATTAAACTg          <  1:4361734/71‑1 (MQ=255)
                                                      cATTGCTGCGGTGCATACGGTTGCCACGGGCGGCTGTTATCTGACGCCGGata                        >  1:3918906/1‑53 (MQ=255)
                                                      cATTGCTGCGGTGCATACGGTTGCCACGGGCGGCTGTTATCTGACGCCGGata                        >  1:3929702/1‑53 (MQ=255)
                                                      cATTGCTGCGGTGCATACGGTTGCCACGGGCGGCTGTTATCTGACGCCGGata                        >  1:2021348/1‑53 (MQ=255)
                                                       aTTGCTGCGGTGCATACGGTTGCCACGGGCGGCTGTTa                                      >  1:3398492/1‑38 (MQ=255)
                                                       aTTGCTGCGGTGCATACGGTTGCCACGGGCGGCTGTTa                                      >  1:1232286/1‑38 (MQ=255)
                                                       aTTGCTGCGGTGCATACGGTTGCCACGGGCGGCTGTTATCTGACGCCGGATATTGCCATTAAACTGGCATc     >  1:5090720/1‑71 (MQ=255)
                                                       aTTGCTGCGGTGAATACGGTTGCCACGGGCGGCTGTTATCTGa                                 >  1:1115233/1‑43 (MQ=255)
                                                       aTTGCTGCGGTGAATACGGTTGCCACGGGCGGCTGTTATCTGACGCCGGATATTGCCATTAAACTGGCAt      >  1:3550722/1‑70 (MQ=255)
                                                          gctgcGGTGCATACGGTTGCCACGGGCGGCTGTTATCTGACGCCGGATATTCCCATTAAACTGGCATCCgg  >  1:1523905/1‑71 (MQ=255)
                                                          gctgcGGTGAATACGGTTGCCACGGGCGGCTGTTATCTGACGCCGGATATTGCCATTAAACTGGCATCCgg  >  1:1141720/1‑71 (MQ=255)
                                                                gTGCATACGGTTGCCACGGGCGGCTGTTATCTGACGCCGGATATTGCCATTAAACTGGCATCCg   <  1:2003314/64‑1 (MQ=255)
                                                                   cATACGGTTGCCACGGGCGGCTGTTATCTGACGCCGGATATTGCCATTa               >  1:3340847/1‑49 (MQ=255)
                                                                   |                                                             
GCTTAACGCGGGGGCACGCGGCTTTCTTTCCAAACGCTGTAGCCCGGATGAACTCATTGCTGCGGTGCATACGGTTGCCACGGGCGGCTGTTATCTGACGCCGGATATTGCCATTAAACTGGCATCCGG  >  W3110S.gb/3789960‑3790088

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: