Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,811,228 A→G 7.1% V506V (GTT→GTC uvrA ← ATPase and DNA damage recognition protein of nucleotide excision repair excinuclease UvrABC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,811,2280AG 7.1% 168.8 / 3.9 56V506V (GTT→GTCuvrAATPase and DNA damage recognition protein of nucleotide excision repair excinuclease UvrABC
Reads supporting (aligned to +/- strand):  ref base A (29/23);  new base G (2/2);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

TACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAA  >  minE/2811158‑2811286
                                                                      |                                                          
tACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCa                                                            <  1:470869/71‑1 (MQ=255)
   ccAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAAcc                                                         <  1:42827/71‑1 (MQ=255)
              tCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCAc                                               >  1:2320166/1‑70 (MQ=255)
                gTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCa                                             <  1:3226517/70‑1 (MQ=255)
                 ttATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCa                                             <  1:3212720/69‑1 (MQ=255)
                    tcgcGCTGGTGCAGGCCGATCGACGGTTCGTCGAGCACGTACATGACGCCGGTCAGGCCCGAACCGATCt                                         <  1:2440271/67‑1 (MQ=255)
                    tcgcGCTGGTGCAGGCCGATCGACGGTTCGTCGAGCACGTACATGACGCCGGTCAGGCCCGAACCGATCt                                         <  1:3387150/67‑1 (MQ=255)
                      aCGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTgg                                       <  1:1199231/70‑1 (MQ=255)
                      aCGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTgg                                       <  1:723193/70‑1 (MQ=255)
                         cTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCg                                   <  1:1128334/71‑1 (MQ=255)
                         cTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCg                                   <  1:1273995/71‑1 (MQ=255)
                         cTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCg                                   <  1:3408779/71‑1 (MQ=255)
                            gTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCa                                <  1:98131/71‑1 (MQ=255)
                            gTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCa                                <  1:804620/71‑1 (MQ=255)
                              gCAGGCCGCTCGATGGCTCGTCCAGCACGTAGACCACCCCGAACAGGCTCGATCGCACCTGGGTCGCCAgg                              >  1:2624099/1‑70 (MQ=255)
                              gCAGGCCGCTCGATGGCTCGTCCAGCACGTAGACCACCCCGAACAGGCTCGATCGCACCTGGGTCGCCAg                               >  1:1541802/1‑70 (MQ=255)
                               cAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCt                                         >  1:1241938/1‑59 (MQ=255)
                                aGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACg                            <  1:2650674/71‑1 (MQ=255)
                                 ggCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACgg                           <  1:3034376/71‑1 (MQ=255)
                                        aGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCt                    >  1:3172595/1‑71 (MQ=255)
                                         gACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCtg                   >  1:815522/1‑71 (MQ=255)
                                         gACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCtg                   <  1:2722402/71‑1 (MQ=255)
                                         gACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCtg                   >  1:2748485/1‑71 (MQ=255)
                                          aCGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTc                                    >  1:1874861/1‑53 (MQ=255)
                                          aCGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTc                                    >  1:146977/1‑53 (MQ=255)
                                          aCGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTc                                    >  1:510876/1‑53 (MQ=255)
                                          aCGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCtgt                  >  1:2401672/1‑71 (MQ=255)
                                          aCGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCtgt                  >  1:1998867/1‑71 (MQ=255)
                                           cGGCTCGTCCAGCACGTACATAACGCCAACCAG‑‑‑‑GCACCAATCTGGCTCGCCAGACGGATACGCTGTGCtt              >  1:1068363/1‑70 (MQ=255)
                                                  tCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCAcc          <  1:647019/71‑1 (MQ=255)
                                                   ccAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCg         >  1:2308063/1‑71 (MQ=255)
                                                   ccAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCg         >  1:1860411/1‑71 (MQ=255)
                                                   ccAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCg         >  1:382028/1‑71 (MQ=255)
                                                    cAGCACGTACATAACGCCAACCAGGCCCGCACCAATc                                          >  1:3277001/1‑37 (MQ=255)
                                                    cAGCACGTACATAACGCCAACCAGGCCCGCACCAATc                                          >  1:2424079/1‑37 (MQ=255)
                                                    cAGCACGTACATAACGCCAACCAGGCCCGCACCAATc                                          >  1:524372/1‑37 (MQ=255)
                                                    cAGCACGTACATAACGCCAACCAGGCCCGCACCAATc                                          >  1:253769/1‑37 (MQ=255)
                                                    cAGCACGTACATAACGCCAACCAGGCCCGCACCAATc                                          >  1:2440669/1‑37 (MQ=255)
                                                    cAGCACGTACATAACGCCAACCAGGCCCGCACCAATc                                          >  1:1864661/1‑37 (MQ=255)
                                                    cAGCACGTACATAACGCCAACCAGGCCCGCACCAATc                                          >  1:145387/1‑37 (MQ=255)
                                                    cAGCACGTACATAACGCCAACCAGGCCCGCACCAATc                                          >  1:1100617/1‑37 (MQ=255)
                                                    cAGCACGTACATAACGCCAACCAGGCCCGCACCAATc                                          >  1:1093219/1‑37 (MQ=255)
                                                    cAGCACGTACATAACGCCAACCAGGCCCGCACCAATc                                          >  1:1022847/1‑37 (MQ=255)
                                                    cAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCa            >  1:2887689/1‑67 (MQ=255)
                                                    cAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGc        >  1:1870083/1‑71 (MQ=255)
                                                         cGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGaa   <  1:1104057/71‑1 (MQ=255)
                                                          gTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGaaa  <  1:2532710/71‑1 (MQ=255)
                                                           tACATAACGCCAACCAGGCCCGCACCAATCTGGCTCg                                   <  1:99193/37‑1 (MQ=255)
                                                           tACATAACGCCAACCAGGCCCGCACCAATCTGGCTCg                                   <  1:444979/37‑1 (MQ=255)
                                                           tACATAACGCCAACCAGGCCCGCACCAATCTGGCTCg                                   <  1:2130111/37‑1 (MQ=255)
                                                           tACATAACGCCAACCAGGCCCGCACCAATCTGGCTCg                                   <  1:1083039/37‑1 (MQ=255)
                                                             cATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGa                          <  1:2541613/44‑1 (MQ=255)
                                                             cATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGa                          <  1:1914330/44‑1 (MQ=255)
                                                               tAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCtgtg                 <  1:1037684/51‑1 (MQ=255)
                                                                 aCGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCt               >  1:3333134/1‑51 (MQ=255)
                                                                 aCGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCt               >  1:2970702/1‑51 (MQ=255)
                                                                 aCGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCt               >  1:775096/1‑51 (MQ=255)
                                                                      |                                                          
TACCCAACAGGCGCTCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAA  >  minE/2811158‑2811286

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: