Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 36,749 C→T 36.5% R364H (CGC→CAC)  caiC ← predicted crotonobetaine CoA ligase:carnitine CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb36,7490CT36.5% 53.8 / 50.4 52R364H (CGC→CAC) caiCpredicted crotonobetaine CoA ligase:carnitine CoA ligase
Reads supporting (aligned to +/- strand):  ref base C (23/10);  new base T (2/17);  total (25/27)
Fisher's exact test for biased strand distribution p-value = 3.99e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTT  >  W3110S.gb/36682‑36818
                                                                   |                                                                     
tttGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGAt                                                                    >  1:4155221/1‑71 (MQ=255)
tttGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGAt                                                                    >  1:4822803/1‑71 (MQ=255)
  tGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGTGa                                                                     <  1:1325787/68‑1 (MQ=255)
  tGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGTGa                                                                     <  1:3631872/68‑1 (MQ=255)
  tGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGTGAt                                                                    <  1:4054462/69‑1 (MQ=255)
  tGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGTGAt                                                                    <  1:4048731/69‑1 (MQ=255)
  tGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATtg                                                                  >  1:3227957/1‑71 (MQ=255)
    aaGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGTGATtgt                                                                 <  1:3564702/70‑1 (MQ=255)
     aGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGTGATtgtg                                                                <  1:1616868/70‑1 (MQ=255)
            tttCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGTGa                                                                     >  1:3352432/1‑58 (MQ=255)
             ttCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTcgcg                                                       >  1:2503534/1‑71 (MQ=255)
               cccAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGTGATTGTGATCGTCGCgg                                                      <  1:877960/70‑1 (MQ=255)
               cccAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGTGATTGTGATCGTCGCgg                                                      <  1:3130025/70‑1 (MQ=255)
               cccAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGTGATTGTGATCGTCGCgg                                                      <  1:1880695/70‑1 (MQ=255)
                         ccTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTc                                           >  1:690586/1‑71 (MQ=255)
                           tttAATGCAGATTTCACCGATCTCACCAGGCGGGAGCGGGTGATtgtg                                                                <  1:3544899/48‑1 (MQ=255)
                           tttAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGTGATtgtg                                                                <  1:2722255/48‑1 (MQ=255)
                             tAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGTGATTGTGATCGTCGCGGATCTCCGCTTCGt                                         <  1:338041/69‑1 (MQ=255)
                             tAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGTGATTGTGATCGTCGCGGATCTCCGCTTCGTAg                                       >  1:3988047/1‑71 (MQ=255)
                             tAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGAt                                                    >  1:2233156/1‑58 (MQ=255)
                             tAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGAt                                                    >  1:3231905/1‑58 (MQ=255)
                             tAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGAt                                                    >  1:3292659/1‑58 (MQ=255)
                             tAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAg                                       >  1:241483/1‑71 (MQ=255)
                             tAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAg                                       >  1:5133534/1‑71 (MQ=255)
                             tAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAg                                       >  1:3712601/1‑71 (MQ=255)
                             tAATGCAGATTTAACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAg                                       >  1:1634099/1‑71 (MQ=255)
                                tGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATtgtg                                                                <  1:4365989/43‑1 (MQ=255)
                                tGCAGATTTCACCGATCTCACCAGCCGGAGCGGGc                                                                       <  1:4635501/35‑1 (MQ=37)
                                 gCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGAtcgt                                                           >  1:590981/1‑47 (MQ=255)
                                  cAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCaaaa                                  >  1:1700406/1‑71 (MQ=255)
                                   aGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAAc                                 >  1:3015578/1‑71 (MQ=255)
                                   aGATTTCACCGATCTCACCAGCCGGGAGCGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAAcc                                <  1:4841889/71‑1 (MQ=255)
                                     aTTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGc                                      <  1:3834538/64‑1 (MQ=255)
                                      tttCACCGATCTCACCAGCCGGGAGCTGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAAcccc                              <  1:3508146/71‑1 (MQ=255)
                                       ttCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCa                             >  1:971627/1‑71 (MQ=255)
                                       ttCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCa                             >  1:4916672/1‑71 (MQ=255)
                                        tCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCAc                            >  1:4734701/1‑71 (MQ=255)
                                             gATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCaa                                    <  1:4692069/58‑1 (MQ=255)
                                              aTCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCg                                          <  1:3603166/51‑1 (MQ=255)
                                              aTCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCg                                          <  1:2095866/51‑1 (MQ=255)
                                                ctcACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCaa                    >  1:2371544/1‑71 (MQ=255)
                                                  cACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATc                  >  1:2089636/1‑71 (MQ=255)
                                                   aCCAGCCGGGAGCGGGTGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCa                     <  1:3525528/67‑1 (MQ=255)
                                                   aCCAGCCGGGAGCGGGTGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCg                 <  1:3597137/71‑1 (MQ=255)
                                                   aCCAGCCGGGAGCGGGTGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCg                 <  1:1199460/71‑1 (MQ=255)
                                                   aCCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATc                  <  1:283918/70‑1 (MQ=255)
                                                   aCCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACAAATc                  >  1:4456091/1‑70 (MQ=255)
                                                    ccAGCCGGGAGCGGGTGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGa                <  1:2921833/71‑1 (MQ=255)
                                                    ccAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGa                >  1:4461670/1‑71 (MQ=255)
                                                       gCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACgg             >  1:1405538/1‑71 (MQ=255)
                                                        ccGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGc            <  1:147674/71‑1 (MQ=255)
                                                         cGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGcc           <  1:2414105/71‑1 (MQ=255)
                                                             aGCGGGTGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGcg       <  1:1008507/71‑1 (MQ=255)
                                                                   cGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGtt  >  1:2350701/1‑70 (MQ=255)
                                                                   |                                                                     
TTTGAAGATGGTTTTCCCAGGTATGCCTTTAATGCAGATTTCACCGATCTCACCAGCCGGGAGCGGGCGATTGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTT  >  W3110S.gb/36682‑36818

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