Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,042,971 G→A 85.2% A39T (GCC→ACC)  livK → leucine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,042,9710GA85.2% 63.3 / 5.3 27A39T (GCC→ACC) livKleucine transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (8/15);  total (10/17)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

TCACACACCGCTATGGCTGACGATATTAAAGTCGCCGTTGTCGGCGCGATGTCCGGCCCGATTGCCCAGTGGGGCGATATGGAATTTAACGGCGCGCGTCAGGCAATTAAAGACATTAATGCCAA  >  W3110S.gb/4042908‑4043032
                                                               |                                                             
tCACACACCGCTATGGCTGACGATATTAAAGTCGCCGTTGTCGGCGCGATGTCCGGCCCGATTACCCAGTg                                                        >  1:1829632/1‑71 (MQ=255)
      acCGCTATGGCTGACGATATTAAAGTCGCCGTTGTCGGCGCGATGTCCGGCCCGATTACCCAGTGGGGCGa                                                  <  1:349719/71‑1 (MQ=255)
             tGGCTGACGATATTAAAGTCGCCGTTGTCGGCGCGATGTCCGGCCCGATTACCCAGTGGGGCGATATgg                                             <  1:1937848/69‑1 (MQ=255)
               gCTGACGATATTAAAGTCGCCGTTGTCGGCGCGATGTCCGGCCCGATTACCCAGTGGGGCGATATGGAAtt                                         >  1:2510539/1‑71 (MQ=255)
                cTGACGATATTAAAGTTGCCGTTGTCGGCGCGATGTCCGGCCCGATTACCCAGTGGGGCGATATGGAAttt                                        <  1:284588/71‑1 (MQ=255)
                    cGATATTAAAGTCGCCGTTGTCGGCGCGATGTCCGGCCCGATTACCCAGTGGGGCGATATgg                                             <  1:2268330/62‑1 (MQ=255)
                        atTAAAGTCGCCGTTGTCGGCGCGATGTCCGGCCCGATTACCCAGTGGGGCGATATGGAATTTaa                                      >  1:2649601/1‑65 (MQ=255)
                           aaaGTCGCCGTTGTCGGCGCGATGTCCGGCCCGATTACCCAGTGGGGCGATATGGAAtt                                         >  1:1452034/1‑59 (MQ=255)
                           aaaGTCGCCGTTGTCGGCGCGATGTCCGGCCCGATTACCCAGTGGGGCGATATGGAATTTAACGGcgcgcg                             >  1:1162656/1‑71 (MQ=255)
                           aaaGTCGCCGTTGTCGGCGCGATGTCCGGCCCGATTACCCAGTGGGGCGATATGGAATTTAACGGcgcgcg                             >  1:1903984/1‑71 (MQ=255)
                               tCGCCGTTGTCGGCGCGATGTCCGGCCCGATTACCCAGTGGGGCGATATGGAAttt                                        <  1:794082/56‑1 (MQ=255)
                                 gCCGTTGTCGGCGCGATGTCCGGCCCGATTACCCAGTGGGGCGATATGGAATTTAACGGCGCGCGTCAGGc                       <  1:1579807/71‑1 (MQ=255)
                                  ccGTTGTCGGCGCGATGTCCGGCCCGATTACCCAGTGGGGCGATATGGAATTTAACGGCGCGCGTCAGGc                       <  1:2300967/70‑1 (MQ=255)
                                     ttGTCGGCGCGATGTCCGGCCCGATTGCCCAGTGGGGCGata                                                <  1:2408220/42‑1 (MQ=255)
                                     ttGTCGGCGCGATGTCCGGCCCGATTACCCAGTGGGGCGa                                                  <  1:229519/40‑1 (MQ=255)
                                       gTCGGCGCGATGTCCGGCCCGATTGCCCAGTGGGGCGATATg                                              >  1:858901/1‑42 (MQ=255)
                                       gTCGGCGCGATGTCCGGCCCGATTACCCAGTGGGGCGATATGGAATTTAACGGCGCGCGTCAGGCAATTa                  >  1:2784449/1‑70 (MQ=255)
                                         ttgcgcGATGTCCGGCCCGATTACCCAGTGGGGCGATATGGAAttt                                        <  1:1679638/44‑1 (MQ=255)
                                          ggCGCGATGTCCGGCCCGATTGCCCAGTGGGGCGATATGGAATTTAACGGCGCGCGTCAGGCAATTAAAGa              >  1:85006/1‑71 (MQ=255)
                                              cgATGTCCGGCCCGATTACCCAGTGGGGCGATATGGAATTTAACGGCGCGCGTCAGGCAATTAAAGACAtt          <  1:1595081/71‑1 (MQ=255)
                                                aTGTCCGGCCCGATTACCCAGTGGGGCGATATGGAATTTAACGGCGCGCGTCAGGCAATTAAAGACAtt          <  1:2534710/69‑1 (MQ=255)
                                                aTGTCCGGCCCGATTACCCAGTGGGGCGATATGGAATTTAACGGCGCGCGTCAGGCAATTAAAGACAtt          <  1:1438957/69‑1 (MQ=255)
                                                aTGTCCGGCCCGATTACCCAGTGGGGCGATATGGAATTTAACGGCGCGCGTCAGGCAATTAAAGACAtt          <  1:1379580/69‑1 (MQ=255)
                                                     cGGCCCGATTACCCAGTGGGGCGATATGGAATTTAACGGCGCGCGTCAGGCAATTAAAGACATTAATGcc    >  1:512917/1‑70 (MQ=255)
                                                      ggCCCGATTGCCCAGTGGGGCGATATGGAATTTAACGGCGCGCGTCAGGCAATTAAAGACATTAATGCCaa  <  1:1853873/71‑1 (MQ=255)
                                                      ggCCCGATTACCCAGTGGGGCGATATGGAATTTAACGGCGCGCGTCAGGCAATTAAAGACATTAATGCCaa  <  1:2554735/71‑1 (MQ=255)
                                                      ggCCCGATTACCCAGTGGGGCGATATGGAATTTAACGGCGCGCGTCAGGCAATTAAAGACATTAATGCCaa  <  1:1312968/71‑1 (MQ=255)
                                                               |                                                             
TCACACACCGCTATGGCTGACGATATTAAAGTCGCCGTTGTCGGCGCGATGTCCGGCCCGATTGCCCAGTGGGGCGATATGGAATTTAACGGCGCGCGTCAGGCAATTAAAGACATTAATGCCAA  >  W3110S.gb/4042908‑4043032

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: