Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,026,083 G→C 46.2% L263L (CTC→CTG malI ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,026,0830GC46.2% 14.6 / 92.8 80L263L (CTC→CTGmalIDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (4/39);  new base C (15/22);  total (19/61)
Fisher's exact test for biased strand distribution p-value = 1.41e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATG  >  minE/1026015‑1026143
                                                                    |                                                            
gccgccAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCTCCCTGCTTTg                                                              >  1:3303347/1‑69 (MQ=255)
gccgccAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTg                                                              >  1:75721/1‑69 (MQ=255)
gccgccAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTg                                                              <  1:1600038/69‑1 (MQ=255)
          aaaCCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGcaaaccaa                                                    <  1:755598/69‑1 (MQ=255)
          aaaCCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGcaaaccaa                                                    <  1:3213800/69‑1 (MQ=255)
          aaaCCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGcaaaccaa                                                    <  1:198734/69‑1 (MQ=255)
          aaaCCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGcaaaccaa                                                    <  1:44896/69‑1 (MQ=255)
          aaaCCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGcaaaccaa                                                    <  1:3107372/69‑1 (MQ=255)
          aaaCCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGcaaaccaa                                                    <  1:789576/69‑1 (MQ=255)
          aaaCCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGcaaaccaa                                                    <  1:970617/69‑1 (MQ=255)
          aaaCCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGcaaaccaa                                                    <  1:940835/69‑1 (MQ=255)
           aaCCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGcaaaccaa                                                    <  1:1598474/68‑1 (MQ=255)
           aaCCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCCCCTGCTTTGAGcaaaccaaac                                                  <  1:1375301/69‑1 (MQ=255)
           aaCCTGTTGCGCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGcaaaccaa                                                    <  1:96283/68‑1 (MQ=255)
            aCCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGcaaaccaaac                                                  <  1:1699981/69‑1 (MQ=255)
            aCCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGcaaaccaaac                                                  <  1:1762204/69‑1 (MQ=255)
            aCCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGcaaaccaaac                                                  <  1:2909621/69‑1 (MQ=255)
                          gTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGcccc                                          >  1:1309378/1‑63 (MQ=255)
                          gTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATca                                     >  1:1527524/1‑67 (MQ=255)
                          gTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGc                                    >  1:3159497/1‑69 (MQ=255)
                          gTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGc                                    <  1:2381270/69‑1 (MQ=255)
                          gTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGc                                    >  1:2893652/1‑69 (MQ=255)
                               gATCGACTCCGCTTTCCCCGCTTTGCCTTCCTGCTTTCAGTAAACCAAACCATGCCCCCATCGCAATa                                <  1:392223/68‑1 (MQ=255)
                               gATCGACTCCGCTTTCCCCGCTTTGCCTTCCTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATa                                <  1:1073414/68‑1 (MQ=255)
                               gATCGACTCCGCTTTCCCCGCTTTGCCTTCCTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATa                                <  1:2805462/68‑1 (MQ=255)
                               gATCGACTCCGCTTTCCCCGCTTTGCCTTCCTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATa                                <  1:2907401/68‑1 (MQ=255)
                               gATCGACTCCGCTTTCCCCGCTTTGCCTTCCTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATa                                <  1:3157857/68‑1 (MQ=255)
                               gATCGACTCCGCTTTCCCCGCTTTGCCTTCCTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATa                                <  1:3161975/68‑1 (MQ=255)
                               gATCGACTCCGCTTTCCCCGCTTTGCCTTCCTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATa                                <  1:3221807/68‑1 (MQ=255)
                               gATCGACTCCGCTTTCCCCGCTTTGCCTTCCTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATa                                <  1:1993022/68‑1 (MQ=255)
                               gATCGACTCCGCTTTCCCCGCTTTGCCTTCCTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATa                                <  1:1823538/68‑1 (MQ=255)
                               gATCGACTCCGCTTTCCCCGCTTTGCCTTCCTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATa                                <  1:548890/68‑1 (MQ=255)
                               gATCGACTCCGCTTTCCCCGCTTTGCCTTCCTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATa                                <  1:562424/68‑1 (MQ=255)
                               gATCGACTCCGCTTTCCCCGCTTTGCCTTCCTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATa                                <  1:78993/68‑1 (MQ=255)
                               gATCGACTCCGCTTTCCCCGCTTTGCCTTCCTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATa                                <  1:1255186/68‑1 (MQ=255)
                               gATCGACTCCGCTTTCCCCGCTTTGCCTTCCTGCTTTCACCAAACCAAACCATGCCCCCATCGCAATa                                <  1:2108215/68‑1 (MQ=255)
                               gATCGACTCCGCTTTCCCCGCTTTCCCTTCCTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATa                                <  1:3375850/68‑1 (MQ=255)
                                aTCGACTCCGCTTTCCCCGCTTTGCCTTCCTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATa                                <  1:226383/67‑1 (MQ=255)
                                aTCGACTCCGCTTTCCCCGCTTTGCCTTCCTGCGTTCAGCAAACCAAACCATGCCCCCATCGCAATa                                <  1:2254521/67‑1 (MQ=255)
                                     cTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGt                         <  1:385688/69‑1 (MQ=255)
                                     cTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGt                         <  1:3186431/69‑1 (MQ=255)
                                     cTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGt                         <  1:400363/69‑1 (MQ=255)
                                     cTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGt                         <  1:996935/69‑1 (MQ=255)
                                     cTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGt                         <  1:2836950/69‑1 (MQ=255)
                                     cTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGt                         <  1:1161603/69‑1 (MQ=255)
                                     cTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGt                         <  1:2639506/69‑1 (MQ=255)
                                     cTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGt                         <  1:2268787/69‑1 (MQ=255)
                                     cTCCGCTTTACCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGt                         <  1:963835/69‑1 (MQ=255)
                                     cTCCGCTTACCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGt                         <  1:2183461/69‑1 (MQ=255)
                                            ttCCCCGCTTTGCCTTCCTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCa                  <  1:1068112/69‑1 (MQ=255)
                                            ttCCCCGCTTTGCCTTCCTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCa                  <  1:979939/69‑1 (MQ=255)
                                            ttCCCCGCTTTGCCTTCCTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCa                  <  1:908169/69‑1 (MQ=255)
                                            ttCCCCGCTTTGCCTTCCTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCa                  <  1:1981830/69‑1 (MQ=255)
                                            ttCCCCGCTTTGCCTTCCTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCa                  <  1:2300000/69‑1 (MQ=255)
                                                ccGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAAc               <  1:1786162/68‑1 (MQ=255)
                                                ccGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAAc               <  1:3385371/68‑1 (MQ=255)
                                                     ttGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGc         <  1:1588512/69‑1 (MQ=255)
                                                     ttGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGc         <  1:463722/69‑1 (MQ=255)
                                                          gccCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTg     <  1:3414051/68‑1 (MQ=255)
                                                            ccTGCTTTCAGCAAACCACACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCGCTGAt   >  1:1893846/1‑68 (MQ=255)
                                                            ccTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAAcc              >  1:1242187/1‑57 (MQ=255)
                                                            ccTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCGCTGAt   >  1:3362006/1‑68 (MQ=255)
                                                            ccTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCGCTGATg  >  1:2611376/1‑69 (MQ=255)
                                                            ccTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCGCTGATg  >  1:351723/1‑69 (MQ=255)
                                                            ccTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCGCTGATg  >  1:1807203/1‑69 (MQ=255)
                                                            ccTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCGCTGATg  >  1:2423142/1‑69 (MQ=255)
                                                            ccTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCGCTGATg  >  1:930032/1‑69 (MQ=255)
                                                            ccTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCGCTGATg  >  1:2658988/1‑69 (MQ=255)
                                                            ccTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCGCTGATg  >  1:729160/1‑69 (MQ=255)
                                                            ccTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCGCTGATg  >  1:2609861/1‑69 (MQ=255)
                                                            ccTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCGCTGATg  >  1:1644499/1‑69 (MQ=255)
                                                            ccTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCGCTGATg  >  1:1365884/1‑69 (MQ=255)
                                                            ccTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCGCTGATg  >  1:2871137/1‑69 (MQ=255)
                                                            ccTGCTTTCAGCAAACCAAACCATGCCCCCATCGCAATAGTTTAGTTATAGCAAACCACGGCGCTGATg  >  1:3093334/1‑69 (MQ=255)
                                                            acTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAAt                                 <  1:1632187/37‑1 (MQ=255)
                                                              tGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAAc               <  1:1246849/54‑1 (MQ=255)
                                                               gCTTTGAGCAAACCAAACCATGCCCCCATCGCAAt                                 <  1:3085584/35‑1 (MQ=255)
                                                               gCTTTGAGCAAACCAAACCATGCCCCCATCGCAAt                                 <  1:2811007/35‑1 (MQ=255)
                                                               gCTTTGAGCAAACCAAACCATGCCCCCATCGCAAt                                 <  1:724183/35‑1 (MQ=255)
                                                               gCTTTGAGCAAACCAAACCATGCCCCCATCGCAAt                                 <  1:744525/35‑1 (MQ=255)
                                                               gCTTTGAGCAAACCAAACCATGCCCCCATCGCAAt                                 <  1:2407023/35‑1 (MQ=255)
                                                               gCTTTGAGCAAACCAAACCATGCCCCCATCGCAAt                                 <  1:1665346/35‑1 (MQ=255)
                                                               gCTTTGAGCAAACCAAACCATGCCCCCATCGCAAt                                 <  1:2419979/35‑1 (MQ=255)
                                                                    |                                                            
GCCGCCAGCGAAACCTGTTGCTCAAAGTAACGATCGACTCCGCTTTCCCCGCTTTGCCGCCCTGCTTTGAGCAAACCAAACCATGCCCCCATCGCAATAGTTTCGTTATAGCAAACCACGGCACTGATG  >  minE/1026015‑1026143

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