Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 3,047,252 T→G 13.7% S111R (AGC→CGC)  CJ019_02862 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,047,2520TG13.7% 91.8 / 2.4 43S111R (AGC→CGC) CJ019_02862hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (12/24);  new base G (5/2);  total (17/26)
Fisher's exact test for biased strand distribution p-value = 9.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTCGCGCGTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGTGGCTTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGCGGCTGGTCTTCTATGAAGCGACCTTGCAGCCATTGACGGCGATAGACGTCGATCAGCACCGGCCGGGGCGT  >  ACB122/3047123‑3047378
                                                                                                                                 |                                                                                                                              
cTTTCGCGCGTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGTGGGGTTTTGGAGGGTTCGCACGAAGTAGTCCTGCCATGCCGGGCGCAGtt                                                                                                                           >  3:77895/1‑135 (MQ=255)
cTTTCGCGCGTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGTGGGGTGTTTGAGGTTTAGGCCGAAGTAGTCCTGGCATGCCGGGCGCAggt                                                                                                                           >  2:208888/1‑133 (MQ=255)
 tttCGCGCGTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGTGGCTTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTtc                                                                                                                          <  3:35796/135‑1 (MQ=255)
 tttCGCGCGTTTTACCGCGGTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGTGGCTTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTtc                                                                                                                          <  3:34350/135‑1 (MQ=255)
      cgcgTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGTGGGTTTTTTCCGGTTACGCACGATGTATTCCTGCCATTACGGGCGAAGTtagccgg                                                                                                                     >  4:13166/1‑134 (MQ=255)
      cgcgTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGTGGCTTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTtcgccgt                                                                                                                     <  1:201649/135‑1 (MQ=255)
      cgcgTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGTGGCGTTTTTGAGGTTACGCACGAAGTAGTCCTGCCATTCCGGGCTCAGTtcgccgt                                                                                                                     >  6:85401/1‑135 (MQ=255)
                cgcgTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGTGGCTTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGtt                                                                                                           <  7:248580/135‑1 (MQ=255)
                cgcgTGCCAGCTTAACCAGCAAGCGGCGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGTGGCTTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGtt                                                                                                           <  1:29487/135‑1 (MQ=255)
                       cAGCTTAACCAGCACACGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGTGGCTTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCGCAGTTCGCCGTCGCCGGTGTTGgccag                                                                                                     <  8:176873/134‑1 (MQ=255)
                        aGCTTAACCAGCAATCGGATCGGCAGTAAAAGAGCAGTGGTTGGTCCAGGTGCCAGGCCGTGGCTTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGgccagcc                                                                                                   <  6:243418/135‑1 (MQ=255)
                                  gCACTCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGTGGCTTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGcc                                                                                         <  6:132060/135‑1 (MQ=255)
                                       cGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGTGGCTTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGa                                                                                    <  3:292429/135‑1 (MQ=255)
                                                 gTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGTGGGGTTTTTGAGGTTTCGCACGCAGTAGTCCAGCCAAGCCCGGCGCAGTTAGCCGTCGCCGGTGTTTGCCCGCCCCAGGCGGCCAGGGGGGGTgggt                                                                            >  8:146077/1‑130 (MQ=255)
                                                 gTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGTGGCTTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATg                                                                          <  1:15310/135‑1 (MQ=255)
                                                  tacaaGACCAGTGGGTGGTCCAGGTGCCCGGCCGTGGCTTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGc                                                                         <  6:47757/132‑1 (MQ=255)
                                                          cAGTGGTTGGTCCAGGTGCCCGGCCGTGGCTTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCGCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGGTGCTATGCGGCTGGtc                                                                 <  6:255236/135‑1 (MQ=255)
                                                                 tggtCCAGGTGCCCGGCCGTGGCTTTTTGCAGGTTACGCAAGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGCGGCTGGTCTTCTATg                                                          <  6:133225/135‑1 (MQ=255)
                                                                  ggtCCAGGTGCCCGGCCGTGGCTTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGCGGCTGGTCTTCTATGa                                                         <  1:107692/135‑1 (MQ=255)
                                                                            cccGGCCGTGGCTGTTTGGAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCGCAGTTCGCCGGCGCCGGGGTTTGCCCGCCACAGGCTGCCCGGGAGGTTGCTGTTCGGCTGGTCTTCTATGCAGCGCCCATg                                               >  6:39798/1‑135 (MQ=255)
                                                                                   gTGGCTTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGCGGCTGGTCTTCTATGAAGCGACCTTGCAGCCAt                                        <  7:230172/135‑1 (MQ=255)
                                                                                              cAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGCGGCTGGTCTTCTATGAAGCGACCTTGCAGCCATTGACGGCg                                >  2:74029/1‑132 (MQ=255)
                                                                                               aGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGCGGCTGGTCTTCTATGAAGCGACCTTGCAGCCATTGACGGCGATAg                            <  3:104750/135‑1 (MQ=255)
                                                                                               aGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCGGCCAGGGAGGTTGCTATGCGGCTGGTCTTCTATGAAGCGACCTTGCAGCCATTGACGGCGATAg                            <  1:208888/135‑1 (MQ=255)
                                                                                                    aCGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGtt                                                                                                           >  5:244716/1‑51 (MQ=255)
                                                                                                    aCGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGtt                                                                                                           <  6:244716/51‑1 (MQ=255)
                                                                                                     cGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGCGGCTGGTCTTCTATGAAGCGACCTTGCAGCCATTGACGGCGATAGACGTc                       <  1:141222/134‑1 (MQ=255)
                                                                                                     cGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGCGGCTGGTCTTCTATGAAGCGACCTTGCAGCCATTGACGGCGATAGACGTc                       <  4:77895/134‑1 (MQ=255)
                                                                                                              tagtCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGCGGCTGGTCTTCTATGAAGCGACCTTGCAGCCATTGACGGCGATAGACGTCGATCAGCAcc             >  7:74518/1‑135 (MQ=255)
                                                                                                              tagtCCTGCCATGCCGGGCTCAGTTCGAAGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGCGGCTGGTCTTCTATGAAGCGACCTTGCAGCCATTGACGGCGATAGACGTCGATCAGCAcc             >  7:129855/1‑135 (MQ=255)
                                                                                                                  ccTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGCGGCTGGTCTTCTATGAAGCGACCTTGCAGCCATTGACGGCGATAGACGTCGATCAGCAccggcc         <  1:74029/135‑1 (MQ=255)
                                                                                                                  ccTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGCGGCTGGTCTTCTATGAAGCGACCTTGCAGCCATTGACGGCGATAGACGTCGATCAGCAccggcc         <  7:146077/135‑1 (MQ=255)
                                                                                                                  ccTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGCGGCTGGTCTTCTATGAAGCGACCTTGCAGCCATTGACGGCGATAGACGTCGATCAGCAccggcc         <  1:161040/135‑1 (MQ=255)
                                                                                                                         tGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCt                                                                                            >  4:56693/1‑45 (MQ=255)
                                                                                                                         tGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCt                                                                                            <  3:56693/45‑1 (MQ=255)
                                                                                                                         tGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCa                                                                                        <  6:46793/49‑1 (MQ=255)
                                                                                                                         tGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCa                                                                                        >  5:46793/1‑49 (MQ=255)
                                                                                                                         tGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGCGGCTg                                                                    >  2:134638/1‑69 (MQ=255)
                                                                                                                         tGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGCGGCTg                                                                    <  1:134638/69‑1 (MQ=255)
                                                                                                                         tGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGCGGCTGGTCTTCTATGAAGCGACCTTGCAGCCATTGACGGCGATAGACGTCGATCAGCACCGGCCGGGGCGt  >  2:290006/1‑135 (MQ=255)
                                                                                                                         tGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGCGGCTGGTCTTCTATGAAGCGACCTTGCAGCCATTGACGGCGATAGACGTCGATCAGCACCGGCCGGGGCGt  >  2:137798/1‑135 (MQ=255)
                                                                                                                         tGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGCGGCTGGTCTTCTATGAAGCGACCTTGCAGCCATTGACGGCGATAGACGTCGATCAGCACCGGCCGGGGCGt  >  7:8423/1‑135 (MQ=255)
                                                                                                                         tGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGACAAAGTCTGCCAGGGATGTTGCTATCCGGCTGGACTTCGATGAAGCGACCTTGCAGCAATTGACGGCGATAGACGTCCATCAGCACCGGCCGGGGCgg  >  5:199614/1‑134 (MQ=255)
                                                                                                                                 |                                                                                                                              
CTTTCGCGCGTTTTACCGCGTTCCAGCTTAACCAGCAATCGGAGCGGCAGTAAAAGACCAGTGGTTGGTCCAGGTGCCCGGCCGTGGCTTTTTGCAGGTTACGCACGAAGTAGTCCTGCCATGCCGGGCTCAGTTCGCCGTCGCCGGTGTTGGCCAGCCACAGGCTGCCAGGGAGGTTGCTATGCGGCTGGTCTTCTATGAAGCGACCTTGCAGCCATTGACGGCGATAGACGTCGATCAGCACCGGCCGGGGCGT  >  ACB122/3047123‑3047378

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