Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,845,620 T→G 51.1% intergenic (+47/+72) apt → / ← recR Adenine phosphoribosyltransferase/Recombination protein RecR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,845,6200TG51.1% 2.0 / 2.3 24intergenic (+47/+72)apt/recRAdenine phosphoribosyltransferase/Recombination protein RecR
Reads supporting (aligned to +/- strand):  ref base T (4/6);  new base G (5/4);  total (10/14)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACTGGACGGGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATG  >  ACB122/4845493‑4845727
                                                                                                                                 |                                                                                                             
gAACTGGACGGGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTTGGGCTGCTTTGCAGCCCATTAGCGGGGTAAACCCGGGCCTag                                                                                                          >  7:463787/1‑133 (MQ=255)
       ccGGGGCGGGTCGGTTGCACGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGGGAGTATTTAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCcgc                                                                                                 <  7:124164/134‑1 (MQ=38)
                    gTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAg                                                                                   <  4:153168/135‑1 (MQ=255)
                     ttGCACGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAgg                                                                                  <  8:384273/135‑1 (MQ=255)
                                  gggTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGTGGCTGCTTTGCAGCCCATTCGCGGGGTAAACCCCGCCGACGCCAGGCGGGGATTGTTTTGTGGCGCGGGGTTTccc                                                                     >  2:419214/1‑135 (MQ=255)
                                  gggTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTTAACCCGCGGCGACGCAAGCCGCCGTTGTTTTTTATGGGCGGGgtttgcc                                                                     >  8:111411/1‑131 (MQ=25)
                                            cTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAAACCGCGCCGACGGAAGCCGCGGTTGTTTTGTAGGCCCGGGTTTTCCCGCGAATTGgc                                                           >  8:357957/1‑135 (MQ=255)
                                                tGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTTATACCCCCGCGCCCCAAGCCGCGTTTGGTTTGTATGCCCGGGGTTAACCGCGCATAGggcgcc                                                       >  8:433208/1‑135 (MQ=255)
                                                    tGACCTAGTTTTCCCGCACCCTGTATGTAGTGGGCTTTTGGCGACGCTTTCAAGCCAATTCGCGGGTAAACCCGCGCCCACCAACGCCGCGTATGTTTGGAAGGGGCGGTTTAACCCGCTAATGGGCCGCCAAGc                                                   <  4:332312/135‑1 (MQ=255)
                                                    tGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCATCCCATTCGCGGGTAAACCCGCTCCTACGAAACCCGCGGATTTTTTTTATTCTCGGGTTTAAACGACAATTGGCCGGCatcc                                                   >  8:455124/1‑132 (MQ=25)
                                                         gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCAACGCAAGCCGCGGATTTTTTTTTGGCGCGGGGTTTCCCGCGCATGGGCCGCCACGCGGccc                                              >  8:354528/1‑135 (MQ=37)
                                                         gAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCCCTCCGACCACATCCGCCGAGGTTGTGGAGGCCCGGGTTTACCCCCGAATGGGCCGCCACGCGGccc                                              >  7:159439/1‑135 (MQ=37)
                                                               ccccTCCGCGTGTTTGGAGTGGGCTTTGGGGGCGGTTTGGCCGCCATTAGGCGGTTAAACCCCCCCCACCTAACGCCGGGTTTGTTTGGTCGGCGGGGGTATACCGGCTAATGGGCCGCCAAGCGGCCCCAATg                                         <  6:124198/133‑1 (MQ=255)
                                                                                  tGGGCTTTTGGGGCTGCTTTGCAGCCCATGCGCGGGTAACCCCGCTCCTCCGAAAGCCGCGTATGTTTTGTCGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCg                    <  7:39044/135‑1 (MQ=255)
                                                                                  ggggCTTTGGGGGCTGCTTGCCAGCCCATTCGCGGGTAAACCCGCGCCAACGAACGCCGCTTTGGTTTTGTCGGCGCGGGTTAACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCg                    <  2:113119/134‑1 (MQ=25)
                                                                                  ggggCTTTGGGGGCGGCTTGGCAGCCCATTCGCGGGAAACCCCGCCCCTCCGAACGCCGCGTTTGTTTGGTCGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCg                    <  4:233178/134‑1 (MQ=37)
                                                                                   gggCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTTAACCCGCGGCGGCGCAATCCGCGGATTTTTTGTTGGCGCGGGTTTAGCCGCGGATGGGGCGCCACCGGGGCCCAAAGAATACGCGCGGAAAGGGCCgc                   >  4:55905/1‑134 (MQ=16)
                                                                                   gggCTTTTGGGGCTGCTTTGAAGCCCATTCGCGGGTTAACCCGCTGCTACGCAAGCCGCGGATGTTTTTTATGCGCGGGGTTAACCCCGGATGGGCCCCCCAGCCGGCCCAATGAATAAGAGCGGAATTGGGCCg                    >  2:123811/1‑135 (MQ=255)
                                                                                          tGGGGCGTGCTTGCCCGCCCAGTCGCGGTAAACCCCGCCCCACCCAACGCCGCGTAGGTTTGGTCGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCt             <  1:148384/135‑1 (MQ=25)
                                                                                           ggggCGGCTTGGCCGCCCTTTCGCGGGAAAACCCGCGCCACCAACCGCCGGGTATTTTTTGTCGGCGCGGTTAACCCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTg            <  5:295709/134‑1 (MQ=255)
                                                                                                     ttGCAGCCCATTCGCGGGTAAACCCGCTCCTCCGAACGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATg  <  2:51520/135‑1 (MQ=255)
                                                                                                     ttGCAGCCAATTCGGGGGAAACCCCGCACCTACCAACGCCGGTTATGTTTTGTCGGCGGGGGTTAACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATg  <  6:409237/135‑1 (MQ=255)
                                                                                                     tgccAGCCAATTCGCGGGTAAACCGGCGCATACGAACGCCGCTTATGTTTTGTCGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATg  <  7:111411/132‑1 (MQ=38)
                                                                                                       gccgCCCATTCGCGGGTAACCCCGCTCCTACGAACGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATg  <  3:108142/130‑1 (MQ=255)
                                                                                                                                 |                                                                                                             
GAACTGGACGGGTCGCGTCGGTTGCAAGCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATG  >  ACB122/4845493‑4845727

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