Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,057,946 A→G 47.8% V257V (GTA→GTG ttdA → L‑tartrate dehydratase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,057,9460AG47.8% ‑2.8 / 38.7 46V257V (GTA→GTGttdAL‑tartrate dehydratase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (14/10);  new base G (3/19);  total (17/29)
Fisher's exact test for biased strand distribution p-value = 2.31e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

AGGACTCAACCGTCTGGGGATTGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTACATATCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGTCTCTACCGGCTGCTGGGCGC  >  minE/2057883‑2058013
                                                               |                                                                   
agGACTCAACCGTCTGGGGATTGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTa                                                                     <  1:725164/64‑1 (MQ=255)
agGACTCAACCGTCTGGGGATTGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTa                                                                     <  1:3079529/64‑1 (MQ=255)
agGACTCAACCGTCTGGGGATTGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTa                                                                     <  1:622608/64‑1 (MQ=255)
agGACTCAACCGTCTGGGGATTGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTACata                                                                 >  1:1043503/1‑68 (MQ=255)
agGACTCAACCGTCTGGGGATTGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTACata                                                                 >  1:2255448/1‑68 (MQ=255)
agGACTCAACCGTCTGGGGATTGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTACata                                                                 >  1:2030258/1‑68 (MQ=255)
agGACTCAACCGTCTGGGGATTGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTACata                                                                 >  1:2587907/1‑68 (MQ=255)
agGACTCAACCGTCTGGGGATTGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTACa                                                                   <  1:469689/66‑1 (MQ=255)
agGACTCAACCGTCTGGGGATTGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTACa                                                                   <  1:1390057/66‑1 (MQ=255)
              tGGGGATTGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTACATATCGAATCTgccgcc                                                   >  1:1198852/1‑68 (MQ=255)
              tGGGGATTGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTACATATCGAATCTGccgccc                                                  >  1:1645533/1‑69 (MQ=255)
              tGGGGATTGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTACATATCGAATCTGccgccc                                                  >  1:3379849/1‑69 (MQ=255)
              tGGGGATTGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTACATATCGAATCTGccgccc                                                  >  1:3318407/1‑69 (MQ=255)
              tGGGGATTGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTACATATCGAATCTGccgccc                                                  >  1:3038221/1‑69 (MQ=255)
              tGGGGATTGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTACATATCGAATCTGccgccc                                                  >  1:1902571/1‑69 (MQ=255)
              tGGGGATTGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTACATATCGAATCTGccgccc                                                  >  1:2970442/1‑69 (MQ=255)
              tGGGGATTGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTACATATCGAATCTGccgccc                                                  >  1:2449759/1‑69 (MQ=255)
              tGGGGATTGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTACATATCGAATCTGccgccc                                                  >  1:95821/1‑69 (MQ=255)
              tGGGGAATGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTACATATCGAATCTGccgccc                                                  >  1:441496/1‑69 (MQ=255)
                     tGGTCCACAAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATcc                                           >  1:445190/1‑69 (MQ=255)
                     tGGTCCACAAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATcc                                           >  1:1068570/1‑69 (MQ=255)
                     tGGTCCACAAGGGCTGACTGGAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATcc                                           >  1:2082076/1‑68 (MQ=255)
                            caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                    <  1:2126044/69‑1 (MQ=255)
                            caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                    <  1:802147/69‑1 (MQ=255)
                            caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                    <  1:756073/69‑1 (MQ=255)
                            caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                    <  1:1329164/69‑1 (MQ=255)
                            caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                    <  1:1484657/69‑1 (MQ=255)
                            caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                    <  1:602792/69‑1 (MQ=255)
                            caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                    <  1:546725/69‑1 (MQ=255)
                            caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                    <  1:1516947/69‑1 (MQ=255)
                            caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                    <  1:1744864/69‑1 (MQ=255)
                            caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                    <  1:3211413/69‑1 (MQ=255)
                            caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                    <  1:3159477/69‑1 (MQ=255)
                            caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                    <  1:3107909/69‑1 (MQ=255)
                            caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                    <  1:2733263/69‑1 (MQ=255)
                            caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                    <  1:1968877/69‑1 (MQ=255)
                            caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                    <  1:2439877/69‑1 (MQ=255)
                            caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                    <  1:2395462/69‑1 (MQ=255)
                            caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                    <  1:2051014/69‑1 (MQ=255)
                                  cTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                    <  1:1964533/63‑1 (MQ=255)
                                        ggCAACAGTTCAGTGATGGGCGTACATATCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCt                        <  1:175529/69‑1 (MQ=255)
                                        ggCAACAGTTCAGTGATGGGCGTACATATCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCt                        <  1:1746616/69‑1 (MQ=255)
                                              aGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                    <  1:3048091/51‑1 (MQ=38)
                                                          ggCGTACATATCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGTCTCTACCGGCTGCTgg      <  1:3046662/69‑1 (MQ=255)
                                                          ggCGTACATATCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGTCTCTACCGGCTGCTgg      <  1:1808226/69‑1 (MQ=255)
                                                          ggCGTACATATCGAATCTGCCGCCCGCCATCCGTCAACCATCGCCGTTGCTGTCTCTACCGGCTGCTgg      <  1:744930/69‑1 (MQ=255)
                                                              tACATATCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGTCTCTACCGGCTGCTGGgcgc  <  1:1597352/69‑1 (MQ=255)
                                                              tACATATCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGTCTCTACCGGCTGCTGGgcgc  <  1:1548813/69‑1 (MQ=255)
                                                              tACATATCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGTCTCTACCGGCTGCTGGgcgc  <  1:1145446/69‑1 (MQ=255)
                                                               aCATATCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGTCTCTACCGGCTGCTGGgcgc  <  1:1963326/68‑1 (MQ=255)
                                                               |                                                                   
AGGACTCAACCGTCTGGGGATTGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTACATATCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGTCTCTACCGGCTGCTGGGCGC  >  minE/2057883‑2058013

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: