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→GGGyjcSPutative 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

AGTTGATGGTCAGGTTGCTATCGGCAGCGTTGTAGGCATTCACCCGAACACCCAGGTAGTCGAAAAACAGGCTGGGTGTCAGGGCCCGCACCATGTCATCGGCGTTGCCGCCGCTACCCGCACCGGCGGCAACGCCTGAGCGCAATTCCTGCGCGCCACTGAGGTAGGCATTGCGCCAGGTGGCGTTCTCACTCTGGTAACCCAGTTGCTCCAGCGCGTCGGCCTGCAGTTCGCGGGCCTCGGTGTTGTCTGGCT  >  ACB122/2318665‑2318919
                                                                                                                          |                                                                                                                                       
aGTTGATGGTCAGGTTGCTATCGGCAGCGTTGTAGGCATTCACCCGAACACCCAGGTAGTCGAAAAACAGGCTGGGTGTCAGGGCCCGCACCATGTCATCGGCCTTTCCCCCCCTACCCCCCCCCGCGGCAACGc                                                                                                                            >  7:133461/1‑135 (MQ=255)
                     cGGCAGCGTTGTAGGCATTCACCCGAACACCCAGGTAGTCGAAAAACAGGCTGGGTGTCAGGGCCCGCACCATGTCATCGGCGTTTCCCCCCCTACCCGCCCCCGCGGCAACGCCAGCGCGCAAATGCAgcgcgc                                                                                                      >  3:251522/1‑135 (MQ=255)
                      agacaggTTGTAGGCATTCACCCGAACACCCTGGTAGTCGAAAAACCGGCTGGGTGTCTGGGCCCGCACCATGTCAGCGGCGTTGCCGCCGCTACCCCCACCGGCGGCAACGCCTGAGCGCAATTCCTGCGCGcc                                                                                                     <  4:295178/129‑1 (MQ=255)
                       gCAGCGTTGTAGGCATTCACCCGAACACCCAGGTAGTCGAAAAACAGGCTGGGTGTCAGGGCCCGCACCATGTCATCGGCGTTGCCGCCGCTACCCGCACCGGCGGAAACGCCTGAGCGCAATTCCTGCGCGcc                                                                                                     <  3:293868/134‑1 (MQ=255)
                                 cggCATCCACCCCAACCCCCTGTGCGGCCAAAAACCGGCGGGGTGCCAGGCCCCGCCCCTTGCCAGCGGCGTGGCCGCCGCTACCCCCCCCGGCGGAAACGCCTGAGCGCAATTCCTGCGCGCCACTGAGGTAg                                                                                          <  8:236868/133‑1 (MQ=255)
                                 aGGCATTCACCCGAACACCCAGGTAGTCGAAAAACAGGCTGGGTGTCAGGGCCCGCACCATGTCATCGGCGTTGCCGCCGCTACCCGCACCGGCGGAAACGCCTGAGCGCAATTCCTGCGCGCCACTGAGGTAg                                                                                          <  7:228321/134‑1 (MQ=255)
                                       acACCCGAACCCCCAGGTCGTCCAAAAACCGGCGGGGGGTCAGGGCCCCCACCATGTCATCGGCGTTGCCGCCGCTACCCGCCCCGGCGGCAACGCCTGAGCGCAATTCCTGCGCGCCACTGAGGTAGGCATTgc                                                                                   <  1:228119/134‑1 (MQ=255)
                                                      ggTAGTCGAAAAACAGGCTGGGTGTCAGGGCCCGCACCATGTCATCGGCGTTTGCCGCGCGGCCCCCCCCCGCGGCCAAGCCCGCGCGGCAATTCCGGCGGCCCCGGAGGGTGGGCATTGGCCAGgggggggtgt                                                                    >  5:378336/1‑129 (MQ=255)
                                                      ggTAGTCGAAAAACAGGCTGGGTGTCAGGACCCGCACCATGTCATCGGCGTTGCCGGCCCGACCCGCACCCGCGGCCACGCCGGCGCGGCATTCCAGCGCGCCCCGGAGGGAGGCCTTGCGCCCGGTgggggtcg                                                                    >  7:372518/1‑129 (MQ=255)
                                                      ggTAGCCCAAAAACCGGCGGGGTGACAGGGCCCGCACCATGTCATCGGCGTTGCCGCCGCTACCCGCCCCGGCGGAAACGCCTGAGCGCAATTCCTGCGCGCCACTGAGGTAGGCATTGCGCCAGGTGGCGTtct                                                                    <  3:141850/135‑1 (MQ=255)
                                                        atgTCCAAAACCGGGCGGGGTTTCAGGGGCCCCACCAGGTCTTGGGCGGGGCCGCCGAACCCCGCACCGGCGCAAACGCCTGAGCGCAATTCCTGCGCGCCACTGAGGTAGGCATTGCGCCAGGTGGCGTtct                                                                    <  1:96369/131‑1 (MQ=255)
                                                              aaaaaCAGGCTGGGTGTCAGGGCCCGCACCATGTCATCGGCATTGCCGCCGCTACCCGCACCGGCGGCAACGCCTGCGCGCAATTCCTGCGCGCCACTGAGGTAGGCATTGCGCCAGGTGGCGTTCTCACTCTg                                                             >  6:88684/1‑134 (MQ=255)
                                                                    aGGCTGGGTGTCAGGGCCCGCACCATGTCATCGGCGTTGCCGCCGCTACCCGCACCGGCGGCAACGCCTGAGCGCAATTCCTGCGCGCCACTGAGGTAGGCATTGCGCCAGGTGGCGTTCTCACTCTGGTAAccc                                                      <  6:345777/135‑1 (MQ=255)
                                                                           gtgtCAGGGCCCGCACCATGTCATCGGCGTTGCCGCCGCTACCCGCCCCGGCGGCAACGCCTGAGCGCAATTCCTGCGCGCCACTGAGGTAGGCATTGCGCCAGGTGGCGTTCTCACTCTGGTAACCCAGTTGCt                                               <  8:69107/135‑1 (MQ=255)
                                                                                ccgcccccCACCATGCCTGGGGGGTGGCCGCCGCCACCCCCCCCGCCGCAACCGCCTGAGCGCAATTCCTGCGCGCCACTGAGGTAGGCATTGCGCCAGGTGGCGTTCTCACTCTGGTAACCCAGTTGCTCCAgc                                          <  1:126869/131‑1 (MQ=255)
                                                                                                 aTCGGCGTGCCCGCCGCTACCCGCACCGGCGGCAACGCCTGAGCGCAATTCCTGCGCGCCACTGAGGTAGGCATTGCGCCAGGTGGCGTTCTCACTCTGGTAACCCAGTTGCTCCAGCGCGTCGGCCTGCAGTTc                         <  7:70163/135‑1 (MQ=255)
                                                                                                    ggCGTTGCCGCCGCTACCCGCACCGGCGGCAACGCCTGAGCGCAATTCCTGCGCGCCACTGAGGTAGGCATTGCGCCAGGTGGCGTTCTCACTCTGGTAACCCAGTTGCTCCAGCGCGTCGGCCTGCAGTTcgcg                      <  4:165980/135‑1 (MQ=255)
                                                                                                      cGTGGCCGCCGCTACCCGCACCGGCGGCAACGCCTGAGCGCAATTCCTGCGCGCCACTGAGGTAGGCATTGCGCCAGGTGGCGTTCTCACTCTGGTAACCCAGTTGCTCCAGCGCGTCGGCCTGCAGTTCGCggg                    <  7:178526/135‑1 (MQ=255)
                                                                                                       gtgggccgccGCTACCCCCCCGGGCGGAAAGGCCTGAGCGCAATTCCTGCGCGCCACTGAGGTAGGCATTGCGCCAGGTGGCGTTCTCACTCTGGTAACCCAGTTGCTCCAGCGCGTCGGCCTGCAGTTCGCggg                    <  7:358099/131‑1 (MQ=255)
                                                                                                                cgcTACCCGCCCCGGCGGCAACGCCTGAGCGCAATTCCTGCGCGCCACTGAGGTAGGCATTGCGCCAGGTGGCGTTCTCACTCTGGTAACCCAGTTGCTCCAGCGCGTCGGCCTGCAGTTCGCGGGCCTCGgtgt           <  5:152856/135‑1 (MQ=255)
                                                                                                                         cACCGGCGGAAACGCCTGAGCGCAATTCCTGCGCGCCACTGAGGTAGGCATTGCGCCAGGTGGCGTTCTCACTCTGGTAACCCAGTTGCTCCAGCGCGTCGGCCTGCAGTTCGCGGGCCTCGGTGTTGTCTGGCt  <  2:63117/135‑1 (MQ=255)
                                                                                                                          |                                                                                                                                       
AGTTGATGGTCAGGTTGCTATCGGCAGCGTTGTAGGCATTCACCCGAACACCCAGGTAGTCGAAAAACAGGCTGGGTGTCAGGGCCCGCACCATGTCATCGGCGTTGCCGCCGCTACCCGCACCGGCGGCAACGCCTGAGCGCAATTCCTGCGCGCCACTGAGGTAGGCATTGCGCCAGGTGGCGTTCTCACTCTGGTAACCCAGTTGCTCCAGCGCGTCGGCCTGCAGTTCGCGGGCCTCGGTGTTGTCTGGCT  >  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: