Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,447,784 G→A 38.2% D64N (GAT→AAT)  atoS → sensory histidine kinase in two‑component regulatory system with AtoC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,447,7840GA38.2% 32.5 / 36.5 60D64N (GAT→AAT) atoSsensory histidine kinase in two‑component regulatory system with AtoC
Reads supporting (aligned to +/- strand):  ref base G (8/29);  new base A (7/16);  total (15/45)
Fisher's exact test for biased strand distribution p-value = 5.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.85e-01

CAGCAGTCTTATCTGAAAAAGAGAAAAAACTTTCTGCCGTGGTCAACCTGCTTAATCAGGCACTAGGCGATCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTATCCGCGCATTAAATGCAGAACT  >  minE/1447716‑1447848
                                                                    |                                                                
cagcagTCTTATCTGAAAAAGAGAAAAAACTTTCTGCCGTGGTCAACCTGCTTAATCAGGCACTAGGCg                                                                  >  1:3425345/1‑69 (MQ=255)
      tCTTATCTGAAAAAGAGAAAAAACTTTCTGCCGTTGTCAACCTGCTTAATCAGGCGCTAGGCAATCGCt                                                            <  1:1399983/69‑1 (MQ=255)
      tCTTATCTGAAAAAGAGAAAAAACTTTCTGCCGTGGTCAACCTGCTTAATCAGGCGCTAGGCAATCGCt                                                            <  1:2055609/69‑1 (MQ=255)
      tCTTATCTGAAAAAGAGAAAAAACTTTCTGCCGTGGTCAACCTGCTTAATCAGGCGCTAGGCAATCGCt                                                            <  1:2551978/69‑1 (MQ=255)
      tCTTATCTGAAAAAGAGAAAAAACTTTCTGCCGTGGTCAACCTGCTTAATCAGGCGCTAGGCAATCGCt                                                            <  1:2779483/69‑1 (MQ=255)
      tCTTATCTGAAAAAGAGAAAAAACTTTCTGCCGTGGTCAACCTGCTTAATCAGGCGCTAGGCAATCGCt                                                            <  1:1806693/69‑1 (MQ=255)
      tCTTATCTGAAAAAGAGAAAAAACTTTCTGCCGTGGTCAACCTGCTTAATCAGGCGCTAGGCAATCGCt                                                            <  1:2478782/69‑1 (MQ=255)
      tCTTATCTGAAAAAGAGAAAAAACTTTCTGCCGTGGTCAACCTGCTTAATCAGGCGCTAGGCAATCGCt                                                            <  1:3351607/69‑1 (MQ=255)
      tCTTATCTGAAAAAGAGAAAAAACTTTCTGCCGTGGTCAACCTGCTTAATCAGGCGCTAGGCAATCGCt                                                            <  1:3354451/69‑1 (MQ=255)
      tCTTATCTGAAAAAGAGAAAAAACTTTCTGCCGTGGTCAACCTGCTTAATCAGGCGCTAGGCAATCGCt                                                            <  1:353731/69‑1 (MQ=255)
      tCTTATCTGAAAAAGAGAAAAAACTTTCTGCCGTGGTCAACCTGCTTAATCAGGCGCTAGGCAATCGCt                                                            <  1:551666/69‑1 (MQ=255)
      tCTTATCTGAAAAAGAGAAAAAACTTTCTGCCGTGGTCAACCTGCTTAATCAGGCGCTAGGCAATCGCt                                                            <  1:1238489/69‑1 (MQ=255)
      tCTTATCTGAAAAAGAGAAAAAACTTTCTGCCGTGGTCAACCTGCTTAATCAGGCGCTAGGCAATCGCt                                                            <  1:2388758/69‑1 (MQ=255)
      tCTTATCTGAAAAAGAGAAAAAACTTTCTGCCGTGGTCAACCTGCTTAATCAGGCGCTAGGCAATCGCt                                                            <  1:216429/69‑1 (MQ=255)
      tCTTATCTGAAAAAGAGAAAAAACTTTCTGCCGTGGTCAACCTGCTTAATCAGGCGCTAGGCAATCGCt                                                            <  1:756819/69‑1 (MQ=255)
      tCTTATCTGAAAAAGAGAAAAAACTTTCTGCCGTGGTCAACCTGCTTAATCAGGCGCTAGGCAATCGCt                                                            <  1:889205/69‑1 (MQ=255)
       cTTATCTGAAAAAGAGAAAAAACTTTCTGCCGTGGTCAACCTGCTTAATCAGGCGCTAGGCAATCGCt                                                            <  1:3314625/68‑1 (MQ=255)
                 aaaGAGAAAAAACTTTCTGCCGTGGTCAACCTGCTTAATCAGGCACTAGGCGATCGCTATGAtctc                                                    <  1:2144073/66‑1 (MQ=255)
                 aaaGAGAAAAAACTTTCTGCCGTGGTCAACCTGCTTAATCAGGCACTAGGCGATCGCTATGAtctc                                                    <  1:2104064/66‑1 (MQ=255)
                 aaaGAGAAAAAACTTTCTGCCGTGGTCAACCTGCTTAATCAGGCACTAGGCGATCGCTATGAtctc                                                    <  1:2184144/66‑1 (MQ=255)
                 aaaGAGAAAAAACTTTCTGCCGTGGTCAACCTGCTTAATCAGGCACTAGGCGATCGCTATGAtctc                                                    <  1:2710889/66‑1 (MQ=255)
                 aaaGAGAAAAAACTTTCTGCCGTGGTCAACCTGCTTAATCAGGCACTAGGCGATCGCTATGAtctc                                                    <  1:1722125/66‑1 (MQ=255)
                 aaaGAGAAAAAACTTTCTGCCGTGGTCAACCTGCTTAATCAGGCACTAGGCGATCGCTATGAtctc                                                    <  1:1182749/66‑1 (MQ=255)
                 aaaGAGAAAAAACTTTCTGCCGTGGTCAACCTGCTTAATCAGGCACTAGGCGATCGCTATGAtctc                                                    <  1:673257/66‑1 (MQ=255)
                          aatcTTTCTGCCGTTGTCAACCTGCTTAATCAGGCACTAGGCGATCGCTATGATCTCTACATCGa                                            <  1:650007/62‑1 (MQ=255)
                            aCTTTCTGCCGTCGTCAACCTGCTTAATCAGGCACTAGGCGATCGCTATGATCTCTACATCGACTTAcc                                      <  1:89210/69‑1 (MQ=255)
                                           cAACCTGTTTAATCAGGCACTAGGCGATCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTa                        <  1:1147749/68‑1 (MQ=255)
                                           cAACCTGCTTAATCAGGCACTAGGCGATCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTa                        <  1:2179066/68‑1 (MQ=255)
                                           cAACCTGCTTAATCAGGCACTAGGCGATCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTa                        <  1:1003270/68‑1 (MQ=255)
                                           cAACCTGCTTAATCAGGCACTAGGCGATCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTa                        <  1:226009/68‑1 (MQ=255)
                                           cAACCTGCTTAATCAGGCACTAGGCGATCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTa                        <  1:360951/68‑1 (MQ=255)
                                           cAACCTGCTTAATCAGGCACTAGGCGATCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTa                        <  1:3388715/68‑1 (MQ=255)
                                           cAACCTGCTTAATCAGGCACTAGGCGATCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTa                        <  1:3199605/68‑1 (MQ=255)
                                           cAACCTGCTTAATCAGGCACTAGGCGATCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTa                        <  1:2694463/68‑1 (MQ=255)
                                           cAACCTGCTTAATCAGGCACTAGGCGATCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTa                        <  1:3130437/68‑1 (MQ=255)
                                            aaCCTGCTTAATCAGGCACTAGGCGATCGCTATGAtctc                                                    <  1:2413016/39‑1 (MQ=255)
                                            aaCCTGCTTAATCAGGCACTAGGCGATCGCTATGATCTCTAc                                                 >  1:2197810/1‑42 (MQ=255)
                                            aaCCTGCTTAATCAGGCACTAGGCGATCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTa                        <  1:781600/67‑1 (MQ=255)
                                            aaCCTGCTTAATCAGGCACTAGGCGATCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTa                        <  1:160689/67‑1 (MQ=255)
                                            aaCCTGCTTAATCAGGCACTAGGCGATCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTATc                      >  1:1019848/1‑69 (MQ=255)
                                            aaCCTGCTTAATCAGGCACTAGGCGATCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTATc                      >  1:1224480/1‑69 (MQ=255)
                                            aaCCTGCTTAATCAGGCACTAGGCGATCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTATc                      >  1:1935666/1‑69 (MQ=255)
                                            aaCCTGCTTAATCAGGCACTAGGCGATCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTATc                      >  1:1709733/1‑69 (MQ=255)
                                            aaCCTGCTTAATCAGGCACTAGGCGATCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTATc                      >  1:2416100/1‑69 (MQ=255)
                                            aaCCTGCTTAATCAGGCACTAGGCGATCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTATc                      >  1:2142192/1‑69 (MQ=255)
                                            aaCCTGCTTAATCAGGCACTAGGCGATCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTATc                      >  1:1313675/1‑69 (MQ=255)
                                                        cAGGCACTAGGCGATCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTATCCGCGCATTaaa           <  1:1262685/68‑1 (MQ=255)
                                                            cACTAGGCGGTCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTATCCGCGCATTAAATGCAg      <  1:2148414/69‑1 (MQ=255)
                                                            cACTAGGCGATCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTATCCGCGCATTAAATGCAg      <  1:1625605/69‑1 (MQ=255)
                                                            cACTAGGCGATCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTATCCGCGCATTAAATGCAg      <  1:369691/69‑1 (MQ=255)
                                                            cACTAGGCGATCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTATCCGCGCATTAAATGCAg      <  1:1821092/69‑1 (MQ=255)
                                                            cACTAGGCGATCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTATCCGCGCATTAAATGCAg      <  1:2051576/69‑1 (MQ=255)
                                                            cACTAGGCGATCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTATCCGCGCATTAAATGCAg      <  1:3079644/69‑1 (MQ=255)
                                                               tAGGCGATCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTa                        <  1:858987/48‑1 (MQ=255)
                                                                aGGCAATCGCTACGATCTCTACATCGACTTACCGCGTGAGGAGCGTATCCGCGCATTAAATGCAGAACt  >  1:2201940/1‑69 (MQ=255)
                                                                aGGCAATCGCTACGATCTCTACATCGACTTACCGCGTGAGGAGCGTATCCGCGCATTAAATGCAGAACt  >  1:2745736/1‑69 (MQ=255)
                                                                aGGCAATCGCTACGATCTCTACATCGACTTACCGCGTGAGGAGCGTATCCGCGCATTAAATGCAGAACt  >  1:1030651/1‑69 (MQ=255)
                                                                aGGCAATCGCTACGATCTCTACATCGACTTACCGCGTGAGGAGCGTATCCGCGCATTAAATGCAGAACt  >  1:386561/1‑69 (MQ=255)
                                                                aGGCAATCGCTACGATCTCTACATCGACTTACCGCGTGAGGAGCGTATCCGCGCATTAAATGCAGAACt  >  1:3277908/1‑69 (MQ=255)
                                                                aGGCAATCGCTACGATCTCTACATCGACTTACCGCGTGAGGAGCGTATCCGCGCATTAAATGCAGAACt  >  1:3076629/1‑69 (MQ=255)
                                                                aGGCAATCGCTACGATCTCTACATCGACTTACCGCGTGAGGAGCGTATCCGCGCATTAAATGCAGAACt  >  1:1916203/1‑69 (MQ=255)
                                                                    |                                                                
CAGCAGTCTTATCTGAAAAAGAGAAAAAACTTTCTGCCGTGGTCAACCTGCTTAATCAGGCACTAGGCGATCGCTATGATCTCTACATCGACTTACCACGTGAGGAGCGTATCCGCGCATTAAATGCAGAACT  >  minE/1447716‑1447848

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