Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 1,126,937 (A)7→8 coding (25/45 nt) pheM ← phenylalanyl‑tRNA synthetase operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,126,9301.A81.5% 73.0 / 13.5 27coding (32/45 nt)pheMphenylalanyl‑tRNA synthetase operon leader peptide
Reads supporting (aligned to +/- strand):  ref base . (0/5);  new base A (0/22);  total (0/27)
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

TGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTT  >  minE/1126878‑1126978
                                                     |                                                 
tGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAaaga                                         <  1:1110287/62‑1 (MQ=255)
tGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAaaga                                         <  1:3075406/62‑1 (MQ=255)
tGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAaaga                                         <  1:1870382/62‑1 (MQ=255)
tGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAaaga                                         <  1:2585200/62‑1 (MQ=255)
                                        cAGGTGCTAAAGTAAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt  <  1:1016843/62‑1 (MQ=255)
                                        cAGGTGCTAAAGTAAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt  <  1:849949/62‑1 (MQ=255)
                                        cAGGTGCTAAAGTAAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt  <  1:708890/62‑1 (MQ=255)
                                        cAGGTGCTAAAGTAAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt  <  1:708833/62‑1 (MQ=255)
                                        cAGGTGCTAAAGTAAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt  <  1:345364/62‑1 (MQ=255)
                                        cAGGTGCTAAAGTAAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt  <  1:340752/62‑1 (MQ=255)
                                        cAGGTGCTAAAGTAAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt  <  1:2941501/62‑1 (MQ=255)
                                        cAGGTGCTAAAGTAAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt  <  1:2716523/62‑1 (MQ=255)
                                        cAGGTGCTAAAGTAAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt  <  1:2655350/62‑1 (MQ=255)
                                        cAGGTGCTAAAGTAAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt  <  1:2638797/62‑1 (MQ=255)
                                        cAGGTGCTAAAGTAAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt  <  1:2278889/62‑1 (MQ=255)
                                        cAGGTGCTAAAGTAAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt  <  1:2209770/62‑1 (MQ=255)
                                        cAGGTGCTAAAGTAAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt  <  1:2056938/62‑1 (MQ=255)
                                        cAGGTGCTAAAGTAAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt  <  1:1761522/62‑1 (MQ=255)
                                        cAGGTGCTAAAGTAAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt  <  1:1455564/62‑1 (MQ=255)
                                        cAGGTGCTAAAGTAAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt  <  1:1445247/62‑1 (MQ=255)
                                        cAGGTGCTAAAGTAAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt  <  1:1337218/62‑1 (MQ=255)
                                        cAGGTGCTAAAGTAAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt  <  1:1127367/62‑1 (MQ=255)
                                        cAGGTGCTAAAGTAAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt  <  1:1109850/62‑1 (MQ=255)
                                        cAGGTGCTAAAGTAAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt  <  1:1043227/62‑1 (MQ=255)
                                        cAGGTGCTAAAGTAAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt  <  1:100574/62‑1 (MQ=255)
                                        cAGGTGCTAAAGTAAAAAAAAGAAGCGGAAAATAGCAGCATTAATTGCTTGCGTTACCtttt  <  1:980482/62‑1 (MQ=255)
                                        cAGGTGCTAAAGTAGAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt  <  1:2875011/62‑1 (MQ=255)
                                                     |                                                 
TGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTT  >  minE/1126878‑1126978

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