Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 223,627 T→A 29.4% 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,6270TA29.4% 85.1 / 42.1 68intergenic (+47/+10)prfH/pepDpredicted peptide chain release factor/aminoacyl‑histidine dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (14/34);  new base A (10/10);  total (24/44)
Fisher's exact test for biased strand distribution p-value = 1.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATG  >  minE/223560‑223687
                                                                   |                                                            
aTGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCata                                                             <  1:2698317/69‑1 (MQ=255)
         aTCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAaataat                                                    >  1:2656521/1‑69 (MQ=255)
         aTCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAaataat                                                    >  1:1125592/1‑69 (MQ=255)
         aTCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAaataat                                                    >  1:2053042/1‑69 (MQ=255)
         aTCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAaataat                                                    >  1:2030019/1‑69 (MQ=255)
         aTCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAaataat                                                    >  1:2026130/1‑69 (MQ=255)
         aTCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAaataat                                                    >  1:1611122/1‑69 (MQ=255)
         aTCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAaataa                                                     >  1:808275/1‑68 (MQ=255)
         aTCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAaataa                                                     >  1:1323228/1‑68 (MQ=255)
         aTCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAaataa                                                     >  1:3263737/1‑68 (MQ=255)
         aTCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAaataa                                                     >  1:1979153/1‑68 (MQ=255)
            gAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCCAATAATTa                                                  <  1:1038117/68‑1 (MQ=39)
            gAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTa                                                  <  1:878786/68‑1 (MQ=255)
            gAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTa                                                  <  1:2636626/68‑1 (MQ=255)
            gAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTa                                                  <  1:2477807/68‑1 (MQ=255)
            gAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTa                                                  <  1:1430513/68‑1 (MQ=255)
            gAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTa                                                  <  1:1359147/68‑1 (MQ=255)
            gAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTa                                                  <  1:515245/68‑1 (MQ=255)
            gAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTa                                                  <  1:1293488/68‑1 (MQ=255)
            gAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGAATAAACGCCATCCGGCAGCTAATCAAATAATTa                                                  <  1:3414537/68‑1 (MQ=255)
             aaGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTa                                                  <  1:1603131/67‑1 (MQ=255)
                  tAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTAc                                                 >  1:368901/1‑63 (MQ=255)
                  tAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTAc                                                 >  1:2990838/1‑63 (MQ=255)
                  tAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTAc                                                 >  1:3231918/1‑63 (MQ=255)
                  tAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTAc                                                 >  1:825590/1‑63 (MQ=255)
                  tAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTAc                                                 >  1:345110/1‑63 (MQ=255)
                       aaGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGaa                                       >  1:2884471/1‑68 (MQ=255)
                           aTGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCata                                                             <  1:1025235/42‑1 (MQ=255)
                                aCTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTcag                             <  1:3307115/69‑1 (MQ=255)
                                aCTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTcag                             <  1:2337587/69‑1 (MQ=255)
                                   gCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTAc                                                 <  1:1932011/46‑1 (MQ=255)
                                   gCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTAc                                                 <  1:832121/46‑1 (MQ=255)
                                   gCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTcagcag                          <  1:992779/69‑1 (MQ=255)
                                   gCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTcagcag                          <  1:851320/69‑1 (MQ=255)
                                   gCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTcagcag                          <  1:2740682/69‑1 (MQ=255)
                                   gCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTcagcag                          <  1:2796/69‑1 (MQ=255)
                                   gCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTcagcag                          <  1:3241507/69‑1 (MQ=255)
                                   gCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTcagcag                          <  1:2568015/69‑1 (MQ=255)
                                   gCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTcagcag                          <  1:233309/69‑1 (MQ=255)
                                   gCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTcagcag                          <  1:2207856/69‑1 (MQ=255)
                                   gCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTcagc                            >  1:260282/1‑67 (MQ=255)
                                   gCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTcagc                            >  1:1017147/1‑67 (MQ=255)
                                   gCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTcagc                            >  1:562458/1‑67 (MQ=255)
                                   gCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTcagc                            >  1:428721/1‑67 (MQ=255)
                                   gCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTcagc                            >  1:1648284/1‑67 (MQ=255)
                                   gCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTcagc                            >  1:1945381/1‑67 (MQ=255)
                                   gCCGGAAGGCGATTAAACGCCATCCGGCAGAATATCAAATAATTACTTCGCCGGAATTTCTTTcagc                            >  1:876988/1‑67 (MQ=255)
                                       gAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtca                      <  1:2715266/69‑1 (MQ=255)
                                       gAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtca                      <  1:955304/69‑1 (MQ=255)
                                       gAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtca                      <  1:953898/69‑1 (MQ=255)
                                       gAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtca                      <  1:1023160/69‑1 (MQ=255)
                                       gAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtca                      <  1:2011449/69‑1 (MQ=255)
                                       gAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtca                      <  1:2258195/69‑1 (MQ=255)
                                       gAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtca                      <  1:3075061/69‑1 (MQ=255)
                                       gAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtca                      <  1:3104816/69‑1 (MQ=255)
                                                tAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAgt             >  1:1008141/1‑69 (MQ=255)
                                                          ccGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAAt   <  1:2002442/69‑1 (MQ=255)
                                                          ccGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAAt   <  1:3004847/69‑1 (MQ=255)
                                                          ccGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAAt   <  1:1492095/69‑1 (MQ=255)
                                                           cGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGt      <  1:1336155/65‑1 (MQ=255)
                                                           cGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGt      <  1:3110107/65‑1 (MQ=255)
                                                           cGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGt      <  1:1287668/65‑1 (MQ=255)
                                                           cGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGt      <  1:480687/65‑1 (MQ=255)
                                                           cGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGt      <  1:1011135/65‑1 (MQ=255)
                                                           cGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGt      <  1:2940267/65‑1 (MQ=255)
                                                           cGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGt      <  1:217442/65‑1 (MQ=255)
                                                           cGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGt      <  1:962682/65‑1 (MQ=255)
                                                           cGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGt      <  1:206034/65‑1 (MQ=255)
                                                           cGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATg  <  1:1285743/69‑1 (MQ=255)
                                                           cGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATg  <  1:3180317/69‑1 (MQ=255)
                                                                   |                                                            
ATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATG  >  minE/223560‑223687

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: