Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,468,051 A→C 35.4% intergenic (‑777/+24) algD ← / ← yaaA GDP‑mannose 6‑dehydrogenase/Peroxide stress resistance protein YaaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,468,0510AC35.4% 3.7 / 2.0 27intergenic (‑777/+24)algD/yaaAGDP‑mannose 6‑dehydrogenase/Peroxide stress resistance protein YaaA
Reads supporting (aligned to +/- strand):  ref base A (10/6);  new base C (2/8);  total (12/15)
Fisher's exact test for biased strand distribution p-value = 5.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.11e-01

TATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTT  >  ACB122/1467923‑1468164
                                                                                                                                 |                                                                                                                 
tATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTAGCGGCGTCAGCcgc                                                                                                              >  7:4568/1‑134 (MQ=255)
          ggCGGTAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGt                                                                                                   >  4:414886/1‑135 (MQ=255)
          ggCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTAGCGGCTTCAGCCGCGGAGAGGGCGt                                                                                                   >  3:48698/1‑135 (MQ=255)
          ggCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTAGCGGCGTCAGCCGCGAAGAGGGCg                                                                                                    >  3:364366/1‑134 (MQ=255)
          ggCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCATAGCCGCTTCAGCCGCGCAGAGGGCg                                                                                                    >  3:236384/1‑134 (MQ=255)
          ggCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCACTTAGCGGCGTGAGCCGCCCAGCGGGCgg                                                                                                   >  3:360479/1‑134 (MQ=255)
          ggCGATAAAGGGAGGTGCGACGGATTTTTGGAGTCGGATTCGTGGCGTGATGTCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCAAATCGCGGGTTAACCCGCGAAGAGGGCg                                                                                                    >  7:73260/1‑134 (MQ=255)
                 aaGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGAGGATCACGGCGGCTGACCTCTTCGCCGGGTAAGCCGCGGAGAGGGCGGGACAGGc                                                                                            >  7:109975/1‑135 (MQ=255)
                      cggggcgcCGGTTTTTGGGCGGCGGATGCGTGGCGTGATGGGGGTTTGCCGCGTGCCCTGTGGCCGGGCAGAGGCGGATCCCGGCGGCTCACCTCTTCGCGGCTGCCGCCGCGAAGAGGGCGTGACAGGCTCACt                                                                                       <  4:364366/131‑1 (MQ=255)
                                          gTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTAGCGGCGTACGCCGCGGAGAGGGCGGGGCCAGCAAAAACCGCGGGGTGGTCgcgccg                                                                    >  8:313188/1‑132 (MQ=255)
                                              gATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCGTCGCGGCGGCCCCCGCGCAGCGGGGGTGGGACGCGCACGCCGCGGGGGGGTGGCGCCGGGacc                                                               >  1:380573/1‑134 (MQ=255)
                                                    tGGCGTGTTGGCGTTTTACAGCGGGACCTGTTCCAGGCCAGCGGGTTAAACCGGGGTATTAACTATTCGCGCCTCCCCCCGCTAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAAt                                                         <  8:73260/135‑1 (MQ=255)
                                                    gggggTGAGGGGGGTTGGCCGCGGGCCCTGTGGCCGGGCTGTGGCGTAGCCGGGCGGCTGACGTCTGCGCGGCTACAGCGGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAAt                                                         <  3:414886/131‑1 (MQ=255)
                                                    gggCGTGAGGGGGGTTCGCAGCGGGCCCTGTTGCAGGGCAGAGGCGCATCACGGCGGCTGACCTGTGCGGGGCTTCCGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAAt                                                         <  8:325449/134‑1 (MQ=255)
                                                            gggCGGTTTGCAGCGGGCCCTGTTGCAGGGCAGAGGCGGAGCCCGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGc                                                 <  4:48698/134‑1 (MQ=255)
                                                                   ttGCCGCGGGCCCGGTGGCTGGGCTGAGGCGGAGCCCGGCGGCTGACCTCTGCGCGGCTCCCGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACt                                          <  7:247580/135‑1 (MQ=255)
                                                                   tggCAGCGTGCCCTGTTGCAGGGCAGCGGCGGATCCGGCGGGCTGACCCCTTCGCGGCTCCCGCGGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACt                                          <  4:461170/133‑1 (MQ=255)
                                                                        gggggCCCTGTGGCTGGGAAGTGGCGGATCACGGCGGCGGACCTCTTCGCGGCGGCCGCGGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCg                                     <  8:4568/131‑1 (MQ=255)
                                                                        gCGGGCCCGGTGGCCGGCCAGAGGCGCCGCACGGCGGCGCACCGCTGCGCGGCTTCCGCCGCCAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCg                                     <  5:209554/135‑1 (MQ=255)
                                                                            gCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCCGCGTCAGCCGCCAAGAGGGCGTGCCAGGCGCCCGCCGCGGGGTGATCGCGCAGGAACCCCAAATGGGCTGGCTTGGCCTGCTCGGCAc                                 >  6:218888/1‑135 (MQ=255)
                                                                            gCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTAGCGGCGTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGGACACCAAATGGTGTGGCTTGGGCTGCTCGGCAc                                 >  4:118247/1‑135 (MQ=255)
                                                                                        gggagagGGGGTTCAGGGCGGCTTACTTGTTCGGGGGTACTGCGGGGATGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTAcc                     <  4:114221/132‑1 (MQ=255)
                                                                                                 tgATCCCGGCGGCTCACCTCTTCGCGGCTTCACCCGCCAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCAt            <  7:289388/134‑1 (MQ=255)
                                                                                                  gaccccggcggCTCACCTCTGCGCGGCTTCACCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTAAGCATc           <  4:364123/130‑1 (MQ=255)
                                                                                                            ggcTTACCTATGCGGGGCTGCCGCGGGTAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCtt  <  3:446734/135‑1 (MQ=255)
                                                                                                            ggcTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCtt  <  4:139490/135‑1 (MQ=255)
                                                                                                            ggcTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCtt  <  7:156475/135‑1 (MQ=255)
                                                                                                                                 |                                                                                                                 
TATGGCGTTTGGCGATAAAGGGAGGTGCGACGGATTTTTGGCGTCGGATTCGTGGCGTGATGGCGGTTTGCAGCGTGCCCTGTTGCAGGGCAGAGGCGGATCACGGCGGCTGACCTCTTCGCGGCTTCAGCCGCGAAGAGGGCGTGACAGGCTCACTCCGCGGGGTGATCGCGCAGGAACACCAAATGGTCTGGCTTGGACTGCTCGGCACTGTAGTAGTACCCCTGAGCATCGAACTGCTT  >  ACB122/1467923‑1468164

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