Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 464,269 A→G 40.9% L42L (CTA→CTG sdhB → succinate dehydrogenase, FeS subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE464,2690AG40.9% 22.2 / 47.9 66L42L (CTA→CTGsdhBsuccinate dehydrogenase, FeS subunit
Reads supporting (aligned to +/- strand):  ref base A (17/22);  new base G (10/17);  total (27/39)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

CAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACGG  >  minE/464204‑464334
                                                                 |                                                                 
cAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTaa                                                                  >  1:974140/1‑67 (MQ=255)
cAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTaa                                                                  >  1:96890/1‑67 (MQ=255)
cAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTaa                                                                  >  1:610379/1‑67 (MQ=255)
cAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTaa                                                                  >  1:431558/1‑67 (MQ=255)
cAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTaa                                                                  >  1:398070/1‑67 (MQ=255)
cAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTaa                                                                  >  1:3377841/1‑67 (MQ=255)
cAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTaa                                                                  >  1:2234041/1‑67 (MQ=255)
cAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTaa                                                                  >  1:1480450/1‑67 (MQ=255)
cAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTa                                                                   >  1:1055595/1‑66 (MQ=255)
cAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGATaa                                                                  >  1:2544083/1‑67 (MQ=255)
cAGGATTACACCCTAGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTaa                                                                  >  1:1262063/1‑67 (MQ=255)
     ttACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAaga                                                              >  1:288879/1‑66 (MQ=255)
     ttACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAaga                                                              >  1:2362094/1‑66 (MQ=255)
     ttACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAaga                                                              >  1:2349463/1‑66 (MQ=255)
     ttACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAaga                                                              >  1:3271268/1‑66 (MQ=255)
     ttACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAaga                                                              >  1:330986/1‑66 (MQ=255)
     ttACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAaga                                                              >  1:1650151/1‑66 (MQ=255)
     ttACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAaga                                                              >  1:1718589/1‑66 (MQ=255)
       acacCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTaa                                                                  >  1:2594964/1‑60 (MQ=255)
       acacCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTaa                                                                  >  1:1830372/1‑60 (MQ=255)
       acacCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTaa                                                                  >  1:1971804/1‑60 (MQ=255)
       acacCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTaa                                                                  >  1:1746848/1‑60 (MQ=255)
       acacCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTaa                                                                  >  1:398401/1‑60 (MQ=255)
          cccTGGAAGCGGAAGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATTCAGCTGAAAGAGGAAGAt                                                       >  1:3375203/1‑68 (MQ=255)
          cccTGGAAGCGGAAGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATTCAGCTGAAAGAGAAAGAt                                                       >  1:2165846/1‑68 (MQ=255)
          cccTGGAAGCGGAAGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATTCAGCTGAAAGAGAAAGAt                                                       >  1:840843/1‑68 (MQ=255)
          cccTGGAAGCGGAAGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATTCAGCTGAAAGAGAAAGAt                                                       >  1:2192077/1‑68 (MQ=255)
          cccTGGAAGCGGAAGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATTCAGCTGAAAGAGAAAGAt                                                       >  1:872885/1‑68 (MQ=255)
          cccTGGAAGCGGAAGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATTCAGCTGAAAGAGAAAGAt                                                       >  1:3191684/1‑68 (MQ=255)
          cccTGGAAGCGGAAGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATTCAGCTGAAAGAGAAAGAt                                                       >  1:959612/1‑68 (MQ=255)
          cccTGGAAGCGGAAGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATTCAGCTGAAAGAGAAAGAt                                                       >  1:315961/1‑68 (MQ=255)
          cccTGGAAGCGGAAGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATTCAGATGAAAGAGAAAGAt                                                       >  1:1921481/1‑68 (MQ=255)
          cccTGCAAGCGGAAGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATTCAGCTGAAAGAGAAAGAt                                                       >  1:2239057/1‑68 (MQ=255)
                 aGCGGAAGATGGTCGCGACATGATGCTGCTGGATGCGCTTATTCAGCTGAAAGTGAAAGATCCCAGc                                                 <  1:3390465/67‑1 (MQ=255)
                 aGCGGAAGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATTCAGCTGAAAGAGAAAGATCCCAGc                                                 <  1:2874834/67‑1 (MQ=255)
                 aGCGGAAGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATTCAGCTGAAAGAGAAAGATCCCAGc                                                 <  1:1037290/67‑1 (MQ=255)
                 aGCGGAAGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATTCAGCTGAAAGAGAAAGATCCCAGc                                                 <  1:2825482/67‑1 (MQ=255)
                 aGCGGAAGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATTCAGCTGAAAGAGAAAGATCCCAGc                                                 <  1:1157106/67‑1 (MQ=255)
                 aGCGGAAGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATTCAGCTGAAAGAGAAAGATCCCAGc                                                 <  1:769529/67‑1 (MQ=255)
                 aGCGGAAGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATTCAGCTGAAAGAGAAAGATCCCAGc                                                 <  1:81480/67‑1 (MQ=255)
                 aGCGGAAGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATTCAGCTGAAAGAGAAAGATCCCAGc                                                 <  1:1097082/67‑1 (MQ=255)
                 aGCGGAAGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATTCAGCTGAAAGAGAAAGATCCCAGc                                                 <  1:2799478/67‑1 (MQ=255)
                 aGCGGAAGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATTCAGCTGAAAGAGAAAGATCCCAGc                                                 <  1:1493266/67‑1 (MQ=255)
                 aGCGGAAGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATTCAGCTGAAAGAGAAAGATCCCAGc                                                 <  1:1570583/67‑1 (MQ=255)
                 aGCGGAAGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATTCAGCTGAAAGAGAAAGATCCCAGc                                                 <  1:2154083/67‑1 (MQ=255)
                 aGCGGAAGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATTCAGCTGAAAGAGAAAGATCCCAGc                                                 <  1:2099582/67‑1 (MQ=255)
                 aGCGGAAGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATTCAGCTGAAAGAGAAAGATCCCAGc                                                 <  1:1756769/67‑1 (MQ=255)
                 aGCGGAAGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATTCAGCTGAAAGAGAAAGATCCCAGc                                                 <  1:2226435/67‑1 (MQ=255)
                  gCGGAAGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATTCAGCTGAAAGAGAAAGATCCCAGc                                                 <  1:342497/66‑1 (MQ=255)
                   cGGAAGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATTCAGCTGAAAGAGAAAGATCCCAGc                                                 <  1:839980/65‑1 (MQ=255)
                          aGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTccg                                      >  1:709383/1‑69 (MQ=255)
                          aGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTccg                                      >  1:2436739/1‑69 (MQ=255)
                          aGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTccg                                      >  1:49810/1‑69 (MQ=255)
                          aGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTccg                                      >  1:3369704/1‑69 (MQ=255)
                          aGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTccg                                      >  1:2107991/1‑69 (MQ=255)
                                    atgatgCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGc                            <  1:60508/69‑1 (MQ=255)
                                    atgatgCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGc                            <  1:791994/69‑1 (MQ=255)
                                    atgatgCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGc                            <  1:558553/69‑1 (MQ=255)
                                    atgatgCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGc                            <  1:1145782/69‑1 (MQ=255)
                                    atgatgCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGc                            <  1:946886/69‑1 (MQ=255)
                                    atgatgCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGc                            <  1:1458258/69‑1 (MQ=255)
                                    atgatgCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGc                            <  1:3319099/69‑1 (MQ=255)
                                    atgatgCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGc                            <  1:3003922/69‑1 (MQ=255)
                                    atgatgCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGc                            <  1:2347163/69‑1 (MQ=255)
                                     tgatgCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGc                            <  1:2516886/68‑1 (MQ=255)
                                       atgCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCg                          >  1:1161396/1‑68 (MQ=255)
                                       atgCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCg                          >  1:2080452/1‑68 (MQ=255)
                                         gctgctGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTg                        <  1:3383054/68‑1 (MQ=255)
                                         gctgctGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGa                       >  1:1107433/1‑69 (MQ=255)
                                         gctgctGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGa                       >  1:1625766/1‑69 (MQ=255)
                                         gctgctGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGa                       >  1:3041074/1‑69 (MQ=255)
                                                   gcgcTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTg                        <  1:578626/58‑1 (MQ=255)
                                                   gcgcTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTg                        <  1:2189278/58‑1 (MQ=255)
                                                      cTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCAGCCGTg                        <  1:3127780/55‑1 (MQ=255)
                                                              gCTAAAAGAGAAAGATCCCAGCCTGTCTTTCCGCCGCTCCTGCCGTTAAGGTGTGTGCGGTTCCTACgg  <  1:2895926/69‑1 (MQ=255)
                                                              gCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACgg  <  1:2611890/69‑1 (MQ=255)
                                                              gCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACgg  <  1:2882030/69‑1 (MQ=255)
                                                              gCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACgg  <  1:2568207/69‑1 (MQ=255)
                                                              gCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACgg  <  1:3417463/69‑1 (MQ=255)
                                                              gCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACgg  <  1:16741/69‑1 (MQ=255)
                                                              gCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACgg  <  1:2166574/69‑1 (MQ=255)
                                                              gCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACgg  <  1:351788/69‑1 (MQ=255)
                                                                 |                                                                 
CAGGATTACACCCTGGAAGCGGATGAAGGTCGCGACATGATGCTGCTGGATGCGCTTATCCAGCTAAAAGAGAAAGATCCCAGCCTGTCGTTCCGCCGCTCCTGCCGTGAAGGTGTGTGCGGTTCCGACGG  >  minE/464204‑464334

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: