Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 148,905 A→G 53.3% L687L (CTA→CTG hrpB → predicted ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE148,9050AG53.3% 6.0 / 50.7 62L687L (CTA→CTGhrpBpredicted ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base A (15/14);  new base G (26/7);  total (41/21)
Fisher's exact test for biased strand distribution p-value = 3.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

TATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCC  >  minE/148848‑148973
                                                         |                                                                    
tatGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTa                                                             >  1:624921/1‑67 (MQ=255)
tatGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTAc                                                            >  1:2909482/1‑68 (MQ=255)
tatGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTAc                                                            >  1:475617/1‑68 (MQ=255)
tatGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTAc                                                            >  1:1513278/1‑68 (MQ=255)
tatGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTAc                                                            >  1:704516/1‑68 (MQ=255)
tatGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTA‑‑CGGATTACTTg                                                         >  1:2262924/1‑69 (MQ=255)
      tGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTg                                                     <  1:295239/69‑1 (MQ=255)
      tGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTg                                                     <  1:2461582/69‑1 (MQ=255)
      tGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTg                                                     <  1:2690422/69‑1 (MQ=255)
      tGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTg                                                     <  1:2212652/69‑1 (MQ=255)
      tGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTg                                                     <  1:2779072/69‑1 (MQ=255)
      tGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTg                                                     <  1:1390568/69‑1 (MQ=255)
      tGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTg                                                     <  1:3226368/69‑1 (MQ=255)
      tGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTg                                                     <  1:1341235/69‑1 (MQ=255)
      tGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTg                                                     <  1:796865/69‑1 (MQ=255)
      tGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTg                                                     <  1:1335083/69‑1 (MQ=255)
       ggCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTg                                                     <  1:2651460/68‑1 (MQ=255)
       ggCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTg                                                     >  1:1503119/1‑68 (MQ=255)
       ggCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTg                                                     <  1:1915262/68‑1 (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:2897989/1‑57 (MQ=255)
                                aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAAc                                       >  1:1466616/1‑57 (MQ=255)
                                aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAAc                                       >  1:673008/1‑57 (MQ=255)
                                aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                      >  1:2785088/1‑58 (MQ=255)
                                aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                      >  1:2796979/1‑58 (MQ=255)
                                aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                      >  1:3427916/1‑58 (MQ=255)
                                aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                      >  1:323742/1‑58 (MQ=255)
                                aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                      >  1:2883623/1‑58 (MQ=255)
                                aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                      >  1:3190253/1‑58 (MQ=255)
                                aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                      >  1:467820/1‑58 (MQ=255)
                                aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                      >  1:1099266/1‑58 (MQ=255)
                                aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                      >  1:1110830/1‑58 (MQ=255)
                                aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                      >  1:1065860/1‑58 (MQ=255)
                                aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                      >  1:907407/1‑58 (MQ=255)
                                aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                      >  1:1182674/1‑58 (MQ=255)
                                aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                      >  1:134512/1‑58 (MQ=255)
                                aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                      >  1:1452912/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:211965/1‑58 (MQ=255)
                                aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                      >  1:2239662/1‑58 (MQ=255)
                                aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                      >  1:2265562/1‑58 (MQ=255)
                                aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                      >  1:447417/1‑58 (MQ=255)
                                aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATACAGCAACg                                      >  1:2790594/1‑58 (MQ=255)
                                      tCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACGTCTGGATAGTGAATTGc                     <  1:444617/69‑1 (MQ=255)
                                      tCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACGTCTGGATAGTGAATTGc                     <  1:992673/69‑1 (MQ=255)
                                      tCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACGTCTGGATAGTGAATTGc                     <  1:141690/69‑1 (MQ=255)
                                      tCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACGTCTGGATAGTGAATTGc                     <  1:2932262/69‑1 (MQ=255)
                                      tCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACGTCTGGATAGTGAATTGc                     <  1:1535134/69‑1 (MQ=255)
                                      tCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACGTCTGGATAGTGAATTGc                     <  1:1871733/69‑1 (MQ=255)
                                      tCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACGTCTGGATAGTGAATTGc                     <  1:2709310/69‑1 (MQ=255)
                                           aTTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTg                           >  1:2303631/1‑58 (MQ=255)
                                              tATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCAt             >  1:3197354/1‑69 (MQ=255)
                                              tATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCAt             >  1:1178586/1‑69 (MQ=255)
                                              tATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCAt             >  1:2592408/1‑69 (MQ=255)
                                              tATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTACCTGCGCAt             >  1:1169119/1‑69 (MQ=255)
                                                      tcTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTg                           <  1:76544/46‑1 (MQ=255)
                                                         aCGCGGATTACTTGGTTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTgcg                >  1:3075882/1‑55 (MQ=255)
                                                         aCGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATa                              >  1:1453328/1‑41 (MQ=255)
                                                         aCGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTgcg                >  1:2720025/1‑55 (MQ=255)
                                                         aCGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTgcg                <  1:322351/55‑1 (MQ=255)
                                                         aCGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTgcg                >  1:736474/1‑55 (MQ=255)
                                                         aCGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTgcg                >  1:79873/1‑55 (MQ=255)
                                                         aCGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTgcg                >  1:92673/1‑55 (MQ=255)
                                                         aCGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTgcg                >  1:1122985/1‑55 (MQ=255)
                                                         aCGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGcc  <  1:3396871/69‑1 (MQ=255)
                                                         aCGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGcc  <  1:3276502/69‑1 (MQ=255)
                                                         aCGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGcc  <  1:2859149/69‑1 (MQ=255)
                                                         |                                                                    
TATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCC  >  minE/148848‑148973

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: