Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 419,576 G→A 20.8% A254A (GCG→GCA brnQ → predicted branched chain amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb419,5760GA20.8% 46.9 / 8.4 24A254A (GCG→GCAbrnQpredicted branched chain amino acid transporter
Reads supporting (aligned to +/- strand):  ref base G (7/12);  new base A (2/3);  total (9/15)
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.16e-01

TGGCGGGTGTTGGTCTGACTCTGCTGTACCTGGCGCTGTTCCGTCTGGGGTCAGACAGCGCGTCGCTGGTCGATCAGTCTGCAAACGGCGCTGCTATTCTGCATGCTTACGTTCAGCACACCTTTGGC  >  W3110S.gb/419515‑419642
                                                             |                                                                  
tGGCGGGTGTTGGTCTGACTCTGCTGTACCTGGCGCTGTTCCGTCTGGGGTCAGACAGCGCATCGCTgg                                                             >  1:813004/1‑69 (MQ=255)
tGGCGGGTGTTGGTCTGACTCTGCTGTACCTGGCGCTGTTCCGTCTGGGGTCAGACAGCGCATCGCTgg                                                             >  1:308157/1‑69 (MQ=255)
     ggTGTTGGTCTGACTCTGCTGTACCTGGCGCTGTTCCGTCTGGGGTCAGACAGCGCGTCGCTCGTCGAt                                                        <  1:480891/69‑1 (MQ=255)
         ttGGTCTGACTCTGCTGTACCTGGCGCTGTTCCGTCTGGGGTCAGACAGCGCGTCGCTGGTCGATCa                                                      >  1:733948/1‑67 (MQ=255)
         ttGGTCTGACTCTGCTGTACCTGGCGCTGTTCCGTCTGGGGTCAGACAGCGCGTCGCTGGTCGATCAGt                                                    >  1:2294610/1‑69 (MQ=255)
         ttGGTCTGACTCTGCTGTACCTGGCGCTGTTCCGTCTGGGGTCAGACAGCGCGTCGCTGGTCGATCAGt                                                    >  1:3680183/1‑69 (MQ=255)
         ttGGTCTGACTCTGCTGTACCTGGCGCTGTTCCGTCTGGGGTAAGACAGCGCGTCGCTGGTCGATCAg                                                     >  1:1591737/1‑68 (MQ=255)
           ggTCTGACTCTGCTGTACCTGGCGCTGTTCCGTCTGGGGTCAGACAGCGCGTCGCTGGTCGATCAGTc                                                   >  1:2261042/1‑68 (MQ=255)
           ggTCTGACTCTGCTGTACCTGGCGCTGTTCCGTCTGGGGTCAGACAGCGCGTCGCTGGTCGATCAGTc                                                   <  1:456673/68‑1 (MQ=255)
                         gTACCTGGCGCTGTTCCGTCTGGGGTCAGACAGCGCGTCGCTGGTCGATCAGTCTGCAAACGGCgctgc                                    <  1:2967887/69‑1 (MQ=255)
                                  gcTGTTCCGTCTGGGGTCAGACAGCGCGTCGCTGGTCGATCAGTCTGCAAACGgcgc                                       <  1:1569456/57‑1 (MQ=255)
                                  gcTGTTCCGTCTGGGGTCAGACAGCGCGTCGCTGGTCGATCAGTCTGCAAACGGCGCTGCTATTCTg                             <  1:1646078/67‑1 (MQ=255)
                                  gcTGTTCCGTCTGGGGTCAGACAGCGCGTCGCTGGTCGATCAGTCTGCAAACGGCGCTGCTATTCTg                             <  1:1893540/67‑1 (MQ=255)
                                  gcTGTTCCGACTGGGGTCAGACAGCGCGTCGCTGGTCGATCAGTCTGCAAACGGCGCTGCTATTCTg                             <  1:1086448/67‑1 (MQ=255)
                                            cTGGGGTCAGACAGCGCGTCACTGGTCGATCAGTCTGCAAACGGCGCTGCTAt                                 >  1:583045/1‑53 (MQ=255)
                                                   cAGACAGCGCGTCGCTTGTCGATCAGTCTGCAAACGGCGCTGCTATTCTGCATGCTTACGTTCAg              >  1:3415002/1‑65 (MQ=255)
                                                    aGACAGCGCGTCGCTGGTCTATCAGTCTGCAAACGGCGCTGCTATTCTg                             <  1:2977284/49‑1 (MQ=255)
                                                      aCAGCGCGTCGCTGGTCGATCAGTCTGCAAACGGCGCTGCTATTCTGCATGCTTACTTTc                <  1:2655582/60‑1 (MQ=255)
                                                       cAGCGCGTCGCTGGTCGATCAGTCTGCAAACGGCGCTGCTATTCTGCATGCTTACGTTCAGCACACCtt      <  1:328917/69‑1 (MQ=255)
                                                       cAGCGCATCGCTGGTCGATCAGTCTGCAAACGGCGCTGCTATTCTGCATGCTTACGTTCAGCACACCtt      <  1:3150049/69‑1 (MQ=255)
                                                       cAGCGCATCGCTGGTCGATCAGTCTGCAAACGGCGCTGCTATTCTGCATGCTTACGTTCAGCACACCtt      <  1:2415901/69‑1 (MQ=255)
                                                       cAGCGCATCGCTGGTCGATCAGTCTGCAAACGGCGCTGCTATTCTGCATGCTTACGTTCAGCACACCtt      <  1:1371051/69‑1 (MQ=255)
                                                        aGCGCGTCGCTGGTCGATCAGTCTGCAAACGGCGCTGCTATTCTGCATGCTTACg                   <  1:260126/55‑1 (MQ=255)
                                                        aGCGCGTCGCTGGTCGATCAGTCTGCAAACGGCGCTGCTATTCTGCATGCTTACg                   <  1:1637782/55‑1 (MQ=255)
                                                           gcgTCGCTGGTCGATCAGTCTGCAAACGGCGCTGCTATTCTGCATGCTTACGTTCAGCACACCTTTggc  >  1:2398171/1‑69 (MQ=255)
                                                             |                                                                  
TGGCGGGTGTTGGTCTGACTCTGCTGTACCTGGCGCTGTTCCGTCTGGGGTCAGACAGCGCGTCGCTGGTCGATCAGTCTGCAAACGGCGCTGCTATTCTGCATGCTTACGTTCAGCACACCTTTGGC  >  W3110S.gb/419515‑419642

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: