Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 950,075 A→G 21.6% I254I (ATT→ATC dos ← cAMP phosphodiesterase, heme‑regulated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE950,0750AG21.6% 92.0 / 18.4 69I254I (ATT→ATCdoscAMP phosphodiesterase, heme‑regulated
Reads supporting (aligned to +/- strand):  ref base A (35/19);  new base G (10/5);  total (45/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

ATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATA  >  minE/950007‑950143
                                                                    |                                                                    
aTTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGa                                                                      <  1:1127451/69‑1 (MQ=255)
aTTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGa                                                                      <  1:3079975/69‑1 (MQ=255)
aTTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGa                                                                      <  1:1179913/69‑1 (MQ=255)
aTTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGa                                                                      <  1:2067486/69‑1 (MQ=255)
  ttCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGaa                                                                     <  1:2236123/68‑1 (MQ=255)
  ttCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGaa                                                                     <  1:3175986/68‑1 (MQ=255)
  ttCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGaa                                                                     <  1:237203/68‑1 (MQ=255)
  ttCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGaa                                                                     <  1:1954271/68‑1 (MQ=255)
  ttCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGaa                                                                     <  1:1944674/68‑1 (MQ=255)
  ttCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGaa                                                                     <  1:1910574/68‑1 (MQ=255)
                     ggCGGGCTGCTGCACATGGCGGCGAGAATGTTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGata                                                 <  1:1254843/69‑1 (MQ=255)
                     ggCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGata                                                 <  1:1853020/69‑1 (MQ=255)
                     ggCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGata                                                 <  1:3375929/69‑1 (MQ=255)
                         ggCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGatat                                                >  1:1091189/1‑66 (MQ=255)
                         ggCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGatat                                                >  1:3254591/1‑66 (MQ=255)
                         ggCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGatat                                                >  1:3434758/1‑66 (MQ=255)
                         ggCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGatat                                                >  1:1594702/1‑66 (MQ=255)
                         ggCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGatat                                                >  1:1622521/1‑66 (MQ=255)
                         ggCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGata                                                 >  1:1381104/1‑65 (MQ=255)
                               tgcACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGaaa                                        >  1:29316/1‑68 (MQ=255)
                               tgcACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGaaa                                        >  1:1469358/1‑68 (MQ=255)
                               tgcACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGaaa                                        >  1:2984494/1‑68 (MQ=255)
                               tgcACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGaaa                                        >  1:1363648/1‑68 (MQ=255)
                                gcACATGGCGGCGTGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGt                                      <  1:3019002/69‑1 (MQ=255)
                                gcACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGt                                      <  1:2937796/69‑1 (MQ=255)
                                gcACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGt                                      <  1:295681/69‑1 (MQ=255)
                                gcACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGt                                      <  1:1381073/69‑1 (MQ=255)
                                gcACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGt                                      <  1:3162779/69‑1 (MQ=255)
                                gcACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGt                                      <  1:1589840/69‑1 (MQ=255)
                                gcACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGt                                      <  1:401562/69‑1 (MQ=255)
                                      ggcggcGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTAc                                <  1:2007924/69‑1 (MQ=255)
                                      ggcggcGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTAc                                <  1:3325624/69‑1 (MQ=255)
                                       gcggcgAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTAc                                <  1:2090764/68‑1 (MQ=255)
                                          gcgAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAgg                            <  1:3231645/69‑1 (MQ=255)
                                          gcgAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAgg                            <  1:417908/69‑1 (MQ=255)
                                          gcgAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAgg                            <  1:2223975/69‑1 (MQ=255)
                                           cgAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAgg                            <  1:305443/68‑1 (MQ=255)
                                           cgAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAgg                            <  1:717137/68‑1 (MQ=255)
                                                   ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTc                         >  1:1166038/1‑63 (MQ=255)
                                                   ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAg                    >  1:2474796/1‑68 (MQ=255)
                                                   ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAg                    >  1:2785020/1‑68 (MQ=255)
                                                   ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGa                   >  1:3427993/1‑69 (MQ=255)
                                                   ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGa                   >  1:1883510/1‑69 (MQ=255)
                                                   ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGa                   >  1:2081042/1‑69 (MQ=255)
                                                   ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGa                   >  1:2876978/1‑69 (MQ=255)
                                                   ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGa                   >  1:705224/1‑69 (MQ=255)
                                                   ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGa                   >  1:2563839/1‑69 (MQ=255)
                                                   ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGa                   >  1:875451/1‑69 (MQ=255)
                                                   ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGa                   >  1:2260567/1‑69 (MQ=255)
                                                   ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGa                   >  1:2110202/1‑69 (MQ=255)
                                                    tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCg                                             >  1:1783806/1‑42 (MQ=37)
                                                    tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGa                   >  1:2201782/1‑68 (MQ=255)
                                                    tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGAt                  >  1:3263500/1‑69 (MQ=255)
                                                    tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGAt                  >  1:2122124/1‑69 (MQ=255)
                                                    tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGAt                  >  1:1283537/1‑69 (MQ=255)
                                                    tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGAt                  >  1:3042303/1‑69 (MQ=255)
                                                    tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGAt                  >  1:2325594/1‑69 (MQ=255)
                                                    tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGAt                  >  1:2806684/1‑69 (MQ=255)
                                                    tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGAt                  >  1:1679364/1‑69 (MQ=255)
                                                    tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGAt                  >  1:1407879/1‑69 (MQ=255)
                                                       cTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGATgc                <  1:2905556/68‑1 (MQ=255)
                                                       cTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGATgc                <  1:2089766/68‑1 (MQ=255)
                                                       cTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGATgc                <  1:2754517/68‑1 (MQ=255)
                                                       cTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGATgc                <  1:1419189/68‑1 (MQ=255)
                                                       cTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGATgc                <  1:2328465/68‑1 (MQ=255)
                                                         tCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGaaact                                      <  1:1855073/44‑3 (MQ=255)
                                                          cAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCa                     >  1:1986410/1‑60 (MQ=255)
                                                          cAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCa                     >  1:1812446/1‑60 (MQ=255)
                                                          cAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCa                     >  1:2574235/1‑60 (MQ=255)
                                                          cAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGcgcg             >  1:970133/1‑68 (MQ=255)
                                                          cAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGcgcg             >  1:1230445/1‑68 (MQ=255)
                                                          cAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGcgcg             >  1:1205740/1‑68 (MQ=255)
                                                          cAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGcgcg             >  1:2587904/1‑68 (MQ=255)
                                                              cTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCAc        >  1:1964826/1‑69 (MQ=255)
                                                              cTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCAc        >  1:2483638/1‑69 (MQ=255)
                                                              cTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCAc        >  1:1706527/1‑69 (MQ=255)
                                                              cTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCAc        >  1:2770031/1‑69 (MQ=255)
                                                              cTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCAc        >  1:394926/1‑69 (MQ=255)
                                                              cTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCAc        >  1:467994/1‑69 (MQ=255)
                                                                    aaTCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATa  <  1:2639793/69‑1 (MQ=255)
                                                                    aaTCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATa  <  1:1477346/69‑1 (MQ=255)
                                                                    aaTCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATa  <  1:1861930/69‑1 (MQ=255)
                                                                    aaTCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATa  <  1:2107135/69‑1 (MQ=255)
                                                                    aaTCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATa  <  1:2838490/69‑1 (MQ=255)
                                                                    aaTCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATa  <  1:1130328/69‑1 (MQ=255)
                                                                    |                                                                    
ATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATA  >  minE/950007‑950143

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: