Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,948,808 C→T 61.7% V4I (GTA→ATA)  gcvH ← glycine cleavage complex lipoylprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,948,8080CT61.7% 53.1 / 72.4 107V4I (GTA→ATA) gcvHglycine cleavage complex lipoylprotein
Reads supporting (aligned to +/- strand):  ref base C (28/13);  new base T (19/47);  total (47/60)
Fisher's exact test for biased strand distribution p-value = 1.13e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

AGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAG  >  minE/1948746‑1948875
                                                              |                                                                   
aGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCt                                                               >  1:3175132/1‑69 (MQ=255)
aGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCt                                                               >  1:2192582/1‑69 (MQ=255)
aGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCt                                                               <  1:1758226/69‑1 (MQ=255)
aGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCt                                                               >  1:151105/1‑69 (MQ=255)
        cGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCatcaatc                                                       <  1:2745741/69‑1 (MQ=255)
        cGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCatcaatc                                                       <  1:2594629/69‑1 (MQ=255)
         ggCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCatcaatca                                                      <  1:1071761/69‑1 (MQ=255)
         ggCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCatcaatca                                                      <  1:976947/69‑1 (MQ=255)
         ggCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCatcaatca                                                      <  1:841425/69‑1 (MQ=255)
         ggCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCatcaatca                                                      <  1:687065/69‑1 (MQ=255)
         ggCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCatcaatca                                                      <  1:3409373/69‑1 (MQ=255)
         ggCTTCTTTACGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCatcaatca                                                      <  1:381974/69‑1 (MQ=255)
         ggCTTCTTTACGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCatcaatca                                                      <  1:214240/69‑1 (MQ=255)
         ggCTTCTTTACGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCatcaatca                                                      <  1:2148475/69‑1 (MQ=255)
         ggCTTCTTTACGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCatcaatca                                                      <  1:2430155/69‑1 (MQ=255)
         ggCTTCTTTACGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCatcaatca                                                      <  1:2553769/69‑1 (MQ=255)
         ggCTTCTTTACGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCatcaatca                                                      <  1:2891626/69‑1 (MQ=255)
         ggCTTCTTTACGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCatcaatca                                                      <  1:3067941/69‑1 (MQ=255)
         ggCTTCTTTACGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCatcaatca                                                      <  1:1983214/69‑1 (MQ=255)
         ggCTTCTTTACGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCatcaatca                                                      <  1:3073641/69‑1 (MQ=255)
         ggCTTCTTTACGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCatcaatca                                                      <  1:3268583/69‑1 (MQ=255)
         ggCTTCTTTACGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCatcaatca                                                      <  1:1815925/69‑1 (MQ=255)
         ggCTTCTTTACGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCatcaatca                                                      <  1:558496/69‑1 (MQ=255)
         ggCTTCTTTACGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCatcaatca                                                      <  1:1118341/69‑1 (MQ=255)
         ggCTTCTTTACGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCatcaatca                                                      <  1:742512/69‑1 (MQ=255)
         ggCTTCTTTACGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCatcaatca                                                      <  1:1655441/69‑1 (MQ=255)
         ggCTTCTTTACGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCatcaatca                                                      <  1:1648490/69‑1 (MQ=255)
         ggCTTCTTTACGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCatcaatca                                                      <  1:1622470/69‑1 (MQ=255)
         ggCTTCTTTACGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCatcaatca                                                      <  1:1605514/69‑1 (MQ=255)
         ggCTTCTTTACGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCatcaatca                                                      <  1:756494/69‑1 (MQ=255)
         ggCTTCTTTACGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCatcaatca                                                      <  1:1091913/69‑1 (MQ=255)
         ggCTTCTTTACGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCatcaatca                                                      <  1:919736/69‑1 (MQ=255)
         ggCTTCTTTACGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCatcaatca                                                      <  1:1386506/69‑1 (MQ=255)
         ggCTTCTTTACGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCatcaatca                                                      <  1:1384481/69‑1 (MQ=255)
           cttcttTACGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCatcaatca                                                      <  1:3017147/67‑1 (MQ=255)
               tttACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAAtctc                                                 >  1:3108293/1‑68 (MQ=255)
               tttACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAAtctc                                                 >  1:3013548/1‑68 (MQ=255)
               tttACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAAtct                                                  >  1:2826558/1‑67 (MQ=255)
                ttACGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCatctatca                                                      <  1:1205875/62‑1 (MQ=255)
                   cGCAACCATTTGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCATCAATCAATCTCCaaa                                             <  1:144375/68‑1 (MQ=255)
                   cGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCATCAATCAATCTCCaaa                                             <  1:2153905/68‑1 (MQ=255)
                   cGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCATCAATCAATCTCCaaa                                             <  1:958451/68‑1 (MQ=255)
                   cGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCATCAATCAATCTCCaaa                                             <  1:155225/68‑1 (MQ=255)
                   cGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCATCAATCAATCTCCaaa                                             <  1:1798653/68‑1 (MQ=255)
                   cGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCATCAATCAATCTCCaaa                                             <  1:3411718/68‑1 (MQ=255)
                   cGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCATCAATCAATCTCCaaa                                             <  1:3310843/68‑1 (MQ=255)
                   cGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCATCAATCAATCTCCaaa                                             <  1:31910/68‑1 (MQ=255)
                   cGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCATCAATCAATCTCCaaa                                             <  1:3151026/68‑1 (MQ=255)
                   cGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCATCAATCAATCTCCaaa                                             <  1:3130573/68‑1 (MQ=255)
                   cGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCATCAATCAATCTCCaaa                                             <  1:1910756/68‑1 (MQ=255)
                   cGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCATCAATCAATCTCCaaa                                             <  1:3046825/68‑1 (MQ=255)
                   cGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCATCAATCAATCTCCaaa                                             <  1:2917462/68‑1 (MQ=255)
                   cGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCATCAATCAATCTCCaaa                                             <  1:2794173/68‑1 (MQ=255)
                   cGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCATCAATCAATCTCCaaa                                             <  1:2141288/68‑1 (MQ=255)
                   cGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCATCAATCAATCTCCaaa                                             <  1:242375/68‑1 (MQ=255)
                   cGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCATCAATCAATCTCCaaa                                             <  1:261134/68‑1 (MQ=255)
                   cGCAACCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTGCTCATCAATCAATCTCCaaa                                             <  1:2140527/68‑1 (MQ=255)
                            tCGTGTTCTTTGCTTTATTTCAGTTCTGCTGGTATGTTg                                                                 >  1:3054055/1‑39 (MQ=38)
                            tCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTg                                                                 >  1:2713695/1‑39 (MQ=255)
                            tCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTg                                                                 >  1:2155647/1‑39 (MQ=255)
                            tCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTg                                                                 >  1:767521/1‑39 (MQ=255)
                            tCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTg                                                                 >  1:1249185/1‑39 (MQ=255)
                            tCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTg                                                                 >  1:2239329/1‑39 (MQ=255)
                            tCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTg                                                                 >  1:374120/1‑39 (MQ=255)
                            tCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTg                                                                 >  1:360643/1‑39 (MQ=255)
                            tCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTg                                                                 >  1:3313110/1‑39 (MQ=255)
                            tCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTg                                                                 >  1:3311261/1‑39 (MQ=255)
                            tCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTg                                                                 >  1:1868376/1‑39 (MQ=255)
                            tCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTg                                                                 >  1:1198834/1‑39 (MQ=255)
                            tCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTg                                                                 >  1:2133992/1‑39 (MQ=255)
                            tCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTg                                                                 >  1:2110562/1‑39 (MQ=255)
                            tCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTg                                                                 >  1:1234594/1‑39 (MQ=255)
                            tCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTg                                                                 >  1:3045324/1‑39 (MQ=255)
                            tCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTg                                                                 >  1:1923910/1‑39 (MQ=255)
                            tCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTATGTTg                                                                 >  1:2726498/1‑39 (MQ=255)
                             cGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCa                                  >  1:670532/1‑69 (MQ=255)
                             cGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCa                                  >  1:2113913/1‑69 (MQ=255)
                             cGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCa                                  >  1:1143205/1‑69 (MQ=255)
                             cGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCa                                  >  1:3294797/1‑69 (MQ=255)
                             cGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCa                                  >  1:2489045/1‑69 (MQ=255)
                             cGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCa                                  >  1:3221104/1‑69 (MQ=255)
                                  tCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGa                             >  1:1273109/1‑69 (MQ=255)
                                  tCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGa                             >  1:2452224/1‑69 (MQ=255)
                                           tATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTaa                                      <  1:1362056/51‑1 (MQ=255)
                                           tATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTaa                                      <  1:2081170/51‑1 (MQ=255)
                                           tATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTaa                                      <  1:807739/51‑1 (MQ=255)
                                            aTTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTaa                                      <  1:889257/50‑1 (MQ=255)
                                                             aCGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCAcgc                               >  1:1356941/1‑40 (MQ=255)
                                                             aCGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCAcgc                               >  1:1033460/1‑40 (MQ=255)
                                                             aCGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCAcgc                               >  1:679648/1‑40 (MQ=255)
                                                             aCGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCAcgc                               >  1:2690327/1‑40 (MQ=255)
                                                             aCGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCAcgc                               >  1:3335177/1‑40 (MQ=255)
                                                             aCGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAg  >  1:1700501/1‑69 (MQ=255)
                                                             aCGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAg  >  1:1677265/1‑69 (MQ=255)
                                                             aCGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAg  >  1:582153/1‑69 (MQ=255)
                                                             aCGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAg  >  1:2857157/1‑69 (MQ=255)
                                                             aCGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAg  >  1:404349/1‑69 (MQ=255)
                                                             aCGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAg  >  1:294384/1‑69 (MQ=255)
                                                             aCGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAg  >  1:203755/1‑69 (MQ=255)
                                                             aCGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAg  >  1:3264515/1‑69 (MQ=255)
                                                              |                                                                   
AGTGCCGTCGGCTTCTTTACGCAGCCATTCGTGTTCTTTGCTGTATTTCAGTTCTGCTGGTACGTTGCTCATCAATCAATCTCCAAAAAAGTAAATCACGCGACGGCTTTGCCGTTACGCACAAAAACAG  >  minE/1948746‑1948875

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: