Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,184,032 G→A 16.7% Y125Y (TAC→TAT potD ← polyamine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,184,0320GA16.7% 70.7 / 7.4 30Y125Y (TAC→TATpotDpolyamine transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (15/10);  new base A (3/2);  total (18/12)
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.48e-01

CAGATTTCGGATCCACCGCATCACCGTTAACACCAATCGCCGTCGCACCCCAGATATACGGAATGGAGTAGTCGTTATTCGGGTCAAAAGGCTTGTTGAGCATGTCTGGATCGAGATTGCTGAAGTTTGTTAACTTC  >  W3110S.gb/1183965‑1184101
                                                                   |                                                                     
cAGATTTCGGATCCACCGCATCACCGTTAACACCAATCGCCGTCGCACCCCAGATATACGGAATGGagtag                                                                    <  1:1439096/71‑1 (MQ=255)
 aGATTTCGGATCCACCGCATCACCGTTAACACCAATCGCCGTCGCACCCCAGATATACGGAATGGagtagt                                                                   >  1:2071992/1‑71 (MQ=255)
      tCGGATCCACCGCATCACCGTTAACACCAATCGCCGTCGCACCCCAGATATACGGAATGGAGTAGTCGtt                                                               <  1:1987577/70‑1 (MQ=255)
       cGGATCCACCGCATCACCGTTAACACCAATCGCCGTCGCACCCCAGATATACGGAATGGAGTAGTCGTTa                                                              <  1:514661/70‑1 (MQ=255)
        ggATCCACCGCATCACCGTTAACACCAATCGCCGTCGCACCCCAGATATACGGAATGGAGTAGTCGTTAt                                                             <  1:1911361/70‑1 (MQ=255)
                 gCATCACCGTTAACACCAATCGCCGTCGCACCCCAGATATACGGAATGGAGTAGTCGTTATTCGGGTCaa                                                    >  1:1164128/1‑70 (MQ=255)
                           tAACACCAATCGCCGTCGCACCCCAGATATACGGAATGGagtagt                                                                   >  1:730424/1‑45 (MQ=255)
                            aaCACCAATCGCCGTCGCACCCCAGATATACGGAATGGAGTAGTCGTTATTCGGGTCAAAAGGCttgttg                                         <  1:1576029/70‑1 (MQ=255)
                             acacCAATCGCCGTCGCACCCCAGATATACGGAATGGAGTAGTCGTTATTCGGGTCAAAAGGCTTGTTGa                                        <  1:1939907/70‑1 (MQ=255)
                              cacCAATCGCCGTCGCACCCCAGATATACGGAATGGAGTAGTCGTTATTCGGGTCAAAAGGCTTGTTGAg                                       >  1:384694/1‑70 (MQ=255)
                               acCAATCGCCGTCGCACCCCAGATATACGGAATGGAGTAGTCGTTATTCGGGTCaaa                                                   >  1:323825/1‑57 (MQ=255)
                                 cAATCGCCGTCGCACCCCAGATATACGGAATGGAGTAGTCGTTATTCGGGTCAAAAGGCTTGTTGAg                                       >  1:1125406/1‑67 (MQ=255)
                                 cAATCGCCGTCGCACCCCAGATATACGGAATGGAGTAGTCGTTATTCGGGTCAAAAGGCTTGTTGAGc                                      >  1:1373431/1‑68 (MQ=255)
                                  aaTCGCCGTCGCACCCCAGATATACGGAATGGAGTAGTCGTTATTCGGGTCaaaa                                                  >  1:808444/1‑55 (MQ=255)
                                  aaTCGCCGTCGCACCCCAGATATACGGAATGGAGTAGTCGTTATTCGGGTCAAAAGGCTTGTTGAGCATGt                                  >  1:1870619/1‑71 (MQ=255)
                                         gTCGCACCCCAGATATACGGAATGGAGTAGTCGTTATTCGGGTCAAAAGGCTTGTTGAGCATGTCTGGAt                            >  1:1831508/1‑70 (MQ=255)
                                         gTCGCACCCCAGATATACGGAATGGAGTAGTCGTTATTCGGGTCAAAAGGCTTGTTGAGCATGTCTGGATc                           >  1:1151389/1‑71 (MQ=255)
                                           cGCACCCCAGATATACGGAATGGAGTAGTCGTTATTCGGGTCAAAAGGCttgtt                                          >  1:824970/1‑54 (MQ=255)
                                             cACCCCAGATATACGGAATGGAATAGTCGTTATTCGGGTCAAACGGCTTGTTGAGCATGTCTGGATCgagg                       <  1:248305/71‑2 (MQ=255)
                                             cACCCCAGATATACGGAATGGAATAGTCGTTATTCGGGTCAAACGGCTTGTTGAGCATGTCTGGATCgagg                       <  1:1255324/71‑2 (MQ=255)
                                               ccccAGATATACGGAATGGAGTAGTCGTTATTCGGGTCAAAAGGCttgttg                                         >  1:562421/1‑51 (MQ=255)
                                               ccccAGATATACGGAATGGAGTAGTCGTTATTCGGGTCAAAAGGCTTGTTGAGCATGTCTGGAt                            >  1:493430/1‑64 (MQ=255)
                                                cccAGATATACGGAATGGAGTAGTCGTTATTCGGGTCaaaa                                                  <  1:576123/41‑1 (MQ=255)
                                                  cAGATATACGGAATGGAATAGTCGTTATTCGGGTCAAACGGCTTGTTGAGCATGTCTGGATCGAGGTTGc                   >  1:749794/1‑70 (MQ=255)
                                                  cAGATATACGGAATGGAATAGTCGTTATTCGGGTCAAACGGCTTGTTGAGCATGTCTGGATCGAGGTTGc                   >  1:136346/1‑70 (MQ=255)
                                                  cAGATATACGGAATGGAATAGTCGTTATTCGGGTCAAACGGCTTGTTGAGAATGTCTGGATCGAGGTTGc                   >  1:273234/1‑70 (MQ=255)
                                                       ataCGGAATGGAGTAGTCGTTATTCGGGTCAAAAGGCTTGTTGAGCATGTCt                                >  1:237988/1‑52 (MQ=255)
                                                        taCGGAATGGAGTAGTCGTTATTCGGGTCAAAAGGCTTGTTGAGCAt                                    <  1:745859/47‑1 (MQ=255)
                                                        taCGGAATGGAGTAGTCGTTATTCGGGTCAAAAGGCTTGTTGAGCATGTCTGGATCGAGATTGCTGAAGtt            >  1:1427228/1‑71 (MQ=255)
                                                          cGGAATGGAGTAGTCGTTATTCGGGTCAAAAGGCTTGTTGAGCATGTCTGGATCGAGATTGCTGAAGTTTg          <  1:838052/71‑1 (MQ=255)
                                                          cGGAATGGAGTAGTCGTTATTCGGGTCAAAAGGCTTGTTGAGCATGTCTGGATCGAGATTGCTGAAGTTTg          <  1:863148/71‑1 (MQ=255)
                                                             aaTGGAGTAGTCGTTATTCGGGTCAAAAGGCTTGTTGAGCATGTCTGGATCGAGATTGCTGAAGTTTGTTa       <  1:2035679/71‑1 (MQ=255)
                                                                 gAGTAGTCGTTATTCGGGTCAAAAGGCTTGTTGAGCATGTCTGGATCGAGATTGCTGaa               >  1:705753/1‑59 (MQ=255)
                                                                   gtagtCGTTATTCGGGTCAAAAGGCTTGTTGAGCATGTCTGGATCGAGATTGCTGAAGTTTGTTAACTTc  >  1:1334836/1‑70 (MQ=255)
                                                                   |                                                                     
CAGATTTCGGATCCACCGCATCACCGTTAACACCAATCGCCGTCGCACCCCAGATATACGGAATGGAGTAGTCGTTATTCGGGTCAAAAGGCTTGTTGAGCATGTCTGGATCGAGATTGCTGAAGTTTGTTAACTTC  >  W3110S.gb/1183965‑1184101

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: