Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,811,239 C→A 6.9% A503S (GCG→TCG)  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,2390CA 6.9% 203.4 / 2.2 58A503S (GCG→TCG) uvrAATPase and DNA damage recognition protein of nucleotide excision repair excinuclease UvrABC
Reads supporting (aligned to +/- strand):  ref base C (31/23);  new base A (2/2);  total (33/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.90e-01

TCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGC  >  minE/2811172‑2811309
                                                                   |                                                                      
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:3408779/71‑1 (MQ=255)
           cTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCg                                                          <  1:1128334/71‑1 (MQ=255)
           cTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCg                                                          <  1:1273995/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:2748485/1‑71 (MQ=255)
                           gACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCtg                                          >  1:815522/1‑71 (MQ=255)
                           gACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCtg                                          <  1:2722402/71‑1 (MQ=255)
                            aCGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTc                                                           >  1:510876/1‑53 (MQ=255)
                            aCGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTc                                                           >  1:1874861/1‑53 (MQ=255)
                            aCGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTc                                                           >  1:146977/1‑53 (MQ=255)
                            aCGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCtgt                                         >  1:1998867/1‑71 (MQ=255)
                            aCGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCtgt                                         >  1:2401672/1‑71 (MQ=255)
                             cGGCTCGTCCAGCACGTACATAACGCCAACCAG‑‑‑‑GCACCAATCTGGCTCGCCAGACGGATACGCTGTGCtt                                     >  1:1068363/1‑70 (MQ=255)
                                    tCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCAcc                                 <  1:647019/71‑1 (MQ=255)
                                     ccAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCg                                >  1:382028/1‑71 (MQ=255)
                                     ccAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCg                                >  1:2308063/1‑71 (MQ=255)
                                     ccAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCg                                >  1:1860411/1‑71 (MQ=255)
                                      cAGCACGTACATAACGCCAACCAGGCCCGCACCAATc                                                                 >  1:253769/1‑37 (MQ=255)
                                      cAGCACGTACATAACGCCAACCAGGCCCGCACCAATc                                                                 >  1:1022847/1‑37 (MQ=255)
                                      cAGCACGTACATAACGCCAACCAGGCCCGCACCAATc                                                                 >  1:1093219/1‑37 (MQ=255)
                                      cAGCACGTACATAACGCCAACCAGGCCCGCACCAATc                                                                 >  1:3277001/1‑37 (MQ=255)
                                      cAGCACGTACATAACGCCAACCAGGCCCGCACCAATc                                                                 >  1:1100617/1‑37 (MQ=255)
                                      cAGCACGTACATAACGCCAACCAGGCCCGCACCAATc                                                                 >  1:145387/1‑37 (MQ=255)
                                      cAGCACGTACATAACGCCAACCAGGCCCGCACCAATc                                                                 >  1:1864661/1‑37 (MQ=255)
                                      cAGCACGTACATAACGCCAACCAGGCCCGCACCAATc                                                                 >  1:2424079/1‑37 (MQ=255)
                                      cAGCACGTACATAACGCCAACCAGGCCCGCACCAATc                                                                 >  1:2440669/1‑37 (MQ=255)
                                      cAGCACGTACATAACGCCAACCAGGCCCGCACCAATc                                                                 >  1:524372/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:444979/37‑1 (MQ=255)
                                             tACATAACGCCAACCAGGCCCGCACCAATCTGGCTCg                                                          <  1:99193/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)
                                                          ccAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAg           >  1:2097285/1‑71 (MQ=255)
                                                                ccGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGaa     <  1:455901/71‑1 (MQ=255)
                                                                 cGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGa      >  1:848252/1‑69 (MQ=255)
                                                                   cACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGc  <  1:2416928/71‑1 (MQ=255)
                                                                   |                                                                      
TCGTTATCACGCTGGTGCAGGCCGATAGACGGCTCGTCCAGCACGTACATAACGCCAACCAGGCCCGCACCAATCTGGCTCGCCAGACGGATACGCTGTGCTTCACCGCCAGAAAGCGTTTCTGCCGAGCGGGAAAGC  >  minE/2811172‑2811309

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: