Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,898,833 C→T 60.7% Q222Q (CAG→CAA recD ← exonuclease V (RecBCD complex), alpha chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,898,8330CT60.7% 25.2 / 40.4 61Q222Q (CAG→CAArecDexonuclease V (RecBCD complex), alpha chain
Reads supporting (aligned to +/- strand):  ref base C (14/10);  new base T (28/9);  total (42/19)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

GGTGCAAAGTGCTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCA  >  minE/1898775‑1898893
                                                          |                                                            
ggTGCAAAGTGCTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAg                                                    <  1:2753446/69‑1 (MQ=255)
ggTGCAAAGTGCTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAg                                                    <  1:744643/69‑1 (MQ=255)
ggTGCAAAGTGCTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAg                                                    <  1:2698583/69‑1 (MQ=255)
ggTGCAAAGTGCTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAg                                                    <  1:2123207/69‑1 (MQ=255)
  tGCAAAGTGCTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCTCAAAGc                                                   <  1:911/68‑1 (MQ=255)
  tGCAAAGTGCTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGc                                                   <  1:3274919/68‑1 (MQ=255)
  tGCAAAGTGCTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGc                                                   <  1:2954325/68‑1 (MQ=255)
  tGCAAAGTGCTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGc                                                   <  1:668408/68‑1 (MQ=255)
  tGCAAAGTGCTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGc                                                   <  1:660176/68‑1 (MQ=255)
  tGCAAAGTGCTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGc                                                   <  1:174238/68‑1 (MQ=255)
        gTGCTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTg                                               >  1:1821496/1‑66 (MQ=255)
           cTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGcc                                             >  1:1705326/1‑65 (MQ=255)
                aTCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGCCGAGAGATTc                                    >  1:1436245/1‑69 (MQ=255)
                 tcctCCGGAATGCGTTTCTTTTGTTCATGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg                                   >  1:2784127/4‑68 (MQ=255)
                 tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTGAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTc                                    >  1:2106968/1‑68 (MQ=255)
                 tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTc                                    >  1:2159918/1‑68 (MQ=255)
                 tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTc                                    >  1:219899/1‑68 (MQ=255)
                 tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTc                                    >  1:2556396/1‑68 (MQ=255)
                 tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTc                                    >  1:593703/1‑68 (MQ=255)
                 tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTc                                    >  1:35772/1‑68 (MQ=255)
                 tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTc                                    >  1:1766312/1‑68 (MQ=255)
                 tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTc                                    >  1:1874844/1‑68 (MQ=255)
                 tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTc                                    >  1:783627/1‑68 (MQ=255)
                 tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg                                   >  1:2564719/1‑69 (MQ=255)
                 tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg                                   >  1:107544/1‑69 (MQ=255)
                 tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg                                   >  1:3057903/1‑69 (MQ=255)
                 tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg                                   >  1:2974366/1‑69 (MQ=255)
                 tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg                                   >  1:115358/1‑69 (MQ=255)
                 tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg                                   >  1:2897097/1‑69 (MQ=255)
                 tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg                                   >  1:1165060/1‑69 (MQ=255)
                 tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg                                   >  1:1555350/1‑69 (MQ=255)
                 tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg                                   >  1:1943514/1‑69 (MQ=255)
                 tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg                                   >  1:1058904/1‑69 (MQ=255)
                 tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg                                   >  1:2437037/1‑69 (MQ=255)
                 tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg                                   >  1:2435756/1‑69 (MQ=255)
                 tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg                                   >  1:2387649/1‑69 (MQ=255)
                 tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCg                                   >  1:845494/1‑69 (MQ=255)
                 tCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGAAATTGTCGCAAAGCCTTGCCGAGAGATTCg                                   >  1:3080494/1‑69 (MQ=255)
                   ttCCGGAATGCGTTTCTTTTGTTCATCGGTTAGCGGTAACTGTCGCAAAGCCTTGCCGAGAGATTCgg                                  >  1:135571/1‑68 (MQ=255)
                   ttCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGCCGAGAGATTCgg                                  >  1:2365755/1‑68 (MQ=255)
                   ttCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGCCGAGAGATTCg                                   >  1:658260/1‑67 (MQ=255)
                                 tCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGCCgagag                                        >  1:2120042/1‑48 (MQ=255)
                                 tCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGCCGAGAGATTCGGTTaa                              >  1:2168670/1‑58 (MQ=255)
                                 tCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAg                   >  1:795182/1‑69 (MQ=255)
                                 tCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAg                   >  1:2443657/1‑69 (MQ=255)
                                 tCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAg                   >  1:2266127/1‑69 (MQ=255)
                                 tCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAg                   >  1:2230096/1‑69 (MQ=255)
                                 tCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAg                   >  1:1966511/1‑69 (MQ=255)
                                 tCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAg                   >  1:1141885/1‑69 (MQ=255)
                                            tCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGtt        >  1:257100/1‑69 (MQ=255)
                                            tCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGtt        >  1:213213/1‑69 (MQ=255)
                                            tCGGTCAGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGtt        >  1:914267/1‑69 (MQ=255)
                                                  aGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCa  <  1:3423337/69‑1 (MQ=255)
                                                  aGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCa  <  1:271726/69‑1 (MQ=255)
                                                  aGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCa  <  1:509073/69‑1 (MQ=255)
                                                  aGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCa  <  1:2704542/69‑1 (MQ=255)
                                                  aGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCa  <  1:2314643/69‑1 (MQ=255)
                                                  aGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCa  <  1:714867/69‑1 (MQ=255)
                                                  aGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCa  <  1:718865/69‑1 (MQ=255)
                                                  aGCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCa  <  1:2044983/69‑1 (MQ=255)
                                                   gCGGTAATTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCa  <  1:3380302/68‑1 (MQ=255)
                                                          |                                                            
GGTGCAAAGTGCTGGCATCTTCCGGAATGCGTTTCTTTTGTTCATCGGTCAGCGGTAACTGTCGCAAAGCCTTGCCGAGAGATTCGGTTAAGCGCGCGGCAGCTTTACCCGTTGGTGCA  >  minE/1898775‑1898893

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: