Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 950,105 A→T 27.0% T244T (ACT→ACA dos ← cAMP phosphodiesterase, heme‑regulated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE950,1050AT27.0% 141.2 / 57.9 100T244T (ACT→ACAdoscAMP phosphodiesterase, heme‑regulated
Reads supporting (aligned to +/- strand):  ref base A (35/38);  new base T (9/18);  total (44/56)
Fisher's exact test for biased strand distribution p-value = 2.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 105 total aligned reads displayed.

TGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGATAGAGGCTTTA  >  minE/950038‑950165
                                                                   |                                                            
tgcACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGaaa                                                              >  1:29316/1‑68 (MQ=255)
tgcACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGaaa                                                              >  1:1363648/1‑68 (MQ=255)
tgcACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGaaa                                                              >  1:2984494/1‑68 (MQ=255)
tgcACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGaaa                                                              >  1:1469358/1‑68 (MQ=255)
 gcACATGGCGGCGTGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGt                                                            <  1:3019002/69‑1 (MQ=255)
 gcACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGt                                                            <  1:401562/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)
       ggcggcGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTAc                                                      <  1:2007924/69‑1 (MQ=255)
       ggcggcGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTAc                                                      <  1:3325624/69‑1 (MQ=255)
        gcggcgAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTAc                                                      <  1:2090764/68‑1 (MQ=255)
           gcgAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAgg                                                  <  1:417908/69‑1 (MQ=255)
           gcgAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAgg                                                  <  1:3231645/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:2785020/1‑68 (MQ=255)
                    ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGa                                         >  1:2563839/1‑69 (MQ=255)
                    ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGa                                         >  1:2260567/1‑69 (MQ=255)
                    ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGa                                         >  1:3427993/1‑69 (MQ=255)
                    ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGa                                         >  1:2876978/1‑69 (MQ=255)
                    ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGa                                         >  1:2110202/1‑69 (MQ=255)
                    ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGa                                         >  1:875451/1‑69 (MQ=255)
                    ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGa                                         >  1:2081042/1‑69 (MQ=255)
                    ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGa                                         >  1:1883510/1‑69 (MQ=255)
                    ttGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGa                                         >  1:705224/1‑69 (MQ=255)
                     tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGa                                         >  1:2201782/1‑68 (MQ=255)
                     tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGAt                                        >  1:2122124/1‑69 (MQ=255)
                     tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGAt                                        >  1:1283537/1‑69 (MQ=255)
                     tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGAt                                        >  1:2806684/1‑69 (MQ=255)
                     tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGAt                                        >  1:1407879/1‑69 (MQ=255)
                     tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGAt                                        >  1:3263500/1‑69 (MQ=255)
                     tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGAt                                        >  1:2325594/1‑69 (MQ=255)
                     tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGAt                                        >  1:3042303/1‑69 (MQ=255)
                     tGCCTTCAAGCTGGCGGATTTGCCTTTCTTCGGTGATATCCGAGAATGTCATTACCAGGTTCTGCAGAt                                        >  1:1679364/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:2574235/1‑60 (MQ=255)
                           cAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCa                                           >  1:1812446/1‑60 (MQ=255)
                           cAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCa                                           >  1:1986410/1‑60 (MQ=255)
                           cAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGcgcg                                   >  1:2587904/1‑68 (MQ=255)
                           cAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGcgcg                                   >  1:1230445/1‑68 (MQ=255)
                           cAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGcgcg                                   >  1:1205740/1‑68 (MQ=255)
                               cTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCAc                              >  1:467994/1‑69 (MQ=255)
                               cTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCAc                              >  1:1706527/1‑69 (MQ=255)
                               cTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCAc                              >  1:2483638/1‑69 (MQ=255)
                               cTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCAc                              >  1:2770031/1‑69 (MQ=255)
                               cTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCAc                              >  1:1964826/1‑69 (MQ=255)
                               cTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCAc                              >  1:394926/1‑69 (MQ=255)
                                     aaTCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATa                        <  1:1130328/69‑1 (MQ=255)
                                     aaTCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATa                        <  1:1477346/69‑1 (MQ=255)
                                     aaTCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATa                        <  1:2838490/69‑1 (MQ=255)
                                     aaTCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATa                        <  1:2639793/69‑1 (MQ=255)
                                     aaTCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATa                        <  1:2107135/69‑1 (MQ=255)
                                      aTCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATaa                       <  1:577141/69‑1 (MQ=255)
                                      aTCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATa                        <  1:2497709/68‑1 (MQ=255)
                                                cttcGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCa                                           >  1:605198/1‑39 (MQ=255)
                                                cttcGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCa                                           >  1:3137654/1‑39 (MQ=255)
                                                cttcGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGa             <  1:2878506/69‑1 (MQ=255)
                                                cttcGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGa             <  1:746894/69‑1 (MQ=255)
                                                cttcGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGa             <  1:449409/69‑1 (MQ=255)
                                                cttcGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGa             <  1:3433465/69‑1 (MQ=255)
                                                cttcGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGa             <  1:2217538/69‑1 (MQ=255)
                                                cttcGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGa             <  1:2863824/69‑1 (MQ=255)
                                                cttcGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGa             <  1:2648888/69‑1 (MQ=255)
                                                 ttcGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGa             <  1:3079010/68‑1 (MQ=255)
                                                      tGATATCCGAGAAAGTCATTACCAGGTTGTGCAGATGCGCGAGCACGTCATAAACCGGGCTGATaga         >  1:2011508/1‑67 (MQ=255)
                                                      tGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGATaga         >  1:2811576/1‑67 (MQ=255)
                                                      tGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGATaga         >  1:2991297/1‑67 (MQ=255)
                                                           tCCGAGAATGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGATa           <  1:376181/60‑1 (MQ=255)
                                                           tCCGAGAATGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGATa           <  1:2902612/60‑1 (MQ=255)
                                                           tCCGAGAATGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGATa           <  1:812218/60‑1 (MQ=255)
                                                           tCCGAGAATGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGATa           <  1:42041/60‑1 (MQ=255)
                                                           tCCGAGAATGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGATa           <  1:832791/60‑1 (MQ=255)
                                                           tCCGAGAATGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGATa           <  1:3030192/60‑1 (MQ=255)
                                                           tCCGAGAATGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGATa           <  1:2639758/60‑1 (MQ=255)
                                                           tCCGAGAATGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGATa           <  1:2137145/60‑1 (MQ=255)
                                                           tCCGAGAATGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGATa           <  1:2371348/60‑1 (MQ=255)
                                                           tCCGAGAATGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGATa           <  1:1006770/60‑1 (MQ=255)
                                                           tCCGAGAATGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGATa           <  1:949424/60‑1 (MQ=255)
                                                           tCCGAGAATGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGATa           <  1:3088796/60‑1 (MQ=255)
                                                           tCCGAGAATGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGATa           <  1:3106306/60‑1 (MQ=255)
                                                           tCCGAGAAAGTCATTACCAGGTTCTGCAGATgcgc                                    >  1:2295411/1‑35 (MQ=255)
                                                           tCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGa             <  1:3101030/58‑1 (MQ=255)
                                                           tCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGa             <  1:964770/58‑1 (MQ=255)
                                                           tCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGa             <  1:1249142/58‑1 (MQ=255)
                                                           tCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGa             <  1:2497465/58‑1 (MQ=255)
                                                           tCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGa             <  1:3099754/58‑1 (MQ=255)
                                                           tCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGATAGAGGCTTTa  >  1:722875/1‑69 (MQ=255)
                                                           tCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGATAGAGGCTTTa  >  1:1634337/1‑69 (MQ=255)
                                                           tCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGATAGAGGCTTTa  >  1:1808731/1‑69 (MQ=255)
                                                           tCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGATAGAGGCTTTa  >  1:3182093/1‑69 (MQ=255)
                                                                   |                                                            
TGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCGAGAAAGTCATTACCAGGTTCTGCAGATGCGCGAGCACGTCATAAACCGGGCTGATAGAGGCTTTA  >  minE/950038‑950165

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: