Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 447,412 G→A 100% R102C (CGC→TGC)  kdpD ← fused sensory histidine kinase in two‑component regulatory system with KdpE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE447,4120GA81.0% 41.7 / 6.2 21R102C (CGC→TGC) kdpDfused sensory histidine kinase in two‑component regulatory system with KdpE
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base A (5/12);  total (6/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGGGATGACGGGAACCTGGCGCATTACTGTGCGCCAGTTCGTCCATTAAGATCAGCGCCGGGCGGCGGGCGAGGGCGGCATCGAGATCAAACTCGCTGATATGCCGCCCACGGTACGCCTGGCGTTTTAACGG  >  minE/447349‑447481
                                                               |                                                                     
tGGGATGACGGGAACCTGGCGCATTACTGTGCGCCAGTTCGTCCATTAAGATCAGCGCCGGGCggcgggc                                                                 <  1:3773906/70‑1 (MQ=255)
     tGACGGGAACCTGGCGCATTACTGTGCGCCAGTTCGTCCATTAAGATCAGCGCCGGgcagcg                                                                    <  1:3209299/62‑1 (MQ=255)
     tGACGGGAACCTGGCGCATTACTGTGCGCCAGTTCGTCCATTAAGATCAGCGCCGGgcagcg                                                                    <  1:3010971/62‑1 (MQ=255)
     tGACGGGAACCTGGCGCATTACTGTGCGCCAGTTCGTCCATTAAGATCAGCGCCGGgcagcg                                                                    <  1:2455633/62‑1 (MQ=255)
            aaCCTGGCGCATTACTGTGCGCCAGTTCGTCCATTAAGATCAGCGCCGGGCAGCGGGCGAGGGCGGCATc                                                     >  1:1282074/1‑70 (MQ=255)
             aCCTGGCGCATTACTGTGCGCCAGTTCGTCCATTAATATCAGCGCCGGGCGGCGGGCGAGggcggc                                                        <  1:2030142/66‑1 (MQ=255)
             aCCTGGCGCATTACTGTGCGCCAGTTCGTCCATTAAGATCAGCGCCGGGCAGCGGGCGAGGGCGGCATc                                                     <  1:2001751/69‑1 (MQ=255)
               cTGGCGCATTACTGTGCGCCAGTTCGTCCATTAAGATCAGCGCCGGGCAGCGGGCGAGGGCGGCATCga                                                   >  1:3545854/1‑69 (MQ=255)
                  gcgcATTACTGTGCGCCAGTTCGTCCATTAAGATCAGCGCAGGGCGGCGGGCGAGGGCGGCATCGAGATCa                                              <  1:983885/71‑1 (MQ=255)
                       ttACTGTGCGCCAGTTCGTCCATTAAGATCAGCGCCGggcagc                                                                     <  1:426626/43‑1 (MQ=255)
                         aCTGTGCGCCAGTTCGTCCATTAAGATCAGCGCCGGGCAGCGGGCGAGGGCGGCATCGAGATCAAACTCGc                                       >  1:1856608/1‑71 (MQ=255)
                             tgCGCCAGTTCGTCCATTAAGATCAGCGCCGGGCGGCGGGCGAGGGCGGCATCGAGATCAAACTCGCTGa                                    <  1:3288517/70‑1 (MQ=255)
                                     ttCGTCCATTAAGATCAGCGCCGGGCGGCGGGCGAGGGCGGCATCGAGATCAAACTCGCTGATATgccgcc                           >  1:364018/1‑71 (MQ=255)
                                        gTCCATTAAGATCAGCGCCGGGCAGCGGGCGAGGGCGGCATCGAGATCAAACTCGCTGATATGCCGCCCa                         <  1:3062706/70‑1 (MQ=255)
                                            aTTAAGATCAGCGCCGGGCAGCGGGCGAGGGCGGCATCGAGATCa                                              <  1:1808996/45‑1 (MQ=255)
                                            aTTAAGATCAGCGCCGGGCAGCGGGCGAGGGCGGCATCGAGATCa                                              <  1:1804088/45‑1 (MQ=255)
                                              tAAGATCAGCGCCGGGCAGCGGGCGAGGGCGTCATc                                                     <  1:709487/36‑1 (MQ=37)
                                              tAAGATCAGCGCCGGGCAGCGGGCGAGGGCGGCATCGAGATCa                                              <  1:1723331/43‑1 (MQ=255)
                                               aaGATCAGCGCCGGGCAGCGGGCGAGGGCGGCATc                                                     >  1:4258932/1‑35 (MQ=255)
                                               aaGATCAGCGCCGGGCAGCGGGCGAGGGCGGCATCGAGATCAAACTCGCTGATATGCCGCCCACGGTACGc                 >  1:3468738/1‑71 (MQ=255)
                                                  aTCAGCGCCGGGCAGCGGGCGAGGGCGGCATCGAGATCAAACTCGCTGATATGCCGCCCACGGTACGCCTg              <  1:2657086/71‑1 (MQ=255)
                                                  aTCAGCGCCGGGCAGCGGGCGAGGGCGGCATCGAGATCAAACTCGCTGATATGCCGCCCACGGTACGCCTg              <  1:1870240/71‑1 (MQ=255)
                                                   gcAGCGCCGGGCAGCGGGCGAGGGCGGCATCGAGATCAAACTCGCTGATATGCCGCCCACGGTACGCCTg              <  1:3066240/69‑1 (MQ=255)
                                                          cGGGCAGCGGGCGAGGGCGGCATCGAGATCAAACTCGCTGATATGCCGCCCACGGTACGCCTGGCGTTTTa      <  1:1974320/71‑1 (MQ=255)
                                                          cGGGCAGCGGGCGAGGGCGGCATCGAGATCAAACTCGCTGATATGCCGCCCACGGTACGCCTGGCGTTTTa      <  1:1751355/71‑1 (MQ=255)
                                                           gggCAGCGGGCGAGGGCGGCATCGAGATCAAACTCGCTGATATGCCGCCCACGGTACGCCTGGCGTTTTa      <  1:2469977/70‑1 (MQ=255)
                                                               ggcgggcgAGGGCGGCATCGAGATCAAACTCGCTGATATGCCGCCCACGGTACGCCTGGCGTTTTAACgg  >  1:1527779/1‑70 (MQ=255)
                                                               |                                                                     
TGGGATGACGGGAACCTGGCGCATTACTGTGCGCCAGTTCGTCCATTAAGATCAGCGCCGGGCGGCGGGCGAGGGCGGCATCGAGATCAAACTCGCTGATATGCCGCCCACGGTACGCCTGGCGTTTTAACGG  >  minE/447349‑447481

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: