Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 767,700 T→C 57.5% V444V (GTA→GTG phoQ ← sensory histidine kinase in two‑compoent regulatory system with PhoP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE767,7000TC57.5% 15.8 / 37.2 54V444V (GTA→GTGphoQsensory histidine kinase in two‑compoent regulatory system with PhoP
Reads supporting (aligned to +/- strand):  ref base T (11/12);  new base C (7/24);  total (18/36)
Fisher's exact test for biased strand distribution p-value = 7.98e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

CTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTC  >  minE/767638‑767758
                                                              |                                                          
ctctctCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCt                                                      >  1:292073/1‑69 (MQ=255)
ctctctCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCt                                                      >  1:2842802/1‑69 (MQ=255)
ctctctCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCt                                                      >  1:1481569/1‑69 (MQ=255)
         gCGACTATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCCACACATTGCCCAggg                                             <  1:1936500/69‑1 (MQ=255)
         gCGACGATTTTACCCTGATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCCACACATTGCCCAggg                                             <  1:3331286/69‑1 (MQ=255)
         gCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCCACACCTTGCCCAggg                                             <  1:1852398/69‑1 (MQ=255)
         gCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCCACACCTTGCCCAggg                                             <  1:2919484/69‑1 (MQ=255)
         gCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCCACACCTTGCCCAggg                                             <  1:2930585/69‑1 (MQ=255)
         gCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCCACACCTTGCCCAggg                                             <  1:2523373/69‑1 (MQ=255)
         gCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCCACACCTTGCCCAggg                                             <  1:2384580/69‑1 (MQ=255)
         gCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCCACACCTTGCCCAggg                                             <  1:3147515/69‑1 (MQ=255)
         gCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCCACACCTTGCCCAggg                                             <  1:2119091/69‑1 (MQ=255)
         gCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCCACACCTTGCCCAggg                                             <  1:1743112/69‑1 (MQ=255)
         gCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCCACACCTTGCCCAggg                                             <  1:923320/69‑1 (MQ=255)
         gCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCCACACCTTGCCCAggg                                             <  1:1179456/69‑1 (MQ=255)
         gCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCCACACCTTCCCAggg                                             <  1:1455970/68‑1 (MQ=255)
         gCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCCACACATTGCCCAggg                                             <  1:817188/69‑1 (MQ=255)
         gCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCCACACATTGCCCAggg                                             <  1:167835/69‑1 (MQ=255)
         gCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCCACACATTGCCCAggg                                             <  1:411460/69‑1 (MQ=255)
         gCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCCACACATTGCCCAggg                                             <  1:1544310/69‑1 (MQ=255)
         gCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCCACACATTGCCCAggg                                             <  1:1932430/69‑1 (MQ=255)
         gCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCCACACATTGCCCAggg                                             <  1:1458078/69‑1 (MQ=255)
         gCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCCACACATTGCCCAggg                                             <  1:647005/69‑1 (MQ=255)
         gCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCCACACATTGCCCAggg                                             <  1:2289915/69‑1 (MQ=255)
         gCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCCACACATTGCCCAggg                                             <  1:838607/69‑1 (MQ=255)
          cgacgaTTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCCACACATTGCCCAggg                                             <  1:1531698/68‑1 (MQ=255)
          cgacgaTTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCCACACATTGCCCAggg                                             <  1:2441146/68‑1 (MQ=255)
             cgaTTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTa                                         <  1:3307582/69‑1 (MQ=255)
             cgaTTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTa                                         <  1:2473306/69‑1 (MQ=255)
             cgaTTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTa                                         <  1:2963296/69‑1 (MQ=255)
             cgaTTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTa                                         <  1:1357135/69‑1 (MQ=255)
             cgaTTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTa                                         <  1:558807/69‑1 (MQ=255)
                ttttACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAg                                      <  1:1651353/69‑1 (MQ=255)
                ttttACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAg                                      <  1:2547857/69‑1 (MQ=255)
                ttttACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAg                                      <  1:2936818/69‑1 (MQ=255)
                     cccTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATc                                  <  1:596688/68‑1 (MQ=255)
                     cccTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATc                                  <  1:1089754/68‑1 (MQ=255)
                     cccTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATc                                  <  1:2867358/68‑1 (MQ=255)
                              tGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTg                         >  1:1130891/1‑68 (MQ=255)
                              tGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTg                         >  1:2287313/1‑68 (MQ=255)
                              tGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTg                         >  1:1716255/1‑68 (MQ=255)
                                         ttCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGaa             <  1:1279072/69‑1 (MQ=255)
                                                gCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACg                   >  1:1195914/1‑56 (MQ=255)
                                                gCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGAcc      >  1:2295688/1‑69 (MQ=255)
                                                gCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGAc       >  1:1074847/1‑68 (MQ=255)
                                                gCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGAc       >  1:478938/1‑68 (MQ=255)
                                                gCTACCGCAAGCCCTACACCTTGCCCAGGGCGTAAAGTATc                                  >  1:887848/1‑41 (MQ=255)
                                                    ccGCCAGCCCCACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACctct   >  1:1602950/1‑68 (MQ=255)
                                                    ccGCCAGCCCCACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACctct   >  1:2434673/1‑68 (MQ=255)
                                                    ccGCCAGCCCCACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCtctc  >  1:3361402/1‑69 (MQ=255)
                                                    ccGCCAGCCCCACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCtctc  >  1:1402122/1‑69 (MQ=255)
                                                    ccGCCAGCCCCACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCtctc  >  1:2156779/1‑69 (MQ=255)
                                                    ccGCCAGCCCCACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCtctc  >  1:2434050/1‑69 (MQ=255)
                                                    ccGCCAGCCCCACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCtctc  >  1:268217/1‑69 (MQ=255)
                                                              |                                                          
CTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTC  >  minE/767638‑767758

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: