Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 148,890 C→T 41.6% D682D (GAC→GAT hrpB → predicted ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE148,8900CT41.6% 32.0 / 59.9 84D682D (GAC→GAThrpBpredicted ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base C (17/32);  new base T (26/9);  total (43/41)
Fisher's exact test for biased strand distribution p-value = 4.14e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGC  >  minE/148825‑148954
                                                                 |                                                                
cTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACAtt                                                               <  1:3388070/69‑1 (MQ=255)
cTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACAtt                                                               <  1:2598280/69‑1 (MQ=255)
cTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACAtt                                                               <  1:2223365/69‑1 (MQ=255)
cTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACAtt                                                               <  1:185983/69‑1 (MQ=255)
cTGGAAACGGGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACAtt                                                               <  1:929772/69‑1 (MQ=255)
 tGGAAACGTGGCTACTGCCGCATATGACAGGCGTACATTCGCTACGCGGCCTGAAATCACTCGATAttt                                                              <  1:2416856/69‑1 (MQ=255)
 tGGAAACGTGGCTACTGCCGCATATGACAGGCGTACATTCGCTACGCGGCCTGAAATCACTCGATAttt                                                              <  1:2757252/69‑1 (MQ=255)
       cGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAgg                                                        <  1:485420/69‑1 (MQ=255)
       cGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAgg                                                        <  1:2602550/69‑1 (MQ=255)
       cGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAgg                                                        <  1:2643014/69‑1 (MQ=255)
       cGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAgg                                                        <  1:2454004/69‑1 (MQ=255)
       cGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAgg                                                        <  1:393385/69‑1 (MQ=255)
       cGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAgg                                                        <  1:2647002/69‑1 (MQ=255)
       cGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAgg                                                        <  1:1262357/69‑1 (MQ=255)
       cGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAgg                                                        <  1:1724245/69‑1 (MQ=255)
       cGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAgg                                                        <  1:653384/69‑1 (MQ=255)
       cGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAgg                                                        <  1:1541273/69‑1 (MQ=255)
       cGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAgg                                                        <  1:819204/69‑1 (MQ=255)
       cGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAgg                                                        <  1:674579/69‑1 (MQ=255)
       cGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAgg                                                        <  1:2853208/69‑1 (MQ=255)
       cGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAgg                                                        <  1:340328/69‑1 (MQ=255)
        gTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGc                                                       >  1:1599397/1‑69 (MQ=255)
        gTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGc                                                       >  1:1858225/1‑69 (MQ=255)
         tGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGc                                                       >  1:2105837/1‑68 (MQ=255)
         tGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGc                                                       >  1:3211148/1‑68 (MQ=255)
         tGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGc                                                       >  1:2245390/1‑68 (MQ=255)
         tGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGc                                                       >  1:857525/1‑68 (MQ=255)
         tGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGc                                                       >  1:1681077/1‑68 (MQ=255)
                       tatGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTa                                          >  1:624921/1‑67 (MQ=255)
                       tatGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTAc                                         >  1:475617/1‑68 (MQ=255)
                       tatGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTAc                                         >  1:1513278/1‑68 (MQ=255)
                       tatGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTAc                                         >  1:704516/1‑68 (MQ=255)
                       tatGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTAc                                         >  1:2909482/1‑68 (MQ=255)
                       tatGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTA‑‑CGGATTACTTg                                      >  1:2262924/1‑69 (MQ=255)
                             tGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTg                                  <  1:2779072/69‑1 (MQ=255)
                             tGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTg                                  <  1:2690422/69‑1 (MQ=255)
                             tGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTg                                  <  1:3226368/69‑1 (MQ=255)
                             tGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTg                                  <  1:295239/69‑1 (MQ=255)
                             tGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTg                                  <  1:1335083/69‑1 (MQ=255)
                             tGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTg                                  <  1:1341235/69‑1 (MQ=255)
                             tGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTg                                  <  1:2461582/69‑1 (MQ=255)
                             tGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTg                                  <  1:1390568/69‑1 (MQ=255)
                             tGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTg                                  <  1:2212652/69‑1 (MQ=255)
                             tGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTg                                  <  1:796865/69‑1 (MQ=255)
                              ggCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTg                                  <  1:2651460/68‑1 (MQ=255)
                              ggCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTg                                  <  1:1915262/68‑1 (MQ=255)
                              ggCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTg                                  >  1:1503119/1‑68 (MQ=255)
                                    cATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGa                              >  1:2904961/1‑66 (MQ=255)
                                       tCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACtt                                       >  1:2875209/1‑54 (MQ=255)
                                       tCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTg                                      >  1:2523818/1‑55 (MQ=255)
                                                ggCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTg               <  1:2707373/69‑1 (MQ=255)
                                                       aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAAc                    >  1:673008/1‑57 (MQ=255)
                                                       aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAAc                    >  1:1466616/1‑57 (MQ=255)
                                                       aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAAc                    >  1:2897989/1‑57 (MQ=255)
                                                       aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                   >  1:1099266/1‑58 (MQ=255)
                                                       aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                   >  1:323742/1‑58 (MQ=255)
                                                       aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                   >  1:467820/1‑58 (MQ=255)
                                                       aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                   >  1:3427916/1‑58 (MQ=255)
                                                       aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                   >  1:447417/1‑58 (MQ=255)
                                                       aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                   >  1:907407/1‑58 (MQ=255)
                                                       aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                   >  1:2239662/1‑58 (MQ=255)
                                                       aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                   >  1:1452912/1‑58 (MQ=255)
                                                       aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                   >  1:134512/1‑58 (MQ=255)
                                                       aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                   >  1:1557250/1‑58 (MQ=255)
                                                       aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                   >  1:1576817/1‑58 (MQ=255)
                                                       aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                   >  1:1753849/1‑58 (MQ=255)
                                                       aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                   >  1:1821939/1‑58 (MQ=255)
                                                       aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                   >  1:1182674/1‑58 (MQ=255)
                                                       aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                   >  1:211965/1‑58 (MQ=255)
                                                       aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                   >  1:1110830/1‑58 (MQ=255)
                                                       aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                   >  1:2265562/1‑58 (MQ=255)
                                                       aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                   >  1:1065860/1‑58 (MQ=255)
                                                       aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                   >  1:2785088/1‑58 (MQ=255)
                                                       aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                   >  1:2796979/1‑58 (MQ=255)
                                                       aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                   >  1:2883623/1‑58 (MQ=255)
                                                       aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                   >  1:3190253/1‑58 (MQ=255)
                                                       aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATACAGCAACg                   >  1:2790594/1‑58 (MQ=255)
                                                             tCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACGTCTGGATAGTGAATTGc  <  1:992673/69‑1 (MQ=255)
                                                             tCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACGTCTGGATAGTGAATTGc  <  1:141690/69‑1 (MQ=255)
                                                             tCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACGTCTGGATAGTGAATTGc  <  1:1535134/69‑1 (MQ=255)
                                                             tCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACGTCTGGATAGTGAATTGc  <  1:1871733/69‑1 (MQ=255)
                                                             tCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACGTCTGGATAGTGAATTGc  <  1:444617/69‑1 (MQ=255)
                                                             tCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACGTCTGGATAGTGAATTGc  <  1:2709310/69‑1 (MQ=255)
                                                             tCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACGTCTGGATAGTGAATTGc  <  1:2932262/69‑1 (MQ=255)
                                                                 |                                                                
CTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGC  >  minE/148825‑148954

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: