Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,974,074 A→C 63.9% Q688P (CAG→CCG)  sasA_19 → Adaptive‑response sensory‑kinase SasA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,974,0740AC63.9% ‑4.7 / 6.2 25Q688P (CAG→CCG) sasA_19Adaptive‑response sensory‑kinase SasA
Reads supporting (aligned to +/- strand):  ref base A (3/4);  new base C (11/4);  total (15/10)
Fisher's exact test for biased strand distribution p-value = 3.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.59e-01

TGGCTTCCTGCAACCGGGCCAGAACGGCTTGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCTGGCCATGCGCCAGGCTTGGCCAGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTACAACGCGCCAGCCCCTGCGGGTTTGCAGCGGATGACTGCGAAGCGTTCCTGCAGCTTTTACGCGAAGCCGCATCATGCGC  >  ACB122/5973946‑5974206
                                                                                                                                  |                                                                                                                                    
tGGCTTCCTGCAACCGGGCCAGAACGGCTTGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCAGGCCCGGCGCCCGGCTTg                                                                                                                                >  7:285786/1‑135 (MQ=255)
                   cAGAACGGCTTGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCCGGGCATGCGCCCGGGGTGGCCCGGGGGCACGCCCgcg                                                                                                             >  7:132521/1‑133 (MQ=255)
                              gggCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCAGGCCCTGCGCCCGGCATGGCCAGGGCGCACGCCCCAGGggttggtgtg                                                                                                  >  3:123162/1‑130 (MQ=255)
                                        cATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCAGGCCAATCGCCCGGCATGGCCCGGGCGCACGCACGACGGGTTGTTGGGGCTCgccacc                                                                                        >  1:45409/1‑133 (MQ=255)
                                                   gATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCGCGCCCTGCGCCCGGCGTGGGCCCGGCGCAAGCCCCAGGGGGTGTTTGCGCGCGCCCGCCCGCGTCAAcc                                                                             >  1:82158/1‑134 (MQ=255)
                                                   gATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCAGGCCCAGCGCCCCGCGTGGGCACGGCGCCCGCCCGAGGGGGTTTTTGCGCTCGCCCCCCACCTACaccc                                                                             >  2:369795/1‑132 (MQ=255)
                                                   gATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCAGGCCATGCGCCAGGCTTGGGCAGGGCGCACGCACCACGTGTTGTTAGAGCTCGCCAGCCAGCTACAAcg                                                                             >  8:215565/1‑135 (MQ=255)
                                                    aTGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCTGGCCCAGCGCCCCGCGTGGGCCGGGGGGACGGACGAGGGGGTGTTTGCGCTCGCGAGCCCGCTAAAAccc                                                                            >  4:253297/1‑133 (MQ=255)
                                                    aTGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCCGGCCCAGCGGCCGGCGTGGGCCCGGCGGACCCCCCCGGTGTTGTTTGGGCTGGCCCGCCCCCGACAAccc                                                                            >  1:241674/1‑133 (MQ=255)
                                                     ggggggcggcTTATCCACCCCCCCCGGGGCACCTCCGCGCTGCGCTGCAGGCAGTAATCCCCCTGGCCATGCGCGCGGCTGGGCCAGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTACAAcgcg                                                                           <  4:109000/131‑1 (MQ=255)
                                                      ggAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCAGGCACGTATCGCCCTGGCCAGGCGCCAGGCTTGGCCAGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTACAAcgcg                                                                           <  8:380497/134‑1 (MQ=255)
                                                            cGCTTAGGCACTCACGCCCGTGCACCCGCGCGCCTCGCGGCATGCACGTATCGCCCTGGCCATGCGCCCGGCTGGGCCAGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTACAACGCGCCAGccc                                                                    <  2:45409/135‑1 (MQ=255)
                                                              cTTATGCACTCCCGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCGGGCCCTGCGGCCCGCCTGGGCACGGCGGACGGCGCCGGGGGTGTTCGCGCTCGCCCGCCCGCGCAAAAGCGCCCGccccg                                                                  >  8:293620/1‑134 (MQ=255)
                                                              cTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCAGGCCCTTCGCCACGCTTGGCCAGGGCGCACGCAACATGTGGTGTTAGAGCTAGCCAGCCCGCGTCAACGCGCCCGCCCCt                                                                  >  8:282149/1‑135 (MQ=255)
                                                                    cACTCCCGCCGGTGCACCAGCGCGCCTGGCGCCTGGACCGTATCGCCGTGGCCATCCGCATGGGTTGGCCAGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTACAACGCGCCAGCCCCTGCGGGt                                                            <  1:306076/135‑1 (MQ=255)
                                                                            cccGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCGGGCCCTGCGCCCGGGGTGGGCCGGGCGCACGCCCGAGGGGGTGGTAGAGGTCGCCCGCCCGCGGCCACCCGCCCGCCCCGGCGGGGTTggctcg                                                    >  4:32526/1‑130 (MQ=255)
                                                                            cccGTGCACCAGCGCGCCTCGCGGCATGCACGTATCGCCCTGGCCATGCGCCAGGCTTGGCCAGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTACAACGCGCCAGCCCCTGCGGGTTTGCAGCg                                                    <  3:206648/135‑1 (MQ=255)
                                                                                   aCCAGCGCGCCTCGCTGCATGCACGTATCGCCCGGGGCAGGCGGCCGGGGTTGGCCGGGCGGTCGCCGCCGGTGGGGTTGGCGCTCGCCCACCCGCGGCAACACGCCCGCCCCCGCGGGTTTGTCGCGGGTgcca                                             >  8:253559/1‑132 (MQ=255)
                                                                                   aCCAGCGCGCCTCGCTGCATGCACGTATCGCCCCGGCCCAGCGCCCGGCCTGGGCCGGGGGCCCGCACGCGGGGGTGGTGGGGCGCGGCCGCCAGCGTCCCCCCGCCCGCCCCGGCGGGGTTGCAGCggggggcg                                             >  6:393713/1‑129 (MQ=255)
                                                                                               ggCTGCAGGCAGGTAGCGCCGTGGCCATGCGCCAGGCTTGGCCCGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTACAACGCGCCAGCCCCTGCGGGTTTGCAGCGGATGACTGCGAAGCGTTcc                                 <  7:85694/134‑1 (MQ=255)
                                                                                                         aCGTCTCGCCCTGGCCAGGCGCCCGGCTGGGCCCGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTACAACGCGCCAGCCCCTGCGGGTTTGCAGCGGATGACTGCGAAGCGTTCCTGCAGCtttt                       <  6:351930/135‑1 (MQ=255)
                                                                                                         aCGTAGCGCCCGGCCCAGCCGCCCGGCTGGCCCCGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTACAACGCGCCAGCCCCTGCGGGTTTGCAGCGGATGACTGCGAAGCGTTCCTGCAGCtttt                       <  1:320942/135‑1 (MQ=255)
                                                                                                            gTAGCGCCCTGCCCATGCGCCAGGCTTGGCCAGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTACAACGCGCCAGCCCCTGCGGGTTTGCAGCGGATGACTGCGAAGCGTTCCTGCAGCTTTTAc                     <  1:54714/135‑1 (MQ=255)
                                                                                                                     tGGCCATGCGCCTGGCTTGGCCTGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTACAACGCGCCAGCCCCTGCGGGTTTGCAGCGGATGACTGCGAAGCGTTCCTGCAGCTTTTACGCGAAGCCg            <  6:241069/135‑1 (MQ=255)
                                                                                                                                ccAGGCTTGGCCAGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTACAACGCGCCAGCCCCTGCGGGTTTGCAGCGGATGACTGCGAAGCGTTCCTGCAGCTCTTACGCGAAGCCGCATCATgcgc  >  7:280058/1‑135 (MQ=255)
                                                                                                                                  |                                                                                                                                    
TGGCTTCCTGCAACCGGGCCAGAACGGCTTGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGCCCTGGCCATGCGCCAGGCTTGGCCAGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTACAACGCGCCAGCCCCTGCGGGTTTGCAGCGGATGACTGCGAAGCGTTCCTGCAGCTTTTACGCGAAGCCGCATCATGCGC  >  ACB122/5973946‑5974206

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: