Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 3,977,481 A→C 44.5% T284P (ACC→CCC)  ilvE → Branched‑chain‑amino‑acid aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,977,4810AC44.5% 10.6 / 3.0 30T284P (ACC→CCC) ilvEBranched‑chain‑amino‑acid aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (11/4);  new base C (11/3);  total (23/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGT  >  ACB122/3977352‑3977614
                                                                                                                                 |                                                                                                                                     
tCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGCCCGTGATCTAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGGCATCGGCCCcgcg                                                                                                                                  >  5:323929/1‑135 (MQ=255)
     aCTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTGCGTCGAAGCCCGTGCCTCCCGCCCCGCGGGCga                                                                                                                             >  5:64549/1‑134 (MQ=255)
               aTCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCTCCGCGGCCGTTCTCAcccc                                                                                                                    >  7:379638/1‑134 (MQ=255)
                 cACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCCCCCCgct                                                                                                                 >  5:87036/1‑133 (MQ=255)
                 cACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGGTATCCCCCCGAt                                                                                                                 >  4:43673/1‑135 (MQ=255)
                 cACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCCGTGCCGGCGGCCCCGCGGCCGTGTTGACCCCCAt                                                                                                                 >  6:340518/1‑135 (MQ=255)
                    ccGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGGTATAAccccccgcgg                                                                                                              >  5:396531/1‑129 (MQ=255)
                          ctcGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCAt                                                                                                        >  3:134953/1‑135 (MQ=255)
                              cTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCGGCGGCACCGCGGCCGGTTTGACCCCCAACggcggccgcga                                                                                                     >  3:216685/1‑129 (MQ=255)
                                            gcAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATAACCCCGCTCGGCGGCATGGCGTACCACGGCCAGCt                                                                                      >  1:121926/1‑135 (MQ=255)
                                                  gcgcCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCCCCGCGGGCGTTATAACCCCGCTCGGGGGGATCGCGTACCACGGCCAGCTGCAccg                                                                                >  2:326410/1‑133 (MQ=255)
                                                  gcgcCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCCCCGCGGCCGTTATCACCCCGCTCGGCGGCCTCGCGTACAACGGCAAGCTGcccctg                                                                               >  8:232262/1‑131 (MQ=255)
                                                                 cGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCCGTGCCTGCGGCCCCCCGGCCGGTATAACCCCGCTCGGCGGCAAAGAGTAAAAAAGCCAGCTGCAAGTGTTACACGCACTgg                                                                  >  3:116136/1‑134 (MQ=255)
                                                                      tCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAgg                                                            >  3:183843/1‑135 (MQ=255)
                                                                      tCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACCGCCAGCTGCACGTGTTCCACCAACTGGAAAAgg                                                            >  2:19251/1‑135 (MQ=255)
                                                                      tCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCCGTTCCGGCGGCACCGCCGCCCGGATAACCCCCCTAGGCGGCCTAGACGAAAACCGCCACCAGCCCCTGTTGCACCACCCGGAAAAggg                                                           >  7:166078/1‑135 (MQ=255)
                                                                        gAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTTCCTGCGGCCCCGCGGCCGGTATCACCCCGCTAGGGGGCATCCAGTACAACGGCCAGCAGCCCCGGTTACACGCCCTGGCAAAggggg                                                         >  7:160513/1‑133 (MQ=255)
                                                                             cGATGTCGAGATCAGCATGCTGCACCGCTTCGTCGAAGCCGGTGCCTGGGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGccc                                                     <  3:291659/135‑1 (MQ=255)
                                                                                  tCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCCCCGCGGCCGTTATCACCCCGCTCGGCGGCATCCAGGACAACCGCCAGCGGCACGTGTTCCAAGCCCGGGAAAAGGTAGGCCCGGGTa                                                >  3:122397/1‑135 (MQ=255)
                                                                                  tCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTAAAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTAc                                               >  7:293760/1‑136 (MQ=255)
                                                                                   cGATATAAGCAAGCGGGACCGCTTCGTCCAAGCGGGGGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTAcc                                              <  8:49845/136‑1 (MQ=255)
                                                                                   cGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCCCCGCGGCCGGTTTAACCCCGATCCGCGGCATCGCGTACAACGGCAACCTGCACCGGTTCCCCGCCCGGGAAAAGGTCGGCCCGGTTTcc                                              >  8:231789/1‑136 (MQ=255)
                                                                                         ccgCAAGCTGCCCGGGTGCGGCGACGCGGGTGCGGGCGCCCCCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGaa                                         <  1:19251/133‑1 (MQ=255)
                                                                                               gCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCCCCGCGGCCGTTTTCCCCCCGCTCGGCGGCCTCGAGTACAACAGCAAGCTGCCCGTGTTCCACGACCGGGAAAAGGTCGGCCCGGTTACCCAGAAGctcg                                    >  3:261580/1‑133 (MQ=255)
                                                                                                   gACCGCTTCGTCGAAGCCGGTGCCTGCAGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGGTACCCAGAAGCTCTACaa                                >  3:202477/1‑134 (MQ=255)
                                                                                                       gCTTCGTCGAAGCCGGTGCCTGCGGCCCCGCGGCCGGTTTCACCCCGCTCGGCGGGCTCGAGGTCAACGGCCAGCTGCACGTGTTGCACCCCCCGGAAAAGGTCGGCCCGGGTACCCAGAACCTCGACAACGAGc                           >  2:29983/1‑135 (MQ=255)
                                                                                                       gCTTCGGCGACGCCGGTGCCGGGGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGc                           <  4:169245/135‑1 (MQ=255)
                                                                                                       gCTGCGCCGAAGCGGGGGCGGGCGGCCCCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGc                           <  3:67213/135‑1 (MQ=255)
                                                                                                          gcggcgACGCGGGGGCCGGCGCCCCCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGa                        <  3:15611/131‑1 (MQ=255)
                                                                                                                                cACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACAGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGt  <  4:261580/135‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                     
TCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGTTCCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGT  >  ACB122/3977352‑3977614

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: