Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,903,353 C→T 10.2% G172D (GGT→GAT)  dgoT_1 ← D‑galactonate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,903,3530CT10.2% 116.5 / 4.2 39G172D (GGT→GAT) dgoT_1D‑galactonate transporter
Reads supporting (aligned to +/- strand):  ref base C (19/16);  new base T (2/2);  total (21/18)
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.10e-01

GCTTACCTTCTTGCAGTTCAGCGGGTCGCGGTACAGCACCAGCCACACGACGCCGAACAGGATGCCGACCAGACCGGTGCTGTAGAACACGTGGCGCCAGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGC  >  ACB122/1903223‑1903471
                                                                                                                                  |                                                                                                                      
gCTTACCTTCTTGCAGTTCAGCGGGTCGCGGTACAGCACCAGCCACACGACGCCGAACAGGATGCCGACCAGACCGGTGCTGTAGAACACGTGGCGCCAGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGt                                                                                                                    <  1:172368/135‑1 (MQ=255)
         cttGCCGTGCCGCGGGGCGCGGTACAGCACCAGCCACACGACGCCGAACAGGATGCCGACCAGACCGGTGCTGTAGAACACGTGGCGCCAGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAgcg                                                                                                            <  5:1220/134‑1 (MQ=255)
            gCAGTTCAGCGGGTCGCGGTACCCCACCAGCCACACGACGCCGAACAGGATGCCGACCAGACCGGTGCTGTAGAACACGTGGCGCCAGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTg                                                                                                       <  3:110878/136‑1 (MQ=255)
             cAGTTCAGCGGGTCGCGGTACAGCACCAGCCACACGACGCCGAACAGGATGCCGACCAGACCGGTGCTGTAGAACACGTGGCGCCAGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGc                                                                                                      >  3:195521/1‑136 (MQ=255)
             cAGTTCAGCGGGTCGCGGTACAGCACCAGCCACACGACGCCGAACAGGATGCCGACCAGACCGGTGCTGTAGAACACGTGGCGCCAGTCATAGGGGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTg                                                                                                       <  4:83629/135‑1 (MQ=255)
                  ccggggggCGCGGTACAGCACCAGCCACACGACGCCGAACAGGATGCCGACCAGACCGGTGCTGTAGAACACGTGGCGCCAGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAg                                                                                                  <  1:139252/131‑1 (MQ=255)
                  cAGCGGGTCGCGGTACAGCACCAGCCACACGACGCCGAACAGGATGCCGACCAGACCGGTGCTGTAGAACACGTGGCGCCAGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAg                                                                                                  <  3:44369/135‑1 (MQ=255)
                        gTCGCGGTACAGCACCAGCCACACGACGCCGAACAGGATGCCGACCAGACCGGTGCTGTAGAACACGTGGCGCCAGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTgg                                                                                           >  3:111753/1‑136 (MQ=255)
                          cgcgGTACACCACCAGCCACACGACGCCGAACAGGATGCCGACCAGACCGGTGCTGTAGAACACGTGGCGCCAGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGcc                                                                                         <  6:42668/136‑1 (MQ=255)
                                    cACCAGCCACACGACGCCGAACAGGATGCCGACCAGACCGGTGCTGTAGAACACGTGGCGCCAGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAg                                                                               <  5:27313/136‑1 (MQ=255)
                                         gCCACACGACGCCGAACAGGATGCCGACCAGACCGGTGCTGTAGAACACGTGGCGCCAGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGAtgct                                                                          >  4:27457/1‑136 (MQ=255)
                                            acacGACGCCGAACAGGATGCCGACCAGACCGGTGCTGTAGAACACGTGGCGCCAGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGTGACCAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATgctgct                                                                       >  7:172377/1‑136 (MQ=255)
                                            acacGACGCCGAACAGGATGCCGACCAGACCGGTGCTGTAGAACACGTGGCGCCAGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATgctgct                                                                       >  2:25674/1‑136 (MQ=255)
                                                 acgCCGAACAGGATGCCGACCAGACCGGTGCTGTAGAACACGTGGCGCCAGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGt                                                                                 >  3:151831/1‑121 (MQ=255)
                                                 acgCCGAACAGGATGCCGACCAGACCGGTGCTGTAGAACACGTGGCGCCAGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGt                                                                                 <  4:151831/121‑1 (MQ=255)
                                                 acgCCGAACAGGATGCCGACCAGACCGGTGCTGTAGAACACGTGGCGCCAGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCa                                                                  >  4:141306/1‑136 (MQ=255)
                                                          aGGATGCCGACCAGCCCGGTGCTGTAGAACACGGGGCGCCAGTCATAGGTGGGCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCt                                                         <  5:6637/136‑1 (MQ=255)
                                                             aTGCCGACCAGACCGGTGCTGTAGAACACGTGGCGCCAGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTcgc                                                      <  4:111753/136‑1 (MQ=255)
                                                                     cAGACCGGTGCTGTAGAACACGTGGCGCCAGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCgggg                                              <  2:111755/136‑1 (MQ=255)
                                                                     cAGACCGGTGCTGTAGAACACGTGGCGCCAGTCATAGGTGGGCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCgggg                                              <  4:203286/136‑1 (MQ=255)
                                                                           ggTGCTGTAGAACACGTGGCGCCAGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCgggg                                              <  6:101585/130‑1 (MQ=255)
                                                                           ggTGCTGTAGAACACGTGGCGCCAGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCgggg                                              >  5:101585/1‑130 (MQ=255)
                                                                                         cGTGGCGCCAGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACg                          >  7:53238/1‑136 (MQ=255)
                                                                                          gTGGCGCCAGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGCACCATACCGTAACCGCACgg                         >  4:15864/1‑136 (MQ=255)
                                                                                          gTGGCGCCAGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACgg                         >  2:65372/1‑136 (MQ=255)
                                                                                               gcCAGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGtt                    >  4:15453/1‑136 (MQ=255)
                                                                                               gcCAGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCCCGCGGGGAACCATACAGTAACCGCACGGCTGtt                    >  5:121035/1‑136 (MQ=255)
                                                                                                                gccagccaCAGCAGCGCATCGGTGAACAGCGCGGTGCCCAGGTACTGGCcac                                                                                       >  3:84029/1‑52 (MQ=255)
                                                                                                                gccagccaCAGCAGCGCATCGGTGAACAGCGCGGTGCCCAGGTACTGGCcac                                                                                       <  4:84029/52‑1 (MQ=255)
                                                                                                                gccagccaCAGCAGCGCATCGGTGAACAGCGCGGTGCCCAGGTACTGGCcac                                                                                       <  4:80719/52‑1 (MQ=255)
                                                                                                                gccagccaCAGCAGCGCATCGGTGAACAGCGCGGTGCCCAGGTACTGGCcac                                                                                       >  3:80719/1‑52 (MQ=255)
                                                                                                                 ccagccaCAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATAc                                      >  5:3918/1‑100 (MQ=255)
                                                                                                                 ccagccaCAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATAc                                      <  6:3918/100‑1 (MQ=255)
                                                                                                                 ccagccaCAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGc  >  2:47860/1‑136 (MQ=255)
                                                                                                                 ccagccaCAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGc  >  4:93882/1‑136 (MQ=255)
                                                                                                                 ccagccaCAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGc  >  6:130456/1‑136 (MQ=255)
                                                                                                                 ccagccaCAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGc  >  8:186139/1‑136 (MQ=255)
                                                                                                                            agcGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGcc                                                                                         <  5:50137/38‑1 (MQ=255)
                                                                                                                            agcGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGcc                                                                                         >  6:50137/1‑38 (MQ=255)
                                                                                                                                  |                                                                                                                      
GCTTACCTTCTTGCAGTTCAGCGGGTCGCGGTACAGCACCAGCCACACGACGCCGAACAGGATGCCGACCAGACCGGTGCTGTAGAACACGTGGCGCCAGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGC  >  ACB122/1903223‑1903471

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