Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 626,915 G→C 65.4% L276L (CTG→CTC ycaQ → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE626,9150GC65.4% 46.2 / 39.7 52L276L (CTG→CTCycaQconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (7/11);  new base C (12/22);  total (19/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTGAAAAATTGGGCAATCTTTGGCTGCATGATGATTTGCTGCCGCTACTCGAGCGAGCGCTGGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGCCTTTTGATCCTGTTGTCTGGGA  >  minE/626852‑626977
                                                               |                                                              
tGTTGAAAAATTGGGCAATCTTTGGCTGCATGATGATTTGCTGCCGCTACTCGAGCGAGCGCTGGCCg                                                            <  1:748250/68‑1 (MQ=255)
tGTTGAAAAATTGGGCAATCTTTGGCTGCATGATGATTTGCTGCCGCTACTCGAGCGAGCGCTGGCCg                                                            <  1:3193305/68‑1 (MQ=255)
tGTTGAAAAATTGGGCAATCTTTGGCTGCATGATGATTTGCTGCCGCTACTCGAGCGAGCGCTGGCCg                                                            <  1:3082670/68‑1 (MQ=255)
tGTTGAAAAATTGGGCAATCTTTGGCTGCATGATGATTTGCTGCCGCTACTCGAGCGAGCGCTGGCCg                                                            <  1:2583547/68‑1 (MQ=255)
  ttGAAAAATTGGGCAATCTTTGGCTGCATGATGATTTGCTGCCGCTACTCGAGCGAGCGCTGGCCGGaa                                                         >  1:2601493/1‑69 (MQ=255)
  ttGAAAAATTGGGCAATCTTTGGCTGCATGATGATTTGCTGCCGCTACTCGAGCGAGCGCTGGCCGGaa                                                         >  1:2646304/1‑69 (MQ=255)
  ttGAAAAATTGGGCAATCTTTGGCTGCATGATGATTTGCTGCCGCTACTCGAGCGAGCGCTGGCCGGaa                                                         >  1:2657764/1‑69 (MQ=255)
  ttGAAAAATTGGGCAATCTTTGGCTGCATGATGATTTGCTGCCGCTACTCGAGCGAGCGCTGGCCGGaa                                                         >  1:1468522/1‑69 (MQ=255)
  ttGAAAAATTGGGCAATCTTTGGCTGCATGATGATTTGCTGCCGCTACTCGAGCGAGCGCTGGCCGGaa                                                         >  1:2809236/1‑69 (MQ=255)
  ttGAAAAATTGGGCAATCTTTGGCTGCATGATGATTTGCTGCCGCTACTCGAGCGAGCGCTGGCCGGa                                                          >  1:27717/1‑68 (MQ=255)
       aaaTTGGGCAATCTTTGGCTGCATGATGATTTGCTGCCGCTACTCGAGCGAGCGCTGGCCGGAAAGCt                                                     <  1:2961995/68‑1 (MQ=255)
       aaaTTGGGCAATCTTTGGCTGCATGATGATTTGCTGCCGCTACTCGAGCGAGCGCTGGCCGGAAAGCt                                                     <  1:1470504/68‑1 (MQ=255)
       aaaTTGGGCAATCTTTGGCTGCATGATGATTTGCTGCCGCTACTCGAGCGAGCGCTGGCCGGAAAGCt                                                     <  1:3237825/68‑1 (MQ=255)
                 aTCTTTGGCTGCATGATGATTTGCTGCCGCTACTCGAGCGAGCGCTGGCCGGAAAGCt                                                     <  1:1395159/58‑1 (MQ=255)
                   cTTTGGCTGCATGATGATTTGCTGCCGCTACTCGAGCGAGCGCTGGCCGGAAAGCTCACTGCAACGcac                                        >  1:2801924/1‑69 (MQ=255)
                                  gatTTGCTGCCGCTACTCGAGCGAGCGCTGGCCGGAAATCTCACTGCAACGCACAGCGCGGTACTTTCg                         <  1:2290381/69‑1 (MQ=255)
                                  gatTTGCTGCCGCTACTCGAGCGAGCGCTGGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCg                         <  1:1886261/69‑1 (MQ=255)
                                  gatTTGCTGCCGCTACTCGAGCGAGCGCTGGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCg                         <  1:1590780/69‑1 (MQ=255)
                                   atTTGCTGCCGCTACTCGAGCGTGCGCTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGc                        >  1:2695084/1‑69 (MQ=255)
                                   atTTGCTGCCGCTACTCGAGCGTGCGCTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGc                        >  1:699224/1‑69 (MQ=255)
                                   atTTGCTGCCGCTACTCGAGCGTGCGCTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGc                        >  1:1157551/1‑69 (MQ=255)
                                   atTTGCTGCCGCTACTCGAGCGTGCGCTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGc                        >  1:2830829/1‑69 (MQ=255)
                                   atTTGCTGCCGCTACTCGAGCGTGCGCTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGc                        >  1:2168998/1‑69 (MQ=255)
                                   atTTGCTGCCGCTACTCGAGCGTGCGCTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGc                        >  1:215915/1‑69 (MQ=255)
                                   atTTGCTGCCGCTACTCGAGCGTGCGCTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGc                        >  1:2151096/1‑69 (MQ=255)
                                   atTTGCTGCCGCTACTCGAGCGTGCGCTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGc                        >  1:75898/1‑69 (MQ=255)
                                   atTTGCTGCCGCTACTCGAGCGTGCGCTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGc                        >  1:766073/1‑69 (MQ=255)
                                   atTTGCTGCCGCTACTCGAGCGTGCGCTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGc                        >  1:828319/1‑69 (MQ=255)
                                   atTTGCTGCCGCTACTCGAGCGTGCGCTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGc                        >  1:1511110/1‑69 (MQ=255)
                                   atTTGCTGCCGCTACTCGAGCGTGCGCTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGc                        >  1:916080/1‑69 (MQ=255)
                                                         tgcgcTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGCCTTTTGATCCTGTTGTCTGGGa  <  1:3321339/68‑1 (MQ=255)
                                                         tgcgcTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGCCTTTTGATCCTGTTGTCTGGGa  <  1:277502/68‑1 (MQ=255)
                                                         tgcgcTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGCCTTTTGATCCTGTTGTCTGGGa  <  1:3348155/68‑1 (MQ=255)
                                                         tgcgcTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGCCTTTTGATCCTGTTGTCTGGGa  <  1:336794/68‑1 (MQ=255)
                                                         tgcgcTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGCCTTTTGATCCTGTTGTCTGGGa  <  1:452399/68‑1 (MQ=255)
                                                         tgcgcTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGCCTTTTGATCCTGTTGTCTGGGa  <  1:533138/68‑1 (MQ=255)
                                                         tgcgcTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGCCTTTTGATCCTGTTGTCTGGGa  <  1:593186/68‑1 (MQ=255)
                                                         tgcgcTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGCCTTTTGATCCTGTTGTCTGGGa  <  1:878160/68‑1 (MQ=255)
                                                         tgcgcTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGCCTTTTGATCCTGTTGTCTGGGa  <  1:3298202/68‑1 (MQ=255)
                                                         tgcgcTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGCCTTTTGATCCTGTTGTCTGGGa  <  1:2918049/68‑1 (MQ=255)
                                                         tgcgcTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGCCTTTTGATCCTGTTGTCTGGGa  <  1:2789082/68‑1 (MQ=255)
                                                         tgcgcTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGCCTTTTGATCCTGTTGTCTGGGa  <  1:2785209/68‑1 (MQ=255)
                                                         tgcgcTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGCCTTTTGATCCTGTTGTCTGGGa  <  1:2407972/68‑1 (MQ=255)
                                                         tgcgcTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGCCTTTTGATCCTGTTGTCTGGGa  <  1:2278347/68‑1 (MQ=255)
                                                         tgcgcTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGCCTTTTGATCCTGTTGTCTGGGa  <  1:2036245/68‑1 (MQ=255)
                                                         tgcgcTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGCCTTTTGATCCTGTTGTCTGGGa  <  1:1818519/68‑1 (MQ=255)
                                                         tgcgcTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGCCTTTTGATCCTGTTGTCTGGGa  <  1:1801980/68‑1 (MQ=255)
                                                         tgcgcTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGCCTTTTGATCCTGTTGTCTGGGa  <  1:1776138/68‑1 (MQ=255)
                                                         tgcgcTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGCCTTTTGATCCTGTTGTCTGGGa  <  1:1375448/68‑1 (MQ=255)
                                                         tgcgcTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGCCTTTTGATCCTGTTGTCTGGGa  <  1:1176303/68‑1 (MQ=255)
                                                          gcgcTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGCCTTTTGATCCTGTTGTCTGGGa  <  1:518048/68‑1 (MQ=255)
                                                          gcgcTCGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGCCTTTTGATCCTGTTGTCTGGGa  <  1:219263/68‑1 (MQ=255)
                                                               |                                                              
TGTTGAAAAATTGGGCAATCTTTGGCTGCATGATGATTTGCTGCCGCTACTCGAGCGAGCGCTGGCCGGAAAGCTCACTGCAACGCACAGCGCGGTACTTTCGCCTTTTGATCCTGTTGTCTGGGA  >  minE/626852‑626977

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: