Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,699,899 G→A 21.2% G169G (GGC→GGT rbsB ← D‑ribose transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,699,8990GA21.2% 67.7 / 13.4 33G169G (GGC→GGTrbsBD‑ribose transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (12/14);  new base A (4/3);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

ATCAAAATCTGCTGGCTGGCTGGCAAGAACATTAAACTTGTGAGCAGCAACGGCCTGCTGGAAGCCTTCGCCACGTTCACGGGCTGCGGATGTACCAGCAATGCCTTGCAGCTCGATAACTTTGGCACCTTC  >  W3110S.gb/3699836‑3699967
                                                               |                                                                    
aTCAAAATCTGCTGGCTGGCTGGCAAGAACATTAAACTTGTGAGCAGCAACGGCCTGCTGGAAGCCTTCGc                                                               <  1:2432585/71‑1 (MQ=255)
  cAAAATCTGCTGGCTGGCTGGCAAGAACATTAAACTTGTGAGCAGCAACGGCCTGCTGGAAGCCTTCGCCa                                                             >  1:4371220/1‑71 (MQ=255)
   aaaaTCTGCTGGCTGGCTGGCAAGAACATTAAACTTGTGAGCAGCAACGGCCTGCTGGAAACCTTCGCCAc                                                            >  1:4022648/1‑71 (MQ=255)
      aTCTGCTGGCTGGCTGGCAAGAACATTAAACTTGTGAGCAGCAACGGCCTGCTGGAAGCCTTCTCCACGtt                                                         <  1:153322/71‑1 (MQ=255)
         tgctgGCTGGCTGGCAAGAACATTAAACTTGTGAGCAGCAACGGCCTGCTGGAAACCTTCGCCACGTTCAc                                                      >  1:1857847/1‑71 (MQ=255)
               ctggctggcAAGAACATTAAACTTGTGAGCAGCAACGGCCTGCTGGAAGCCTTCGcc                                                              <  1:4152719/57‑1 (MQ=255)
                   ctggcAAGAACATTAAACTTGTGAGCAGCAACGGCCTGCTGGAAGCCTTCGc                                                               >  1:4179914/1‑52 (MQ=255)
                   ctggcAAGAACATTAAACTTGTGAGCAGCAACGGCCTGCTGGAAGCCTTCGCCACGTTCACGGGCTGCGGa                                            <  1:1410235/71‑1 (MQ=255)
                        aaGAACATTAAACTTGTGAGCAGCAACGGCCTGCTGGAAGCCTTCGcc                                                              <  1:4303357/48‑1 (MQ=255)
                            aCATTAAACTTGTGAGCAGCAACGGCCTGCTGGAAGCCTTCGCCACGTTCACGGGCTGCGGATGTACCAGc                                   <  1:2032472/71‑1 (MQ=255)
                            aCATTAAACTTGTGAGCAGCAACGGCCTGCTGGAAACCTTCGCCACGTTCACGGGCTGCGGATGTACCAGc                                   <  1:348206/71‑1 (MQ=255)
                            aCATTAAACTTGTGAGCAGCAACGGCCTGCTGGAAACCTTCGCCACGTTCACGGGCTGCGGATGTACCAGc                                   <  1:2272000/71‑1 (MQ=255)
                            aCATTAAAATTGTGAGCAGCAACGGCCTGCTGGAAGCCTTCGCCACGTTCACggg                                                   >  1:1543795/1‑55 (MQ=255)
                               ttAAACTTGTGAGCAGCAACGGCCTGCTGGAAACCTTCGCCACGTTCACGGGCTGCGGATGTACCAGCAAt                                >  1:4243681/1‑71 (MQ=255)
                                tAAACTTGTGAGCAGCAACGGCCTGCTGGAAGCCTTCGCCACGTTCACGGGCTGCGGATGTACCAGCAATg                               >  1:814052/1‑71 (MQ=255)
                                 aaaCTTGTGAGCAGCAACGGCCTGCTGGAAGCCTTCGCCACGTTCa                                                       >  1:3239925/1‑46 (MQ=255)
                                 aaaCTTGTGAGCAGCAACGGCCTGCTGGAAACCTTCg                                                                >  1:3934400/1‑37 (MQ=255)
                                 aaaCTTGTGAGCAGCAACGGCCTGCTGGAAACCTTCGCCACGTTCACGGGCTGCGGATGTACCAGCAATg                               <  1:787786/70‑1 (MQ=255)
                                    cTTGTGAGCAGCAACGGCCTGCTGGAAGCCTTCGCCACGTTCACGGGCt                                                 >  1:2021016/1‑49 (MQ=255)
                                     ttGTGAGCAGCAACGGCCTGCTGGAAGCCTTCGCCACGTTCACg                                                     <  1:3296323/44‑1 (MQ=255)
                                       gtgAGCAGCAACGGCCTGCTGGAAGCCTTCGCCACGTTCACGGGCTGCGGATGTACCAGCaa                                 >  1:4472600/1‑62 (MQ=255)
                                       gtgAGCAGCAACGGCCTGCTGGAAGCCTTCGCCACGTTCACGGGCTGCGGATGTACCAGCAATGCCTTg                          <  1:3599948/69‑1 (MQ=255)
                                       gtgAGCAGCAACGGCCTGCTGGAAGCCTTCGCCACGTTCACGGGCTGCGGATGTACCAGCAATGCCTTg                          <  1:3648915/69‑1 (MQ=255)
                                          agcagcAACGGACTGCTGGAAGCCTTCGCCACGTTCACGGGCTGCGGATGTACCAGCAATGCCTTGCAGc                      >  1:1568897/1‑70 (MQ=255)
                                           gcagcaACGGCCTGCTGGAAGCCTTCGCCACGTTCACGGGCTGCGGATGTACCAGCaa                                 >  1:4073290/1‑58 (MQ=255)
                                            cagcaACGGCCTGCTGGAAGCCTTCGCCACGTTCACGGGCTGCGGATGTACCAGCAATGCCTTGCAGCTc                    >  1:3297107/1‑70 (MQ=255)
                                               caACGGCCTGCTGGAAGCCTTCGCCACGTTCACGGGCTGCGGATGTACCAGCAATGCCTTGCAGCTCGAt                 <  1:3463977/70‑1 (MQ=255)
                                                    gCCTGCTGGAAGCCTTCGCCACGTTCACGGGCTGCGGATGTACCAGCAATGCCTTGCAGCTCGATAACttt           <  1:2809634/71‑1 (MQ=255)
                                                      ctgctgGAAGCCTTCGCCACGTTCACGGGCTGCGGATGTACCAGCAATGCCTTGCAGCTCGATAACTTTgg         <  1:1661860/71‑1 (MQ=255)
                                                         ctgGAAGCCTTCGCCACGTTCACGGGCTGCGGATGTAc                                       >  1:1293106/1‑38 (MQ=255)
                                                           ggAAGCCTTCGCCACGTTCACGGGCTGCGGATGTACCAGCAATGCCTTGCAGCTCGATAACTTTg          >  1:1809888/1‑65 (MQ=255)
                                                            gAAGCCTTCGCCACGTTCACGGGCTGCGGATGTACCAGCAATGCCTTGCAGCTCGATAACTTTGGCAcc     <  1:2569300/69‑1 (MQ=255)
                                                             aaGCCTTCGCCACGTTCACGGGCTGCGGATGTACCAGCAATGCCTTGCAGCTCGATAACTTTGGCACCTTc  <  1:2768274/71‑1 (MQ=255)
                                                               gCCTTCGCCACGTTCACGGGCTGCGGATGTACCAGCAATGCCTTGCAGCTCGATAACTTTGGCACCtt   <  1:1071715/68‑1 (MQ=255)
                                                               |                                                                    
ATCAAAATCTGCTGGCTGGCTGGCAAGAACATTAAACTTGTGAGCAGCAACGGCCTGCTGGAAGCCTTCGCCACGTTCACGGGCTGCGGATGTACCAGCAATGCCTTGCAGCTCGATAACTTTGGCACCTTC  >  W3110S.gb/3699836‑3699967

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: