Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,626,383 C→T 12.1% A33A (GCG→GCA yigL ← predicted hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,626,3830CT12.1% 90.3 / 4.3 33A33A (GCG→GCAyigLpredicted hydrolase
Reads supporting (aligned to +/- strand):  ref base C (13/16);  new base T (2/2);  total (15/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 = 8.75e-01

CCAGATTATCGCGAATTTGCCCCACATCAACGTGGTGACGACCGGTCGCAAACACAAAGTTGATGCCGCGCGCGGTGAGCAGCTTCAGAGTTTCTTTGGCGTAAGGGGATAACGTATGGTCGGGAGAAAGTAA  >  W3110S.gb/3626313‑3626445
                                                                      |                                                              
ccAGATTATCGCGAATTTGCCCCACATCAACGTGGTGACGACCGGTCGCAAACACAAAGTTGATGCcgcgc                                                                >  1:900244/1‑71 (MQ=255)
       aTCGCGAATTTGCCCCACATCAACGTGGTGACGACCGGTCGCAAACACAAAGTTGATGCCGCGCGCGGTg                                                          <  1:1737968/70‑1 (MQ=255)
         cgcgAATTTGCCCCACATCAACGTGGTGACGACCGGTCGCAAACACAAAGTTGATGCCGCGCGCGGt                                                           <  1:1433237/67‑1 (MQ=255)
            gAATTTGCCCCACATCAACGTGGTGACGACCGGTAGCAAACACAAAGTTGATGCCGCGCGCGGTTagcag                                                     <  1:2691428/70‑1 (MQ=255)
                          tCAACGTGGTGACGACCGGTCGCAAACACAAAGTTGATGCCGCGCGCGGTGAGCAGCTTCAGAGTTTCttt                                      >  1:1118026/1‑71 (MQ=255)
                          tCAACGTGGTGACGACCGGTCGCAAACACAAAGTTGATGCCGCGCGCGGTGAGCAGCTTCAGAGTTTCttt                                      >  1:455369/1‑71 (MQ=255)
                            aaCGTGGTGACGACCGGTCGCAAACACAAAGTTGATGCCGCGTGCGGTGAGCAGCTTCAGAGTTTCTTTgg                                    >  1:1758437/1‑71 (MQ=255)
                            aaCGTGGTGACGACCGGTCGCAAACACAAAGTTGATGCCGCGCGCGGTGAGCAGCTTCAGAGTTTCTTTg                                     >  1:2332076/1‑70 (MQ=255)
                            aaCGTGGTGACGACCGGTCGCAAACACAAAGTTGATGCCGCGCGCGGTGAGCAGCTTCAGAGTTTCTTTg                                     >  1:2328994/1‑70 (MQ=255)
                            aaCGTGGTGACGACCGGTCGCAAACACAAAGTTGATGCCGCGCGCGGTGAGCAGCTTCAGAGCTTCTTTgg                                    >  1:619765/1‑71 (MQ=255)
                               gtggtgACGACCGGTCGCAATCACAAAGTTGATGCCGCGCGCGGTGAGCAGCTTCAGAGTTTCTTTGGCg                                  >  1:1017707/1‑70 (MQ=255)
                               gtggtgACGACCGGTCGCAAACACAAAGTTGATGCCGCGCGCGGTGAGCAGCTTCAGAGTTTCtt                                       >  1:1775825/1‑65 (MQ=255)
                                tggttACGACCGGTCGCAAACACAAAGTTGATGCCGCGCGCGGTGAGCAGCTTCAGAGTTTCTTTGGCGTa                                <  1:2048506/71‑1 (MQ=255)
                                tggtgACGACCGGTCGCAAACACAAAGTTGATGCcgcgcgc                                                              <  1:73222/41‑1 (MQ=255)
                                tggtgACGACCGGTCGCAAACACAAAGTTGATGCCGCGCGCGGTGAGCAGCTTCAGAGTTTCTTTGGCGTa                                <  1:159788/71‑1 (MQ=255)
                                tggtgACGACCGGTCGCAAACACAAAGTTGATGCCGCGCGCGGTGAGCAGCTTCAGAGTTTCTTTGGCGTa                                <  1:2129793/71‑1 (MQ=255)
                                  gtgACGACCGGTCGCAAACACAAAGTTGATGCCGCGTGCGGTGAGCAGCTTCAGAGTTTCTTTGGCGTaa                               <  1:524934/70‑1 (MQ=255)
                                   tgACGACCGGTCGCAAACACAAAGTTGATGCCGCGCGCGGTGagca                                                      <  1:1725990/46‑1 (MQ=255)
                                     acgacCGGTCGCAAACAGAAAGTTGATGCCGCGCGCGGTGAGCAGCTTCAGAGTTTCTTTGGCGTAAggg                            >  1:694187/1‑70 (MQ=255)
                                     acgacCGGTCGCAAACACAAAGTTGATGCCGCGCGCGGTGAGCAGCTTCAGAGTTTCTTTGGCGTAAgggg                           <  1:1891079/71‑1 (MQ=255)
                                       gacCGGTCGCAAACACAAAGTTGATGCCGCGTGCGGTGAGCAGCTTCag                                               <  1:2595506/49‑1 (MQ=255)
                                         ccGGTCGCAAACACAAAGTTGATGCCGCGCGCGGTGAGCAGCTTCAGAGTTTCTTTGGCGTAAGGGGATaa                       >  1:697335/1‑71 (MQ=255)
                                           ggTCGCAAACACAAAGTTGATGCCGCGCGCGGTGAGCAGCTTCAGAGTTTCTTTGGCGTAAGGGGATAACg                     <  1:1786404/71‑1 (MQ=255)
                                           ggTCGCAAACACAAAGTTGATGCCGCGCGCGGTGAGCAGCTTCAGAGTTTCTTTGGCGTAAGGGGATAACg                     <  1:627818/71‑1 (MQ=255)
                                           ggTCGCAAACACAAAGTTGATGCCGCGCGCGGTGAGCAGCTTCAGAGTTTCTTTGGCGTAAGGGCATAACg                     <  1:2448842/71‑1 (MQ=255)
                                            gTCGCAAACACAAAGTTGATGCCGCGCGCGGTGAGCAGCTTCAGAGTTTCTTTGGCGTAAGGGGATAACGt                    <  1:2759907/71‑1 (MQ=255)
                                              cGCAAACACAATGTTGATGCCGCGCGCGGTGAGCAGCTTCAGAGTTTCTTTGGCGTa                                <  1:2192031/57‑1 (MQ=255)
                                               gCAAACACAAAGTTGATGCCGCGCGCGGTGAGCAGCTTCAGAGTTTCTTTGGCGTAAg                              >  1:2061262/1‑58 (MQ=255)
                                               gCAAACACAAAGTTGATGCCGCGCGCGGTGAGCAGCTTCAGAGTTTCTTTGGCGTAAGGGGATAACGTATg                 >  1:1698250/1‑71 (MQ=255)
                                                     acaAAGTTGATGCCGCGCGCGGTGAGCAGCTTCAGAGTTTCTTTGGCGTAAGGGGATAACGTATGGTCgg            >  1:2499301/1‑70 (MQ=255)
                                                          gTTGATGCCGCGTGCGGTGAGCAGCTTCAGAGTTTCTTTg                                     >  1:1854259/1‑40 (MQ=255)
                                                             gATGCCGCGCGCGGTGAGCAGCTTCAGAGTTTCTTTGGCGTAAGGGGATAACGTATGGTCGGGAGAAAGt    <  1:550738/70‑1 (MQ=255)
                                                             gATGCCGCGCGCGGTGAGCAGCTTCAGAGTTTCTTTGGCGTAAGGGGATAACGTATGGTCGGGAGAAAGt    <  1:2725898/70‑1 (MQ=255)
                                                              aTGCCGCGCGCGGTGAGCAGTTTCAGAGTTTCTTTGGCGTAAGGGGATAACGTATGGTCGGGAGAAAGTaa  <  1:2393770/71‑1 (MQ=255)
                                                                 ccGCGCGCGGTGAGCAGCTTCAGAGTTTCTTTGGc                                   >  1:584572/1‑35 (MQ=255)
                                                                   gcgcgcgGTGAGCAGCTTCAGAGTTTCTTTGGCGTAAg                              >  1:1359500/1‑38 (MQ=255)
                                                                   gcgcgcgGTGAGCAGCTTCAGAGTTTCTTTGGCGTAAg                              >  1:1515045/1‑38 (MQ=255)
                                                                   gcgcgcgGTGAGCAGCTTCAGAGTTTCTTTGGCGTAAg                              >  1:1985827/1‑38 (MQ=255)
                                                                      |                                                              
CCAGATTATCGCGAATTTGCCCCACATCAACGTGGTGACGACCGGTCGCAAACACAAAGTTGATGCCGCGCGCGGTGAGCAGCTTCAGAGTTTCTTTGGCGTAAGGGGATAACGTATGGTCGGGAGAAAGTAA  >  W3110S.gb/3626313‑3626445

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