Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,902,495 G→C 34.1% V182V (GTC→GTG 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,4950GC34.1% 33.8 / 29.2 41V182V (GTC→GTGrecBexonuclease V (RecBCD complex), beta subunit
Reads supporting (aligned to +/- strand):  ref base G (8/19);  new base C (7/7);  total (15/26)
Fisher's exact test for biased strand distribution p-value = 3.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

GCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGCAGTGGCGAC  >  minE/1902427‑1902535
                                                                    |                                        
gcTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAgacg                                          <  1:474399/69‑1 (MQ=255)
gcTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAgacg                                          <  1:2613601/69‑1 (MQ=255)
gcTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAgacg                                          <  1:1353290/69‑1 (MQ=255)
gcTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAgacg                                          <  1:531135/69‑1 (MQ=255)
gcTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAgacg                                          <  1:2981802/69‑1 (MQ=255)
 cTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAgacga                                         <  1:2510363/69‑1 (MQ=255)
 cTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAgacga                                         <  1:3103627/69‑1 (MQ=255)
 cTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAgacga                                         <  1:1731513/69‑1 (MQ=255)
         tGCAGATACCGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCt                                 <  1:2497673/69‑1 (MQ=255)
         tGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCt                                 <  1:2871571/69‑1 (MQ=255)
            aGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAAt                               <  1:1677939/68‑1 (MQ=255)
            aGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAAt                               <  1:546076/68‑1 (MQ=255)
            aGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAAt                               <  1:1922681/68‑1 (MQ=255)
            aGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAAt                               <  1:3192673/68‑1 (MQ=255)
            aGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAAt                               <  1:3343544/68‑1 (MQ=255)
            aGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAAt                               <  1:415111/68‑1 (MQ=255)
            aGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAAt                               <  1:659324/68‑1 (MQ=255)
              aTAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTc                            >  1:2615652/1‑69 (MQ=255)
              aTAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTc                            >  1:1031002/1‑69 (MQ=255)
              aTAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTc                            >  1:77147/1‑69 (MQ=255)
              aTAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTc                            >  1:2326513/1‑69 (MQ=255)
                      tAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAATTTCACGCGTCa                    >  1:982098/1‑69 (MQ=255)
                      tAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAATTTCACGCGTCa                    >  1:2840337/1‑69 (MQ=255)
                      tAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAATTTCACGCGTCa                    >  1:2207575/1‑69 (MQ=255)
                      tAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAATTTCACGCGTCa                    >  1:2205114/1‑69 (MQ=255)
                      tAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAATTTCACGCGTCa                    >  1:3344204/1‑69 (MQ=255)
                      tAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAATTTCACGCGTCa                    >  1:3191779/1‑69 (MQ=255)
                      tAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACCACCTGGGCAATTTCACGCGTCa                    >  1:3273282/1‑69 (MQ=255)
                               gcagcaACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGc           <  1:3171391/69‑1 (MQ=255)
                               gcagcaACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGc           <  1:523077/69‑1 (MQ=255)
                               gcagcaACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGc           <  1:1422159/69‑1 (MQ=255)
                               gcagcaACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGc           <  1:3414423/69‑1 (MQ=255)
                               gcagcaACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGc           <  1:1513764/69‑1 (MQ=255)
                               gcagcaACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGc           <  1:1888128/69‑1 (MQ=255)
                               gcagcaACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGc           <  1:3098734/69‑1 (MQ=255)
                               gcagcaACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGc           <  1:2924623/69‑1 (MQ=255)
                               gcagcaACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGc           <  1:2789219/69‑1 (MQ=255)
                               gcagcaACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGc           <  1:2216720/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:811735/1‑69 (MQ=255)
                                        ccTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGCAGTGGCGAc  >  1:90501/1‑69 (MQ=255)
                                        ccTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGCAGTGGCGAc  >  1:1225100/1‑69 (MQ=255)
                                        ccTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGCAGTGGCGAc  >  1:1221280/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)
                                                                    |                                        
GCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCCAGGTTTCAAAGACGACCTGGGCTATTTCACGCGGCAGCGGGTAGCAGTGGCGAC  >  minE/1902427‑1902535

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: