Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 6,020,340 A→C 39.3% L58R (CTG→CGG)  ycaD_4 ← putative MFS‑type transporter YcaD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1226,020,3400AC39.3% 11.3 / 7.5 41L58R (CTG→CGG) ycaD_4putative MFS‑type transporter YcaD
Reads supporting (aligned to +/- strand):  ref base A (5/14);  new base C (13/3);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 2.05e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CGATCATCAGCCGCAGCAGCAGCCATAACGGGTAACTAGGCAGCAGCGCGAGCAGGCCAATCGACAGCGCCCCGCCCCACAGGCACAGGCGCATCGCCGAAGGTACGCCAACCCAACCTGCCAGGCGGCTGGCCAGGCTGGCGCCGGCCAGTACGCCCAACGCCGGCATTGCCGCCATCACGCCGATGGCGAAGCTGCCGTAGCCCCAGGCCTCCAGGCGCAGGGACACCAGCGGCATGCTCACGCCAAGCGCCAGCCCGACACTGAGC  >  ACB122/6020206‑6020474
                                                                                                                                      |                                                                                                                                       
cGATCAGCCGCCGCAGCAGCCGCCAAACCGGGTACCTAGGCCGCAGCGCGAGCAGGCCAATCGACAGCGCCCCGCCCCACAGGCACAGGCGCATCGCCGAAGGTACGCCAACCCAACCTGCCAGGCGGCTGGCCa                                                                                                                                         <  7:352396/135‑1 (MQ=255)
 gATCATCAGCCGCAGCAGCAGCCATAACGGGTAACTAGGCAGCAGCGCGAGCAGGCCAATCGACAGCGCCCCGCCCCACAGGCACAGGCGCATCGCCGACGGTACGCCAACCCAACCTGCACTGCGGCTGGGCac                                                                                                                                        >  4:150632/1‑134 (MQ=255)
 gATCATCAGCCGCAGCAGCAGCCATAACGGGTAACTAGGCAGCAGCGCGAGCAGGCCAATCGACAGCGCCCCGCCCCACAGGCACAGGCGCATCGCCGAAGGTACGCCAACCCAACCTGCCAGGCGGCAGGCCac                                                                                                                                        >  2:177241/1‑134 (MQ=255)
          ccGCAGCAGCAGCCATAACGGGTAACTAGGCAGCAGCGCGAGCAGGCCAATCGACAGCGCCCCGCCCCACAGGCACAGGCGCATCGCCGAAGGTACGCCAACCCAACCTGCCAGGCGGGTGGGCCGGCTGGCGcc                                                                                                                               >  6:150812/1‑135 (MQ=255)
          ccGCAGCAGCAGCCATAACGGGTAACTAGGCAGCAGCGCGAGCAGGCCAATCGACAGCGCCCCGCCCCACAGGCACAGGCGCATCGCCGAAGGTACGCCAACCCAACCTGCCAGGCGGGTGGCCCCGCTGGCGcc                                                                                                                               >  1:33043/1‑135 (MQ=255)
          ccGCAGCAGCAGCCATAACGGGTAACTAGGCAGCAGCGCGAGCAGGCCAATCGACAGCGCCCCGCCCCACAGGCACAGGCGCATCGCCGAAGGTACGCCAACCCAACCTGCCAGGCGGCTGGCCAGGCTGGCGcc                                                                                                                               <  6:335829/135‑1 (MQ=255)
                    agcCATACGGGGTAACTCGGCAGCCGCGCGAGCAGCCCAATCGACCGCGCCCCGCCCCACCGCCACCGGCGCATCGCCGAAGGTACGCCAACCCAACCTGCCAGGCGGCGGGCCTGGCTGGCGCCGGCCAGTACg                                                                                                                    <  2:338770/135‑1 (MQ=255)
                                   cTAGGCAGCAGCGCGAGCAGGCCAATCGACAGCGCCCCGCCCCACAGGCACAGGCGCATCGCCGAAGGTACGCCAACCCAACCTGCCAGGCGGCTGGCCACGCAGGCGCCGGCCACTAACCCCAACCCCGGCAtt                                                                                                     >  4:124213/1‑135 (MQ=255)
                                     aGGCAGCAGCGCGAGCAGGCCAATCGACAGCGCCCCGCCCCACAGGCACAGGCGCATCGCCGAAGGTACGCCAACCCAACCTGCCAGGCGGCTGGCCACGCTGGCGCCCGCCCGTACGCCCACCGCCGGccatgc                                                                                                   >  2:235289/1‑130 (MQ=255)
                                           cagcGCGAGCAGGCCAATCGACAGCGCCCCGCCCCACAGGCACAGGCGCATCGCCGAAGGTACGCCAACCCAACCTGCCAGGCGGCTGGCCCGGCTGGCGCCGGCCAGTAAGCCCCACGCCGGCATTgccgccc                                                                                              >  8:265831/1‑133 (MQ=255)
                                            agcGCGAGCAGGCCAATCGACAGCGCCCCGCCCCACAGGCACAGGCGCATCGCCGAAGGTACGCCAACCCAACCTGCCAGGCGGCTGGCCACGCGGGCGCCCGCCCCGACGCCCAACACCCGCATTGGCGCCCTc                                                                                            >  3:52407/1‑135 (MQ=255)
                                                 gAGCAGGCCAATCGACAGCGCCCCGCCCCACAGGCACAGGCGCATCGCCGAAGGTACGCCAACCCAACCTGCCAGGCGGCTGGCCCGGCGGGCGGCGGCCAGTACGCCCCCCCCCGGCCTTGCCGCCCTCCCGcc                                                                                       >  7:103718/1‑135 (MQ=255)
                                                      ggCAAATCGACCGCGCCCCGCCCCACCGCCCCCGGCGCAGCGCCGACGGTACGCCAACCCACCCTCCCCGGCGGCGGGCCAGGCTGGCGCCGGCCAGTACGCCCAACGCCGGCATTGCCGCCATCACGCCGATgg                                                                                  <  5:154278/135‑1 (MQ=255)
                                                              gACCGCCCCCCGCCCCACGGGCACCGGCGCATCGCCGAAGTTACCCCAACCCACCCTCCCGGGCGGGGGGCAGGGGGGGGGCCGGCTTGTACGCCCAACGCCGGCATTGCCGCCATCACGCCGATGGCGAAGCTg                                                                          <  2:82832/135‑1 (MQ=255)
                                                               aCAGCGCCCCGCCCCACAGCCACAGGCGCAGCGCCGAAGGTACCCCAACCCACCCTGCCAGCGGGCTGGCCCGGCTGGCGCCGGCCAGTACGCCCAACGCCGGCATTGCCGCCATCACGCCGATGGCGAAGCTGc                                                                         <  7:4613/135‑1 (MQ=255)
                                                                    gccccgccccACAGGCACAGGCGCATCGCCGAAGGTACGCCAACCCAACCTGCCAGGCGGCTGGCCAGGCTGGTGCCGGCCCGTACGCCCCACGCCGGCCTTGCCGCCCTCACGCCGATGGCGAAGCTGCCGTAg                                                                    >  6:326185/1‑135 (MQ=255)
                                                                    gccccgccccACAGGCACAGGCGCATCGCCGAAGGTACGCCAACCCAACCTGCCAGGCGGCGGGCCCGGCGGGCGCCGGCCCGTACGCCCCACGCCCGCCGTGCCGCCCTCCCCCCGCAGGCCAAGCTTCCGTAg                                                                    >  6:360349/1‑135 (MQ=255)
                                                                        cgccccACAGGCACAGGCGCATCGCCGAAGGTACGCCCACCCAACCTGCCAGGCGGCTGGCCCGGCAGGCGCCGGCCAGTACGCCCAACGCCGGCATTGCCGCCATCACGCCGATGGCGAAGCTGCCGTAGcccc                                                                >  4:264952/1‑135 (MQ=255)
                                                                             cccaggcactggcGCAGCGCCCAAGGAACGCCAACCCAACCTGCCCGGCGGGTGGCCTGGCTGGCGCGGGCCAGTACGCCCAACGCCGGCATTGCCGCCATCACGCCGATGGCGAAGCTGCCGTAGCCCCAGGcc                                                           <  1:74419/133‑1 (MQ=255)
                                                                                  gcacaggcGCATCGCCGAAGGTACGCCAACCCAACCTGCCAGGCGGCGGGGCCGGCGGGCGCCCGCCCGGTCGCCCCACGCCCGCAAGGCCCCCAGCCCGCCGCTGGCGAAGCAGACGGCGCCCCCCGCCGCCac                                                      >  8:170687/1‑134 (MQ=255)
                                                                                     caggcGCATCGCCGAAGGTACGCCAACCCAACCTGCCAGGCGGCGGGCCCCGCGGGCGCCGGCCAGGCAGGCCCCCCCCCGCCAGTCCCCCAACCCCCCCCTGGCGCCGCTGGCGTAGCCCCCGGCCTccccgcc                                                   >  2:338079/1‑130 (MQ=255)
                                                                                      aggcGCATCGCCTAAGGTACGCCAACCCAACCTGCCAGGCGGCTGGCCAGGCTGGCGCCGGCCAGTACGCCCAACGCCGGCATTGCCGCCATCACGCCGATGGCGAAGCTGCCGTAGCCCCAGGCCTCCAGgcgc                                                  <  2:73769/135‑1 (MQ=255)
                                                                                        gcgcATCGCCGAAGGTACGCCAACCCAACCTGCCAGGCGGCGGGCCCCGCTGGCGCCGGCCAGGACGCCCCACCCCCGCATTGCCGCCCACACCCCGAAGGCGAAGCTGCCGTAGCCCCCCGCCTCCACGCGCAg                                                >  6:379985/1‑135 (MQ=255)
                                                                                           aaTCGCCGAAGGTACGCCAACCCAACCTGCCCGGCGGGTGGCCAGGCTGGCGCCGGCCAGTACGCCCAACGCCGGCATTGCCGCCATCACGCCGATGGCGAAGCTGCCGTAGCCCCAGGCCTCCAGGCGCAGGGa                                             <  7:26685/134‑1 (MQ=255)
                                                                                               gCCGAAGGTACGCCAACCCAACCTGCCAGGCGGCTGGCCCGGCTGGCGCCGGCCCGTACGCCCCCCGCCGGCATTGCCGCCCTCACGCCGATGGCGAAGCTGCCGTAGCCCCAGGCCGCCCGGCGCAgggccccc                                         >  3:344259/1‑130 (MQ=255)
                                                                                               gCCGAAGGTACGCCAACCCAACCTGCCAGGCGGCAGGCCAGGCTGGCGCCGGGCAGTACGCCCAACGCCGGCATTTCCGCCATCACCCCGCTGGCGCAGCTGCCGTAGCCCCAGGCCTCCAGGCGCACGGCCAcc                                         >  2:78032/1‑135 (MQ=255)
                                                                                                 cGATGGTACGCCACCCCACCCTGCCCGGCGGCTGGCCTGGCTGGCGCCGGCCAGAACGCCCAACGCCGGCATTGCCGCCATCACGCCGATGGCGAAGCTGCCGTAGCCCCAGGCCTCCAGGCGCAGGGACACCAg                                       <  8:178358/135‑1 (MQ=255)
                                                                                                 cGACGGTCCCCCACCCCACCCGGCCCGGCGGCGGGCCCGGGGGGCGCCGCCCAGTACGCCCAACGCCGGCATTGCCGCCATCACGCCGATGGCGAAGCTGCCGTAGCCCCAGGCCTCCAGGCGCAGGGACACCAg                                       <  3:287724/135‑1 (MQ=255)
                                                                                                 cGAAGGTACGCCAACCCACCCGGCCAGGCGGCGGGCCAGGCTGGCGCCGGCCAGTACGCCCAACGCCGGCATTGCCGCCATCACGCCGATGGCGAAGCTGCCGTAGCCCCAGGCCTCCAGGCGCAGGGACACCAg                                       <  3:314990/135‑1 (MQ=255)
                                                                                                  gAAGGTACGCCAACCCACCGTGCCTGGCGGCTGGCCAGGCTGGCGCCGCCCAGTACGCCCAACGCCGGCATTGCCGCCATCACGCCGATGGCGAAGCTGCCGTAGCCCCAGGCCTCCAGGCGCAGGGACACCAGc                                      <  7:234436/135‑1 (MQ=255)
                                                                                                  gAAGGTACGCCAACCCAACCTGCCAGGCGGCTGGCCAGGCTGGCGCCGGCCAGTACGCCCAACGCCGGCATTGCCGCCATCACGCCGATGGCGAAGCTGCCGTAGCCCCAGGCCTCCAGGCGCAGGGACACCAGc                                      <  4:52407/135‑1 (MQ=255)
                                                                                                       aaCGCCACCCCAACCTGCCAGGCGGCTGGCCAGGGGGGCGCCGGCCAGTACGCCCAACGCCGGCATTGCCGCCATCACGCCGATGGCGAAGCTGCCGTAGCCCCAGGCCTCCAGGCGCAGGGACACCAGCGGCAt                                 <  2:344490/134‑1 (MQ=255)
                                                                                                                  aaccTGCCAGGCGGCTGGCCCCGCTGGCGCCGGCCCGTACGCCCCACCCCCGCATTGCCGCCCTCACGCCGATGGCGAAGCTGCCGTAGCCCCCGGCCGCCAGGCGCCGGGCCCCCCGCGGCCTGcgcccgccca                      >  7:215969/1‑133 (MQ=255)
                                                                                                                   cccTGCCCGGGGGGGGGCCAGGCTGGCGCCGGCAAGTACCCCCAACGCCGGCATTGCCGCCATCACGCCGATGGCGAAGCTGCCGTAGCCCCAGGCCTCCAGGCGCAGGGACACCAGCGGCATGCTCACGCCAAg                     <  5:360349/134‑1 (MQ=255)
                                                                                                                   accTGCCCGGCGGCTGGCCAGGCGGGCGCCGGCCAGTACCCCCAACGCCGGCATTGCCGCCATCACGCCGATGGCGAAGCTGCCGTAGCCCCAGGCCTCCAGGCGCAGGGACACCAGCGGCATGCTCACGCCAAg                     <  1:50617/135‑1 (MQ=255)
                                                                                                                   accTGCCAGGGGGCGGGCCAGGCTGGGGCCGGCCAGTACGCCCAACGCCGGCATTGCCGCCATCACGCCGATGGCGAAGCTGCCGTAGCCCCAGGCCTCCAGGCGCAGGGACACCAGCGGCATGCTCACGCCAAg                     <  2:33043/135‑1 (MQ=255)
                                                                                                                   accTGCCAGGCGGCTGGCCAGGCTGGCGCGGGCCAGTACGCCCAACGCCGGCATTGCCGCCATCACGCCGATGGCGAAGCTGCCGTAGCCCCAGGCCTCCAGGCGCAGGGACACCAGCGGCATGCTCACGCCAAg                     <  7:265831/135‑1 (MQ=255)
                                                                                                                        cccggcggcTGGCCAGGCGGGCGCCGGCCAGTACGCCCAACGCCGGCATTGCCGCCATCACGCCGATGGCGAAGCTGCCGTAGCCCCAGGCCTCCAGGCGCAGGGACACCAGCGGCATGCTCACGCCAAGCGCCa                <  3:189825/132‑1 (MQ=255)
                                                                                                                          aGGCGGGTGGCCTGGCGTGGGCCGGCAAGTACGCCCAACGCCGGCATTGCCGCCATCACGCCGATGGCTAAGCTTGCGTAGCCCCAGGCCTCAAGGCGCAGGGTCACCAGCGTCATGCTCACGCCATGCTCCAGc              <  8:103078/135‑1 (MQ=255)
                                                                                                                             gggcTGGCCTGGGTGGCGCCGGCCAGAACGCCCAACGCCGGCATTGCCGCCATCACGCCGATGGCGAAGCTGCCGTAGCCCCAGGCCTCCAGGCGCAGGGACACCAGCGGCATGCTCACGCCAAGCGCCAGCCCg           <  7:345857/134‑1 (MQ=255)
                                                                                                                             cggcggGCCAGGCTGGCGCCGGCCAGTACGCCCAACGCCGGCATTGCCGCCATCACGCCGATGGCGAAGCTGCCGTAGCCCCAGGCCTCCCGGCGCAGGGACACCAGCGGCATGCTCACGCCAAGCGCCAGCCCg           <  7:24795/135‑1 (MQ=255)
                                                                                                                             cggcTGGCCCGGCTGGCGCCGGCCAGGTCGCCCCACGCCCGCATTGCCGCCCTCACGCCGATGGCGAAGCTGCCGTAGCCCCCGGCCGCCCGGCGCAGGGACCCCCGCCGCATGCTCACGCCAAGCCCCCCCCCg           >  6:65868/1‑135 (MQ=255)
                                                                                                                             cggcTGGCCCGGCTGGCGCCCGCCCGTACCCCCCACCCCCGCATTTCCCCCAGCACGCCCAAGGCGAAGCTGCCGTAGCCCCCCGCCTCCAGGCGCACGGCCACCAGCCGCCAGCTCACGCCAACCGCCAGcccc           >  4:336976/1‑134 (MQ=255)
                                                                                                                             cggcTGGCCCGGCGGGCGCCGCCCAGTACGCCCAACGCCTGCATTGCCGCCAGCAGGCCGATGGCGAAGCTGCCGTAGCCCCAGGCCTACAGGCGCAGGGACACCAGCGGCATGCTCACGCCAAGCGCCAGCCCg           <  3:295659/135‑1 (MQ=255)
                                                                                                                                      aGGCGGGCGCCGGCCAGTACGCCCAACGCCGGCATTGCCGCCATCACGCCGATGGCGAAGCTGCCGTAGCCCCAGGCCTCCAGGCGCAGGGACACCAGCGGCATGCTCACGCCAAGCGCCAGCCCGACACTGAGc  <  5:142506/135‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                       
CGATCATCAGCCGCAGCAGCAGCCATAACGGGTAACTAGGCAGCAGCGCGAGCAGGCCAATCGACAGCGCCCCGCCCCACAGGCACAGGCGCATCGCCGAAGGTACGCCAACCCAACCTGCCAGGCGGCTGGCCAGGCTGGCGCCGGCCAGTACGCCCAACGCCGGCATTGCCGCCATCACGCCGATGGCGAAGCTGCCGTAGCCCCAGGCCTCCAGGCGCAGGGACACCAGCGGCATGCTCACGCCAAGCGCCAGCCCGACACTGAGC  >  ACB122/6020206‑6020474

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