Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 847,104 A→G 17.2% 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,1040AG17.2% 63.9 / 8.2 29R146R (CGT→CGCglnPglutamine transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (6/18);  new base G (3/2);  total (9/20)
Fisher's exact test for biased strand distribution p-value = 2.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

GCCGATCACAATAAACAGCGAGGTGTCTTTAATGCTGATGATCCACTGGTTACCCAGCGGCGGCAGCATACGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAACGTGAAAGACCGAGCGCC  >  W3110S.gb/847035‑847167
                                                                     |                                                               
gCCGATCACAATAAACAGCGAGGTGTCTTTAATGCTGATGATCCACTGGTTACCCAGCGGCGGCAGCATac                                                                <  1:3519336/71‑1 (MQ=255)
gCCGATCACAATAAACAGCGAGGTGTCTTTAATGCTGATGATCCACTGGTTACCCAGCGGCGGCAGCATac                                                                <  1:3347838/71‑1 (MQ=255)
gCCGATCACAATAAACAGCGAGGTGTCTTTAATGCTGATGATCCACTGGTTACCCAGCGGCGGCAGCATac                                                                <  1:2065574/71‑1 (MQ=255)
     tCACAATAAACAGCGAGGTGTCTTTAATGCTGATGATCCACTGGTTACCCAGCGGCGGCAGCATACGACGc                                                           <  1:1185064/71‑1 (MQ=255)
      cacaATAAACAGCGAGGTGTCTTTAATGCTGATGATCCACTGGTTACCCAGCGGCGGCAGCATacgacg                                                            >  1:755934/1‑69 (MQ=255)
         aaTAAACAGCGAGGTGTCTTTAATGCTGATGATCCACTGGTTACCCAGCGGCGGCAGCATACGACGCAGt                                                        <  1:2028307/70‑1 (MQ=255)
           tAAACAGCGAGGTGTCTTTAATGCTGATGATCCACTGGTTACCCAGCGGCGGCAGCATGCGACGCAGTGc                                                      >  1:2635483/1‑70 (MQ=255)
             aaCAGCGAGGTGTCTTTAATGCTGATGATCCACTGGTTACCCAGCGGCGGCAGCATACGACGCAGTGCCAg                                                   >  1:2943578/1‑71 (MQ=255)
               cAGCGAGGTGTCTTTAATGCTGATGATCCACTGGTTACCCAGCGGCGGCAGCATGCGACGCAGt                                                        <  1:2320899/64‑1 (MQ=255)
                 gCGAGGTGTCTTTAATGCTGATGATCCACTGGTTACCCAGCGGCGGCAGCATACGACGCAGTGCCAGCgg                                                >  1:3117676/1‑70 (MQ=255)
                   gAGGTGTCTTTAATGCTGATGATCCACTGGTTACCCAGCGGCGGCAGCATACGACGCAGTGCCAGCGGTaa                                             <  1:2460854/71‑1 (MQ=255)
                       tgtCTTTAATGCTGATGATCCACTGGTTACCCAGCGGCGGCAGCATGCGACGCAGTGCCAGCGGTAAAATg                                         <  1:2306782/71‑1 (MQ=255)
                       tgtCTTTAATGCTGATGATCCACTGGTTACCCAGCGGCGGCAGCATACGACGCAGTGCCAGCGGTaaa                                            <  1:3922457/68‑1 (MQ=255)
                       tgtCTTTAATGCTGATGATCCACTGGTTACCCAGCGGCGGCAGCATACGACGCAGTGCCAGCGGTaaa                                            <  1:2450983/68‑1 (MQ=255)
                         tCTTTAATGCTGATGATCCACTGGTTACCCAGTGGCGGCAGCATACGACGCAGt                                                        <  1:3571513/54‑1 (MQ=255)
                           tttAATGCTGATGATCCACTGGTTACCCAGCGGCGGCAGCATACGACGCAGTGCCAGCGGTaa                                             <  1:3869136/63‑1 (MQ=255)
                             tAATGCTGATGATCCACTGGTTACCCAGCGGCGGCAGCATACGACGCAGTGCCAGCGGTAAAATGACg                                      <  1:1639829/68‑1 (MQ=255)
                             tAATGCTGATGATCCACTGGTTACCCAGCGGCGGCAGCATACGACGCAGTGCCAGCGGTAAAATGACGTAg                                   >  1:1583600/1‑71 (MQ=255)
                               aTGCTGATGATCCACTGGTTACCCAGCGGCGGCAGCATACGACGCAGTGCCAGCGGTAAAATGACGTAGCg                                 <  1:1273341/71‑1 (MQ=255)
                                  cTGATGATCCACTGGTTACCCAGCGGCGGCAGCATACGACGCAGTGCCAGCGGTAAAATGACGTAGCGaa                               <  1:2666467/70‑1 (MQ=255)
                                  cTGATGATCCACTGGTTACCCAGCGGCGGCAGCATACGACGCAGTGCCAGCGGTAAAATGACGTAGCGaa                               <  1:2450672/70‑1 (MQ=255)
                                  cTGATGATCCACTGGTTACCCAGCGGCGGCAGCATACGACGCAGTGCCAGCGGTAAAATGACGTAGCGaa                               <  1:2082060/70‑1 (MQ=255)
                                       gatCCACTGGTTACCCAGCGGCGGCAGCATACGACGCAGTGCCAGCGGTAAAATGACGTAGCg                                 <  1:2376524/63‑1 (MQ=255)
                                             cTGGTTACCCAGCGGCGGCAGCATACGACGCAGTGCCAGCGGTaaa                                            <  1:3305067/46‑1 (MQ=255)
                                                   aCCCAGCGGCGGCAGCATGCGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTcc                       >  1:3837165/1‑61 (MQ=255)
                                                   aCCCAGCGGCGGCAGCATGCGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAACGTGa               >  1:1928398/1‑69 (MQ=255)
                                                    cccAGCGGCGGCAGCATACGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAACGTGaaa             <  1:1303659/70‑1 (MQ=255)
                                                           gcggcAGCATACGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAACGTGAAAGACCGAgc     >  1:3224158/1‑71 (MQ=255)
                                                           gcggcAGCATACGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAACGTGAAAGACCGAgc     >  1:555755/1‑71 (MQ=255)
                                                           gcggcAGCATACGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAACGTGAAAGACCGAgc     >  1:74023/1‑71 (MQ=255)
                                                                    tACGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAACGTGAAAGACCGAGCGcc  <  1:3093154/65‑1 (MQ=255)
                                                                    tACGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAACGTGAAAGACCGAGCGcc  <  1:3577110/65‑1 (MQ=255)
                                                                    tACGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAACGTGAAAGACCGAGCGcc  <  1:661739/65‑1 (MQ=255)
                                                                     acgacgCAGTGCCAGCGGTAAAATGACGTAGCGaa                               >  1:3018714/1‑35 (MQ=255)
                                                                     |                                                               
GCCGATCACAATAAACAGCGAGGTGTCTTTAATGCTGATGATCCACTGGTTACCCAGCGGCGGCAGCATACGACGCAGTGCCAGCGGTAAAATGACGTAGCGAATGGTTTCCCAACGTGAAAGACCGAGCGCC  >  W3110S.gb/847035‑847167

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: