Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,088,571 T→C 34.4% E282E (GAA→GAG yhaJ ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,088,5710TC34.4% 37.8 / 31.2 58E282E (GAA→GAGyhaJpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (21/17);  new base C (5/15);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 5.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTCCGGGCTGAC  >  minE/2088506‑2088636
                                                                 |                                                                 
aaCAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTc                                                                 <  1:1406345/68‑1 (MQ=255)
aaCAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTc                                                                 <  1:2240388/68‑1 (MQ=255)
    ataataAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCa                                                            <  1:2677064/69‑1 (MQ=255)
    ataataAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTCCCCCCa                                                            <  1:3285322/69‑1 (MQ=255)
     taataaATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCAt                                                           <  1:2204017/69‑1 (MQ=255)
          aataTTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCAt                                                           <  1:3260228/64‑1 (MQ=255)
          aataTTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCAt                                                           <  1:469123/64‑1 (MQ=255)
          aataTTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCAt                                                           <  1:549974/64‑1 (MQ=255)
          aataTTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCAt                                                           <  1:235389/64‑1 (MQ=255)
          aataTTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCAt                                                           <  1:1947124/64‑1 (MQ=255)
               tattTTCCGCTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCCTcc                                                                >  1:658538/1‑54 (MQ=255)
               tattTTCCGCTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCCTCCCCCATACTGTCacga                                                 >  1:1046123/1‑69 (MQ=255)
               tattTTCCGCTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCCTCCCCCATACTGTCacga                                                 >  1:2882405/1‑69 (MQ=255)
               tattTTCCGCTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCCTCCCCCATACTGTCacga                                                 >  1:2811477/1‑69 (MQ=255)
               tattTTCCGCTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCCTCCCCCATACTGTCacga                                                 >  1:2307818/1‑69 (MQ=255)
                   ttCCGTTAAAAAGTTTGGGAATTTCCCGCAGCCACCAGGATTTCGCTTCCCCCATACTGTCACGACGcc                                             >  1:2954089/1‑69 (MQ=255)
                   ttCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGcc                                             >  1:2221668/1‑69 (MQ=255)
                      cGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAg                                           <  1:3155856/68‑1 (MQ=255)
                        ttAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCAt                                                           <  1:128825/50‑1 (MQ=255)
                                tttGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCataataat                                <  1:401178/69‑1 (MQ=255)
                                tttGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCataataat                                <  1:1625545/69‑1 (MQ=255)
                                tttGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCataataat                                <  1:1736290/69‑1 (MQ=255)
                                tttGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCataataat                                <  1:2222402/69‑1 (MQ=255)
                                          taccGCAGACACCAGGATTTCGCTTCCCCCATACTGTCac                                                   >  1:3170164/3‑40 (MQ=255)
                                          tCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCac                                                   >  1:588717/1‑40 (MQ=255)
                                          tCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCac                                                   >  1:1067121/1‑40 (MQ=255)
                                          tCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCac                                                   >  1:1194561/1‑40 (MQ=255)
                                          tCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCac                                                   >  1:3226804/1‑40 (MQ=255)
                                          tCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCac                                                   >  1:2160647/1‑40 (MQ=255)
                                          tCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCac                                                   >  1:1261924/1‑40 (MQ=255)
                                                     ccAGGATTTCGCCTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTc           <  1:928809/69‑1 (MQ=255)
                                                     ccAGGATTTCGCCTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTc           <  1:1833889/69‑1 (MQ=255)
                                                     ccAGGATTTCGCCTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTc           <  1:3434742/69‑1 (MQ=255)
                                                     ccAGGATTTCGCCTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTc           <  1:419762/69‑1 (MQ=255)
                                                     ccAGGATTTCGCCTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTc           <  1:3322071/69‑1 (MQ=255)
                                                     ccAGGATTTCGCCTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTc           <  1:1249249/69‑1 (MQ=255)
                                                     ccAGGATTTCGCCTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTc           <  1:275929/69‑1 (MQ=255)
                                                     ccAGGATTTCGCCTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTc           <  1:576746/69‑1 (MQ=255)
                                                     ccAGGATTTCGCCTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTc           <  1:1457285/69‑1 (MQ=255)
                                                     ccAGGATTTCGCCTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTc           <  1:618471/69‑1 (MQ=255)
                                                     ccAGGATTTCGCCTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTc           <  1:2429112/69‑1 (MQ=255)
                                                     ccAGGATTTCGCCTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTc           <  1:1569604/69‑1 (MQ=255)
                                                     ccAGGATTTCGCCTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTc           <  1:1758081/69‑1 (MQ=255)
                                                     ccAGGATTTCGCCTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGGCGATTc           <  1:2940193/69‑1 (MQ=255)
                                                           tttCGCTTCCCCCATACTGTCACGACGCCAGGCCat                                      >  1:700575/1‑36 (MQ=255)
                                                           tttCGCTTCCCCCATACTGTCACGACGCCAGGCCat                                      >  1:432661/1‑36 (MQ=255)
                                                           tttCGCTTCCCCCATACTGTCACGACGCCAGGCCat                                      >  1:275387/1‑36 (MQ=255)
                                                           tttCGCTTCCCCCATACTGTCACGACGCCAGGCCat                                      >  1:3425695/1‑36 (MQ=255)
                                                           tttCGCTTCCCCCATACTGTCACGACGCCAGGCCat                                      >  1:3255428/1‑36 (MQ=255)
                                                           tttCGCTTCCCCCATACTGTCACGACGCCAGGCCat                                      >  1:3146354/1‑36 (MQ=255)
                                                           tttCGCTTCCCCCATACTGTCACGACGCCAGGCCat                                      >  1:2474310/1‑36 (MQ=255)
                                                           tttCGCTTCCCCCATACTGTCACGACGCCAGGCCat                                      >  1:2465629/1‑36 (MQ=255)
                                                           tttCGCTTCCCCCATACTGTCACGACGCCAGGCCat                                      >  1:1959992/1‑36 (MQ=255)
                                                           tttCGCTTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGAtt            >  1:970043/1‑62 (MQ=255)
                                                           tttCGCTTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTCCGGGCt     >  1:1782006/1‑69 (MQ=255)
                                                           tttCGCCTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGGCGATTc           <  1:3283334/63‑1 (MQ=255)
                                                            ttCGCTTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTCCGGGCTg    >  1:1038972/1‑69 (MQ=255)
                                                               gCTTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTCCGGGCTgac  <  1:1825883/68‑1 (MQ=255)
                                                                 |                                                                 
AACAATAATAAATATTATTTTCCGTTAAAAAGTTTGGGAATTTCCCGCAGACACCAGGATTTCGCTTCCCCCATACTGTCACGACGCCAGGCCATAATAATATCGATCTCGCTGGTCGATTCCGGGCTGAC  >  minE/2088506‑2088636

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: