Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 5,563 G→A 15.4% intergenic (+33/+120) yaaX → / ← yaaA hypothetical protein/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported5,5630GA15.4% 89.2 / 8.2 39intergenic (+33/+120)yaaX/yaaAhypothetical protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (20/13);  new base A (4/2);  total (24/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

ACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCT  >  Exported/5496‑5632
                                                                   |                                                                     
aCGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAgcgc                                                                    >  1:4793265/1‑71 (MQ=255)
         aTGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGa                                                           >  1:168099/1‑71 (MQ=255)
          tGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGAt                                                               >  1:3303570/1‑66 (MQ=255)
          tGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGAc                                                          >  1:5011463/1‑71 (MQ=255)
             tCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGAt                                                               >  1:1736092/1‑63 (MQ=255)
               cAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTgg                                                     >  1:4560764/1‑71 (MQ=255)
               cAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTg                                                      <  1:4487782/70‑1 (MQ=255)
               cAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTg                                                      <  1:1934310/70‑1 (MQ=255)
                           aCCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCa                                         <  1:82975/71‑1 (MQ=255)
                             cGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCt                                                                  <  1:681862/44‑1 (MQ=255)
                                 aaaTGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTa                                   >  1:4630571/1‑71 (MQ=255)
                                     gACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGt                                >  1:3674851/1‑70 (MQ=255)
                                      aCAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGtt                               <  1:3212405/70‑1 (MQ=255)
                                      aCAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTa                              <  1:3452523/71‑1 (MQ=255)
                                       cAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTaa                             >  1:4630314/1‑71 (MQ=255)
                                       cAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTaa                             <  1:4412898/71‑1 (MQ=255)
                                            gCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTg                        >  1:1866947/1‑71 (MQ=255)
                                                ggTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTg                        <  1:53032/67‑1 (MQ=255)
                                                ggTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTg                        <  1:910727/67‑1 (MQ=255)
                                                 gTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGTCGCGTCTTAt                                           >  1:2487018/1‑47 (MQ=255)
                                                 gTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTAt                                           >  1:3831613/1‑47 (MQ=255)
                                                  tAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGt                                                >  1:5568036/1‑41 (MQ=255)
                                                    cgcatccGGTATGTAACGCCTGATGCGACGCTGACGCGTCTTATCTGGCCTACacgctgcgattttgtagg                >  1:5540228‑M1/4‑53 (MQ=255)
                                                    aCAATCAGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAatatat                >  1:740861/1‑71 (MQ=255)
                                                     cAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATAtt               <  1:2741661/71‑1 (MQ=255)
                                                       aTCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTg              >  1:5468280/1‑70 (MQ=255)
                                                       aTCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGa             <  1:4792640/71‑1 (MQ=255)
                                                          cGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAgg                                       >  1:340860/1‑42 (MQ=255)
                                                          cGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATc          >  1:4198018/1‑71 (MQ=255)
                                                          cGGCATTCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggccggat          <  1:3323539‑M1/71‑25 (MQ=255)
                                                           ggCATTCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACgcgctgcgattttgtaggccggat          >  1:3154301‑M1/1‑46 (MQ=255)
                                                            gCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTg        >  1:1033395/1‑71 (MQ=255)
                                                             cATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTg        >  1:3429588/1‑70 (MQ=255)
                                                             cATTCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggccggataag       <  1:577759‑M1/71‑28 (MQ=255)
                                                               ttCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCt         <  1:4489936/67‑1 (MQ=255)
                                                               ttCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggccggataagc      >  1:264260‑M1/1‑42 (MQ=255)
                                                               ttCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggccgg            >  1:5316084‑M1/1‑42 (MQ=255)
                                                                tCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCt                                    <  1:5617861/39‑1 (MQ=255)
                                                                tCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCAt     >  1:1566045/1‑70 (MQ=255)
                                                                  aGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCt  >  1:626255/1‑71 (MQ=255)
                                                                   |                                                                     
ACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCT  >  Exported/5496‑5632

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: