Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 223,627 T→A 74.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
*Exported223,6270TA74.4% 64.0 / 23.3 39intergenic (+47/+10)prfH/pepDpredicted peptide chain release factor/aminoacyl‑histidine dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (3/7);  new base A (12/17);  total (15/24)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

ATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGAC  >  Exported/223560‑223689
                                                                   |                                                              
aTGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCatat                                                              >  1:513021/1‑70 (MQ=255)
aTGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCaaat                                                              >  1:1405190/1‑70 (MQ=255)
   gCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCaaa                                                          <  1:660179/71‑1 (MQ=255)
       ttATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAaataat                                                      >  1:4802136/1‑71 (MQ=255)
       ttATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAaataa                                                       <  1:1876211/70‑1 (MQ=255)
        tATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAaataa                                                       >  1:3172538/1‑69 (MQ=255)
        tATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAaataat                                                      <  1:5349290/70‑1 (MQ=255)
          tCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTAc                                                   <  1:1142462/71‑1 (MQ=255)
                  tAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATc                                                             <  1:2138658/53‑1 (MQ=255)
                    aTGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGa                                          <  1:1167490/70‑1 (MQ=255)
                       aaGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAttt                                      <  1:4453716/71‑1 (MQ=255)
                       aaGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAttt                                      >  1:2846230/1‑71 (MQ=255)
                        aGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAttt                                      <  1:525247/70‑1 (MQ=255)
                          gATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAt                                        <  1:435002/66‑1 (MQ=255)
                                    ccGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGt                           <  1:5091499/69‑1 (MQ=255)
                                    ccGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAAtt                                                     <  1:2889224/43‑1 (MQ=255)
                                    ccGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGtt                          <  1:803871/70‑1 (MQ=255)
                                     cGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtca                        <  1:1241154/71‑1 (MQ=255)
                                     cGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtca                        <  1:1785425/71‑1 (MQ=255)
                                     cGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtca                        <  1:5318717/71‑1 (MQ=255)
                                      ggAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGtt                          <  1:1319656/68‑1 (MQ=255)
                                       gAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagt                      >  1:5414332/1‑71 (MQ=255)
                                         tggCGTTTAATCGCCTTCCGGCAGTAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagt                      >  1:3658445/2‑69 (MQ=25)
                                         tggCGTTTAATCGCCTTCCGGCAGTAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcag                       >  1:5426807/2‑68 (MQ=25)
                                         aGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagtca                    <  1:4434603/71‑1 (MQ=255)
                                         aGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGaa                                         >  1:4945021/1‑50 (MQ=255)
                                           gCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTcagc                  <  1:2478385/71‑1 (MQ=255)
                                            cGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTcagca                 <  1:3647297/71‑1 (MQ=255)
                                            cGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTcagca                 <  1:1178868/71‑1 (MQ=255)
                                              aTTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTcagc                  <  1:1453631/68‑1 (MQ=255)
                                                tAAACGCCATCCGGCAGCAAATCAAATAATTACTTc                                                <  1:4733738/36‑1 (MQ=255)
                                                tAAACGCCATCCGGCAGCAAATCAAATAATTACTTc                                                <  1:874660/36‑1 (MQ=255)
                                                tAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagca                             >  1:5023277/1‑55 (MQ=255)
                                                tAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAgtgt             >  1:499984/1‑71 (MQ=255)
                                                  aaCGCCATCCGGCAGCAAATCAAATAATTACTTCGc                                              >  1:2138884/1‑36 (MQ=255)
                                                  aaCGCCATCCGGCAGCAAATCAAATAATTACTTCGc                                              >  1:5603102/1‑36 (MQ=255)
                                                  aaCGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTc            >  1:5139722/1‑70 (MQ=255)
                                                           cGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTcagcag                >  1:2450013/1‑57 (MQ=255)
                                                           cGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGAc  >  1:476510/1‑71 (MQ=255)
                                                            ggCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagt                      >  1:568237/1‑50 (MQ=255)
                                                                gcaAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagtcag                   <  1:3017110/49‑1 (MQ=255)
                                                                   |                                                              
ATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGAC  >  Exported/223560‑223689

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: