Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 161,956 T→C 46.7% D223D (GAT→GAC cdaR → DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE161,9560TC46.7% 4.0 / 63.7 75D223D (GAT→GACcdaRDNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (24/16);  new base C (17/18);  total (41/34)
Fisher's exact test for biased strand distribution p-value = 3.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CTCGCTAACCGAAATGGTGGTGTTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGATCATCGTAAGCGAGTTGAACAACTGATTACCCGCATGAAAGAGTACGGCCAGCTGC  >  minE/161899‑162020
                                                         |                                                                
ctcGCTAACCGAAATGGTGGTGTTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCagaag                                                           >  1:1016952/1‑65 (MQ=255)
ctcGCTAACCGAAATGGTGGTGTTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCagaag                                                           >  1:3144395/1‑65 (MQ=255)
ctcGCTAACCGAAATGGTGGTGTTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCagaag                                                           >  1:2137870/1‑65 (MQ=255)
ctcGCTAACCGAAATGGTGGTGTTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCagaag                                                           >  1:817931/1‑65 (MQ=255)
ctcCCTAACCGAAATGGTGGTGTTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCagaag                                                           >  1:2550868/1‑65 (MQ=255)
 tcGCTAACCGAAATGGTGGTGTTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGatcat                                                      >  1:3208446/1‑69 (MQ=255)
 tcGCTAACCGAAATGGTGGTGTTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGatcat                                                      >  1:3399258/1‑69 (MQ=255)
 tcGCTAACCGAAATGGTGGTGTTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGatcat                                                      >  1:1857336/1‑69 (MQ=255)
 tcGCTAACCGAAATGGTGGTGTTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGatcat                                                      >  1:787374/1‑69 (MQ=255)
 tcGCTAACCGAAATGGTGGTGTTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGatcat                                                      >  1:3104148/1‑69 (MQ=255)
 tcGCTAACCGAAATGGTGGTGTTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGatcat                                                      >  1:3100105/1‑69 (MQ=255)
 tcGCTAACCGAAATGGTGGTGTTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGatcat                                                      >  1:2408084/1‑69 (MQ=255)
 tcGCTAACCGAAATGGTGGTGTTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGatcat                                                      >  1:3399572/1‑69 (MQ=255)
 tcGCTAACCGAAATGGTGGTGTTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGatcat                                                      >  1:1387265/1‑69 (MQ=255)
 tcGCTAACCGAAATGGTGGTGTTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGatcat                                                      >  1:609861/1‑69 (MQ=255)
 tcGCTAACCGAAATGGTGGTGTTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGatcat                                                      >  1:2456012/1‑69 (MQ=255)
         cGAAATGGTGGTGTTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGATCATCGTAAGCg                                              <  1:1307921/69‑1 (MQ=255)
              tggtggTGTTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGATCATCGTAAGCGAGTTg                                         <  1:2901445/69‑1 (MQ=255)
              tggtggTGTTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGATCATCGTAAGCGAGTTg                                         <  1:1567778/69‑1 (MQ=255)
              tggtggTGTTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGATCATCGTAAGCGAGTTg                                         <  1:1563641/69‑1 (MQ=255)
              tggtggTGTTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGATCATCGTAAGCGAGTTg                                         <  1:470370/69‑1 (MQ=255)
               ggtggtGTTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGATCATCGTAAGCGAGTTg                                         <  1:1239581/68‑1 (MQ=255)
                  ggtgTTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGATCATCGTAAGCGAGTTGaac                                      <  1:1784634/68‑1 (MQ=255)
                  ggtgTTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGATCATCGTAAGCGAGTTGaac                                      <  1:2290957/68‑1 (MQ=255)
                   gtgtTGAAACCGGCGTTGAACTCTTTTGGGCGCTTGGATGCAGAAGATCATCGTAAGCGAGTTGaacaa                                    <  1:1341022/69‑1 (MQ=255)
                   gtgtTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGATCATCGTAAGCGAGTTGaacaa                                    <  1:3381441/69‑1 (MQ=255)
                   gtgtTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGATCATCGTAAGCGAGTTGaacaa                                    <  1:3113960/69‑1 (MQ=255)
                   gtgtTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGATCATCGTAAGCGAGTTGaacaa                                    <  1:1680750/69‑1 (MQ=255)
                    tgtTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGATCATCGTAAGCGAGTTGaacaa                                    <  1:1491589/68‑1 (MQ=255)
                     gtTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGaaca                                     <  1:34036/66‑1 (MQ=255)
                     gtTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGaaca                                     <  1:1182458/66‑1 (MQ=255)
                     gtTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGaaca                                     <  1:1189785/66‑1 (MQ=255)
                     gtTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGaaca                                     <  1:1222957/66‑1 (MQ=255)
                     gtTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGaaca                                     <  1:875442/66‑1 (MQ=255)
                     gtTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGaaca                                     <  1:1337875/66‑1 (MQ=255)
                     gtTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGaaca                                     <  1:2871501/66‑1 (MQ=255)
                     gtTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGaaca                                     <  1:1997506/66‑1 (MQ=255)
                     gtTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGaaca                                     <  1:395956/66‑1 (MQ=255)
                     gtTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGaaca                                     <  1:469266/66‑1 (MQ=255)
                     gtTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGaaca                                     <  1:2469862/66‑1 (MQ=255)
                     gtTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGaaca                                     <  1:204532/66‑1 (MQ=255)
                     gtTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGaaca                                     <  1:2335527/66‑1 (MQ=255)
                     gtTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGaaca                                     <  1:2377756/66‑1 (MQ=255)
                     gtTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGaaca                                     <  1:639488/66‑1 (MQ=255)
                     gtTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGaaca                                     <  1:817154/66‑1 (MQ=255)
                     gtTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGaaca                                     <  1:1149298/66‑1 (MQ=255)
                      ttGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGaaca                                     <  1:3274851/65‑1 (MQ=255)
                        gAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGATCATCGTAAGCGAGTTGAACAACTGAt                               >  1:2916485/1‑69 (MQ=255)
                        gAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGATCATCGTAAGCGAGTTGAACAACTGAt                               >  1:1356623/1‑69 (MQ=255)
                        gAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGATCATCGTAAGCGAGTTGAACAACTGAt                               >  1:699764/1‑69 (MQ=255)
                                    gAACTCTTTTGGGCGCTGGGATGCAGAAGATCATCGTAAGCGAGTTGAACAACTGAt                               >  1:1683283/1‑57 (MQ=255)
                                    gAACTCTTTTGGGCGCTGGGATGCAGAAGATCATCGTAAGCGAGTTGAACAACTGAt                               >  1:1491112/1‑57 (MQ=255)
                                    gAACTCTTTTGGGCGCTGGGATGCAGAAGATCATCGTAAGCGAGTTGAACAACTGAt                               >  1:3351232/1‑57 (MQ=255)
                                    gAACTCTTTTGGGCGCTGGGATGCAGAAGATCATCGTAAGCGAGTTGAACAACTGAt                               >  1:3355222/1‑57 (MQ=255)
                                    gAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGAACAACTGATTGCCCGCa                       >  1:3107527/1‑65 (MQ=255)
                                    gAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGAACAACTGATTACCCGCa                       >  1:1334856/1‑65 (MQ=255)
                                    gAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGAACAACTGATTACCCGCa                       >  1:714168/1‑65 (MQ=255)
                                    gAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGAACAACTGATTACCCGCAt                      >  1:918250/1‑66 (MQ=255)
                                    gAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGAACAACTGATTACCCGCAt                      >  1:897752/1‑66 (MQ=255)
                                    gAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGAACAACTGATTACCCGCAt                      >  1:3005618/1‑66 (MQ=255)
                                    gAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGAACAACTGATTACCCGCAt                      >  1:2754532/1‑66 (MQ=255)
                                    gAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGAACAACTGATTACCCGCAt                      >  1:2699703/1‑66 (MQ=255)
                                    gAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGAACAACTGATTACCCGCAt                      >  1:2684178/1‑66 (MQ=255)
                                    gAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGAACAACTGATTACCCGCAt                      >  1:2498174/1‑66 (MQ=255)
                                    gAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGAACAACTGATTACCCGCAt                      >  1:2452618/1‑66 (MQ=255)
                                    gAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGAACAACTGATTACCCGCAt                      >  1:2435753/1‑66 (MQ=255)
                                    gAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGAACAACTGATTACCCGCAt                      >  1:2285833/1‑66 (MQ=255)
                                    gAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGAACAACTGATTACCCGCAt                      >  1:1977264/1‑66 (MQ=255)
                                    gAACTCTTTTGGGCGCTGGGACGCAGAAGATCATCGTAAGCGAGTTGAACAACTGATTACCCGCAt                      >  1:1549731/1‑66 (MQ=255)
                                    gAACTCTTTTGGGCGCTGGGACGCAGAAGATAATCGTAAGCGAGTTGAACAACTGATTACCCGCAt                      >  1:3131685/1‑66 (MQ=255)
                                    gAACTCTTTTGGGCGCTGGACGCAGAAGATCATCGTAAGCGAGTTGAACAACTGATTACCCGCATGa                    >  1:1265994/1‑67 (MQ=255)
                                        tctTTTGGGCGCTGGGATGCAGAAGATCATCGTAAGCGAGTTGAACAACTGATTACCCGCATGAAagag               <  1:1056495/69‑1 (MQ=255)
                                           tttGGGCGCTGGGATGCAGAAGATCATCGTAAGCGAGTTGAACAACTGATTACCCGCATGAAagag               <  1:1607882/66‑1 (MQ=255)
                                             tGGGCGCTGGGATGCAGAAGATCATCGTAAGCGAGTTg                                         <  1:1739761/38‑1 (MQ=255)
                                                     gggATGCAGAAGATCATCGTAAGCGAGTTGAACAACTGATTACCCGCATGAAAGAGTACGGCCAGCTGc  >  1:1128168/1‑69 (MQ=255)
                                                         |                                                                
CTCGCTAACCGAAATGGTGGTGTTGAAACCGGCGTTGAACTCTTTTGGGCGCTGGGATGCAGAAGATCATCGTAAGCGAGTTGAACAACTGATTACCCGCATGAAAGAGTACGGCCAGCTGC  >  minE/161899‑162020

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: