Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 148,925 G→T 58.9% G694V (GGA→GTA)  hrpB → predicted ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE148,9250GT58.9% 16.9 / 54.7 56G694V (GGA→GTA) hrpBpredicted ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (13/10);  new base T (26/7);  total (39/17)
Fisher's exact test for biased strand distribution p-value = 8.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATC  >  minE/148861‑148989
                                                                |                                                                
cATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGa                                                                 >  1:2904961/1‑66 (MQ=255)
            ggCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTg                                                  <  1:2707373/69‑1 (MQ=255)
                   aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAAc                                                       >  1:673008/1‑57 (MQ=255)
                   aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAAc                                                       >  1:2897989/1‑57 (MQ=255)
                   aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAAc                                                       >  1:1466616/1‑57 (MQ=255)
                   aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                                      >  1:211965/1‑58 (MQ=255)
                   aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                                      >  1:1099266/1‑58 (MQ=255)
                   aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                                      >  1:2239662/1‑58 (MQ=255)
                   aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                                      >  1:2265562/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)
                   aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                                      >  1:323742/1‑58 (MQ=255)
                   aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                                      >  1:3427916/1‑58 (MQ=255)
                   aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                                      >  1:447417/1‑58 (MQ=255)
                   aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                                      >  1:467820/1‑58 (MQ=255)
                   aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                                      >  1:907407/1‑58 (MQ=255)
                   aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                                      >  1:1821939/1‑58 (MQ=255)
                   aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                                      >  1:1065860/1‑58 (MQ=255)
                   aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                                      >  1:1753849/1‑58 (MQ=255)
                   aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                                      >  1:1110830/1‑58 (MQ=255)
                   aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                                      >  1:1182674/1‑58 (MQ=255)
                   aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                                      >  1:1576817/1‑58 (MQ=255)
                   aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                                      >  1:134512/1‑58 (MQ=255)
                   aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                                      >  1:1557250/1‑58 (MQ=255)
                   aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACg                                                      >  1:1452912/1‑58 (MQ=255)
                   aaTCACTCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATACAGCAACg                                                      >  1:2790594/1‑58 (MQ=255)
                         tCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACGTCTGGATAGTGAATTGc                                     <  1:992673/69‑1 (MQ=255)
                         tCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACGTCTGGATAGTGAATTGc                                     <  1:1871733/69‑1 (MQ=255)
                         tCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACGTCTGGATAGTGAATTGc                                     <  1:444617/69‑1 (MQ=255)
                         tCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACGTCTGGATAGTGAATTGc                                     <  1:2932262/69‑1 (MQ=255)
                         tCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACGTCTGGATAGTGAATTGc                                     <  1:141690/69‑1 (MQ=255)
                         tCGATATTTATCAGGCACTGCGTGGATTACTTGATTGGGTAATGCAGCAACGTCTGGATAGTGAATTGc                                     <  1:1535134/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:2592408/1‑69 (MQ=255)
                                 tATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCAt                             >  1:1178586/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:79873/1‑55 (MQ=255)
                                            aCGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTgcg                                >  1:92673/1‑55 (MQ=255)
                                            aCGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTgcg                                <  1:322351/55‑1 (MQ=255)
                                            aCGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTgcg                                >  1:2720025/1‑55 (MQ=255)
                                            aCGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTgcg                                >  1:1122985/1‑55 (MQ=255)
                                            aCGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTgcg                                >  1:736474/1‑55 (MQ=255)
                                            aCGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGcc                  <  1:2859149/69‑1 (MQ=255)
                                            aCGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGcc                  <  1:3396871/69‑1 (MQ=255)
                                            aCGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGcc                  <  1:3276502/69‑1 (MQ=255)
                                             cgcgGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTgcg                                <  1:845935/54‑1 (MQ=255)
                                                           ttGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCgg     <  1:612233/67‑1 (MQ=255)
                                                           ttGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCgg     <  1:2695823/67‑1 (MQ=255)
                                                           ttGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCgg     <  1:2089727/67‑1 (MQ=255)
                                                             ggggAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATc  >  1:583318/1‑68 (MQ=255)
                                                             ggggAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATc  >  1:1809096/1‑68 (MQ=255)
                                                                |                                                                
CATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATC  >  minE/148861‑148989

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: