Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 223,934 G→A 60.7% G387G (GGC→GGT pepD ← aminoacyl‑histidine dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE223,9340GA60.7% 22.2 / 36.6 56G387G (GGC→GGTpepDaminoacyl‑histidine dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (10/12);  new base A (16/18);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

TGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGAATCCAGCATGCTCACCACGTAGTCTTTACCGCTGTCGAT  >  minE/223869‑223993
                                                                 |                                                           
tGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCgcca                                                          >  1:3127103/1‑69 (MQ=255)
tGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCgcc                                                           >  1:1815211/1‑68 (MQ=255)
         ggAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTa                                                 <  1:2565878/69‑1 (MQ=255)
                 aGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGa                                         <  1:893389/69‑1 (MQ=255)
                 aGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGa                                         >  1:421904/1‑69 (MQ=255)
                 aGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGa                                         >  1:3369392/1‑69 (MQ=255)
                 aGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGa                                         >  1:1953049/1‑69 (MQ=255)
                 aGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGa                                         >  1:1918986/1‑69 (MQ=255)
                 aGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGa                                         >  1:1890195/1‑69 (MQ=255)
                 aGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGa                                         >  1:1833307/1‑69 (MQ=255)
                 aGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGa                                         >  1:1697279/1‑69 (MQ=255)
                 aGCGTCCGGCTGACAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCTGTTTACCCAGCGa                                         >  1:107000/1‑69 (MQ=255)
                       cGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGAATCCAg                                   <  1:2630829/69‑1 (MQ=255)
                       cGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGAATCCAg                                   <  1:3388135/69‑1 (MQ=255)
                           tGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGAATCCAGCATg                               <  1:286331/69‑1 (MQ=255)
                           tGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGAATCCAGCATg                               <  1:1025663/69‑1 (MQ=255)
                           tGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGAATCCAGCATg                               <  1:319156/69‑1 (MQ=255)
                            gccagccaGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGAATCCAGCATg                               <  1:1450951/68‑1 (MQ=255)
                               agccagGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGAATCCAGCATGCTca                           <  1:2182369/69‑1 (MQ=255)
                               agccagGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGAATCCAGCATGCTca                           <  1:2178518/69‑1 (MQ=255)
                                              cTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGAATCCAGCATGCTCACCACGTAGTCTTTAc            >  1:2592032/1‑69 (MQ=255)
                                              cTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGAATCCAGCATGCTCACCACGTAGTCTTTAc            >  1:2558478/1‑69 (MQ=255)
                                                       ctgTTTTCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:2992747/3‑68 (MQ=255)
                                                       ctgTTTTCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:3319319/3‑68 (MQ=255)
                                                       ctgTTTTCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:2866231/3‑68 (MQ=255)
                                                       ctgTTTTCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:3391807/3‑68 (MQ=255)
                                                       ctgTTTTCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:3435200/3‑68 (MQ=255)
                                                       ctgTTTTCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:414526/3‑68 (MQ=255)
                                                       ctgTTTTCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:1572751/3‑68 (MQ=255)
                                                       ctgTTTTCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:2546679/3‑68 (MQ=255)
                                                       ctgTTTTCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:2495377/3‑68 (MQ=255)
                                                       ctgTTTTCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:870780/3‑68 (MQ=255)
                                                       ctgTTTTCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:1292815/3‑68 (MQ=255)
                                                       ctgTTTTCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:193818/3‑68 (MQ=255)
                                                       ctgTTTTCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:1310975/3‑68 (MQ=255)
                                                       ctgTTTTCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:183462/3‑68 (MQ=255)
                                                       ctgTTTTCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:1322440/3‑68 (MQ=255)
                                                       ctgTTTTCGCACCTGACAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:2328473/3‑68 (MQ=255)
                                                         gTTTTCGCGCCAGCCAGTTTACCCAGCGAATCCAGCATGCTCACCACGTAGTCTTTACCGCTGTCGAt  <  1:762993/68‑1 (MQ=255)
                                                         gTTTTCGCGCCAGCCAGTTTACCCAGCGAATCCAGCATGCTCACCACGTAGTCTTTACCGCTGTCGAt  <  1:603652/68‑1 (MQ=255)
                                                          ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:1004121/40‑1 (MQ=37)
                                                          ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:728164/40‑1 (MQ=37)
                                                          ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:1008461/40‑1 (MQ=37)
                                                          ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:1595574/40‑1 (MQ=37)
                                                          ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:875990/40‑1 (MQ=37)
                                                          ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:423930/40‑1 (MQ=37)
                                                          ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:501357/40‑1 (MQ=37)
                                                          ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:635133/40‑1 (MQ=37)
                                                          ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:389418/40‑1 (MQ=37)
                                                          ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:331761/40‑1 (MQ=37)
                                                          ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:3163756/40‑1 (MQ=37)
                                                          ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:2912181/40‑1 (MQ=37)
                                                          ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:2904239/40‑1 (MQ=37)
                                                          ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:1251312/40‑1 (MQ=37)
                                                          ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:2562617/40‑1 (MQ=37)
                                                          ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:2120147/40‑1 (MQ=37)
                                                          ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:1893660/40‑1 (MQ=37)
                                                          ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:1536991/40‑1 (MQ=37)
                                                                 |                                                           
TGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGAATCCAGCATGCTCACCACGTAGTCTTTACCGCTGTCGAT  >  minE/223869‑223993

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: