Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,811,009 C→T 31.8% A579A (GCG→GCA 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,0090CT31.8% 39.3 / 19.8 47A579A (GCG→GCAuvrAATPase and DNA damage recognition protein of nucleotide excision repair excinuclease UvrABC
Reads supporting (aligned to +/- strand):  ref base C (5/27);  new base T (8/7);  total (13/34)
Fisher's exact test for biased strand distribution p-value = 1.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCTTGCGTTTGCCGCTCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCAC  >  minE/2810953‑2811077
                                                        |                                                                    
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:113150/69‑1 (MQ=255)
                 cATGTACTGTCCGGTCAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTc                                         <  1:2548668/69‑1 (MQ=255)
                 cATGTACTGTCCGGTCAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTc                                         <  1:1179268/69‑1 (MQ=255)
                 cATGTACTGTCCGGTCAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTc                                         <  1:2266888/69‑1 (MQ=255)
                 cATGTACTGTCCGGTCAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTc                                         <  1:2265202/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:2680108/67‑1 (MQ=255)
                            cGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTg                                <  1:1175498/67‑1 (MQ=255)
                            cGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTg                                <  1:1911007/67‑1 (MQ=255)
                            cGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTg                                <  1:1579193/67‑1 (MQ=255)
                            cGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTg                                <  1:445040/67‑1 (MQ=255)
                            cGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTg                                <  1:371973/67‑1 (MQ=255)
                            cGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTg                                <  1:1242613/67‑1 (MQ=255)
                            cGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTg                                <  1:3176029/67‑1 (MQ=255)
                            cGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTg                                <  1:3408890/67‑1 (MQ=255)
                                cAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTCACCGCCGTGAacac                           >  1:1237747/1‑68 (MQ=255)
                                cAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTCACCGCCGTGAACAcc                          >  1:1775243/1‑69 (MQ=255)
                                cAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTCACCGCCGTGAACAcc                          >  1:833300/1‑69 (MQ=255)
                                cAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTCACCGCCGTGAACAcc                          >  1:2518271/1‑69 (MQ=255)
                                cAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTCACCGCCGTGAACAcc                          >  1:2170489/1‑69 (MQ=255)
                                cAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTCACCGCCGTGAACAcc                          >  1:610811/1‑69 (MQ=255)
                                cAACGACTCCGGCACCGCCATAATTGCTTCCAGTGGACCTTCCGCGACCACTTCACCGCCGTGAACAcc                          >  1:1830465/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:2133438/1‑69 (MQ=255)
                                    gACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTgcg                      >  1:1224031/1‑69 (MQ=255)
                                    gACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTgcg                      >  1:1129511/1‑69 (MQ=255)
                                    gACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTgcg                      >  1:1055770/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:2412864/37‑1 (MQ=255)
                                                      aTCGCTTCCAGCGGACCTTCTGCGACCACTTCAccgc                                    <  1:2171759/37‑1 (MQ=255)
                                                      aTCGCTTCCAGCGGACCTTCTGCGACCACTTCAccgc                                    <  1:1066778/37‑1 (MQ=255)
                                                        cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATcac  >  1:3249791/1‑69 (MQ=255)
                                                        cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATcac  >  1:2775427/1‑69 (MQ=255)
                                                        cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATcac  >  1:3381611/1‑69 (MQ=255)
                                                        cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATcac  >  1:359974/1‑69 (MQ=255)
                                                        cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATcac  >  1:2040968/1‑69 (MQ=255)
                                                        cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATcac  >  1:598942/1‑69 (MQ=255)
                                                        cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATcac  >  1:1848275/1‑69 (MQ=255)
                                                        cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATcac  >  1:1645318/1‑69 (MQ=255)
                                                        cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATcac  >  1:1335812/1‑69 (MQ=255)
                                                        cGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATca   >  1:1178604/1‑68 (MQ=255)
                                                        |                                                                    
ATCTTGCGTTTGCCGCTCATGTACTGCCCGGTCAACGACTCCGGCACCGCCATAATCGCTTCCAGCGGACCTTCTGCGACCACTTCACCGCCGTGAACACCTGCGCCCGGGCCAATGTCGATCAC  >  minE/2810953‑2811077

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: