Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,902,504 T→A 30.4% I179I (ATA→ATT recB ← exonuclease V (RecBCD complex), beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,902,5040TA30.4% 49.9 / 28.7 46I179I (ATA→ATTrecBexonuclease V (RecBCD complex), beta subunit
Reads supporting (aligned to +/- strand):  ref base T (8/24);  new base A (7/7);  total (15/31)
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.27e-01

TGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGA  >  minE/1902436‑1902572
                                                                    |                                                                    
tGCAGATACCGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCt                                                                      <  1:2497673/69‑1 (MQ=255)
tGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCt                                                                      <  1:2871571/69‑1 (MQ=255)
   aGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAAt                                                                    <  1:3343544/68‑1 (MQ=255)
   aGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAAt                                                                    <  1:3192673/68‑1 (MQ=255)
   aGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAAt                                                                    <  1:415111/68‑1 (MQ=255)
   aGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAAt                                                                    <  1:659324/68‑1 (MQ=255)
   aGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAAt                                                                    <  1:1677939/68‑1 (MQ=255)
   aGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAAt                                                                    <  1:546076/68‑1 (MQ=255)
   aGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAAt                                                                    <  1:1922681/68‑1 (MQ=255)
     aTAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTc                                                                 >  1:2326513/1‑69 (MQ=255)
     aTAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTc                                                                 >  1:2615652/1‑69 (MQ=255)
     aTAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTc                                                                 >  1:1031002/1‑69 (MQ=255)
     aTAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTc                                                                 >  1:77147/1‑69 (MQ=255)
             tAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAATTTCACGCGTCa                                                         >  1:2840337/1‑69 (MQ=255)
             tAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAATTTCACGCGTCa                                                         >  1:3344204/1‑69 (MQ=255)
             tAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAATTTCACGCGTCa                                                         >  1:3273282/1‑69 (MQ=255)
             tAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAATTTCACGCGTCa                                                         >  1:982098/1‑69 (MQ=255)
             tAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAATTTCACGCGTCa                                                         >  1:3191779/1‑69 (MQ=255)
             tAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAATTTCACGCGTCa                                                         >  1:2207575/1‑69 (MQ=255)
             tAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAATTTCACGCGTCa                                                         >  1:2205114/1‑69 (MQ=255)
                      gcagcaACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGc                                                <  1:1422159/69‑1 (MQ=255)
                      gcagcaACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGc                                                <  1:2789219/69‑1 (MQ=255)
                      gcagcaACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGc                                                <  1:1888128/69‑1 (MQ=255)
                      gcagcaACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGc                                                <  1:1513764/69‑1 (MQ=255)
                      gcagcaACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGc                                                <  1:2924623/69‑1 (MQ=255)
                      gcagcaACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGc                                                <  1:3098734/69‑1 (MQ=255)
                      gcagcaACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGc                                                <  1:3414423/69‑1 (MQ=255)
                      gcagcaACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGc                                                <  1:3171391/69‑1 (MQ=255)
                      gcagcaACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGc                                                <  1:2216720/69‑1 (MQ=255)
                      gcagcaACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGc                                                <  1:523077/69‑1 (MQ=255)
                       cagcaACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGc                                                <  1:748325/68‑1 (MQ=255)
                       cagcaACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGc                                                <  1:991288/68‑1 (MQ=255)
                        agcaACGCCTGCGGCCCTTT‑‑AGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGc                                                <  1:242933/65‑1 (MQ=255)
                          caACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGc                                                <  1:3432739/65‑1 (MQ=255)
                               ccTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGCAGTGGCGAc                                       >  1:1221280/1‑69 (MQ=255)
                               ccTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGCAGTGGCGAc                                       >  1:90501/1‑69 (MQ=255)
                               ccTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGCAGTGGCGAc                                       >  1:1225100/1‑69 (MQ=255)
                               ccTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGCAGTGGCGAc                                       >  1:811735/1‑69 (MQ=255)
                                           tCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGc                                                <  1:650533/48‑1 (MQ=255)
                                           tCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGc                                                <  1:1461794/48‑1 (MQ=255)
                                                     aaaGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGc                                                <  1:2130527/38‑1 (MQ=255)
                                                       cgacgacCTGGGCTATTTCACGCGGCAGCGGGTAGc                                                <  1:850194/35‑1 (MQ=255)
                                                        gacgacCTGGGCTATTTCACGCGGCAGCGGGTAGc                                                <  1:3153232/35‑1 (MQ=255)
                                                         acgacCTGGGCTATTTCACGCGGCAGCGGGTAGCAGTGGc                                          <  1:86270/40‑1 (MQ=255)
                                                            acCTGGGCTATTTCACGCGGCATCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGgtagcgt          <  1:503994/69‑1 (MQ=255)
                                                            acCTGGGCTATTTCACGCGGCAGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGgtagcgt          <  1:1537736/69‑1 (MQ=255)
                                                             ccTGGGCTATTTCACGCGGCAGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGgtagcgt          <  1:1951621/68‑1 (MQ=255)
                                                                gggCTATTTCACGCGGCAGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCa      <  1:3287224/69‑1 (MQ=255)
                                                                    tATTTCACGCGGCAGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAgaga  <  1:1105933/69‑1 (MQ=255)
                                                                    |                                                                    
TGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGCAGTGGCGACGCCAGAAATCGGCGCAGGCCTGGTAGCGTAGCAGAGA  >  minE/1902436‑1902572

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: