Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 147,876 C→T 29.6% A344A (GCC→GCT hrpB → predicted ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE147,8760CT29.6% 62.7 / 26.9 64A344A (GCC→GCThrpBpredicted ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base C (22/23);  new base T (10/9);  total (32/32)
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.87e-01

TCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATGGAAT  >  minE/147822‑147934
                                                      |                                                          
tCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGa                                              <  1:196224/69‑1 (MQ=255)
 cTGGAGCCGGGTATCTGCCTTCATTTAATCGCCAAAGAACAAGCAGAACGCGCTGCCGCGCAAAGTg                                               >  1:1105008/1‑67 (MQ=255)
 cTGGAGCCGGGTATCTGCCTTCATTTAATCGCCAAAGAACAAGCAGAACGCGCTGCCGCGCAAAGTg                                               >  1:2817785/1‑67 (MQ=255)
 cTGGAGCCGGGTATCTGCCTTCATTTAATCGCCAAAGAACAAGCAGAACGCGCTGCCGCGCAAAGTg                                               >  1:1576222/1‑67 (MQ=255)
 cTGGAGCCGGGTATCTGCCTTCATTTAATCGCCAAAGAACAAGCAGAACGCGCTGCCGCGCAAAGTg                                               >  1:662494/1‑67 (MQ=255)
 cTGGAGCCGGGTATCTGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCTGCCGCGCAAAGTg                                               >  1:3088000/1‑67 (MQ=255)
 cTGGAGCCGGGTATCTGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCTGCCGCGCAAAGTg                                               >  1:200945/1‑67 (MQ=255)
 cTGGAGCCGGGTATCTGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCTGCCGCGCAAAGTg                                               >  1:22536/1‑67 (MQ=255)
 cTGGAGCCGGGTATCTGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCTGCCGCGCAAAGTg                                               >  1:3055889/1‑67 (MQ=255)
 cTGGAGCCGGGTATCTGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCTGCCGCGCAAAGTg                                               >  1:1418568/1‑67 (MQ=255)
 cTGGAGCCGGGTATCTGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCTGCCGCGCAAAGTg                                               >  1:1203729/1‑67 (MQ=255)
        cGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGgaga                                      <  1:616759/69‑1 (MQ=255)
        cGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGgaga                                      <  1:2874394/69‑1 (MQ=255)
        cGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGgaga                                      <  1:984877/69‑1 (MQ=255)
                aGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCTCAAAGTGAACCGGAGATCtt                                  >  1:579749/1‑65 (MQ=255)
                aGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCtt                                  >  1:7675/1‑65 (MQ=255)
                aGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCtt                                  >  1:2493771/1‑65 (MQ=255)
                aGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCtt                                  >  1:17961/1‑65 (MQ=255)
                aGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACaa                              >  1:1183566/1‑69 (MQ=255)
                aGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACaa                              >  1:2210875/1‑69 (MQ=255)
                aGCCTGAATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCg                                          >  1:6911/1‑57 (MQ=255)
                    tGCATTTAATCGCCAAAGAACAAGCAGAACGCGCTGCCGCGCAAAGTGAACCGGAGATCTTACAAAg                            <  1:2544047/67‑1 (MQ=255)
                    tGCATTTAATCGCCAAAGAACAAGCAGAACGCGCTGCCGCGCAAAGTGAACCGGAGATCTTACAAAg                            <  1:861589/67‑1 (MQ=255)
                    tGCATTTAATCGCCAAAGAACAAGCAGAACGCGCTGCCGCGCAAAGTGAACCGGAGATCTTACAAAg                            <  1:2142046/67‑1 (MQ=255)
                    tGCATTTAATCGCCAAAGAACAAGCAGAACGCGCTGCCGCGCAAAGTGAACCGGAGATCTTACAAAg                            <  1:2023367/67‑1 (MQ=255)
                    tGCATTTAATCGCCAAAGAACAAGCAGAACGCGCTGCCGCGCAAAGTGAACCGGAGATCTTACAAAg                            <  1:11256/67‑1 (MQ=255)
                    tGCATTTAATCGCCAAAGAACAAGCAGAACGCGCTGCCGCGCAAAGTGAACCGGAGATCTTACAAAg                            <  1:1322127/67‑1 (MQ=255)
                    tGCATTTAATCGCCAAAGAACAAGCAGAACGCGCTGCCGCGCAAAGTGAACCGGAGATCTTACAAAg                            <  1:1386991/67‑1 (MQ=255)
                    tGCATTTAATCGCCAAAGAACAAGCAGAACGCGCTGCCGCGAAAAGTGAACCGGAGATCTTACAAAg                            <  1:687153/67‑1 (MQ=255)
                     gCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCg                          <  1:3098256/68‑1 (MQ=255)
                      cATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGAt                        >  1:1990316/1‑69 (MQ=255)
                      cATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGAt                        >  1:893442/1‑69 (MQ=255)
                      cATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGAt                        >  1:1938597/1‑69 (MQ=255)
                      cATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGAt                        >  1:2600774/1‑69 (MQ=255)
                      cATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGAt                        >  1:2654182/1‑69 (MQ=255)
                      cATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGAt                        >  1:1423584/1‑69 (MQ=255)
                      cATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGAt                        >  1:1739089/1‑69 (MQ=255)
                      cATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGAt                        >  1:1736260/1‑69 (MQ=255)
                          tAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCttt                    >  1:1104293/1‑69 (MQ=255)
                          tAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCttt                    >  1:459864/1‑69 (MQ=255)
                              cGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCgg                <  1:202868/69‑1 (MQ=255)
                              cGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCgg                <  1:2192783/69‑1 (MQ=255)
                              cGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCgg                <  1:879705/69‑1 (MQ=255)
                               gCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCgg                >  1:2913918/1‑68 (MQ=255)
                                    aGAACAAGCAGAACGCGCTGCCGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCt          <  1:3218178/69‑1 (MQ=255)
                                    aGAACAACCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAACCGATCTTTCCGGTTTGCt          <  1:2707830/69‑1 (MQ=255)
                                     gAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTg         >  1:2656772/1‑69 (MQ=255)
                                     gAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTg         >  1:1512417/1‑69 (MQ=255)
                                          aGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATgg     <  1:972994/68‑1 (MQ=255)
                                            cAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATgg     <  1:2691489/66‑1 (MQ=255)
                                            cAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATgg     <  1:3297983/66‑1 (MQ=255)
                                            cAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATgg     <  1:598981/66‑1 (MQ=255)
                                            cAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATgg     <  1:3107266/66‑1 (MQ=255)
                                            cAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATgg     >  1:1742102/1‑66 (MQ=255)
                                            cAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATgg     <  1:1375237/66‑1 (MQ=255)
                                            cAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATgg     <  1:736424/66‑1 (MQ=255)
                                            cAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATgg     <  1:2471930/66‑1 (MQ=255)
                                            cAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATgg     >  1:1966119/1‑66 (MQ=255)
                                            cAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATgg     <  1:1533246/66‑1 (MQ=255)
                                            cAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATgg     <  1:1588845/66‑1 (MQ=255)
                                            cAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATGGAAt  <  1:577806/69‑1 (MQ=255)
                                            cAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATGGAAt  <  1:1035440/69‑1 (MQ=255)
                                            cAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATGGAAt  <  1:925294/69‑1 (MQ=255)
                                            cAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATGGAAt  <  1:2943021/69‑1 (MQ=255)
                                                      |                                                          
TCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATGGAAT  >  minE/147822‑147934

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: