Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,394,275 T→C 100% T120T (ACT→ACC yigI → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,394,2750TC76.6% 17.3 / 0.8 17T120T (ACT→ACCyigIconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base C (13/0);  total (14/3)
Fisher's exact test for biased strand distribution p-value = 5.88e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGGGACCATTGATCTTCGCGTTGATTATCTGCGCCCAGGCAGGGGCGAGCGTTTTACTGCTACTAGTAGCCTGTTGCGTGCAGGCAATAAAGTCGCCGTCGCCCGCGTTGAATTACACAATGAAGAACAG  >  minE/2394212‑2394341
                                                               |                                                                  
ggggACCATTGATCTTCGCGTTGATTATCTGCGCCCAGGCAGGGGCGAGCGTTTTACTGCTACtagtag                                                               >  1:839114/1‑69 (MQ=255)
ggggACCATTGATCTTCGCGTTGATTATCTGCGCCCAGGCAGGGGCGAGCGTTTTACTGCTACtagtag                                                               >  1:2596000/1‑69 (MQ=255)
          gATCTTCGCGTTGATTATCTGCGCCCAGGCAGGGGCGAGCGTTTTACTGCTACTAGTAGCCTGTTGCg                                                      >  1:3060145/1‑68 (MQ=255)
            tCTTCGCGTTGATTATCTGCGCCCAGGCAGGGGCGAGCGTTTTACTGCTACTAGTAGCCTGTTGCGTg                                                    <  1:566576/68‑1 (MQ=255)
            tCTTCGCGTTGATTATCTGCGCCCAGGCAGGGGCGAGCGTTTTACTGCTACTAGTAGCCTGTTGCGTg                                                    <  1:2005835/68‑1 (MQ=255)
                                         ggggCGAGCGTTTTACTGCTACTAGTAGCCTGTTGCGTGCAGGCAATAAAGTCGCCGTCGCCCGCGTTg                      <  1:997114/69‑1 (MQ=255)
                                                      tACTGCTACCAGTAGCCTGTTACGTGCAGGCAATAAAGTCGCCGTCGCCCGCGTTGAATTACACaa            >  1:2743944/1‑66 (MQ=255)
                                                      tACTGCTACCAGTAGCCTGTTACGTGCAGGCAATAAAGTCGCCGTCGCCCGCGTTGAATTACACaa            >  1:852917/1‑66 (MQ=255)
                                                      tACTGCTACCAGTAGCCTGTTACGTGCAGGCAATAAAGTCGCCGTCGCCCGCGTTGAATTACACaa            >  1:634851/1‑66 (MQ=255)
                                                      tACTGCTACCAGTAGCCTGTTACGTGCAGGCAATAAAGTCGCCGTCGCCCGCGTTGAATTACACaa            >  1:32581/1‑66 (MQ=255)
                                                      tACTGCTACCAGTAGCCTGTTACGTGCAGGCAATAAAGTCGCCGTCGCCCGCGTTGAATTACACaa            >  1:312267/1‑66 (MQ=255)
                                                      tACTGCTACCAGTAGCCTGTTACGTGCAGGCAATAAAGTCGCCGTCGCCCGCGTTGAATTACACaa            >  1:2972406/1‑66 (MQ=255)
                                                      tACTGCTACCAGTAGCCTGTTACGTGCAGGCAATAAAGTCGCCGTCGCCCGCGTTGAATTACACaa            >  1:2483347/1‑66 (MQ=255)
                                                      tACTGCTACCAGTAGCCTGTTACGTGCAGGCAATAAAGTCGCCGTCGCCCGCGTTGAATTACACaa            >  1:1229607/1‑66 (MQ=255)
                                                      tACTGCTACCAGTAGCCTGTTACGTGCAGGCAATAAAGTCGCCGTCGCCCGCGTTGAATTACACaa            >  1:1233507/1‑66 (MQ=255)
                                                      tACTGCTACCAGTAGCCTGTTACGTGCAGGCAATAAAGTCGCCGTCGCCCGCGTTGAATTACACaa            >  1:2159238/1‑66 (MQ=255)
                                                      tACTGCTACCAGTAGCCTGTTACGTGCAGGCAATAAAGTCGCCGTCGCCCGCGTTGAATTACACaa            >  1:138585/1‑66 (MQ=255)
                                                      tACTGCTACCAGTAGCCTGTTACGTGCAGGCAATAAAGTCGCCGTCGCCCGCGTTGAATTACACaa            >  1:1518260/1‑66 (MQ=255)
                                                      tACTGCTACCAGTAGCCTCTTACGTGCAGGCAATAAAGTCGCCGTCGCCCGCGTTGAATTACACaa            >  1:545158/1‑66 (MQ=255)
                                                             actaGTAGCCTGTTGCGTGCAGGCAATAAAGTCGCGGTCGCCCGCGTTGAATTACACAATGAAGAACa   <  1:2308597/68‑1 (MQ=255)
                                                             actaGTAGCCTGTTGCGTGCAGGCAATAAAGTCGCCGTCGCCCGCGTTGAATTACACAATGAAGAagag  <  1:989419/69‑1 (MQ=255)
                                                             actaGTAGCCTGTTGCGTGCAGGCAATAAAGTCGCCGTCGCCCGCGTTGAATTACACAATGAAGAACa   <  1:118217/68‑1 (MQ=255)
                                                             actaGTAGCCTGTTGCGTGCAGGCAATAAAGTCGCCGTCGCCCGCGTTGAATTACACAATGAAGAACa   <  1:3153127/68‑1 (MQ=255)
                                                             actaGTAGCCTGTTGCGTGCAGGCAATAAAGTCGCCGTCGCCCGCGTTGAATTACACAATGAAGAACa   <  1:1059557/68‑1 (MQ=255)
                                                             actaGTAGCCTGTTGCGTGCAGGCAATAAAGTCGCCGTCGCCCGCGTTGAATTACACAATGAAGAACAg  <  1:1700482/69‑1 (MQ=255)
                                                             actaGTAGCCTGTTGCGTGCAGGCAATAAAGTCGCCGTCGCCCGCGTTGAATTACACAATGAAGAACAg  <  1:2009711/69‑1 (MQ=255)
                                                             actaGTAGCCTGTTGCGTGCAGGCAATAAAGTCGCCGTCGCCCGCGTTGAATTACACAATGAAGAACAg  <  1:2769746/69‑1 (MQ=255)
                                                             actaGTAGCCTGTTGCGTGCAGGCAATAAAGTCGCCGTCGCCCGCGTTGAATTACACAATGAAGAACAg  <  1:2743705/69‑1 (MQ=255)
                                                             actaGTAGCCTGTTGCGTACAGGCAATAAAGTCGCCGTCGCCCGCGTTGAATTACACAATGAAGAACAg  <  1:1590863/69‑1 (MQ=255)
                                                             actaGTAGCCTGTTGCGTACAGGCAATAAAGTCGCCGTCGCCCGCGTTGAATTACACAATGAAGAACAg  <  1:2849396/69‑1 (MQ=255)
                                                               |                                                                  
GGGGACCATTGATCTTCGCGTTGATTATCTGCGCCCAGGCAGGGGCGAGCGTTTTACTGCTACTAGTAGCCTGTTGCGTGCAGGCAATAAAGTCGCCGTCGCCCGCGTTGAATTACACAATGAAGAACAG  >  minE/2394212‑2394341

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: