Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 847,104 A→G 34.6% R146R (CGT→CGC glnP ← glutamine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb847,1040AG34.6% 26.2 / 22.0 28R146R (CGT→CGCglnPglutamine transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (8/10);  new base G (5/4);  total (13/15)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCTTTAATGCTGATGATCCACTGGTTACCCAGCGGCGGCAGCATACGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAACGTGAAAGACCGAGCGC  >  W3110S.gb/847059‑847166
                                             |                                                              
gtCTTTAATGCTGATGATCCACTGGTTACCCAGCGGCGGCAGCATACGACGCAGTGCCAGCGGTAAAAt                                         >  1:2118545/1‑69 (MQ=255)
   tttAATGCTGATGATCCACTGGTTACCCAGCGGCGGCAGCATGCGACGCAGTGCCAGCGGTAAAATGa                                       >  1:2607391/1‑68 (MQ=255)
   tttAATGCTGATGATCCACTGGTTACCCAGCGGCGGCAGCATGCGACGCAGTGCCAGCGGTAAAATGa                                       >  1:1321699/1‑68 (MQ=255)
   tttAATGCTGATGATCCACTGGTTACCCAGCGGCGGCAGCATGCGACGCAGTGCCAGCGGTAAAATGa                                       >  1:1744755/1‑68 (MQ=255)
   tttAATGCTGATGATCCACTGGTTACCCAGCGGCGGCAGCATACGACGCAGTGCCAGCGGTAAAATGa                                       >  1:1309247/1‑68 (MQ=255)
             atgatCCACTGGTTACCCAGCGGCGGCAGCATACGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATg                            >  1:1392819/1‑69 (MQ=255)
                 tCCACTGGTTACCCAGCGGCGGCAGCATGCGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGttt                        >  1:43647/1‑69 (MQ=255)
                 tCCACTGGTTACCCAGCGGCGGCAGCATGCGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGttt                        >  1:2477123/1‑69 (MQ=255)
                  ccACTGGTTACCCAGCGGCGGCAGCATACGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGttt                        >  1:1270231/1‑68 (MQ=255)
                  ccACTGGTTACCCAGCGGCGGCAGCATACGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGttt                        >  1:3118829/1‑68 (MQ=255)
                  ccACTGGTTACCCAGCGGCGGCAGCATACGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGttt                        >  1:1861745/1‑68 (MQ=255)
                  ccACTGGTTACCCAGCGGCGGCAGCATACGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGttt                        >  1:386124/1‑68 (MQ=255)
                       ggTTACCCAGCGGCGGCAGCATACGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAAc                  >  1:2595432/1‑69 (MQ=255)
                           aCCCAGCGGCGGCAGCATACGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAACGTGa              <  1:2177584/69‑1 (MQ=255)
                           aCCCAGCGGCGGCAGCATACGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAACGTGa              <  1:2467843/69‑1 (MQ=255)
                                atggcggcAGCATACGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAACGTg               <  1:1656221/61‑1 (MQ=255)
                                 cggcggCAGCATGCGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAACGTGAAAGAcc        <  1:156944/69‑1 (MQ=255)
                                 cggcggCAGCATGCGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAACGTGAAAGAcc        <  1:2836105/69‑1 (MQ=255)
                                 cggcggCAGCATCCGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAACGTGAAAGAcc        <  1:2398931/69‑1 (MQ=255)
                                  ggcggcAGCATGCGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAACGTGAAAGAcc        <  1:1107031/68‑1 (MQ=255)
                                  ggcggcAGCATACGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAACGTGAAAGACCg       <  1:2898092/69‑1 (MQ=255)
                                   gcggcAGCATACGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAACGTGAAAGACCGa      <  1:3029179/69‑1 (MQ=255)
                                   gcggcAGCATACGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAACGTGAAAGACCGa      <  1:2104926/69‑1 (MQ=255)
                                   gcggcAGCATACGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAACGTGAAAGACCGa      <  1:84074/69‑1 (MQ=255)
                                      gcagcaTACGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAACGTGAAAGACCGAgcg   <  1:248826/69‑1 (MQ=255)
                                      gcagcaTACGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAACGTGAAAGACCGAgcg   <  1:1502838/69‑1 (MQ=255)
                                       cagcaTACGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAACGTGAAAGACCGAgcgc  <  1:3675012/69‑1 (MQ=255)
                                          caTGCGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAACGTGAAAGAcc        <  1:1369572/60‑1 (MQ=255)
                                             |                                                              
GTCTTTAATGCTGATGATCCACTGGTTACCCAGCGGCGGCAGCATACGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAACGTGAAAGACCGAGCGC  >  W3110S.gb/847059‑847166

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: