Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,811,018 C→T 24.0% P576P (CCG→CCA 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,0180CT24.0% 77.2 / 18.5 62P576P (CCG→CCAuvrAATPase and DNA damage recognition protein of nucleotide excision repair excinuclease UvrABC
Reads supporting (aligned to +/- strand):  ref base C (15/32);  new base T (8/7);  total (23/39)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

ATCTTGCGTTTGCCGCTCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCACATGGTC  >  minE/2810953‑2811083
                                                                 |                                                                 
aTCTTGCGTTTGCCGCTCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCgg                                                                 <  1:3354421/68‑1 (MQ=255)
  cTTGCGTTTGCCGCTCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGAc                                                               >  1:1936213/1‑68 (MQ=255)
        tttGCCGCTCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGc                                                        <  1:2576438/69‑1 (MQ=255)
            ccGCTCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGAc                                                     <  1:3327521/68‑1 (MQ=255)
              gCTCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCAc                                                  <  1:1092708/69‑1 (MQ=255)
               cTCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCgg                                                                 <  1:2504129/53‑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:2548668/69‑1 (MQ=255)
                 cATGTACTGTCCGGTCAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTc                                               <  1:1179268/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)
                 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:3176029/67‑1 (MQ=255)
                            cGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTg                                      <  1:1911007/67‑1 (MQ=255)
                            cGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTg                                      <  1:3408890/67‑1 (MQ=255)
                            cGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTg                                      <  1:1579193/67‑1 (MQ=255)
                            cGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTg                                      <  1:445040/67‑1 (MQ=255)
                            cGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTg                                      <  1:1242613/67‑1 (MQ=255)
                            cGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTg                                      <  1:1175498/67‑1 (MQ=255)
                            cGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTg                                      <  1:2680108/67‑1 (MQ=255)
                                cAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTCACCGCCGTGAacac                                 >  1:1237747/1‑68 (MQ=255)
                                cAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTCACCGCCGTGAACAcc                                >  1:2518271/1‑69 (MQ=255)
                                cAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTCACCGCCGTGAACAcc                                >  1:1830465/1‑69 (MQ=255)
                                cAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTCACCGCCGTGAACAcc                                >  1:2170489/1‑69 (MQ=255)
                                cAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTCACCGCCGTGAACAcc                                >  1:1775243/1‑69 (MQ=255)
                                cAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTCACCGCCGTGAACAcc                                >  1:610811/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:2824844/69‑1 (MQ=255)
                                            cACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGcc                    <  1:2561721/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:2775427/1‑69 (MQ=255)
                                                        cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATcac        >  1:598942/1‑69 (MQ=255)
                                                        cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATcac        >  1:1335812/1‑69 (MQ=255)
                                                        cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATcac        >  1:1645318/1‑69 (MQ=255)
                                                        cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATcac        >  1:3381611/1‑69 (MQ=255)
                                                        cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATcac        >  1:1848275/1‑69 (MQ=255)
                                                        cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATcac        >  1:2040968/1‑69 (MQ=255)
                                                        cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATcac        >  1:3249791/1‑69 (MQ=255)
                                                        cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATca         >  1:1178604/1‑68 (MQ=255)
                                                               aGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCACATGGTc  <  1:1570395/68‑1 (MQ=255)
                                                               aGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCACATGGTc  <  1:3018886/68‑1 (MQ=255)
                                                               aGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCACATGGTc  <  1:1901479/68‑1 (MQ=255)
                                                               aGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCACATGGTc  <  1:3014955/68‑1 (MQ=255)
                                                               aGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCACATGGTc  <  1:2823480/68‑1 (MQ=255)
                                                                 |                                                                 
ATCTTGCGTTTGCCGCTCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCACATGGTC  >  minE/2810953‑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: