Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,057,964 C→T 36.0% A263A (GCC→GCT ttdA → L‑tartrate dehydratase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,057,9640CT36.0% 40.1 / 33.6 61A263A (GCC→GCTttdAL‑tartrate dehydratase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (14/25);  new base T (3/19);  total (17/44)
Fisher's exact test for biased strand distribution p-value = 7.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGGGATTGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTACATATCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGTCTCTACCGGCTGCTGGGCGCATCGTCGCGG  >  minE/2057897‑2058023
                                                                   |                                                           
tGGGGATTGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTACATATCGAATCTgccgcc                                                             >  1:1198852/1‑68 (MQ=255)
tGGGGATTGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTACATATCGAATCTGccgccc                                                            >  1:95821/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:2970442/1‑69 (MQ=255)
tGGGGATTGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTACATATCGAATCTGccgccc                                                            >  1:2449759/1‑69 (MQ=255)
tGGGGATTGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTACATATCGAATCTGccgccc                                                            >  1:1902571/1‑69 (MQ=255)
tGGGGATTGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTACATATCGAATCTGccgccc                                                            >  1:1645533/1‑69 (MQ=255)
tGGGGAATGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTACATATCGAATCTGccgccc                                                            >  1:441496/1‑69 (MQ=255)
       tGGTCCACAAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATcc                                                     >  1:1068570/1‑69 (MQ=255)
       tGGTCCACAAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATcc                                                     >  1:445190/1‑69 (MQ=255)
       tGGTCCACAAGGGCTGACTGGAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATcc                                                     >  1:2082076/1‑68 (MQ=255)
              caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                              <  1:602792/69‑1 (MQ=255)
              caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                              <  1:2126044/69‑1 (MQ=255)
              caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                              <  1:1329164/69‑1 (MQ=255)
              caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                              <  1:2051014/69‑1 (MQ=255)
              caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                              <  1:3211413/69‑1 (MQ=255)
              caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                              <  1:1968877/69‑1 (MQ=255)
              caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                              <  1:2733263/69‑1 (MQ=255)
              caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                              <  1:546725/69‑1 (MQ=255)
              caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                              <  1:1484657/69‑1 (MQ=255)
              caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                              <  1:2395462/69‑1 (MQ=255)
              caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                              <  1:2439877/69‑1 (MQ=255)
              caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                              <  1:3159477/69‑1 (MQ=255)
              caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                              <  1:1744864/69‑1 (MQ=255)
              caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                              <  1:756073/69‑1 (MQ=255)
              caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                              <  1:1516947/69‑1 (MQ=255)
              caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                              <  1:802147/69‑1 (MQ=255)
              caAGGGCTGACTGGCAACAGTTCAGTGATGGGCGTGCATATCGAATCTGCCGCTCGCCATCCGTCAAcc                                              <  1:3107909/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:1145446/69‑1 (MQ=255)
                                                tACATATCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGTCTCTACCGGCTGCTGGgcgc            <  1:1597352/69‑1 (MQ=255)
                                                tACATATCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGTCTCTACCGGCTGCTGGgcgc            <  1:1548813/69‑1 (MQ=255)
                                                 aCATATCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGTCTCTACCGGCTGCTGGgcgc            <  1:1963326/68‑1 (MQ=255)
                                                     atCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGtctc                             <  1:2541438/47‑1 (MQ=255)
                                                     atCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGtctc                             <  1:2946525/47‑1 (MQ=255)
                                                     atCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGtctc                             <  1:3197523/47‑1 (MQ=255)
                                                     atCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGtctc                             <  1:530806/47‑1 (MQ=255)
                                                     atCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGtctc                             <  1:702090/47‑1 (MQ=255)
                                                     atCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGtctc                             <  1:2086797/47‑1 (MQ=255)
                                                     atCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGtctc                             <  1:1598758/47‑1 (MQ=255)
                                                     atCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGtctc                             <  1:1731203/47‑1 (MQ=255)
                                                     atCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGtctc                             <  1:2784236/47‑1 (MQ=255)
                                                     atCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGtctc                             <  1:1787143/47‑1 (MQ=255)
                                                     atCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGtctc                             <  1:2481488/47‑1 (MQ=255)
                                                     atCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGtctc                             <  1:1928408/47‑1 (MQ=255)
                                                     atCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGtctc                             <  1:1444467/47‑1 (MQ=255)
                                                     atCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGtctc                             <  1:2872006/47‑1 (MQ=255)
                                                      tCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCt                                  >  1:3201020/1‑41 (MQ=255)
                                                      tCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGtctc                             <  1:532275/46‑1 (MQ=255)
                                                      tCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGTCTCTACCGGCTGCTGGGCGCAtcgtc      >  1:1183067/1‑69 (MQ=255)
                                                      tCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGTCTCTACCGGCTGCTGGGCGCAtcgtc      >  1:3141471/1‑69 (MQ=255)
                                                      tCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGTCTCTACCGGCTGCTGGGCGCAtcgtc      >  1:1066989/1‑69 (MQ=255)
                                                        gAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGTCTCTACCGGCTGCTGGGCGCATCGTcgc    >  1:1892725/1‑69 (MQ=255)
                                                        gAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGTCTCTACCGGCTGCTGGGCGCATCGTcgc    >  1:1840650/1‑69 (MQ=255)
                                                        gAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGTCTCTACCGGCTGCTGGGCGCATCGTcgc    >  1:2225108/1‑69 (MQ=255)
                                                        gAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGTCTCTACCGGCTGCTGGGCGCATCGTcgc    >  1:201074/1‑69 (MQ=255)
                                                        gAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGTCTCTACCGGCTGCTGGGCGCATCGTcgc    >  1:1168086/1‑69 (MQ=255)
                                                        gAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGTCTCTACCGGCTGCTGGGCGCATCGTcgc    >  1:2334850/1‑69 (MQ=255)
                                                        gAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGTCTCTACCGGCTGCTGGGCGCATCGTcgc    >  1:1336036/1‑69 (MQ=255)
                                                        gAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGTCTCTACCGGCTGCTGGGCGCATCGTcgc    >  1:299456/1‑69 (MQ=255)
                                                        gAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGTCTCTACCCgctg                    >  1:2249358/1‑53 (MQ=255)
                                                          aTCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGTCTCTACCGGCTGCTGGGCGCATCGTCGCgg  >  1:1356779/1‑69 (MQ=255)
                                                           tCTGCCGCCCGCCATCCGTCAACCTTCGGCGTTGCTGtctc                             <  1:2660152/41‑1 (MQ=255)
                                                                   |                                                           
TGGGGATTGGTCCACAAGGGCTGACCGGCAACAGTTCAGTGATGGGCGTACATATCGAATCTGCCGCCCGCCATCCGTCAACCATCGGCGTTGCTGTCTCTACCGGCTGCTGGGCGCATCGTCGCGG  >  minE/2057897‑2058023

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: