Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,126,684 T→C 61.7% intergenic (‑51/+233) pheS ← / ← pheM phenylalanine tRNA synthetase, alpha subunit/phenylalanyl‑tRNA synthetase operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,126,6840TC61.7% 32.5 / 40.9 66intergenic (‑51/+233)pheS/pheMphenylalanine tRNA synthetase, alpha subunit/phenylalanyl‑tRNA synthetase operon leader peptide
Reads supporting (aligned to +/- strand):  ref base T (13/12);  new base C (20/20);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

CTGCGAGATGTGACATGGTTTTCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCTATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGCTGTGCTCGAAATTCTCACGTACTTAAA  >  minE/1126618‑1126751
                                                                  |                                                                   
cTGCGAGATGTGACATGGTTTTCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCTAt                                                                   >  1:94148/1‑69 (MQ=255)
cTGCGAGATGTGACATGGTTTTCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCTAt                                                                   >  1:1876717/1‑69 (MQ=255)
cTGCGAGATGTGACATGGTTTTCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCTAt                                                                   >  1:634862/1‑69 (MQ=255)
cTGCGAGATGTGACATGGTTTTCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCTAt                                                                   >  1:730374/1‑69 (MQ=255)
         gtgACATGTTTTTCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCCATCCCa                                                               >  1:569376/1‑64 (MQ=255)
         gtgACATGTTTTTCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCCATCCCATCa                                                            >  1:729704/1‑67 (MQ=255)
         gtgACATGTTTTTCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCCATCCCATCAg                                                           >  1:2604364/1‑68 (MQ=255)
         gtgACATGTTTTTCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCCATCCCATCAg                                                           >  1:2649392/1‑68 (MQ=255)
         gtgACATGTTTTTCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCCATCCCATCAg                                                           >  1:3301393/1‑68 (MQ=255)
         gtgACATGTTTTTCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCCATCCCATCAg                                                           >  1:2389929/1‑68 (MQ=255)
         gtgACATGTTTTTCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCCATCCCATCAg                                                           >  1:33084/1‑68 (MQ=255)
         gtgACATGTTTTTCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCCATCCCATCAg                                                           >  1:3338686/1‑68 (MQ=255)
         gtgACATGTTTTTCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCCATCCCATCAg                                                           >  1:2446376/1‑68 (MQ=255)
         gtgACATGTTTTTCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCCATCCCATCAg                                                           >  1:2068314/1‑68 (MQ=255)
         gtgACATGTTTTTCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCCATCCCATCAg                                                           >  1:1949877/1‑68 (MQ=255)
         gtgACATGTTTTTCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCCATCCCATCAg                                                           >  1:1917444/1‑68 (MQ=255)
         gtgACATGTTTTTCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCCATCCCATCAg                                                           >  1:246411/1‑68 (MQ=255)
         gtgACATGTTTTTCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCCATCCCATCAg                                                           >  1:1653566/1‑68 (MQ=255)
         gtgACATGTTTTTCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCCATCCCATCAg                                                           >  1:659385/1‑68 (MQ=255)
         gtgACATGTTTTTCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCCATCCCATCAg                                                           >  1:1436878/1‑68 (MQ=255)
         gtgACATGTTTTTCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCCATCCCATCAg                                                           >  1:3401095/1‑68 (MQ=255)
         gtgACATGTTTTTCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCCATCCCATCAg                                                           >  1:778476/1‑68 (MQ=255)
         gtgACATGTTTTTCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCCATCCCATCAg                                                           >  1:901405/1‑68 (MQ=255)
         gtgACATGTTTTTCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCAATCCCATCa                                                            >  1:1680570/1‑67 (MQ=255)
         gtgACATGTTTTTACTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCCATCCCATCAg                                                           >  1:3169225/1‑68 (MQ=255)
                    ttCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCTATCCCATCAGGCTATTTTACtt                                               <  1:56929/69‑1 (MQ=255)
                    ttCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCTATCCCATCAGGCTATTTTACtt                                               <  1:3297621/69‑1 (MQ=255)
                             gtgtCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCTATCCCATCAGGCTATTTTACTTGCCATTTTg                                      >  1:1841144/1‑69 (MQ=255)
                             gtgtCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCTATCCCATCAGGCTATTTTACTTGCCATTTTg                                      >  1:1749946/1‑69 (MQ=255)
                             gtgtCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCTATCCCATCAGGCTATTTTACTTGCCATTTTg                                      >  1:2527057/1‑69 (MQ=255)
                             gtgtCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCTATCCCATCAGGCTATTTTACTTGCCATTTTg                                      >  1:1273722/1‑69 (MQ=255)
                                  aGTGGTGACACTGGTTCGTTGGACTTAGAGCCTATCCCATCAGGCTATTTTACtt                                               <  1:3394714/55‑1 (MQ=255)
                                              ggTTCGTTGGACTTAGAGCCTATCCCATCAGGCTATTTTACTTGCCa                                           >  1:1349301/1‑47 (MQ=255)
                                              ggTTCGTTGGACTTAGAGCCTATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGCTGTGCTc                      >  1:828802/1‑68 (MQ=255)
                                              ggTTCGTTGGACTTAGAGCCTATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGCTGTGCTc                      >  1:2522544/1‑68 (MQ=255)
                                              ggTTCGTTGGACTTAGAGCCTATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGCTGTGCTc                      >  1:2230837/1‑68 (MQ=255)
                                                 tCGTTGGACTTAGAGCCTATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGCTGTGCTCGaa                   <  1:1590797/68‑1 (MQ=255)
                                                 tCGTTGGACTTAGAGCCTATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGCTGTGCTCGaa                   <  1:1159677/68‑1 (MQ=255)
                                                     tGGACTTAGAGCCCATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGc                             <  1:3127276/54‑1 (MQ=255)
                                                     tGGACTTAGAGCCCATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGc                             <  1:1657384/54‑1 (MQ=255)
                                                     tGGACTTAGAGCCCATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGc                             <  1:486680/54‑1 (MQ=255)
                                                     tGGACTTAGAGCCCATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGc                             <  1:1590505/54‑1 (MQ=255)
                                                     tGGACTTAGAGCCCATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGc                             <  1:365739/54‑1 (MQ=255)
                                                     tGGACTTAGAGCCCATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGc                             <  1:1785311/54‑1 (MQ=255)
                                                     tGGACTTAGAGCCCATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGc                             <  1:2049410/54‑1 (MQ=255)
                                                     tGGACTTAGAGCCCATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGc                             <  1:3369662/54‑1 (MQ=255)
                                                     tGGACTTAGAGCCCATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGc                             <  1:1465690/54‑1 (MQ=255)
                                                     tGGACTTAGAGCCCATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGc                             <  1:3309025/54‑1 (MQ=255)
                                                     tGGACTTAGAGCCCATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGc                             <  1:2206569/54‑1 (MQ=255)
                                                     tGGACTTAGAGCCCATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGc                             <  1:3067842/54‑1 (MQ=255)
                                                     tGGACTTAGAGCCCATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGc                             <  1:1317729/54‑1 (MQ=255)
                                                     tGGACTTAGAGCCCATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGc                             <  1:1283126/54‑1 (MQ=255)
                                                     tGGACTTAGAGCCCATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGc                             <  1:2981705/54‑1 (MQ=255)
                                                     tGGACTTAGAGCCCATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGc                             <  1:2887428/54‑1 (MQ=255)
                                                     tGGACTTAGAGCCCATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGc                             <  1:2608152/54‑1 (MQ=255)
                                                     tGGACTTAGAGCCCATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGc                             <  1:2556928/54‑1 (MQ=255)
                                                     tGGACTTAGAGCCCATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGc                             <  1:1215570/54‑1 (MQ=255)
                                                      ggACTTAGAGCCTATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGCTGTGCTCGAAATtctc             >  1:891414/1‑69 (MQ=255)
                                                      ggACTTAGAGCCCATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGc                             <  1:214126/53‑1 (MQ=255)
                                                               gCCTATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGCTGTGCTCGAAATTCTCACGTACTTa    <  1:3082689/69‑1 (MQ=255)
                                                               gCCTATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGCTGTGCTCGAAATTCTCACGTACTTa    <  1:1362020/69‑1 (MQ=255)
                                                               gCCTATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGCTGTGCTCGAAATTCTCACGTACTTa    <  1:2557357/69‑1 (MQ=255)
                                                               gCCTATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGCTGTGCTCGAAATTCTCACGTACTTa    <  1:2455183/69‑1 (MQ=255)
                                                                ccTATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGCTGTGCTCGAAATTCTCACGTACTTa    <  1:2337937/68‑1 (MQ=255)
                                                                 cTATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGCTGTGCTCGAAATTCTCACGTACTTaaa  <  1:2515198/69‑1 (MQ=255)
                                                                 cTATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGCTGTGCTCGAAATTCTCACGTACTTaaa  <  1:23273/69‑1 (MQ=255)
                                                                  tATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGCTGTGCTCGa                    <  1:478920/50‑1 (MQ=255)
                                                                  tATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGCTGTGCTCGa                    <  1:2121010/50‑1 (MQ=255)
                                                                  tATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGCTGTGCTCGa                    <  1:1375799/50‑1 (MQ=255)
                                                                  tATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGCTGTGCTCGa                    <  1:3062236/50‑1 (MQ=255)
                                                                  tATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGCTGTGCTCGa                    <  1:3034888/50‑1 (MQ=255)
                                                                  tATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGCTGTGCTCGa                    <  1:896067/50‑1 (MQ=255)
                                                                  tATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGCTGTGCTCGa                    <  1:2452599/50‑1 (MQ=255)
                                                                  tATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGCTGTGCTCGa                    <  1:1214939/50‑1 (MQ=255)
                                                                  tATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGCTGTGCTCGAAATTCTCACGTACTTaaa  <  1:2204293/68‑1 (MQ=255)
                                                                  |                                                                   
CTGCGAGATGTGACATGGTTTTCCTCATTGTGTCAGTGGTGACACTGGTTCGTTGGACTTAGAGCCTATCCCATCAGGCTATTTTACTTGCCATTTTGGTCCCGGGCTGTGCTCGAAATTCTCACGTACTTAAA  >  minE/1126618‑1126751

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: