Predicted mutation
evidence seq id position mutation annotation gene description
RA ACB143 1,669,326 A→G S25G (AGC→GGC)  thrC_2 → Threonine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1431,669,3260AG97.6% 143.9 / ‑4.0 41S25G (AGC→GGC) thrC_2Threonine synthase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (30/10);  total (31/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

TGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGCCCTACCACGAGCTGGCCTTCCGCGTGATGCGCCCGTTCGTCGAAGGCAGCATCGCC  >  ACB143/1669202‑1669457
                                                                                                                            |                                                                                                                                   
tGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGCCggcggc                                                                                                                          >  3:364135/1‑136 (MQ=255)
  ccGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCt                                                                                                                        >  6:310768/1‑136 (MQ=255)
    gACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGAAGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTg                                                                                                                       >  5:119375/1‑135 (MQ=255)
                          ggcgccAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTAGCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACg                                                                                                <  1:228645/136‑1 (MQ=255)
                                   aTCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCAcc                                                                                         >  7:138256/1‑134 (MQ=255)
                                        aGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGgaag                                                                                  <  4:375260/136‑1 (MQ=255)
                                                aCCCATGCGCTATATCAGTACCCGCGGCCAGGCACAGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGc                                                                           >  6:124469/1‑135 (MQ=255)
                                                    aTGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGt                                                                      <  8:73915/136‑1 (MQ=255)
                                                        gcTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGc                                                                  >  8:33323/1‑136 (MQ=255)
                                                        gcTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGc                                                                  >  8:296786/1‑136 (MQ=255)
                                                         cTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGcc                                                                 <  7:296786/136‑1 (MQ=255)
                                                              tCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCtt                                                            <  1:366290/136‑1 (MQ=255)
                                                              tCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCtt                                                            <  5:345590/136‑1 (MQ=255)
                                                               cAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTg                                                           >  4:204180/1‑136 (MQ=255)
                                                               cAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTg                                                           >  7:98698/1‑136 (MQ=255)
                                                                  tACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGccc                                                        >  5:45473/1‑136 (MQ=255)
                                                                  tACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGccc                                                        >  2:312019/1‑136 (MQ=255)
                                                                    cccGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGCCCTa                                                      >  6:113419/1‑136 (MQ=255)
                                                                    cccGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGCCCTa                                                      >  6:335322/1‑136 (MQ=255)
                                                                     ccGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGACACCCTGCCACGATTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGCCCTa                                                      >  2:129388/1‑135 (MQ=255)
                                                                     ccGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGCCCTAc                                                     >  8:276397/1‑136 (MQ=255)
                                                                     ccGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGCCCTAc                                                     >  2:326716/1‑136 (MQ=255)
                                                                     ccGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGCCCTAc                                                     >  1:129716/1‑136 (MQ=255)
                                                                     ccGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGCCCTAc                                                     >  3:144259/1‑136 (MQ=255)
                                                                      cgcgGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGCCCTAcc                                                    <  3:372541/136‑1 (MQ=255)
                                                                       gctgCCAGGCACCGGCCCTGATTTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACATGCAACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGCCCTAcc                                                    >  7:111110/4‑135 (MQ=255)
                                                                       gcgGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGCCCTACCa                                                   >  1:178991/1‑136 (MQ=255)
                                                                       gcgGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGCCCTACCa                                                   >  7:345320/1‑136 (MQ=255)
                                                                         ggCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGCCCTACCAc                                                  >  3:143229/1‑135 (MQ=255)
                                                                          gCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGCCCTCCCAc                                                  >  2:298368/1‑134 (MQ=255)
                                                                           ccAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGCCCTCCCACGAg                                               >  8:240321/1‑136 (MQ=255)
                                                                           ccAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGCCCTACCACGAg                                               >  8:74911/1‑136 (MQ=255)
                                                                             aGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGCCCTACCACGAGCt                                             <  2:75359/136‑1 (MQ=255)
                                                                               gCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGCCCTACCACGAGCTgg                                           <  6:16174/136‑1 (MQ=255)
                                                                               gCACCGGCCCTGAATTTCGAAGACGTGCGGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGCCCTACCACGAGCTgg                                           <  6:20484/136‑1 (MQ=255)
                                                                                 aCCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGCCCTACCACGAGCTGGcc                                         >  5:178659/1‑136 (MQ=255)
                                                                                      cccTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGCCCTACCACGAGCTGGCCTTCcg                                    >  2:236746/1‑136 (MQ=255)
                                                                                        cTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGCCCTACCACGAGCTGGCCTTCcgcg                                  >  5:161792/1‑136 (MQ=255)
                                                                                        cTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGCCCTACCACGAGCTGGCCTTCcgcg                                  >  4:50465/1‑136 (MQ=255)
                                                                                                        tgctgcTGGCTGGCCTGGCCGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGCCCTACCACGAGCTGGCCTTCCGCGTGATGCGCCCGTtcgt                  >  5:81162/1‑136 (MQ=255)
                                                                                                                        ggccGGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGCCCTACCACGAGCTGGCCTTCCGCGTGATGCGCCCGTTCGTCGAAGGCAGCATcgcc  >  6:310878/1‑136 (MQ=255)
                                                                                                                            |                                                                                                                                   
TGCCGACGGGCCGCCCGCGCGGCCCCGGCCCCAGCATCGAAGGTTTGCACCCATGCGCTATATCAGTACCCGCGGCCAGGCACCGGCCCTGAATTTCGAAGACGTGCTGCTGGCTGGCCTGGCCAGCGACGGCGGCCTGTACGTGCCGGAAAACCTGCCACGCTTCACCCAGGAAGAAATCGCCTCGTGGGCCGGCTTGCCCTACCACGAGCTGGCCTTCCGCGTGATGCGCCCGTTCGTCGAAGGCAGCATCGCC  >  ACB143/1669202‑1669457

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