Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,325,781 T→C 36.3% T410T (ACA→ACG wcaK ← predicted pyruvyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,325,7810TC36.3% 23.5 / 24.3 44T410T (ACA→ACGwcaKpredicted pyruvyl transferase
Reads supporting (aligned to +/- strand):  ref base T (16/12);  new base C (9/7);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

AGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTA  >  minE/1325721‑1325847
                                                            |                                                                  
agCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACg                                                            >  1:2888094/1‑69 (MQ=255)
agCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACg                                                            >  1:2346645/1‑69 (MQ=255)
agCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACg                                                            >  1:1699070/1‑69 (MQ=255)
       cTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGa                                                     >  1:1100955/1‑69 (MQ=255)
       cTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGa                                                     >  1:2431438/1‑69 (MQ=255)
       cTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGa                                                     >  1:3140069/1‑69 (MQ=255)
             ttCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGAc                                                <  1:1681056/68‑1 (MQ=255)
             ttCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGAc                                                <  1:637596/68‑1 (MQ=255)
             ttCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGAc                                                <  1:3001321/68‑1 (MQ=255)
             ttCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGAc                                                <  1:3052545/68‑1 (MQ=255)
             ttCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGAc                                                <  1:1950095/68‑1 (MQ=255)
             ttCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGAc                                                <  1:3425260/68‑1 (MQ=255)
             ttCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGAc                                                <  1:1526288/68‑1 (MQ=255)
              tCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCt                                                          <  1:2280132/57‑1 (MQ=255)
              tCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCt                                                          <  1:136255/57‑1 (MQ=255)
                  cTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTCGCGACTGACAGCTTc                                          <  1:3288939/69‑1 (MQ=255)
                  cTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTCGCGACTGACAGCTTc                                          <  1:551730/69‑1 (MQ=255)
                  cTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTCGCGACTGACAGCTTc                                          <  1:2101050/69‑1 (MQ=255)
                  cTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTCGCGACTGACAGCTTc                                          <  1:385639/69‑1 (MQ=255)
                  cTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTCGCGACTGACAGCTTc                                          <  1:2904518/69‑1 (MQ=255)
                  cTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTCGCGACTGACAGCTTc                                          <  1:2507002/69‑1 (MQ=255)
                                           tGCACCATCTGCATTCCCGTCTGACGTTCGCGACTGACAGCTTc                                          <  1:1553656/44‑1 (MQ=255)
                                                    tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                             >  1:3065589/1‑48 (MQ=255)
                                                    tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                             >  1:2455628/1‑48 (MQ=255)
                                                    tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                             >  1:2757351/1‑48 (MQ=255)
                                                    tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                             >  1:1681624/1‑48 (MQ=255)
                                                    tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                             >  1:1804963/1‑48 (MQ=255)
                                                    tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                             >  1:2038858/1‑48 (MQ=255)
                                                    tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                             >  1:1809017/1‑48 (MQ=255)
                                                      cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa      >  1:2812535/1‑69 (MQ=255)
                                                      cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa      >  1:706214/1‑69 (MQ=255)
                                                      cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa      >  1:1127694/1‑69 (MQ=255)
                                                      cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa      >  1:3408956/1‑69 (MQ=255)
                                                      cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa      >  1:3279818/1‑69 (MQ=255)
                                                      cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa      >  1:1876401/1‑69 (MQ=255)
                                                      cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa      >  1:219926/1‑69 (MQ=255)
                                                      cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa      >  1:1029790/1‑69 (MQ=255)
                                                      cATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAGCGCCGGAAGCTGGCCTaa      >  1:2850444/1‑69 (MQ=255)
                                                        ttCCTGTCTGACGCTCACTACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCt        <  1:2404167/65‑1 (MQ=255)
                                                          ccTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTa  >  1:3178922/1‑69 (MQ=255)
                                                          ccTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTa  >  1:3304678/1‑69 (MQ=255)
                                                          ccTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTa  <  1:1925704/69‑1 (MQ=255)
                                                          ccTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTa  <  1:859209/69‑1 (MQ=255)
                                                          ccTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGACGGAAGCTGGCCTAAGGTa  >  1:1520830/1‑69 (MQ=255)
                                                            |                                                                  
AGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTA  >  minE/1325721‑1325847

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: