Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,950,906 T→C 63.3% L135L (TTA→CTA)  bglJ → DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,950,9060TC63.3% 21.1 / 19.4 38L135L (TTA→CTA) bglJDNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (6/8);  new base C (8/16);  total (14/24)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAAT  >  minE/2950842‑2950958
                                                                |                                                    
aaaTGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAg                                                  >  1:1735234/1‑69 (MQ=255)
   tGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGccc                                               >  1:346775/1‑69 (MQ=255)
       gtACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGCTAAGCCCGa                                             >  1:953660/1‑67 (MQ=255)
       gtACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGCTAAGCCCGa                                             >  1:637548/1‑67 (MQ=255)
       gtACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGCTAAGCCCGa                                             >  1:3430611/1‑67 (MQ=255)
       gtACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGCTAAGCCCGa                                             >  1:1797554/1‑67 (MQ=255)
       gtACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGCTAAGCCCGa                                             >  1:1883098/1‑67 (MQ=255)
       gtACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGCTAAGCCCGa                                             >  1:2814926/1‑67 (MQ=255)
       gtACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGCTAAGCCCGa                                             >  1:2474903/1‑67 (MQ=255)
       gtACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACAAAAGCCGGACGCTAAGCCCGa                                             >  1:2737551/1‑67 (MQ=255)
          cGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGga                                         <  1:65946/68‑1 (MQ=255)
             cAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGCTAAGCCCGACGgaggga                                     <  1:3198932/69‑1 (MQ=255)
             cAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGCTAAGCCCGACGgagaga                                     <  1:1345029/69‑1 (MQ=255)
             cAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGCTAAGCCCGACGgagaga                                     <  1:2926989/69‑1 (MQ=255)
             cAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGCTAAGCCCGACGgagaga                                     <  1:1410383/69‑1 (MQ=255)
             cAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACAAAAGCCGGACGCTAAGCCCGACGgagaga                                     <  1:1704595/69‑1 (MQ=255)
              aGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGgagaga                                     >  1:1946829/1‑68 (MQ=255)
              aGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGgagaga                                     >  1:2518630/1‑68 (MQ=255)
              aGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGCTAAGCCCGACGgtgaga                                     <  1:2015459/68‑1 (MQ=255)
                                            aCATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGaaa                                 <  1:693197/42‑1 (MQ=255)
                                            aCATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGaaa                                 <  1:1788326/42‑1 (MQ=255)
                                            aCATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGaaa                                 <  1:1762095/42‑1 (MQ=255)
                                             cATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTc     >  1:3235170/1‑69 (MQ=255)
                                               ttAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCTCGTGGc          <  1:1542948/62‑1 (MQ=255)
                                                tAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAAt  <  1:2588130/69‑1 (MQ=255)
                                                tAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAAt  <  1:1680683/69‑1 (MQ=255)
                                                 aaCCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAAt  <  1:2383424/68‑1 (MQ=255)
                                                 aaCCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAAt  >  1:2732632/1‑68 (MQ=255)
                                                    cAAAGCCGGACGCTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGt             <  1:1076246/54‑1 (MQ=255)
                                                    cAAAGCCGGACGCTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGt             <  1:445273/54‑1 (MQ=255)
                                                    cAAAGCCGGACGCTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGt             <  1:3376645/54‑1 (MQ=255)
                                                    cAAAGCCGGACGCTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGt             <  1:190879/54‑1 (MQ=255)
                                                    cAAAGCCGGACGCTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGt             <  1:1921280/54‑1 (MQ=255)
                                                    cAAAGCCGGACGCTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGt             <  1:2622999/54‑1 (MQ=255)
                                                    cAAAGCCGGACGCTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGt             <  1:2592548/54‑1 (MQ=255)
                                                    cAAAGCCGGACGCTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGt             <  1:2455381/54‑1 (MQ=255)
                                                    cAAAGCCGGACGCTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGt             <  1:2121187/54‑1 (MQ=255)
                                                    cAAAGCCCGACGCTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGt             <  1:78447/54‑1 (MQ=255)
                                                                |                                                    
AAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAAT  >  minE/2950842‑2950958

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: