Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 21,574 C→T 45.2% V502V (GTC→GTT ileS → isoleucyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE21,5740CT45.2% 7.4 / 41.0 53V502V (GTC→GTTileSisoleucyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (10/19);  new base T (9/15);  total (19/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGACACGGAAGAGCTGCATCCGCGTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTA  >  minE/21506‑21639
                                                                    |                                                                 
aaaGACACGGAAGAGCTGCATCCGCGTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTc                                                                   <  1:2593601/69‑1 (MQ=255)
aaaGACACGGAAGAGCTGCATCCGCGTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTc                                                                   <  1:760468/69‑1 (MQ=255)
aaaGACACGGAAGAGCTGCATCCGCGTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTc                                                                   <  1:1441852/69‑1 (MQ=255)
aaaGACACGGAAGAGCTGCATCCGCGTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTc                                                                   <  1:2625284/69‑1 (MQ=255)
aaaGACACGGAAGAGCTGCATCCGCGTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTc                                                                   <  1:2905767/69‑1 (MQ=255)
aaaGACACGGAAGAGCTGCATCCGCGTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTc                                                                   <  1:3196549/69‑1 (MQ=255)
 aaGACACGGAAGAGCTGCATCCGCGTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTc                                                                   <  1:1679322/68‑1 (MQ=255)
  aGACACAGAAGAACTGCATCCGCGTACTCTCGAACTGATGGAAGAAGTGGCAAAACGCTTTGAAGTTGa                                                                 >  1:551225/1‑69 (MQ=255)
  aGACACAGAAGAACTGCATCCGCGTACTCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTGa                                                                 >  1:2053686/1‑69 (MQ=255)
  aGACACAGAAGAACTGCATCCGCGTACTCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTGa                                                                 >  1:2390660/1‑69 (MQ=255)
  aGACACAGAAGAACTGCATCCGCGTACTCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTGa                                                                 >  1:2378297/1‑69 (MQ=255)
  aGACACAGAAGAACTGCATCCGCGTACTCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTGa                                                                 >  1:2645525/1‑69 (MQ=255)
  aGACACAGAAGAACTGCATCCGCGTACTCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTGa                                                                 >  1:2957925/1‑69 (MQ=255)
  aGACACAGAAGAACTGCATCCGCGTACTCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTGa                                                                 >  1:3340560/1‑69 (MQ=255)
  aGACACAGAAGAACTGCATCCGCGTACTCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTGa                                                                 >  1:2455100/1‑69 (MQ=255)
  aGACACAGAAGAACTGCATCCGCGTACTCTCGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTTGa                                                                 >  1:665643/1‑69 (MQ=255)
                tGCATCCGCGTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCgt                                                   <  1:3365331/69‑1 (MQ=255)
                tGCATCCGCGTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCgt                                                   <  1:745130/69‑1 (MQ=255)
                      cgcgTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCa                                                        >  1:199037/1‑58 (MQ=255)
                      cgcgTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGtggtgg                                              >  1:332539/1‑68 (MQ=255)
                      cgcgTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGtggtgg                                              >  1:2382572/1‑68 (MQ=255)
                      cgcgTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTggg                                             <  1:3336667/69‑1 (MQ=255)
                      cgcgTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTggg                                             >  1:1371620/1‑69 (MQ=255)
                      cgcgTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTggg                                             <  1:1003257/69‑1 (MQ=255)
                           aCCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGAtctc                                        <  1:2542559/69‑1 (MQ=255)
                           aCCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGAtctc                                        <  1:886584/69‑1 (MQ=255)
                                  aaCTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGc                                   >  1:3365479/1‑67 (MQ=255)
                                  aaCTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGc                                   >  1:2612913/1‑67 (MQ=255)
                                      gATGGAAGAAGTGGCAAAACGCGTTGAAGTTGACGGCATTCAGGCGTGGTGTGATCTCGATGCGAAaga                             <  1:3314634/69‑1 (MQ=255)
                                      gATGGAAGAAGTGGCAAAACGCGTTGAAGTTGACGGCATTCAGGCGTGGTGGGATCTCGATGCGAAaga                             <  1:2510958/69‑1 (MQ=255)
                                      gATGGAAGAAGTGGCAAAACGCGTTGAAGTTGACGGCATTCAGGCGTGGTGGGATCTCGATGCGAAaga                             <  1:1317559/69‑1 (MQ=255)
                                      gATGGAAGAAGTGGCAAAACGCGTTGAAGTTGACGGCATTCAGGCGTGGTGGGATCTCGATGCGAAaga                             <  1:1340887/69‑1 (MQ=255)
                                      gATGGAAGAAGTGGCAAAACGCGTTGAAGTTGACGGCATTCAGGCGTGGTGGGATCTCGATGCGAAaga                             <  1:1454261/69‑1 (MQ=255)
                                      gATGGAAGAAGTGGCAAAACGCGTTGAAGTTGACGGCATTCAGGCGTGGTGGGATCTCGATGCGAAaga                             <  1:605578/69‑1 (MQ=255)
                                      gATGGAAGAAGTGGCAAAACGCGTTGAAGTTGACGGCATTCAGGCGTGGTGGGATCTCGATGCGAAaga                             <  1:1655350/69‑1 (MQ=255)
                                      gATGGAAGAAGTGGCAAAACGCGTTGAAGTTGACGGCATTCAGGCGTGGTGGGATCTCGATGCGAAaga                             <  1:1958253/69‑1 (MQ=255)
                                      gATGGAAGAAGTGGCAAAACGCGTTGAAGTTGACGGCATTCAGGCGTGGTGGGATCTCGATGCGAAaga                             <  1:3362794/69‑1 (MQ=255)
                                      gATGGAAGAAGTGGCAAAACGCGTTGAAGTTGACGGCATTCAGGCGTGGTGGGATCTCGATGCGAAaga                             <  1:2236105/69‑1 (MQ=255)
                                      gATGGAAGAAGTGGCAAAACGCGTTGAAGTTGACGGCATTCAGGCGTGGTGGGATCTCGATGCGAAaga                             <  1:226772/69‑1 (MQ=255)
                                      gATGGAAGAAGTGGCAAAACGCGTTGAAGTTGACGGCATTCAGGCGTGGTGGGATCTCGATGCGAAaga                             <  1:2419886/69‑1 (MQ=255)
                                      gATGGAAGAAGTGGCAAAACGCGTTGAAGTTGACGGCATTCAGGCGTGGTGGGATCTCGATGCGAAaga                             <  1:2308676/69‑1 (MQ=255)
                                      gATGGAAGAAGTGGCAAAACGCGTTGAAGTTGACGGCATTCAGGCGTGGTGGGATCTCGATGCGAAaga                             <  1:2923600/69‑1 (MQ=255)
                                      gATGGAAGAAGTGGCAAAACGCGTTGAAGTTGACGGCATTCAGGCGTGGTGGGATCTCGATGCGAAaga                             <  1:2358817/69‑1 (MQ=255)
                                         ggAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTTGTGGGATCTCGATGCGAAAGAGAt                          >  1:2968643/1‑69 (MQ=255)
                                         ggAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGAt                          >  1:470472/1‑69 (MQ=255)
                                          gaagaaGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGtggtggt                                             <  1:3194857/49‑2 (MQ=255)
                                          gaagaaGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTggg                                             <  1:2287758/49‑1 (MQ=255)
                                          gaagaaGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTggg                                             <  1:685185/49‑1 (MQ=255)
                                          gaagaaGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTggg                                             <  1:2371954/49‑1 (MQ=255)
                                          gaagaaGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTggg                                             <  1:23164/49‑1 (MQ=255)
                                          gaagaaGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTggg                                             <  1:471436/49‑1 (MQ=255)
                                          gaagaaGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTggg                                             <  1:3435657/49‑1 (MQ=255)
                                          gaagaaGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTggg                                             <  1:1859321/49‑1 (MQ=255)
                                          gaagaaGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTggg                                             <  1:3021836/49‑1 (MQ=255)
                                          gaagaaGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTggg                                             <  1:3175608/49‑1 (MQ=255)
                                           aagaagTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTggg                                             <  1:2557749/48‑1 (MQ=255)
                                                  ggCAAAACGCGTTGAAGTCGATGGCATCCAGGCgt                                                   <  1:208019/35‑1 (MQ=255)
                                                        aCGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAg            >  1:55307/1‑68 (MQ=255)
                                                        aCGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAg            >  1:1341681/1‑68 (MQ=255)
                                                                  gTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTa  <  1:721391/68‑1 (MQ=255)
                                                                    |                                                                 
AAAGACACGGAAGAGCTGCATCCGCGTACCCTTGAACTGATGGAAGAAGTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTA  >  minE/21506‑21639

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: