Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,046,323 G→A 46.5% G38D (GGC→GAC)  livG → leucine/isoleucine/valine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,046,3230GA46.5% 8.0 / 56.6 43G38D (GGC→GAC) livGleucine/isoleucine/valine transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (7/16);  new base A (12/8);  total (19/24)
Fisher's exact test for biased strand distribution p-value = 6.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-01

GGCCTGCTGGCGGTGAACAACGTCAATCTTGAACTGTACCCGCAGGAGATCGTCTCGTTAATCGGCCCTAACGGTGCCGGAAAAACCACGGTTTTTAACTGTCTGACCGGATTCTACAAACCCACC  >  W3110S.gb/4046259‑4046384
                                                                |                                                             
ggcCTGCTGGCGGTGAACAACGTCAATCTTGAACTGTACCCGCAGGAGATCGTCTCGTTAATCGGCCCTaa                                                         <  1:2131217/71‑1 (MQ=255)
 gcCTGCTGGCGGTGAACAACGTCAATCTTGAACTGTACCCGCAGGAGATCGTCTCGTTAATCGGCCCTaa                                                         <  1:3952231/70‑1 (MQ=255)
   ctgctgGCGGTGAACAACGTCAATCTTGAACTGTACCCGCAGGAGATCGTCTCGTTAATCGGCCCt                                                           <  1:1482255/66‑1 (MQ=255)
   ctgctgGCGGTGAACAACGTCAATCTTGAACTGTACCCGCAGGAGATCGTCTCGTTAATCGGCCCTaa                                                         >  1:3479878/1‑68 (MQ=255)
           ggTGAACAACGTCAATCTTGAACTGTACCCGCAGGAGATCGTCTCGTTAATCGACCCTAACGGTGCCgg                                                >  1:1249365/1‑69 (MQ=255)
            gTGAACAACGTCAATCTTGAACTGTACCCGCAGGAGATCGTCTCGTTAATCGACCCTAACGGTGCCGGaaa                                             <  1:2293734/71‑1 (MQ=255)
            gTGAACAACGTCAATCTTGAACTGTACCCGCAGGAGATCGTCTCGTTAATCGACCCTAACGGTGCCGGaaa                                             <  1:2056839/71‑1 (MQ=255)
              gAACAACGTCAATCTTGAACTGTACCCGCAGGAGATCGTCTCGTTAATCGACCCTAACGGTGCCGGaaaaa                                           >  1:2578873/1‑71 (MQ=255)
                 caacGTCAATCTTGAACTGTACCCGCAGGAGATCGTCTCGTTAATCGACCCTaa                                                         >  1:3457890/1‑54 (MQ=255)
                 caacGTCAATCTTGAACTGTACCCGCAGGAGATCGTCTCGTTAATCGACCCTAACGGTGCCGGAAAAACCa                                        >  1:3037387/1‑71 (MQ=255)
                 caacGTCAATCTTGAACTGTACCCGCAGGAGATCGTCTCGTTAATCGACCCTAACGGTGCCGGAAAAACCa                                        >  1:2924482/1‑71 (MQ=255)
                  aacGTCAATCTTGAACTGTACCCGCAGGAGATCGTCTCGTTAATCGGCCCTAACGGTGCCGGAAAAAc                                          <  1:1785982/68‑1 (MQ=255)
                  aacGTCAATCTTGAACTGTACCCGCAGGAGATCGTCTCGTTAATCGGCCCTAACGGTGCCGGAAAAAc                                          <  1:1926956/68‑1 (MQ=255)
                      tCAATCTTGAACTGTACCCGCAGGAGATCGTCTCGTTAATCGGACCTAAc                                                        >  1:3435095/1‑50 (MQ=255)
                      tCAATCTTGAACTGTACCCGCAGGAGATCGTCTCGTTAATCGACCCTAAc                                                        <  1:3513726/50‑1 (MQ=255)
                      tCAATCTTGAACTGTACCCGCAGGAGATCGTCTCGTTAATCGACCCTAAc                                                        <  1:2652634/50‑1 (MQ=255)
                        aaTCTTGAACTGTACCCGCAGGAGATCGTCTCGTTAATCGACCCTAACGGTGCCGGAAAAACCACGGtttt                                 >  1:1849013/1‑71 (MQ=255)
                         aTCTTGAACTGTACCCGCAGGAGATCGTCTCGTTAATCGGCCCTAACGGTGCCGGAAAAACCACGGttttt                                <  1:155644/71‑1 (MQ=255)
                            ttGAACTGTACCCGCAGGAGATCGTCTCGTTAATCGGCCCTAACGGTGCCGGaaaa                                            <  1:2127343/56‑1 (MQ=255)
                            ttGAACTGTACCCGCAGGAGATCGTCTCGTTAATCGGCCCTAACGGTGCCGGaaa                                             <  1:778928/55‑1 (MQ=255)
                            ttGAACTGTACCCGCAGGAGATCGTCTCGTTAATCGGCCCTAACGGTGCCGGaaa                                             <  1:2550317/55‑1 (MQ=255)
                            ttGAACTGTACCCGCAGGAGATCGTCTCGTTAATCGGACCTAACGGTGCCGGaaa                                             <  1:2562402/55‑1 (MQ=255)
                            ttGAACTGTACCCGCAGGAGATCGACTCGTTAATCGGCCCTAACGGTGCCGGAAAAACCACGGTTTTTaa                              >  1:1177538/1‑70 (MQ=255)
                              gAACTGTACCCGCAGGAGATCGTCTCGTTAATCGACCCTAACGGTGCCGGAAAAACCACGGTTTTTa                               <  1:1309788/67‑1 (MQ=255)
                                 cTGTACCCGCAGGAGATCGTCTCGTTAATCGGCCCt                                                           <  1:2115228/36‑1 (MQ=255)
                                 cTGTACCCGCAGGAGATCGTCTCGTTAATCGGCCCt                                                           <  1:540220/36‑1 (MQ=255)
                                 cTGTACCCGCAGGAGATCGTCTCGTTAATCGACCCt                                                           <  1:3673674/36‑1 (MQ=255)
                                  tGTACCCGCAGGAGATCGTCTCGTTAATCGACCCTAACGGTGCCGGAAAAACCACGGTTTTTAACTGTc                         >  1:2871212/1‑69 (MQ=255)
                                   gTACCCGCAGGAGATCGTCTCGTTAATCGGCCCTaa                                                         >  1:2416125/1‑36 (MQ=255)
                                   gTACCCGCAGGAGATCGTCTCGTTAATCGGCCCTaa                                                         >  1:3925625/1‑36 (MQ=255)
                                   gTACCCGCAGGAGATCGTCTCGTTAATCGGCCCTaa                                                         >  1:3517036/1‑36 (MQ=255)
                                      cccGCAGGAGATCGTCTCGTTAATCGGCCCTAACGGTGCCGGAAAAACCACGGTTTTTAACTGTCTGAcag                   <  1:3905408/71‑3 (MQ=255)
                                      cccGCAGGAGATCGTCTCGTTAATCGGCCCTAACGGTGCCGGAAAAACCACGGTTTTTAACTGTCTGACCg                   <  1:2656609/71‑1 (MQ=255)
                                           aGGAGATCGTCTCGTTAATCGGCCCTAACGGTGCCGGAAAAACCACGGTTTTTAACTGTCTGACCGGATTc              <  1:1298018/71‑1 (MQ=255)
                                               gaTCGTCTCGTTAATCGACCCTAACGGTGCCGGaaaaa                                           >  1:1152106/1‑38 (MQ=255)
                                               gaTCGTCTCGTTAATCGACCCTAACGGTGCCGGAAAAACCa                                        >  1:2971139/1‑41 (MQ=255)
                                               gaTCGTCTCGTTAATCGACCCTAACGGTGCCGGAAAAACCa                                        >  1:1919024/1‑41 (MQ=255)
                                               gaTCGTCTCGTTAATCGACCCTAACGGTGCCGGAAAAACCa                                        >  1:3964600/1‑41 (MQ=255)
                                                aTCGTCTCGTTAATCGACCCTAACGGTGCCGGAAAAACCACGGTTTTTAACTGTCTGACCGGATTCTACaa         <  1:3503700/71‑1 (MQ=255)
                                                 tCGTCTCGTTAATCGACCCTAACGGTGCCGGAAAAACCACGGTTTTTAACTGTCTGACCGGATTCTACaa         <  1:930828/70‑1 (MQ=255)
                                                    tctcGTTAATCGGCCCTAACGGTGCCGGAAAAACCACGGTTTTTAACTGTCTGACCGGATTCTACAAAccc     <  1:1839596/71‑1 (MQ=255)
                                                        gTTAATCGGCCCTAACGGTGCCGGAAAAACCACGGTTTTTAACTGTCTGACCGGATTCTACAAACCCAcc  >  1:1005348/1‑70 (MQ=255)
                                                           aaTCGACCCTAACGGTGCCGGAAAAACCACGGTTTTTAAc                             >  1:483402/1‑40 (MQ=255)
                                                                |                                                             
GGCCTGCTGGCGGTGAACAACGTCAATCTTGAACTGTACCCGCAGGAGATCGTCTCGTTAATCGGCCCTAACGGTGCCGGAAAAACCACGGTTTTTAACTGTCTGACCGGATTCTACAAACCCACC  >  W3110S.gb/4046259‑4046384

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: