Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 2,318,786 A→C 45.4% G525G (GGT→GGG)  yjcS ← Putative alkyl/aryl‑sulfatase YjcS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1222,318,7860AC45.4% 7.6 / 2.4 21G525G (GGT→GGG) yjcSPutative alkyl/aryl‑sulfatase YjcS
Reads supporting (aligned to +/- strand):  ref base A (2/9);  new base C (3/7);  total (5/16)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTTGATGGTCAGGTTGCTATCGGCAGCGTTGTAGGCATTCACCCGAACACCCAGGTAGTCGAAAAACAGGCTGGGTGTCAGGGCCCGCACCATGTCATCGGCG‑TTGCCGCCGCTACCCGCACCGGCGGCAACGCCTGAGC‑GCAATTCCTGCGCGCCACTGA‑GGTAGGCATTGCGCCAGGTGGCGTTCTCACTCTGGTAACCCAGTTGCTCCAGCGCGTCGGCCTGCAGTTCGCGGGCCTCGGTGTTGTCTGGCT  >  ACB122/2318665‑2318919
                                                                                                                          |                                                                                                                                       
aGTTGATGGTCAGGTTGCTATCGGCAGCGTTGTAGGCATTCACCCGAACACCCAGGTAGTCGAAAAACAGGCTGGGTGTCAGGGCCCGCACCATGTCATCGGCC‑TTTCCCCCCCTACCCCCCCCCGCGGCAACGc                                                                                                                            >  7:133461/1‑135 (MQ=255)
                     cGGCAGCGTTGTAGGCATTCACCCGAACACCCAGGTAGTCGAAAAACAGGCTGGGTGTCAGGGCCCGCACCATGTCATCGGCG‑TTTCCCCCCCTACCCGCCCCCGCGGCAACGCCAGCGC‑GCAAATGCAgcgcgc                                                                                                      >  3:251522/1‑135 (MQ=255)
                      agacaggTTGTAGGCATTCACCCGAACACCCTGGTAGTCGAAAAACCGGCTGGGTGTCTGGGCCCGCACCATGTCAGCGGCG‑TTGCCGCCGCTACCCCCACCGGCGGCAACGCCTGAGC‑GCAATTCCTGCGCGcc                                                                                                     <  4:295178/129‑1 (MQ=255)
                       gCAGCGTTGTAGGCATTCACCCGAACACCCAGGTAGTCGAAAAACAGGCTGGGTGTCAGGGCCCGCACCATGTCATCGGCG‑TTGCCGCCGCTACCCGCACCGGCGGAAACGCCTGAGC‑GCAATTCCTGCGCGcc                                                                                                     <  3:293868/134‑1 (MQ=255)
                                 cggCATCCACCCCAACCCCCTGTGCGGCCAAAAACCGGCGGGGTGCCAGGCCCCGCCCCTTGCCAGCGGCG‑TGGCCGCCGCTACCCCCCCCGGCGGAAACGCCTGAGC‑GCAATTCCTGCGCGCCACTGA‑GGTAg                                                                                          <  8:236868/133‑1 (MQ=255)
                                 aGGCATTCACCCGAACACCCAGGTAGTCGAAAAACAGGCTGGGTGTCAGGGCCCGCACCATGTCATCGGCG‑TTGCCGCCGCTACCCGCACCGGCGGAAACGCCTGAGC‑GCAATTCCTGCGCGCCACTGA‑GGTAg                                                                                          <  7:228321/134‑1 (MQ=255)
                                       acACCCGAACCCCCAGGTCGTCCAAAAACCGGCGGGGGGTCAGGGCCCCCACCATGTCATCGGCG‑TTGCCGCCGCTACCCGCCCCGGCGGCAACGCCTGAGC‑GCAATTCCTGCGCGCCACTGA‑GGTAGGCATTgc                                                                                   <  1:228119/134‑1 (MQ=255)
                                                      ggTAGTCGAAAAACAGGCTGGGTGTCAGGGCCCGCACCATGTCATCGGCGTTTGCCG‑CGCGGCCCCCCCCCGCGGCCAAGCCCGCGCGGCAATTCCGGCG‑GCCCCGGAGGGTGGGCATTG‑GCCAGgggggggtgt                                                                    >  5:378336/1‑129 (MQ=255)
                                                      ggTAGTCGAAAAACAGGCTGGGTGTCAGGACCCGCACCATGTCATCGGCG‑TTGCCGGCCCGACCCGCACCCGCGGCCACGCCGGCGC‑GGCATTCCAGCGCGCCCCGGA‑GGGAGGCCTTGCGCCCGGTgggggtcg                                                                    >  7:372518/1‑129 (MQ=255)
                                                      ggTAGCCCAAAAACCGGCGGGGTGACAGGGCCCGCACCATGTCATCGGCG‑TTGCCGCCGCTACCCGCCCCGGCGGAAACGCCTGAGC‑GCAATTCCTGCGCGCCACTGA‑GGTAGGCATTGCGCCAGGTGGCGTtct                                                                    <  3:141850/135‑1 (MQ=255)
                                                        atgTCCAAAACCGGGCGGGGTTTCAGGGGCCCCACCAGGTCTTGGGCG‑GGGCCGCCGAACCCCGCACCGGCGCAAACGCCTGAGC‑GCAATTCCTGCGCGCCACTGA‑GGTAGGCATTGCGCCAGGTGGCGTtct                                                                    <  1:96369/131‑1 (MQ=255)
                                                              aaaaaCAGGCTGGGTGTCAGGGCCCGCACCATGTCATCGGCA‑TTGCCGCCGCTACCCGCACCGGCGGCAACGCCTGCGC‑GCAATTCCTGCGCGCCACTGA‑GGTAGGCATTGCGCCAGGTGGCGTTCTCACTCTg                                                             >  6:88684/1‑134 (MQ=255)
                                                                    aGGCTGGGTGTCAGGGCCCGCACCATGTCATCGGCG‑TTGCCGCCGCTACCCGCACCGGCGGCAACGCCTGAGC‑GCAATTCCTGCGCGCCACTGA‑GGTAGGCATTGCGCCAGGTGGCGTTCTCACTCTGGTAAccc                                                      <  6:345777/135‑1 (MQ=255)
                                                                           gtgtCAGGGCCCGCACCATGTCATCGGCG‑TTGCCGCCGCTACCCGCCCCGGCGGCAACGCCTGAGC‑GCAATTCCTGCGCGCCACTGA‑GGTAGGCATTGCGCCAGGTGGCGTTCTCACTCTGGTAACCCAGTTGCt                                               <  8:69107/135‑1 (MQ=255)
                                                                                ccgcccccCACCATGCCTGGGGGG‑TGGCCGCCGCCACCCCCCCCGCCGCAACCGCCTGAGC‑GCAATTCCTGCGCGCCACTGA‑GGTAGGCATTGCGCCAGGTGGCGTTCTCACTCTGGTAACCCAGTTGCTCCAgc                                          <  1:126869/131‑1 (MQ=255)
                                                                                                 aTCGGCG‑TGCCCGCCGCTACCCGCACCGGCGGCAACGCCTGAGC‑GCAATTCCTGCGCGCCACTGA‑GGTAGGCATTGCGCCAGGTGGCGTTCTCACTCTGGTAACCCAGTTGCTCCAGCGCGTCGGCCTGCAGTTc                         <  7:70163/135‑1 (MQ=255)
                                                                                                    ggCG‑TTGCCGCCGCTACCCGCACCGGCGGCAACGCCTGAGC‑GCAATTCCTGCGCGCCACTGA‑GGTAGGCATTGCGCCAGGTGGCGTTCTCACTCTGGTAACCCAGTTGCTCCAGCGCGTCGGCCTGCAGTTcgcg                      <  4:165980/135‑1 (MQ=255)
                                                                                                      cG‑TGGCCGCCGCTACCCGCACCGGCGGCAACGCCTGAGC‑GCAATTCCTGCGCGCCACTGA‑GGTAGGCATTGCGCCAGGTGGCGTTCTCACTCTGGTAACCCAGTTGCTCCAGCGCGTCGGCCTGCAGTTCGCggg                    <  7:178526/135‑1 (MQ=255)
                                                                                                       gtgggccgccGCTACCCCCCCGGGCGGAAAGGCCTGAGC‑GCAATTCCTGCGCGCCACTGA‑GGTAGGCATTGCGCCAGGTGGCGTTCTCACTCTGGTAACCCAGTTGCTCCAGCGCGTCGGCCTGCAGTTCGCggg                    <  7:358099/131‑1 (MQ=255)
                                                                                                                cgcTACCCGCCCCGGCGGCAACGCCTGAGC‑GCAATTCCTGCGCGCCACTGA‑GGTAGGCATTGCGCCAGGTGGCGTTCTCACTCTGGTAACCCAGTTGCTCCAGCGCGTCGGCCTGCAGTTCGCGGGCCTCGgtgt           <  5:152856/135‑1 (MQ=255)
                                                                                                                         cACCGGCGGAAACGCCTGAGC‑GCAATTCCTGCGCGCCACTGA‑GGTAGGCATTGCGCCAGGTGGCGTTCTCACTCTGGTAACCCAGTTGCTCCAGCGCGTCGGCCTGCAGTTCGCGGGCCTCGGTGTTGTCTGGCt  <  2:63117/135‑1 (MQ=255)
                                                                                                                          |                                                                                                                                       
AGTTGATGGTCAGGTTGCTATCGGCAGCGTTGTAGGCATTCACCCGAACACCCAGGTAGTCGAAAAACAGGCTGGGTGTCAGGGCCCGCACCATGTCATCGGCG‑TTGCCGCCGCTACCCGCACCGGCGGCAACGCCTGAGC‑GCAATTCCTGCGCGCCACTGA‑GGTAGGCATTGCGCCAGGTGGCGTTCTCACTCTGGTAACCCAGTTGCTCCAGCGCGTCGGCCTGCAGTTCGCGGGCCTCGGTGTTGTCTGGCT  >  ACB122/2318665‑2318919

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