Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 223,952 C→T 61.9% S381S (TCG→TCA pepD ← aminoacyl‑histidine dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE223,9520CT61.9% 24.3 / 34.1 55S381S (TCG→TCApepDaminoacyl‑histidine dipeptidase
Reads supporting (aligned to +/- strand):  ref base C (10/11);  new base T (16/18);  total (26/29)
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.14e-01

AGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGAATCCAGCATGCTCACCACGTAGTCTTTACCGCTGTCGAT  >  minE/223886‑223993
                                                                  |                                         
aGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGa                                         <  1:893389/69‑1 (MQ=255)
aGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGa                                         >  1:1918986/1‑69 (MQ=255)
aGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGa                                         >  1:1890195/1‑69 (MQ=255)
aGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGa                                         >  1:1953049/1‑69 (MQ=255)
aGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGa                                         >  1:1833307/1‑69 (MQ=255)
aGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGa                                         >  1:1697279/1‑69 (MQ=255)
aGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGa                                         >  1:421904/1‑69 (MQ=255)
aGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGa                                         >  1:3369392/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:319156/69‑1 (MQ=255)
          tGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGAATCCAGCATg                               <  1:286331/69‑1 (MQ=255)
          tGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGAATCCAGCATg                               <  1:1025663/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:2558478/1‑69 (MQ=255)
                             cTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGAATCCAGCATGCTCACCACGTAGTCTTTAc            >  1:2592032/1‑69 (MQ=255)
                                      ctgTTTTCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:870780/3‑68 (MQ=255)
                                      ctgTTTTCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:3435200/3‑68 (MQ=255)
                                      ctgTTTTCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:3391807/3‑68 (MQ=255)
                                      ctgTTTTCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:1292815/3‑68 (MQ=255)
                                      ctgTTTTCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:3319319/3‑68 (MQ=255)
                                      ctgTTTTCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:414526/3‑68 (MQ=255)
                                      ctgTTTTCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:1310975/3‑68 (MQ=255)
                                      ctgTTTTCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:1322440/3‑68 (MQ=255)
                                      ctgTTTTCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:2992747/3‑68 (MQ=255)
                                      ctgTTTTCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:2866231/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:183462/3‑68 (MQ=255)
                                      ctgTTTTCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:193818/3‑68 (MQ=255)
                                      ctgTTTTCGCACCTGACAGTTTACCCAGTGAATCCAGCATGCTCACTACGTAGTCTTTACCGCTGTCg    >  1:2328473/3‑68 (MQ=255)
                                        gTTTTCGCGCCAGCCAGTTTACCCAGCGAATCCAGCATGCTCACCACGTAGTCTTTACCGCTGTCGAt  <  1:603652/68‑1 (MQ=255)
                                        gTTTTCGCGCCAGCCAGTTTACCCAGCGAATCCAGCATGCTCACCACGTAGTCTTTACCGCTGTCGAt  <  1:762993/68‑1 (MQ=255)
                                         ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:728164/40‑1 (MQ=37)
                                         ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:501357/40‑1 (MQ=37)
                                         ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:423930/40‑1 (MQ=37)
                                         ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:2562617/40‑1 (MQ=37)
                                         ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:635133/40‑1 (MQ=37)
                                         ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:875990/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:1004121/40‑1 (MQ=37)
                                         ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:2120147/40‑1 (MQ=37)
                                         ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:1893660/40‑1 (MQ=37)
                                         ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:1595574/40‑1 (MQ=37)
                                         ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:1536991/40‑1 (MQ=37)
                                         ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:1251312/40‑1 (MQ=37)
                                         ttttCGCACCTGCCAGTTTACCCAGTGAATCCAGCATGCt                             <  1:1008461/40‑1 (MQ=37)
                                                                  |                                         
AGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCCAGCGAATCCAGCATGCTCACCACGTAGTCTTTACCGCTGTCGAT  >  minE/223886‑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: