Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 247,548 T→G 23.4% A397A (GCT→GCG phoR → sensory histidine kinase in two‑component regulatory system with PhoB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE247,5480TG23.4% 114.6 / 38.7 81A397A (GCT→GCGphoRsensory histidine kinase in two‑component regulatory system with PhoB
Reads supporting (aligned to +/- strand):  ref base T (22/40);  new base G (9/10);  total (31/50)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

AAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTT  >  minE/247490‑247614
                                                          |                                                                  
aaGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAatcatc                                                          >  1:1091678/1‑69 (MQ=255)
aaGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAatcatc                                                          >  1:1177246/1‑69 (MQ=255)
aaGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAatcatc                                                          >  1:72935/1‑69 (MQ=255)
aaGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAatcatc                                                          >  1:121115/1‑69 (MQ=255)
aaGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAatcatc                                                          >  1:1229424/1‑69 (MQ=255)
aaGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAatcatc                                                          >  1:335155/1‑69 (MQ=255)
aaGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAatcatc                                                          >  1:68342/1‑69 (MQ=255)
aaGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAatcatc                                                          >  1:1534427/1‑69 (MQ=255)
aaGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAatcatc                                                          >  1:651211/1‑69 (MQ=255)
   cgcgTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGa                                                                >  1:2247935/1‑60 (MQ=255)
   cgcgTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGa                                                                >  1:3080572/1‑60 (MQ=255)
   cgcgTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGa                                                                >  1:185692/1‑60 (MQ=255)
   cgcgTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGa                                                                <  1:1736618/60‑1 (MQ=255)
   cgcgTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGa                                                                <  1:1600041/60‑1 (MQ=255)
   cgcgTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGa                                                                >  1:452087/1‑60 (MQ=255)
   cgcgTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGa                                                                <  1:678683/60‑1 (MQ=255)
   cgcgTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGa                                                                <  1:470123/60‑1 (MQ=255)
           cGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCg                                                <  1:2748778/68‑1 (MQ=255)
                aaCCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGaa                                          <  1:1344675/69‑1 (MQ=255)
                aaCCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGaa                                          <  1:2223415/69‑1 (MQ=255)
                aaCCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGaa                                          <  1:2099995/69‑1 (MQ=255)
                aaCCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGaa                                          <  1:2702616/69‑1 (MQ=255)
                aaCCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGaa                                          <  1:1148175/69‑1 (MQ=255)
                aaCCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGaa                                          <  1:223491/69‑1 (MQ=255)
                aaCCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGaa                                          <  1:1782381/69‑1 (MQ=255)
                aaCCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGaa                                          <  1:3304018/69‑1 (MQ=255)
                aaCCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGaa                                          <  1:1237740/69‑1 (MQ=255)
                aaCCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGaa                                          <  1:1081732/69‑1 (MQ=255)
                aaCCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGaa                                          <  1:3406169/69‑1 (MQ=255)
                 aCCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGaa                                          <  1:781519/68‑1 (MQ=255)
                         tAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGAATATTGagag                                 <  1:3125502/69‑1 (MQ=255)
                         tAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGAATATTGagag                                 <  1:1083263/69‑1 (MQ=255)
                         tAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGAATATTGagag                                 <  1:2519969/69‑1 (MQ=255)
                           gCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTa                               <  1:1870308/69‑1 (MQ=255)
                             ggATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTAc                              <  1:3192294/68‑1 (MQ=255)
                             ggATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTAc                              <  1:706118/68‑1 (MQ=255)
                             ggATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTAc                              <  1:2443968/68‑1 (MQ=255)
                             ggATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTAc                              <  1:1263906/68‑1 (MQ=255)
                               aTTAGGGTTAGCGATCGTTAAACATGCTGTGAATCATCACGAAACTCg                                                <  1:1352355/48‑1 (MQ=255)
                               aTTAGGGTTAGCGATCGTGAAACATGCGGTGAATCATCACGa                                                      <  1:392979/42‑1 (MQ=255)
                               aTTAGGGTTAGCGATCGTGAAACATGCGGTGAATCATCACGa                                                      <  1:1175882/42‑1 (MQ=255)
                               aTTAGGGTTAGCGATCGTGAAACATGCGGTGAATCATCACGa                                                      <  1:2240051/42‑1 (MQ=255)
                               aTTAGGGTTAGCGATCGTGAAACATGCGGTGAATCATCACGa                                                      <  1:1515650/42‑1 (MQ=255)
                               aTTAGGGTTAGCGATCGTGAAACATGCGGTGAATCATCACGa                                                      <  1:160789/42‑1 (MQ=255)
                               aTTAGGGTTAGCGATCGTGAAACATGCGGTGAATCATCACGa                                                      <  1:2299513/42‑1 (MQ=255)
                               aTTAGGGTTAGCGATCGTGAAACATGCGGTGAATCATCACGa                                                      <  1:555159/42‑1 (MQ=255)
                               aTTAGGGTTAGCGATCGTGAAACATGCGGTGAATCATCACGa                                                      <  1:1920006/42‑1 (MQ=255)
                               aTTAGGGTTAGCGATCGTGAAACATGCGGTGAATCATCACGa                                                      <  1:3062487/42‑1 (MQ=255)
                               aTTAGGGTTAGCGATCGTGAAACATGCGGTGAATCATCACGa                                                      <  1:3017395/42‑1 (MQ=255)
                                      ttAGCGATCGTGAAACATGCGGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGaaaa                    >  1:230926/1‑69 (MQ=255)
                                      ttAGCGATCGTGAAACATGCGGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGaaaa                    >  1:1368935/1‑69 (MQ=255)
                                      ttAGCGATCGTGAAACATGCGGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGaaaa                    >  1:1923869/1‑69 (MQ=255)
                                      ttAGCGATCGTGAAACATGCGGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGaaaa                    >  1:2918230/1‑69 (MQ=255)
                                      ttAGCGATCGTGAAACATGCGGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGaaaa                    >  1:2826802/1‑69 (MQ=255)
                                      ttAGCGATCGTGAAACATGCGGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGaaaa                    >  1:2005075/1‑69 (MQ=255)
                                      ttAGCGATCGTGAAACATGCGGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGaaaa                    >  1:2672444/1‑69 (MQ=255)
                                      ttAGCGATCGTGAAACATGCGGTGAATCATCACGAAAGTCGCCAGAATATTGAGAGTACAGTAGGaaaa                    >  1:375550/1‑69 (MQ=255)
                                      ttAGCGATCGTGAAACATGCGGTGAATAATCACGAAAGACGCCTGAATATTGAGAGTAAAGTAGGaaaa                    >  1:2766684/1‑69 (MQ=255)
                                        aGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAgg                  >  1:2988446/1‑69 (MQ=255)
                                        aGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAgg                  >  1:130788/1‑69 (MQ=255)
                                        aGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAgg                  >  1:2782868/1‑69 (MQ=255)
                                        aGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAgg                  >  1:1290462/1‑69 (MQ=255)
                                        aGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAgg                  >  1:1181706/1‑69 (MQ=255)
                                        aGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAgg                  >  1:74472/1‑69 (MQ=255)
                                              cGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACg            >  1:1823054/1‑69 (MQ=255)
                                              cGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACg            >  1:2673840/1‑69 (MQ=255)
                                                       tGCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGttt   >  1:3271325/1‑69 (MQ=255)
                                                        gCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTATGAAAAGGAACACGTTTCAGtttt  <  1:3178167/69‑1 (MQ=255)
                                                        gCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGtttt  <  1:1521954/69‑1 (MQ=255)
                                                        gCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGtttt  <  1:1588332/69‑1 (MQ=255)
                                                        gCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGtttt  <  1:60445/69‑1 (MQ=255)
                                                        gCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGtttt  <  1:302519/69‑1 (MQ=255)
                                                        gCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGtttt  <  1:2856928/69‑1 (MQ=255)
                                                        gCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGtttt  <  1:2097492/69‑1 (MQ=255)
                                                        gCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGtttt  <  1:717027/69‑1 (MQ=255)
                                                        gCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGtttt  <  1:2573222/69‑1 (MQ=255)
                                                        gCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGtttt  <  1:2232091/69‑1 (MQ=255)
                                                        gCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGtttt  <  1:2422236/69‑1 (MQ=255)
                                                        gCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGtttt  <  1:820743/69‑1 (MQ=255)
                                                        gCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGtttt  <  1:821444/69‑1 (MQ=255)
                                                        gCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGtttt  <  1:899335/69‑1 (MQ=255)
                                                          |                                                                  
AAGCGCGTTCCCGGCAAACCGGCGGTAGCGGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTT  >  minE/247490‑247614

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