Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 223,627 T→A 51.9% intergenic (+47/+10) prfH → / ← pepD predicted peptide chain release factor/aminoacyl‑histidine dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE223,6270TA51.9% ‑2.6 / 37.3 27intergenic (+47/+10)prfH/pepDpredicted peptide chain release factor/aminoacyl‑histidine dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (7/6);  new base A (6/8);  total (13/14)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

TGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACG  >  minE/223561‑223694
                                                                  |                                                                   
tGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATc                                                                  >  1:533248/1‑70 (MQ=255)
      ttATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAaataat                                                           >  1:2515293/1‑71 (MQ=255)
      ttATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAaataat                                                           >  1:2070993/1‑71 (MQ=255)
         tCGAAGGATAATGAAGGAAGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTAc                                                        >  1:3207112/1‑71 (MQ=255)
                      aaGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAAttt                                           >  1:2837535/1‑71 (MQ=255)
                         gATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAAtttctt                                        <  1:1569891/71‑1 (MQ=255)
                            gAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttcttt                                       <  1:599184/69‑1 (MQ=255)
                             aaaCTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcag                                    <  1:880056/71‑1 (MQ=255)
                               aCTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagca                                  <  1:1731403/71‑1 (MQ=255)
                                cTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagcag                                 >  1:170806/1‑71 (MQ=255)
                                   ccGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGt                                <  1:1481138/69‑1 (MQ=255)
                                             aTTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttctt                                        <  1:1763287/51‑1 (MQ=255)
                                              ttAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAgtg                   >  1:3270125/1‑71 (MQ=255)
                                               tAAACGCCATCCGGCAGCATATGAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAgtg                   <  1:3448158/70‑1 (MQ=255)
                                               tAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAAtttcttt                                       >  1:1564573/1‑50 (MQ=255)
                                               tAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagt                           >  1:342379/1‑62 (MQ=255)
                                               tAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGaa                                              >  1:1234902/1‑43 (MQ=255)
                                               tAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttctttc                                      >  1:188391/1‑51 (MQ=255)
                                                 aaCGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTc                 <  1:1196873/70‑1 (MQ=255)
                                                      cATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttctt                                        <  1:1595159/42‑1 (MQ=255)
                                                       aTCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCa               >  1:3441972/1‑66 (MQ=255)
                                                        tCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagtcag                        <  1:2413048/56‑1 (MQ=255)
                                                            gcagcaAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCt     <  1:2855291/71‑1 (MQ=255)
                                                               gcaTATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACg  >  1:3441905/1‑71 (MQ=255)
                                                               gcaAATCAAATAATTACTTCGCCGGAATTTCTTTcagc                                   <  1:135424/38‑1 (MQ=255)
                                                               gcaAATCAAATAATTACTTCGCCGGAATTTCTTTcagc                                   <  1:961532/38‑1 (MQ=255)
                                                                caTATCAAATAATTACTTCGCGGGAATTTCTTTcag                                    <  1:3392643/36‑1 (MQ=255)
                                                                  |                                                                   
TGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACG  >  minE/223561‑223694

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