Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,502,952 C→A 69.6% intergenic (‑325/‑595) lrhA ← / → yfbQ DNA‑binding transcriptional repressor/predicted aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,502,9520CA69.6% 68.6 / 46.3 69intergenic (‑325/‑595)lrhA/yfbQDNA‑binding transcriptional repressor/predicted aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (5/16);  new base A (15/33);  total (20/49)
Fisher's exact test for biased strand distribution p-value = 5.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATTTATCCGAGCATTTTACCCCAAAAA  >  minE/1502891‑1503019
                                                             |                                                                   
tgtaaatgtGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGc                                                               <  1:105159/68‑1 (MQ=255)
tgtaaatgtGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGc                                                               <  1:2637792/68‑1 (MQ=255)
tgtaaatgtGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGc                                                               <  1:1778542/68‑1 (MQ=255)
    aatgtGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATg                                                           <  1:2185477/68‑1 (MQ=255)
    aatgtGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATg                                                           <  1:1336625/68‑1 (MQ=255)
          attCCTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggc                                                    <  1:3285466/67‑1 (MQ=255)
          aGTCCTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggc                                                    <  1:999851/69‑1 (MQ=255)
          aGTCCTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggc                                                    <  1:2128752/69‑1 (MQ=255)
          aGTCCTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggc                                                    <  1:2297487/69‑1 (MQ=255)
          aGTCCTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggc                                                    <  1:233715/69‑1 (MQ=255)
          aGTCCTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggc                                                    <  1:2571485/69‑1 (MQ=255)
          aGTCCTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggc                                                    <  1:2685445/69‑1 (MQ=255)
          aGTCCTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggc                                                    <  1:1752861/69‑1 (MQ=255)
          aGTCCTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggc                                                    <  1:2824856/69‑1 (MQ=255)
          aGTCCTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggc                                                    <  1:1623513/69‑1 (MQ=255)
          aGTCCTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggc                                                    <  1:773265/69‑1 (MQ=255)
          aGTCCTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggc                                                    <  1:1121190/69‑1 (MQ=255)
          aGTCCTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggc                                                    <  1:1530624/69‑1 (MQ=255)
          aGTCCTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggc                                                    <  1:2895999/69‑1 (MQ=255)
          aGTCCTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggc                                                    <  1:1124778/69‑1 (MQ=255)
          aGTCCTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggc                                                    <  1:1496837/69‑1 (MQ=255)
          aGTCCTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggc                                                    <  1:690439/69‑1 (MQ=255)
          aGTCCTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggc                                                    <  1:1359896/69‑1 (MQ=255)
          aGTCCTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggc                                                    <  1:508096/69‑1 (MQ=255)
           gTCCTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggc                                                    <  1:2842402/68‑1 (MQ=255)
            tCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAggcgg                                                  >  1:2881782/1‑69 (MQ=255)
              cTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAg                                                                >  1:1674584/1‑53 (MQ=255)
              cTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggcggc                                                 >  1:2607486/1‑68 (MQ=255)
              cTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggcggc                                                 >  1:2609234/1‑68 (MQ=255)
              cTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggcggc                                                 >  1:2437623/1‑68 (MQ=255)
              cTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggcggc                                                 >  1:291894/1‑68 (MQ=255)
              cTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggcggc                                                 >  1:2710972/1‑68 (MQ=255)
              cTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggcggc                                                 >  1:2856106/1‑68 (MQ=255)
              cTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggcggc                                                 >  1:2243848/1‑68 (MQ=255)
              cTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggcggc                                                 >  1:1118410/1‑68 (MQ=255)
              cTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggcggc                                                 >  1:1147388/1‑68 (MQ=255)
              cTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggcggc                                                 >  1:2062396/1‑68 (MQ=255)
              cTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggcggc                                                 >  1:2218484/1‑68 (MQ=255)
              cTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggcggc                                                 >  1:1582428/1‑68 (MQ=255)
              cTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggcggc                                                 >  1:1873550/1‑68 (MQ=255)
              cTTGTTCCACTCTCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAggcggc                                                 >  1:1718598/1‑68 (MQ=255)
                       ctctCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAGGCGGCATCAAGAGc                                        <  1:2634398/68‑1 (MQ=255)
                       ctctCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAGGCGGCATCAAGAGc                                        <  1:2959506/68‑1 (MQ=255)
                       ctctCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAGGCGGCATCAAGAGc                                        <  1:1475228/68‑1 (MQ=255)
                       ctctCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAGGCGGCATCAAGAGc                                        <  1:1621562/68‑1 (MQ=255)
                       ctctCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAGGCGGCATCAAGAGc                                        <  1:321172/68‑1 (MQ=255)
                       ctctCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAGGCGGCATCAAGAGc                                        <  1:1804243/68‑1 (MQ=255)
                       ctctCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAGGCGGCATCAAGAGc                                        <  1:1853669/68‑1 (MQ=255)
                       ctctCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAGGCGGCATCAAGAGc                                        <  1:134154/68‑1 (MQ=255)
                       ctctCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAGGCGGCATCAAGAGc                                        <  1:600284/68‑1 (MQ=255)
                       ctctCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAGGCGGCATCAAGAGc                                        <  1:608009/68‑1 (MQ=255)
                       ctctCGTACAGCATCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAGGCGGCATCAAGAGc                                        <  1:2259692/68‑1 (MQ=255)
                            gTGCAGCATCGCTGGTCATACGCGAACACGTACCAACCGCAATg                                                           <  1:441313/44‑1 (MQ=255)
                              gcagcaTCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGa                                   <  1:1530161/66‑1 (MQ=255)
                              gcagcaTCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGa                                   <  1:1090424/66‑1 (MQ=255)
                              gcagcaTCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGa                                   <  1:16582/66‑1 (MQ=255)
                                    tCGCTGGTCATACGCGAACACGTACAAACAGCAATGGTGAGGCGGCATCAAGAGc                                        <  1:2165429/55‑1 (MQ=255)
                                         ggTCATACGCGAACACGTACAAACAGCAATGGTGAggc                                                    <  1:1980600/38‑1 (MQ=255)
                                               aCGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATTTATCCGAGCa                <  1:650014/68‑1 (MQ=255)
                                                cgcgAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATcc                                >  1:2033820/1‑51 (MQ=255)
                                                cgcgAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATTTATCCGAGCAt               >  1:3051981/1‑68 (MQ=255)
                                                cgcgAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATTTATCCGAGCAt               >  1:3124887/1‑68 (MQ=255)
                                                cgcgAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATTTATCCGAGCAt               >  1:3374537/1‑68 (MQ=255)
                                                    aaCACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATTTATCCGAGCATTTTAc          <  1:2305548/69‑1 (MQ=255)
                                                    aaCACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATTTATCCGAGCATTTCAc          <  1:1464508/69‑1 (MQ=255)
                                                      cacGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATTTATCCGAGCATTTTAc          <  1:1951858/67‑1 (MQ=255)
                                                           aCCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATTTATCCGAGCATTTTACCCCaaaa   <  1:210314/69‑1 (MQ=255)
                                                           aCCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATTTATCCGAGCATTTTACCCCaaaa   <  1:932603/69‑1 (MQ=255)
                                                           aCCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATTTATCCGAGCATTTTACCACaaaa   <  1:1756776/69‑1 (MQ=255)
                                                            ccAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATTTATCCGAGCATTTTACCCCaaaaa  >  1:3112837/1‑69 (MQ=255)
                                                             |                                                                   
TGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATTTATCCGAGCATTTTACCCCAAAAA  >  minE/1502891‑1503019

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: