Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,670,690 A→C 29.0% intergenic (+27/‑104) thrC_2 → / → CJ019_01569 Threonine synthase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,670,6900AC29.0% 44.2 / 3.3 50intergenic (+27/‑104)thrC_2/CJ019_01569Threonine synthase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/4);  new base C (9/7);  total (37/13)
Fisher's exact test for biased strand distribution p-value = 2.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGATTCA  >  ACB122/1670559‑1670819
                                                                                                                                     |                                                                                                                                 
aGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTGATTGCATGGCGTTTGGGGCTGCAGCGCAGcc                                                                                                                                >  2:176413/1‑135 (MQ=255)
                   aaGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCACAGCCGTTTc                                                                                                             >  4:156369/1‑135 (MQ=255)
                    agagCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGGTGCCGCGCCGCCCCCACGTCTCAGCCgtttgc                                                                                                            >  5:248838/1‑133 (MQ=255)
                      agCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCCGCCGTTTCcct                                                                                                          >  2:45401/1‑135 (MQ=255)
                      agCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCcct                                                                                                          >  4:239533/1‑135 (MQ=255)
                           tGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTTGGGCTGCAGCGCAGCCCCCACGTCGCAGCCGTTTAcctcct                                                                                                       >  3:77573/1‑133 (MQ=255)
                              aCGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCAGCCCCCACGTCTCACCCGTTTACCTCCGGTCt                                                                                                   >  6:12111/1‑134 (MQ=255)
                              aCGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCt                                                                                                   >  7:160731/1‑134 (MQ=255)
                              aCGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTa                                                                                                  >  6:179613/1‑135 (MQ=255)
                                 gTCCTGGCCACCTCCCTGATGGGTGTACTGGGTTTTGTCTGCCCGCGCGGAAACGGGGGAAAGCCGGTGTAATGCATTGGGTTTGGGGGGTGAGGGGCTGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCa                                                                                               <  6:99490/135‑1 (MQ=255)
                                     tGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCCAGCCGCTGTAATTGCATGGCGTTTTGGGCTGCAGAGCAGCCCCAACGTGTCAGCCTAGTCACTCCTGTCTATAAGaa                                                                                            >  6:86141/1‑134 (MQ=255)
                                     tGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCAGCCGCCCACCCCCCACGAATCCCCCTTTTCCCACCTGTCTATCACaa                                                                                            >  2:392221/1‑134 (MQ=255)
                                          aaCTACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCCACGGCGCCGCCGTTTGCCGCCTGGCGATGATCAACATc                                                                                       >  3:39456/1‑134 (MQ=255)
                                          aaCGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCACGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCCGCCCCCACGTCGCAGCCTTTTCCCTCCGGTCTATCAGAAACCTc                                                                                       >  3:320478/1‑134 (MQ=255)
                                          aaCGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCCGCCCCCAAGTCTCCGCCGTTTCCCTCCTGTCTATCAGAAACATc                                                                                       >  5:129021/1‑134 (MQ=255)
                                               ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGGAGCGCCGCCCCCACCGCTCAGCCTTTTGCCTCCAGTCTATGAGGAACATCCGATTc                                                                                 >  1:197760/1‑135 (MQ=255)
                                               ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCAGCCCCAACGTCTCCCCCGTTTACCTCCTGTCTATCAGCAACCTCCTAtta                                                                                 >  4:464259/1‑134 (MQ=255)
                                               ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCAGTCGATCAGAAACATCCTATTc                                                                                 >  4:159996/1‑135 (MQ=255)
                                               ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCAGCCCCAACGTCTCAGCCGTTTCCCTCCTGTCTATCAGCAACATCCTAt                                                                                   >  1:79835/1‑133 (MQ=255)
                                               ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCCGCCCCCACGTCTCAGCCTTTTCCCTCCTGTCTATCAGCAAAATCCTATTc                                                                                 >  1:367192/1‑135 (MQ=255)
                                               ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCCGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCGATCAGAAACATCCTAt                                                                                   >  2:100947/1‑133 (MQ=255)
                                               ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCCCCCCCCACCTCTCAGCCTATTCGCTCCTGTCAATCAGAAACATGCGAtta                                                                                 >  6:419643/1‑134 (MQ=255)
                                               ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCCACGTCTCAGCCATTTCCCTCCTGTCTATCAGAAACATCCTAt                                                                                   >  4:325475/1‑133 (MQ=255)
                                               ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTAt                                                                                   >  8:214134/1‑133 (MQ=255)
                                               ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCGATCAGAAACATCCTATTc                                                                                 >  8:51717/1‑135 (MQ=255)
                                               ccTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACCTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCGAtta                                                                                 >  3:433282/1‑134 (MQ=255)
                                               ccTGAAGGCTGTACAGGGCTTTGACAGCCACCCCGCCACCCGCGCCACGCCGCTGTAATTGAATGGCTTTGGGGCCTGCCGCGAAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTc                                                                                 <  4:458867/135‑1 (MQ=255)
                                                     ggCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTg                                                                           >  8:446043/1‑135 (MQ=255)
                                                         gtccAGGGCTTGGCCACCCCCCCCGCCAACCGCGCCAAGCCGCTGAATACCCTGGGCGTGGGGGGCCGCAGCCCCCCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCgg                                                                       <  3:395654/132‑1 (MQ=255)
                                                           cccgggCTTTGCCAGCCACCCCTGCACCCGCGCCACCCCGCTTTAATGGCGGGGGGTTGGGGGCCCCCGCCCCGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGgc                                                                     <  6:6357/132‑1 (MQ=255)
                                                                gCTTTGTCCGCCAGCACGGCAACCGCGGCACGCCGCTGTAATTGCAGGGGGTTGGGGGCGGCCGCGCCGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCa                                                                 <  4:43932/134‑1 (MQ=255)
                                                                gCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCCCGCCACCCCAACGTCTCACCCTTTTCCCACCTGGCTATCAGAACCATCCTATTCAGCCTGCCGGGCGCCaa                                                                >  2:389573/1‑135 (MQ=255)
                                                                  tttGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCCGCCCCCCCGGGGCCGCCGTTTGCCGGCGGTGTTTTAGAAACATCCGATTAAGCCGGCCCGgcggccac                                                               >  6:41450/1‑129 (MQ=255)
                                                                  tttGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCCGCGCAGCCCCAACGGCGCCGCCTTGTACCGCCGGGATATCAGAAACACACGATTCAGCCTGCCGGGCGcccac                                                               >  7:242882/1‑131 (MQ=255)
                                                                  tttGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCATCGCAGCCCCAACGTCTCAGCCGTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGc                                                              >  5:74662/1‑135 (MQ=255)
                                                                  tttGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGGCCCCACGTCTCAGCCGTTTCCCTCCTGTCTATGAGAAACATGCTATTCAGCCTGCCGGGCGCCAAGc                                                              >  3:409597/1‑135 (MQ=255)
                                                                            agcACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGc                                                    >  8:342958/1‑135 (MQ=255)
                                                                              cACGGCAACCGCGGCAAGCCGCTGAAATTGCTGGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCg                                                  <  3:340477/135‑1 (MQ=255)
                                                                              cACGCCCACCGCGGCAAGCCGCTGTATTGGCCTGGGGGTGGGGGCGGCCGCGCAGCCCCCACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCg                                                  <  1:162351/135‑1 (MQ=255)
                                                                                ccccacccGCGCCACGCCGGTTTATTTGCTGGGGGTTGGGGGGTGCCGCCCCCCCCCCACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGt                                                <  1:146061/129‑1 (MQ=255)
                                                                                     aCCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGtcctc                                            >  4:394237/1‑134 (MQ=255)
                                                                                            catgCCGCTGTTTTTGCAGGTGGTTTGGGGCGGCTGCGCTGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAgg                                    <  4:73861/132‑1 (MQ=255)
                                                                                                  gCTGTAATTGCAGGGGGTTGGGGGCTCCAGCGCCCCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGc                              <  7:238345/135‑1 (MQ=255)
                                                                                                       aaTTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCGATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATGAGGAGCGGCTCAt                          >  7:246564/1‑134 (MQ=255)
                                                                                                               gggggTTGGGGGCTGCCGCGCCGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGc                  <  8:327856/131‑1 (MQ=255)
                                                                                                                  cGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCCGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCCTGCACCCGcccg               >  5:212258/1‑134 (MQ=255)
                                                                                                                  cGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCt               >  4:337760/1‑135 (MQ=255)
                                                                                                                    tttGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTAt             >  3:158247/1‑135 (MQ=255)
                                                                                                                         gggCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATgtgt         <  4:158247/135‑1 (MQ=255)
                                                                                                                         gggCGGCCGCCCCGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATgtgt         <  7:263823/135‑1 (MQ=255)
                                                                                                                                agcgccgcCCCAACGTCTCAGCCTTTTACCTACTGTGTATAATAAACATCATATTAAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCCTCAGGAGCGGCTCATGCACCAGCCCGATGTGTTTATTCa  >  8:28424/1‑135 (MQ=255)
                                                                                                                                agcgcagcCCCAACGTCTCAGCCGTTTACCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCAAAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGCTTCa  >  8:304482/1‑135 (MQ=255)
                                                                                                                                agcccagcccCAACGTCTCAGCCATTTACCAACTGTCTATAAGAAACAAACAATTCAGCCTGCCGGGCGCCCAGCATAAAATATCCGGTCCTCCCATCACGAGCGGCTCATTCACCACCCCAAAGTGTTGATTCa  >  2:27004/1‑135 (MQ=255)
                                                                                                                                     |                                                                                                                                 
AGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGCCAGCACGGCAACCGCGGCAAGCCGCTGTAATTGCATGGCGTTTGGGGCTGCAGCGCAGCCCCAACGTCTCAGCCTTTTCCCTCCTGTCTATCAGAAACATCCTATTCAGCCTGCCGGGCGCCAAGCATAAATTAGCCGGTCCTCCCATCAGGAGCGGCTCATGCACCAGCCCTATGTGTTGATTCA  >  ACB122/1670559‑1670819

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: