Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,335,299 G→A 49.3% Q916Q (CAG→CAA fdoG → formate dehydrogenase‑O, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,335,2990GA49.3% ‑0.3 / 70.2 79Q916Q (CAG→CAAfdoGformate dehydrogenase‑O, large subunit
Reads supporting (aligned to +/- strand):  ref base G (22/18);  new base A (22/17);  total (44/35)
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

TCGTCTGACCGAGCACTTCCACTACTGGACCAAACACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAGGGCGA  >  minE/2335239‑2335361
                                                            |                                                              
tCGTCTGACCGAGCACTTCCACTACTGGACCAAACACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCa                                                        >  1:1792563/1‑69 (MQ=255)
tCGTCTGACCGAGCACTTCCACTACTGGACCAAACACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCa                                                        >  1:921522/1‑69 (MQ=255)
tCGTCTGACCGAGCACTTCCACTACTGGACCAAACACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCa                                                        >  1:1807821/1‑69 (MQ=255)
tCGTCTGACCGAGCACTTCCACTACTGGACCAAACACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCa                                                        >  1:1234007/1‑69 (MQ=255)
tCGTCTGACCGAGCACTTCCACTACTGGACCAAACACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCa                                                        >  1:2026246/1‑69 (MQ=255)
tCGTCTGACCGAGCACTTCCACTACTGGACCAAACACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCa                                                        >  1:2688975/1‑69 (MQ=255)
              aCTTCCACTACTGGACCAAACACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCg                                          <  1:1340126/69‑1 (MQ=255)
              aCTTCCACTACTGGACCAAACACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCg                                          <  1:1828089/69‑1 (MQ=255)
              aCTTCCACTACTGGACCAAACACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCg                                          <  1:2290441/69‑1 (MQ=255)
                 tCCACTACTGGACCAAACACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCgggg                                       >  1:2616722/1‑69 (MQ=255)
                 tCCACTACTGGACCAAACACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCgggg                                       >  1:2601880/1‑69 (MQ=255)
                 tCCACTACTGGACCAAACACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCgggg                                       >  1:2074525/1‑69 (MQ=255)
                 tCCACTACTGGACCAAACACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCgggg                                       >  1:3430148/1‑69 (MQ=255)
                 tCCACTACTGGACCAAACACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCgggg                                       >  1:1969458/1‑69 (MQ=255)
                 tCCACTACTGGACCAAACACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCgggg                                       >  1:378846/1‑69 (MQ=255)
                 tCCACTACTGGACCAAACACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCgggg                                       >  1:70653/1‑69 (MQ=255)
                 tCCACTACTGGACCAAACACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCgggg                                       >  1:992016/1‑69 (MQ=255)
                 tCCACTACTGGACCAAACACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCgggg                                       >  1:1195966/1‑69 (MQ=255)
                 tCCACTACTGGACCAAACACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCgggg                                       >  1:1199859/1‑69 (MQ=255)
                  ccACTACTGGACCAAACACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCgggg                                       >  1:802759/1‑68 (MQ=255)
                       actGGACCAAACACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCGGGGAGTCGc                                 >  1:2428774/1‑69 (MQ=255)
                       actGGACCAAACACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCGGGGAGTCGc                                 >  1:2011032/1‑69 (MQ=255)
                       actGGACCAAACACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCGGGGAGTCGc                                 >  1:533127/1‑69 (MQ=255)
                          ggACCAAACACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTg                               >  1:705809/1‑68 (MQ=255)
                          ggACCAAACACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTg                               >  1:475443/1‑68 (MQ=255)
                                aaCACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATa                        <  1:1054824/69‑1 (MQ=255)
                                aaCACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATa                        <  1:1128772/69‑1 (MQ=255)
                                 acacGCGCTGTTGAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATaa                       <  1:3282903/69‑1 (MQ=255)
                                  cacGCGCTGTTGAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATaaa                      <  1:212495/69‑1 (MQ=255)
                                  cacGCGCTGTTGAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATaaa                      <  1:420498/69‑1 (MQ=255)
                                       tatGTTGAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTgg                 <  1:390526/67‑1 (MQ=255)
                                             gAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTgg                 <  1:55935/63‑1 (MQ=255)
                                             gAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTgg                 <  1:697869/63‑1 (MQ=255)
                                             gAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTgg                 <  1:3261340/63‑1 (MQ=255)
                                             gAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTgg                 <  1:210128/63‑1 (MQ=255)
                                             gAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTgg                 <  1:3280983/63‑1 (MQ=255)
                                             gAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTgg                 <  1:3344491/63‑1 (MQ=255)
                                             gAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTgg                 <  1:3159465/63‑1 (MQ=255)
                                             gAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTgg                 <  1:3322502/63‑1 (MQ=255)
                                             gAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCAAATAAACTTgg                 <  1:2969834/63‑1 (MQ=255)
                                                  cgATTTTGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAgg      >  1:432337/1‑69 (MQ=255)
                                                  cgATTTTGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAgg      >  1:3308452/1‑69 (MQ=255)
                                                  cgATTTTGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAgg      >  1:1315847/1‑69 (MQ=255)
                                                  cgATTTTGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAgg      >  1:2923897/1‑69 (MQ=255)
                                                  cgATTTTGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAgg      >  1:836397/1‑69 (MQ=255)
                                                  cgATTTTGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAgg      >  1:1128214/1‑69 (MQ=255)
                                                  cgATTTTGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAgg      >  1:238771/1‑69 (MQ=255)
                                                  cgATTTTGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAgg      >  1:2199455/1‑69 (MQ=255)
                                                  cgATTTTGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAgg      >  1:2158933/1‑69 (MQ=255)
                                                  cgATTTTGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAgg      >  1:2117298/1‑69 (MQ=255)
                                                  cgATTTTGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAgg      >  1:1376482/1‑69 (MQ=255)
                                                  cgATTTTGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAgg      >  1:2293880/1‑69 (MQ=255)
                                                  cgATTTTGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAgg      >  1:1401238/1‑69 (MQ=255)
                                                  cgATTTTGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAgg      >  1:1509529/1‑69 (MQ=255)
                                                  cgATTTTGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAgg      >  1:2369621/1‑69 (MQ=255)
                                                  cgATTTTGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAgg      >  1:1359966/1‑69 (MQ=255)
                                                  cgATTTTGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAgg      >  1:2620673/1‑69 (MQ=255)
                                                  cgATTTTGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAgg      >  1:2718578/1‑69 (MQ=255)
                                                  cgATTTTGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAgg      >  1:2719081/1‑69 (MQ=255)
                                                  cgATTTTGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAgg      >  1:2757666/1‑69 (MQ=255)
                                                  cgATTTTGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAg       >  1:389551/1‑68 (MQ=255)
                                                  cgATTTTGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAg       >  1:1536717/1‑68 (MQ=255)
                                                      tttGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAGGGCGa  <  1:2017892/69‑1 (MQ=255)
                                                      tttGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAGGGCGa  <  1:1616487/69‑1 (MQ=255)
                                                      tttGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAGGGCGa  <  1:820893/69‑1 (MQ=255)
                                                      tttGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAGGGCGa  <  1:1978955/69‑1 (MQ=255)
                                                      tttGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAGGGCGa  <  1:974962/69‑1 (MQ=255)
                                                      tttGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAGGGCGa  <  1:546737/69‑1 (MQ=255)
                                                      tttGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAGGGCGa  <  1:530773/69‑1 (MQ=255)
                                                      tttGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAGGGCGa  <  1:507126/69‑1 (MQ=255)
                                                      tttGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAGGGCGa  <  1:2257445/69‑1 (MQ=255)
                                                      tttGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAGGGCGa  <  1:2316447/69‑1 (MQ=255)
                                                      tttGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAGGGCGa  <  1:2338096/69‑1 (MQ=255)
                                                      tttGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAGGGCGa  <  1:2492969/69‑1 (MQ=255)
                                                      tttGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAGGGCGa  <  1:2626137/69‑1 (MQ=255)
                                                      tttGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAGGGCGa  <  1:1323216/69‑1 (MQ=255)
                                                      tttGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAGGGCGa  <  1:2964385/69‑1 (MQ=255)
                                                      tttGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAGGGCGa  <  1:2839711/69‑1 (MQ=255)
                                                      tttGCAACCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGAATTGCCCAGGGCGa  <  1:1701298/69‑1 (MQ=255)
                                                            |                                                              
TCGTCTGACCGAGCACTTCCACTACTGGACCAAACACGCGCTGTTGAACGCGATTTTGCAGCCAGAGCAGTTTGTGGAAATCGGGGAGTCGCTGGCGAATAAACTTGGCATTGCCCAGGGCGA  >  minE/2335239‑2335361

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: