Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,956,370 C→T 39.0% intergenic (+313/+76) ygfA → / ← serA predicted ligase/D‑3‑phosphoglycerate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,956,3700CT39.0% 31.8 / 40.1 64intergenic (+313/+76)ygfA/serApredicted ligase/D‑3‑phosphoglycerate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (14/25);  new base T (12/13);  total (26/38)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

AACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAG  >  minE/1956305‑1956437
                                                                 |                                                                   
aaCACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAAcgc                                                                   >  1:3010305/1‑68 (MQ=255)
          aTTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAATGCGGCAACGgtgt                                                        <  1:2975225/69‑1 (MQ=255)
          aTTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAATGCGGCAACGgtgt                                                        <  1:1139366/69‑1 (MQ=255)
          aTTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAATGCGGCAACGgtgt                                                        <  1:898893/69‑1 (MQ=255)
          aTTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAATGCGGCAACGgtgt                                                        <  1:2796267/69‑1 (MQ=255)
          aTTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAATGCGGCAACGgtgt                                                        <  1:334285/69‑1 (MQ=255)
          aTTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAATGCGGCAACGgtgt                                                        <  1:882493/69‑1 (MQ=255)
          aTTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAATGCGGCAACGgtgt                                                        <  1:1646551/69‑1 (MQ=255)
          aTTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAATGCGGCAACGgtgt                                                        <  1:2131134/69‑1 (MQ=255)
          aTTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAATGCGGCAACGgtgt                                                        <  1:2729953/69‑1 (MQ=255)
          aTTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAATGCGGCAACGgtgt                                                        <  1:2117385/69‑1 (MQ=255)
          aTTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAATGCGGCAACGgtgt                                                        <  1:1867170/69‑1 (MQ=255)
          aTTGATCTGACTTTGATTTATTTACTGGAGCAGACTCGCAAAGTAGAATGCGCAATGCGGCAACGgtgt                                                        <  1:1369656/69‑1 (MQ=255)
           ttGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAATGCGGCAACGgtgt                                                        <  1:2393865/68‑1 (MQ=255)
             gATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGtgtg                                                       >  1:1078006/1‑67 (MQ=255)
             gATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTgg                                                      >  1:2194217/1‑68 (MQ=255)
             gATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTgg                                                      >  1:3094117/1‑68 (MQ=255)
             gATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTgg                                                      >  1:3181185/1‑68 (MQ=255)
             gATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTgg                                                      >  1:2280904/1‑68 (MQ=255)
             gATCTGACTTTGATTTATTTTCTGGAGCAGAATCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTgg                                                      >  1:542384/1‑68 (MQ=255)
                           ttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATaaaaa                                        <  1:2632099/68‑1 (MQ=255)
                           ttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATaaaaa                                        <  1:1365390/68‑1 (MQ=255)
                                          gACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTg                        <  1:2642915/69‑1 (MQ=255)
                                          gACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTg                        <  1:2938822/69‑1 (MQ=255)
                                          gACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTg                        <  1:3063563/69‑1 (MQ=255)
                                          gACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTg                        <  1:279615/69‑1 (MQ=255)
                                          gACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTg                        <  1:3252674/69‑1 (MQ=255)
                                          gACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTg                        <  1:2244310/69‑1 (MQ=255)
                                          gACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTg                        <  1:2024406/69‑1 (MQ=255)
                                          gACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTg                        <  1:1714953/69‑1 (MQ=255)
                                          gACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTg                        <  1:1649372/69‑1 (MQ=255)
                                          gACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTg                        <  1:1342986/69‑1 (MQ=255)
                                          gACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTg                        <  1:1322169/69‑1 (MQ=255)
                                          gACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTg                        <  1:1104974/69‑1 (MQ=255)
                                                   aGTAGAATGCGCAATGCGGCAACGGTGTGGAGAAGGGATaaaaaaa                                      >  1:1549790/1‑46 (MQ=255)
                                                   aGTAGAATGCGCAATGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGccc                >  1:3286784/1‑68 (MQ=255)
                                                   aGTAGAATGCGCAATGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGccc                >  1:3330322/1‑68 (MQ=255)
                                                   aGTAGAATGCGCAATGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGccc                >  1:783650/1‑68 (MQ=255)
                                                   aGTAGAATGCGCAATGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGccc                >  1:87054/1‑68 (MQ=255)
                                                   aGTAGAATGCGCAATGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGccc                >  1:2684132/1‑68 (MQ=255)
                                                   aGTAGAATGCGCAATGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGccc                >  1:2775569/1‑68 (MQ=255)
                                                   aGTAGAATGCGCAATGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGccc                >  1:1674936/1‑68 (MQ=255)
                                                   aGTAGAATGCGCAATGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGccc                >  1:254868/1‑68 (MQ=255)
                                                   aGTAGAATGCGCAATGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGccc                >  1:2264034/1‑68 (MQ=255)
                                                   aGTAGAATGCGCAATGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGcc                 >  1:1597776/1‑67 (MQ=255)
                                                   aGTAGAATGCGCAATGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGcc                 >  1:2234603/1‑67 (MQ=255)
                                                          tGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGt                             <  1:2395614/48‑1 (MQ=255)
                                                          tGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGt                             <  1:694304/48‑1 (MQ=255)
                                                          tGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGt                             <  1:644736/48‑1 (MQ=255)
                                                          tGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGt                             <  1:874923/48‑1 (MQ=255)
                                                          tGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGt                             <  1:115640/48‑1 (MQ=255)
                                                          tGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGt                             <  1:2412692/48‑1 (MQ=255)
                                                          tGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGt                             <  1:2426341/48‑1 (MQ=255)
                                                          tGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGt                             <  1:288261/48‑1 (MQ=255)
                                                          tGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGt                             <  1:2834635/48‑1 (MQ=255)
                                                          tGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGt                             <  1:2939294/48‑1 (MQ=255)
                                                          tGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAtt         >  1:90395/1‑68 (MQ=255)
                                                          tGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAtt         >  1:1775089/1‑68 (MQ=255)
                                                          tGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAtt         >  1:2120459/1‑68 (MQ=255)
                                                          tGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAtt         >  1:456972/1‑68 (MQ=255)
                                                            cgcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTg                        >  1:2736408/1‑51 (MQ=255)
                                                            cgcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAtt         >  1:8931/1‑66 (MQ=255)
                                                            cgcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAtt         >  1:3294252/1‑66 (MQ=255)
                                                            cgcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAt          >  1:786846/1‑65 (MQ=255)
                                                                aCGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCagag  <  1:371018/69‑1 (MQ=255)
                                                                 |                                                                   
AACACAACGCATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAG  >  minE/1956305‑1956437

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: