Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,811,027 T→C 25.3% A573A (GCA→GCG 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,0270TC25.3% 64.9 / 20.2 59A573A (GCA→GCGuvrAATPase and DNA damage recognition protein of nucleotide excision repair excinuclease UvrABC
Reads supporting (aligned to +/- strand):  ref base T (14/30);  new base C (8/7);  total (22/37)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

TTTGCCGCTCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCACATGGTC  >  minE/2810961‑2811083
                                                                  |                                                        
tttGCCGCTCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGc                                                        <  1:2576438/69‑1 (MQ=255)
    ccGCTCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGAc                                                     <  1:3327521/68‑1 (MQ=255)
      gCTCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCAc                                                  <  1:1092708/69‑1 (MQ=255)
        tCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACtt                                                <  1:3287930/69‑1 (MQ=255)
        tCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACtt                                                <  1:2318147/69‑1 (MQ=255)
         cATGTACTGTCCGGTCAACGACTCCGGCACCGGCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTc                                               <  1:6413/69‑1 (MQ=255)
         cATGTACTGTCCGGTCAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTc                                               <  1:2265202/69‑1 (MQ=255)
         cATGTACTGTCCGGTCAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTc                                               <  1:113150/69‑1 (MQ=255)
         cATGTACTGTCCGGTCAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTc                                               <  1:2266888/69‑1 (MQ=255)
         cATGTACTGTCCGGTCAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTc                                               <  1:1179268/69‑1 (MQ=255)
         cATGTACTGTCCGGTCAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTc                                               <  1:2548668/69‑1 (MQ=255)
         cATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACtt                                                <  1:283982/68‑1 (MQ=255)
                   ccGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGc                                                        <  1:3041261/50‑1 (MQ=255)
                    cgttCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTg                                      <  1:1213546/64‑1 (MQ=255)
                    cGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTg                                      <  1:371973/67‑1 (MQ=255)
                    cGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTg                                      <  1:2680108/67‑1 (MQ=255)
                    cGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTg                                      <  1:3176029/67‑1 (MQ=255)
                    cGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTg                                      <  1:3408890/67‑1 (MQ=255)
                    cGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTg                                      <  1:1911007/67‑1 (MQ=255)
                    cGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTg                                      <  1:445040/67‑1 (MQ=255)
                    cGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTg                                      <  1:1579193/67‑1 (MQ=255)
                    cGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTg                                      <  1:1242613/67‑1 (MQ=255)
                    cGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTg                                      <  1:1175498/67‑1 (MQ=255)
                        cAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTCACCGCCGTGAacac                                 >  1:1237747/1‑68 (MQ=255)
                        cAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTCACCGCCGTGAACAcc                                >  1:1830465/1‑69 (MQ=255)
                        cAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTCACCGCCGTGAACAcc                                >  1:1775243/1‑69 (MQ=255)
                        cAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTCACCGCCGTGAACAcc                                >  1:2170489/1‑69 (MQ=255)
                        cAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTCACCGCCGTGAACAcc                                >  1:610811/1‑69 (MQ=255)
                        cAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTCACCGCCGTGAACAcc                                >  1:2518271/1‑69 (MQ=255)
                        cAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTCACCGCCGTGAACAcc                                >  1:833300/1‑69 (MQ=255)
                        cAACGACTCCGGCACCGCCATAATTGCTCCCAGTGGACCTTCCGCGACCACTTCACCGCCGTGAACAcc                                >  1:2855219/1‑69 (MQ=255)
                        cAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCAc                                                  <  1:1054167/51‑1 (MQ=255)
                            gACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTc                                               <  1:626308/50‑1 (MQ=255)
                            gACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTgcg                            >  1:1055770/1‑69 (MQ=255)
                            gACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTgcg                            >  1:1129511/1‑69 (MQ=255)
                            gACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTgcg                            >  1:1224031/1‑69 (MQ=255)
                            gACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTgcg                            >  1:2133438/1‑69 (MQ=255)
                                    cACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGcc                    <  1:2561721/69‑1 (MQ=255)
                                    cACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGcc                    <  1:2824844/69‑1 (MQ=255)
                                    cACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGcc                    <  1:3350900/69‑1 (MQ=255)
                                              aTCGCTTCCAGCGGACCTTCTGCGACCACTTCAccgc                                          <  1:978447/37‑1 (MQ=255)
                                              aTCGCTTCCAGCGGACCTTCTGCGACCACTTCAccgc                                          <  1:1066778/37‑1 (MQ=255)
                                              aTCGCTTCCAGCGGACCTTCTGCGACCACTTCAccgc                                          <  1:2412864/37‑1 (MQ=255)
                                              aTCGCTTCCAGCGGACCTTCTGCGACCACTTCAccgc                                          <  1:2171759/37‑1 (MQ=255)
                                                cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATcac        >  1:359974/1‑69 (MQ=255)
                                                cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATcac        >  1:1335812/1‑69 (MQ=255)
                                                cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATcac        >  1:1645318/1‑69 (MQ=255)
                                                cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATcac        >  1:598942/1‑69 (MQ=255)
                                                cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATcac        >  1:1848275/1‑69 (MQ=255)
                                                cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATcac        >  1:2040968/1‑69 (MQ=255)
                                                cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATcac        >  1:3381611/1‑69 (MQ=255)
                                                cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATcac        >  1:2775427/1‑69 (MQ=255)
                                                cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATcac        >  1:3249791/1‑69 (MQ=255)
                                                cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATca         >  1:1178604/1‑68 (MQ=255)
                                                       aGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCACATGGTc  <  1:3018886/68‑1 (MQ=255)
                                                       aGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCACATGGTc  <  1:1901479/68‑1 (MQ=255)
                                                       aGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCACATGGTc  <  1:1570395/68‑1 (MQ=255)
                                                       aGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCACATGGTc  <  1:2823480/68‑1 (MQ=255)
                                                       aGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCACATGGTc  <  1:3014955/68‑1 (MQ=255)
                                                                  |                                                        
TTTGCCGCTCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCACATGGTC  >  minE/2810961‑2811083

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: