Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,431,447 G→C 49.3% A250A (GCC→GCG rho ← transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,431,4470GC49.3% ‑1.5 / 88.5 71A250A (GCC→GCGrhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base G (26/10);  new base C (13/22);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 4.20e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-04

CGCGCGCCAGACGAGTGATGGAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTGGCCTTCTCGATCACCATTTCCGCAACCTGAACGTGGCGAGATGCGGGTTCGTCAAAGGTAGAAGCA  >  minE/2431380‑2431513
                                                                   |                                                                  
cgcgcgcCAGACGAGTGATGGAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTg                                                                    <  1:2285803/68‑1 (MQ=255)
cgcgcgcCAGACGAGTGATGGAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTg                                                                    <  1:819680/68‑1 (MQ=255)
cgcgcgcCAGACGAGTGATGGAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTg                                                                    <  1:111217/68‑1 (MQ=255)
cgcgcgcCAGACGAGTGATGGAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTg                                                                    <  1:2764129/68‑1 (MQ=255)
cgcgcgcCAGACGAGTGATGGAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTg                                                                    <  1:1946368/68‑1 (MQ=255)
cgcgcgcCAGACGAGTGATGGAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTg                                                                    <  1:3281599/68‑1 (MQ=255)
cgcgcgcCAGACGAGTGATGGAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTg                                                                    <  1:161956/68‑1 (MQ=255)
cgcgcgcCAGACGAGTGATGGAGTCGAGCAGAATGAAAACGTCTTTCTTGTGCTCAACCAGGCGTTTg                                                                    <  1:898733/68‑1 (MQ=255)
 gcgcgcCAGACGAGTGATGGAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTGGc                                                                  <  1:403137/69‑1 (MQ=255)
  cgcgcCAGACGAGTGATGGAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTCGcc                                                                 <  1:983704/69‑1 (MQ=255)
  cgcgcCAGACGAGTGATGGAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTCGcc                                                                 <  1:2468170/69‑1 (MQ=255)
  cgcgcCAGACGAGTGATGGAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTCGcc                                                                 <  1:2350951/69‑1 (MQ=255)
  cgcgcCAGACGAGTGATGGAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTCGcc                                                                 <  1:2702986/69‑1 (MQ=255)
  cgcgcCAGACGAGTGATGGAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTCGcc                                                                 <  1:2734469/69‑1 (MQ=255)
  cgcgcCAGACGAGTGATGGAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTCGcc                                                                 <  1:3012656/69‑1 (MQ=255)
  cgcgcCAGACGAGTGATGGAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTCGcc                                                                 <  1:84072/69‑1 (MQ=255)
  cgcgcCAGACGAGTGATGGAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTCGcc                                                                 <  1:560695/69‑1 (MQ=255)
   gcgcCAGACGAGTGATGGAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTCGcc                                                                 <  1:53454/68‑1 (MQ=255)
   gcgcCAGACGAGTGATGGAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTCGcc                                                                 <  1:640872/68‑1 (MQ=255)
   gcgcCAGACGAGTGATGGAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTCGcc                                                                 <  1:2119389/68‑1 (MQ=255)
   gcgcCAGACGAGTGATGGAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTCGcc                                                                 <  1:3418831/68‑1 (MQ=255)
   gcgcCAGACGAGTGATGGAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTCGcc                                                                 <  1:1768078/68‑1 (MQ=255)
   gcgcCAGACGAGTGATGGAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTCGcc                                                                 <  1:1634499/68‑1 (MQ=255)
    cgcCAGACGAGTGATGGAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTCGcc                                                                 <  1:1281360/67‑1 (MQ=255)
        aGACGAGTGATGGAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTGGCCTTCTCg                                                           <  1:354956/69‑1 (MQ=255)
               tGATGGAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTg                                                                    <  1:1972342/53‑1 (MQ=255)
                   ggAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTCGCCTTCTCGATCACCAttt                                                 <  1:1501961/68‑1 (MQ=255)
                   ggAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTCGCCTTCTCGATCACCAttt                                                 <  1:510919/68‑1 (MQ=255)
                   ggAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTCGCCTTCTCGATCACCAttt                                                 <  1:1916951/68‑1 (MQ=255)
                   ggAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTCGCCTTCTCGATCACCAttt                                                 <  1:1940405/68‑1 (MQ=255)
                   ggAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTCGCCTTCTCGATCACCAttt                                                 <  1:841179/68‑1 (MQ=255)
                   ggAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTCGCCTTCTCGATCACCAttt                                                 <  1:2043450/68‑1 (MQ=255)
                    gAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTCGCCTTCTCGATCACCAttt                                                 <  1:505362/67‑1 (MQ=255)
                           gCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTGGCCTTCTCGATCACCATTTCCGCa                                            <  1:98820/65‑1 (MQ=255)
                           gCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTGGCCTTCTCGATCACCATTTCCGCa                                            <  1:2661326/65‑1 (MQ=255)
                           gCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTGGCCTTCTCGATCACCATTTCCGCa                                            <  1:2799510/65‑1 (MQ=255)
                           gCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTGGCCTTCTCGATCACCATTTCCGCa                                            <  1:2825073/65‑1 (MQ=255)
                           gCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTGGCCTTCTCGATCACCATTTCCGCa                                            <  1:3293255/65‑1 (MQ=255)
                           gCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTGGCCTTCTCGATCACCATTTCCGCa                                            <  1:2151712/65‑1 (MQ=255)
                                 tGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTGGCCt                                                                >  1:212695/1‑39 (MQ=255)
                                 tGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTGGCCTtctc                                                            >  1:2951337/1‑43 (MQ=255)
                                 tGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTGGCCTTCTCg                                                           >  1:1991852/1‑44 (MQ=255)
                                 tGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTGGCCTTCTCg                                                           >  1:1825315/1‑44 (MQ=255)
                                 tGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTGGCCTTCTCg                                                           >  1:461701/1‑44 (MQ=255)
                                 tGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTGGCCTTCTCg                                                           >  1:1490704/1‑44 (MQ=255)
                                 tGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTGGCCTTCTCg                                                           >  1:2736259/1‑44 (MQ=255)
                                 tGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTGGCCTTCTCg                                                           >  1:2732191/1‑44 (MQ=255)
                                 tGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTGGCCTTCTCg                                                           >  1:2188951/1‑44 (MQ=255)
                                 tGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTGGCCTTCTCg                                                           >  1:1097089/1‑44 (MQ=255)
                                 tGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTGGCCTTCTCGATCACCATTTCCGCa                                            >  1:80142/1‑59 (MQ=255)
                                 tGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTGGCCTTCTCGATCACCATTTCCGCa                                            >  1:776902/1‑59 (MQ=255)
                                 tGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTGGCCTTCTCGATCACCATTTCCGCa                                            >  1:1224893/1‑59 (MQ=255)
                                 tGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTGGCCTTCTCGATCACCATTTCCGCAACCTGAACgg                                  <  1:1081593/69‑2 (MQ=255)
                                 tGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTGGCCTTCTCGATCACCATTTCCGCAACCTGAACGt                                  <  1:1502553/69‑1 (MQ=255)
                                         tctttcttGTGCTCAACCAGGCGTTTGGCCTTCTCGATCACCATTTCCGCAACCTGAACGTGGCGAGAt                          >  1:1090296/1‑69 (MQ=255)
                                           tttcttGTGCTCAACCAGGCGTTTGGCCTTCTCGATCACCATTTCCGCAACCTGAACGTGGCGAGATGc                        >  1:1827653/1‑69 (MQ=255)
                                           tttcttGTGCTCAACCAGGCGTTTGGCCTTCTCGATCACCATTTCCGCAACCTGAACGTGGCGAGATGc                        >  1:1339923/1‑69 (MQ=255)
                                               ttGTGCTCAACCAGGCGTTTCGCCTTCTCGATCACCATTTc                                                >  1:1630143/1‑41 (MQ=255)
                                               ttGTGCTCAACCAGGCGTTTCGCCTTCTCGATCACCATTTCCg                                              >  1:2880897/1‑43 (MQ=255)
                                               ttGTGCTCAACCAGGCGTTTCGCCTTCTCGATCACCATTTCCGCa                                            >  1:3296389/1‑45 (MQ=255)
                                               ttGTGCTCAACCAGGCGTTTCGCCTTCTCGATCACCATTTCCGCa                                            >  1:1809698/1‑45 (MQ=255)
                                               ttGTGCTCAACCAGGCGTTTCGCCTTCTCGATCACCATTTCCGCa                                            >  1:2593440/1‑45 (MQ=255)
                                               ttGTGCTCAACCAGGCGTTTCGCCTTCTCGATCACCATTTCCGCa                                            >  1:393965/1‑45 (MQ=255)
                                               ttGTGCTCAACCAGGCGTTTCGCCTTCTCGATCACCATTTCCGCa                                            >  1:419428/1‑45 (MQ=255)
                                               ttGTGCTCAACCAGGCGTTTCGCCTTCTCGATCACCATTTCCGCa                                            >  1:1951261/1‑45 (MQ=255)
                                               ttGTGCTCAACCAGGCGTTTCGCCTTCTCGATCACCATTTCCGCa                                            >  1:1162314/1‑45 (MQ=255)
                                               ttGTGCTCAACCAGGCGTTTCGCCTTCTCGATCACCATTTCCGCa                                            >  1:664222/1‑45 (MQ=255)
                                               ttGTGCTCAACCAGGCGTTTCGCCTTCTCGATCACCATTTCCGCa                                            >  1:1036849/1‑45 (MQ=255)
                                               ttGTGCTCAACCAGGCGTTTCGCCTTCTCGATCACCATTTCCGCa                                            >  1:2770737/1‑45 (MQ=255)
                                               ttGTGCTCAACCAGGCGTTTCGCCTTCTCGATCACCATTTCCGCa                                            >  1:2309911/1‑45 (MQ=255)
                                                                  tGGCCTTCTCGATCACCATTTCCGCAACCTGAACGTGGCGAGATGCGGGTTCGTCAAAGGTAGAAGca  >  1:1753035/1‑68 (MQ=255)
                                                                  tGGCCTTCTCGATCACCATTTCCGCAACCTGAACGTGGCGAGATGCGGGTTCGTCAAAGGTAGAAGca  >  1:896062/1‑68 (MQ=255)
                                                                  tGGCCTTCTCGATCACCATTTCCGCAACCTGAACGTGGCGAGATGCGGGTTCGTCAAAGGTAGAAGca  >  1:919439/1‑68 (MQ=255)
                                                                  tGGCCTTCTCGATCACCATTTCCGCAACCTGAACGTGGCGAGATGCGGGTTCGTCAAAGGTAGAAGca  >  1:3196072/1‑68 (MQ=255)
                                                                  tGGCCTTCTCGATCACCATTTCCGCAACCTGAACGTGGCGAGATGCGGGTTCGTCAAAGGTAGAAGca  >  1:3105279/1‑68 (MQ=255)
                                                                  tGGCCTTCTCGATCACCATTTCCGCAACCTGAACGTGGCGAGATGCGGGTTCGTCAAAGGTAGAAGca  >  1:3014489/1‑68 (MQ=255)
                                                                  tGGCCTTCTCGATCACCATTTCCGCAACCTGAACGTGGCGAGATGCGGGTTCGTCAAAGGTAGAAGca  >  1:1910222/1‑68 (MQ=255)
                                                                  tGGCCTTCTCGATCACCATTTCCGCAACCTGAACGTGGCGAGATGCGGGTTCGTCAAAGGTAGAAGca  >  1:2781633/1‑68 (MQ=255)
                                                                  tGGCCTTCTCGATAACCATTTCCGCAACCTGAACGTGGCGAGATGCGGGTTCGTCAAAGGTAGAAGca  >  1:1092729/1‑68 (MQ=255)
                                                                  tGGACTTCTCGATCACCATTTCCGCAACCTGAACGTGGCGAGATGCGGGTTCGTCAAAGGTAGAAGca  >  1:1028204/1‑68 (MQ=255)
                                                                   |                                                                  
CGCGCGCCAGACGAGTGATGGAGTCGAGCAGAATGATAACGTCTTTCTTGTGCTCAACCAGGCGTTTGGCCTTCTCGATCACCATTTCCGCAACCTGAACGTGGCGAGATGCGGGTTCGTCAAAGGTAGAAGCA  >  minE/2431380‑2431513

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: