Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 2,041,924 T→G 35.4% intergenic (+106/+33) yfcG_2 → / ← CJ019_01936 Disulfide‑bond oxidoreductase YfcG/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1432,041,9240TG35.4% 23.1 / 6.5 40intergenic (+106/+33)yfcG_2/CJ019_01936Disulfide‑bond oxidoreductase YfcG/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (11/14);  new base G (4/10);  total (15/25)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACAC‑GGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCTTGCCGGCTTCTTCAATGCGCAGGCATGGCACT  >  ACB143/2041793‑2042056
                                                                                                                                    |                                                                                                                                    
cgcgGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACAC‑GGGCGCCGCCCGGGCCAACCACCGCGGGGGCTGCGTCGCCGGCGCGcc                                                                                                                                  >  3:265144/1‑136 (MQ=37)
   ggCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGCACAC‑CCGCCCCCCCCCTGCCCACCACCCCGGTGCCAGCTTCGCCGGCTCGccccc                                                                                                                               >  2:226097/1‑134 (MQ=255)
                      cTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACAC‑GGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTAGCGGGCGCGCCCCCGCCCACAAGGGATCTGCAg                                                                                                            >  1:215988/1‑136 (MQ=255)
                       tGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGCACAC‑GGGCGCAGCCCGTGCCAACCCCCGCGGTGCCTGCTTCGCGGGCGCGCCCGCGCCCACAAGGGATCTggcgc                                                                                                           >  5:204350/1‑132 (MQ=37)
                                       ttCCAGGCGCTTGCGCAATCGTTTTGGCAGCGGGCGTGCCCGCTAACAC‑GGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATggcg                                                                                            <  6:324274/135‑1 (MQ=255)
                                       tgcccggCGCTTGCGAATGCGTGGGGGGCGCGGGCGGCCCCGCAAACAC‑GGGCGCAGCCCGGGCCAACCCCCGCGGTGCCTGCTTCGCGGGCGCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATggcg                                                                                            <  8:164196/130‑1 (MQ=21)
                                       tccccggCGCTGGCCAAAGCGTTGGGGGAGGGGGCGGCCCGGCAACCAC‑GGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCGCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATggcg                                                                                            <  3:353169/130‑1 (MQ=25)
                                                               tggccgcgggcggGCCCGCTAACAG‑GGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTGCGCGGGCTCGCCCGCGCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCga                                                                   <  8:204102/131‑1 (MQ=255)
                                                               ggggCGCGGGGGTGCCCGCTAACAC‑GGTCGCAGCCGGTGCCAACCACCGGGGTGCCTGCTTGGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCga                                                                   <  2:1755/135‑1 (MQ=38)
                                                               ggggCGCGGGCGGCCCCGGGAACAC‑GGGCGCAGCCCGTGCAAACCCCCGCGGGGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCga                                                                   <  5:326002/135‑1 (MQ=255)
                                                                gggCGCGGGGGGGCCCGCCAACAC‑GGGCGCAGCCCGTGCCAACCACCGCGGTGCCGGCTTCGCGGGCTCGCCCGCCCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCgat                                                                  <  6:263196/136‑1 (MQ=255)
                                                                     ggggcgggcCGGCTAACAC‑GGGCGCAGCCCGTGCCAACCACCGCGGGGCCTGCTTCGCGGGCGCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTg                                                             <  4:197023/135‑1 (MQ=38)
                                                                         ggTGCCCGCGAACAC‑GGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCtt                                                         <  8:2076/135‑1 (MQ=255)
                                                                         cttccccGCGAACAC‑GGGCGCAGCCCGTGCCAACCACCGAGGTGCCTCGTTCGCGTGCTCGCCCTCTCCCACAAGGTATCTGCAGCGTCAGTCAGATGGCGGAAAACCGCTCATCCAGGTAAGCGATGATTGCCtt                                                         <  2:277731/132‑1 (MQ=255)
                                                                         cgggCCCGCTAACAG‑GGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCtt                                                         <  8:47086/133‑1 (MQ=255)
                                                                         cggccccGCCAACAC‑GGGCGCAGCCCGTGCCAACCACCGCGGGGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCtt                                                         <  8:199015/132‑1 (MQ=255)
                                                                          ggccccGCCAACACGGGGCCCAGCCCGTGCCAACCACGGCGGGGCCGGCTTCGCGGGCGCGCCCGCGCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCtt                                                         <  4:193123/133‑1 (MQ=255)
                                                                          ggccccGCAAACCG‑GGGCGCCGCCGGTGCCAACCACCGCGGGGCCGGCTTCGCGGGCCCGCCCGCCCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTg                                                        <  7:172594/133‑1 (MQ=255)
                                                                              ccGCTAACAC‑GGGCGCAGCCGGTGACAACCACCGCGGGGCCTGCTTCGCGGGCGCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACt                                                    <  2:98818/136‑1 (MQ=39)
                                                                                gcgAACAC‑GGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCGCCCCCAAGGGATCTGCAGCGTCAGTCAGATGGCGGGAAAACGCTCATCCAGGGTAGCGATGATTGGCTTggggtgc                                                  >  1:223571/1‑131 (MQ=38)
                                                                                gcgAACAC‑GGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCGCGCCCGCTCCCACAAGGGATCTGCCGCGGCACTCAGATGGCGGCAAACCCCTCCTCCAGGGAAGCGATGATTGCCTTGGACTc                                                   >  5:126726/1‑135 (MQ=38)
                                                                                 gtaaCAC‑GGGCCCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCAt                                                 <  7:20237/134‑1 (MQ=255)
                                                                                       cggggCGCCCCGCGTGCCACCCCCGGCGGGGCGTCCTGCGGGGGCGCGCCCGCGCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATc                                            <  3:136393/134‑1 (MQ=255)
                                                                                         gggCGCAGCCGGTGCCACCCACCGCGGTGCCTCCTGCGCGGGGGCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCa                                          <  7:22046/136‑1 (MQ=255)
                                                                                                  ccGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCtt                                 <  6:57663/136‑1 (MQ=255)
                                                                                                    gTGCCAACCACCGCGGTGCCTGCTTCGCGCGCTCGCCCGCTCCCACAAGGGAACTGCACCGTCAGTCAGATCGCGGCCAACCGCTCAACCAGTTAAGCGATGATTTCCTTGGACTCATACATCCAGGTCACAtttc                               >  3:350777/1‑134 (MQ=255)
                                                                                                     tGCCAACCACCGCGGTGCCGGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCTTGcc                              <  7:203378/136‑1 (MQ=255)
                                                                                                        cAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCTTGCCGGc                           >  3:361036/1‑136 (MQ=255)
                                                                                                        cAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAACGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGCTGATTTCCTTGGACTCATACATCCAGGTCCACTTGCAGGc                           >  3:302467/1‑136 (MQ=255)
                                                                                                        cAACCACCCCGGTGCCTGCTTCGCGGGCTCGCCCCCTCCCACAAGCGAACTGCAGCGTCAGTCAGATGGCGGCAAACAGCTCATCAAGGAAAGCGAGGATTGCCTTGGACTCATAAATCCAGGGCACATTGCCGGc                           >  4:50498/1‑136 (MQ=255)
                                                                                                                 cgGTGCCTGGTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCTTGCCGGCTTCTTCAAt                  <  6:275758/136‑1 (MQ=255)
                                                                                                                 cgGTGCCTGCTCCGCGGGCGCGCCCGCCCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCTTGCCGGCTTCTTCAAt                  <  3:124379/136‑1 (MQ=255)
                                                                                                                 cgGGGCCTGCTTGGCGGGCGCGCCCTCCCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCTTGCCGGCTTCTTCAAt                  <  1:140507/136‑1 (MQ=255)
                                                                                                                  ggTGCCTGCTTCGCTGGCGCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCTTGCCGGCTTCTTCAATg                 <  5:324431/136‑1 (MQ=255)
                                                                                                                  ggTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCTTGCCGGCTTCTTCAATg                 >  5:122189/1‑136 (MQ=255)
                                                                                                                  ggTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCTTGCCGGCTTCTTCAATg                 >  3:292423/1‑136 (MQ=255)
                                                                                                                  ggTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCTTGCCGGCTTCTTCAATg                 <  1:214365/136‑1 (MQ=255)
                                                                                                                    tGCCTGCTTCGCGGGCTCGCCCGCTCCCATAAGGGATCTGCAGCGGCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCTTGCCGGCTTCTTCAATgcg               >  1:299778/1‑136 (MQ=255)
                                                                                                                              gcgGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCAGCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCTTGCCGGCTTCTTCAATGCGCAGGCATGGc     >  1:403824/1‑136 (MQ=255)
                                                                                                                                 ggCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCTTGCCGGCTTCTTCAATGCGCAGGCATGGCACt  >  3:222724/1‑136 (MQ=255)
                                                                                                                                    |                                                                                                                                    
CGCGGCATGTGCCGATGGTGGGCTGAAAACCGACCTGTTTTCCAGGCGCTTGCGCAATCCTTGTGGGAGCGGGCGTGCCCGCGAACAC‑GGGCGCAGCCCGTGCCAACCACCGCGGTGCCTGCTTCGCGGGCTCGCCCGCTCCCACAAGGGATCTGCAGCGTCAGTCAGATGGCGGCAAACCGCTCATCCAGGTAAGCGATGATTGCCTTGGACTCATACATCCAGGTCACCTTGCCGGCTTCTTCAATGCGCAGGCATGGCACT  >  ACB143/2041793‑2042056

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: ‑