Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,739,599 C→A 13.4% S14S (TCC→TCA purR → DNA‑binding transcriptional repressor, hypoxanthine‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,739,5990CA13.4% 105.2 / 3.4 37S14S (TCC→TCApurRDNA‑binding transcriptional repressor, hypoxanthine‑binding
Reads supporting (aligned to +/- strand):  ref base C (18/14);  new base A (3/2);  total (21/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

GAAGTTAGGGTCTGGAGTGAAATGGAATGGCAACAATAAAAGATGTAGCGAAACGAGCAAACGTTTCCACTACAACTGTGTCACACGTGATCAACAAAACACGTTTCGTCGCTGAAGAAACGCGCAACGCCGTGTG  >  W3110S.gb/1739532‑1739667
                                                                   |                                                                    
gAAGTTAGGGTCTTGAGTGAAATGGAATGGCAACAATAAAAGATGTAGCGAAACGAGCAAACGTTTCCAc                                                                    >  1:2510302/1‑70 (MQ=255)
gAAGTTAGGGTCTGGAGTGAAATGGAATGGCAACAATAAAAGATGTAGCGAAACGAGCAAACGTTTCCAc                                                                    >  1:1068069/1‑70 (MQ=255)
gAAGTTAGGGTCTGGAGTGAAATGGAATGGCAACAATAAAAGATGTAGCGAAACGAGCAAACGTTTCCAc                                                                    >  1:2228509/1‑70 (MQ=255)
gAAGTTAGGGTCTGGAGTGAAATGGAATGGCAACAATAAAAGATGTAGCGAAACGAGCAAACGTTTCCAc                                                                    >  1:153581/1‑70 (MQ=255)
gAAGTTAGGGTCTGGAGTGAAATGGAATGGCAACAATAAAAGATGTAGCGAAACGAGCAAACGTTTCCAc                                                                    >  1:423057/1‑70 (MQ=255)
            tGGAGTGAAATGGAATGGCAACAATAAAAGATGTAGCGAAACGAGCAAACGTTTCCACTACAACTGTGTca                                                       >  1:2080933/1‑71 (MQ=255)
              gAGTGAAATGGAATGGCAACAATAAAAGATGTAGCGAAACGAGCAAACGTTTCCACTACAACTGTGTcaca                                                     >  1:3326793/1‑71 (MQ=255)
               aGTGAAATGGAATGGCAACAATAAAAGATGTAGCGAAACGAGCAAACGTTTCAACTACAACTGTGTCacac                                                    <  1:420422/71‑1 (MQ=255)
               aGTGAAATGGAATGGCAACAATAAAAGATGTAGCGAAACGAGCAAACGTTTCAACTACAACTGTGTCacac                                                    <  1:206928/71‑1 (MQ=255)
                          atggCAACAATAAAAGATGTAGCGAAACGAGCAAACGTTTCCACTACAACTGTGTCACACGTGATcaacaa                                         <  1:2154695/71‑1 (MQ=255)
                          atggCAACAATAAAAGATGTAGCGAAACGAGCAAACGTTTCCACTACAACTGTGTCACACGTGATcaacaa                                         <  1:75638/71‑1 (MQ=255)
                             gCAACAATAAAAGATGTAGCGAAACGAGCAAACGTTTCCACTACAACTGTGTCACACGTGATCAACaaaa                                       >  1:2400994/1‑70 (MQ=255)
                             gCAACAATAAAAGATGTAGCGAAACGAGCAAACGTTTCCACTACAACTGTGTCACACGTGATCAACaaaa                                       >  1:2624992/1‑70 (MQ=255)
                                  aaTAAAAGATGTAGCGAAACGAGCAAACGTTTCCACTACAACTGTGTc                                                        >  1:1844104/1‑48 (MQ=255)
                                     aaaaGATGTAGCGAAACGAGCAAACGTTTCCACTACAAc                                                              >  1:1247728/1‑39 (MQ=255)
                                        aGATGTAGCGAAACGAGCAAACGTTTCCACTACAACTGTGTCACACGTGATCAACAAAACACGTTtcgtcg                           <  1:748520/71‑1 (MQ=255)
                                         gATGTAGCGAAACGAGCAAACGTTTCCACTACAACTGTGTCACACGTGATCAACAAAACACGTTTCGTCGc                          <  1:2792965/71‑1 (MQ=255)
                                           tGTAGCGAAACGAGCAAACGTTTCCACTACAACTATGTCACACGTGATCAACAAAACACGTTTCGTCGCTg                        <  1:3759795/71‑1 (MQ=255)
                                             tAGCGAAACGAGCAAACGTTTCCACTACAACTGTGTCACACGTGATCAACAAAACACGTTTCGTCGCTga                       <  1:319092/70‑1 (MQ=255)
                                               gCGAAACGAGCAAACGTTTCCACTACAACTGTGTCACACGTGATCAACAAAACACGTTTCGTCGCTgaaga                    <  1:1106458/71‑1 (MQ=255)
                                               gCGAAACGAGCAAACGTTTCCACTACAACTGTGTCACACGTGATCAACAAAACACGTTTCGTCGCTgaaga                    <  1:3430651/71‑1 (MQ=255)
                                                     cGAGCAAACGTTTCCACTACAACTGTGTCACACGTGATCAACAAAaca                                     >  1:225731/1‑48 (MQ=255)
                                                     cGAGCAAACGTTTCCACTACAACTGTGTCACACGTGATCAACAAAACACGTTTCGTCGCTGAAGAAAcgcg              <  1:205667/71‑1 (MQ=255)
                                                      gAGCAAACGTTTCCACTACAACTGTGTCACACGTGATCAACAAAACACGTTTCGTCGCTGAAGAAACgcgc             <  1:2236462/71‑1 (MQ=255)
                                                            aCGTTTCCACTACAACTGTGTCACACGTGATCAACAAAACACGt                                  >  1:1234116/1‑44 (MQ=255)
                                                            aCGTTTCCACTACAACTGTGTCACACGTGATCAACAAAACACGt                                  >  1:3580568/1‑44 (MQ=255)
                                                            aCGTTTCCACTACAACTGTGTCACACGTGATCAACAAAACACGTTTCGTCGCTGAAGAAACGCGCAACGcc       <  1:3620363/71‑1 (MQ=255)
                                                             cGTTTCAACTACAACTGTGTCACACGTGATCAACAAAACACGTTTCGTCGCTGAAGAAACGCGCAATGCCg      >  1:3706926/1‑71 (MQ=255)
                                                             cGTTTCAACTACAACTGTGTCACACGTGATCAACAAAACACGTTTCGTCGCTGAAGAAACGCGCAATGCCg      >  1:2035711/1‑71 (MQ=255)
                                                             cGTTTCAACTACAACTGTGTCACACGTGATCAACAAAACACGTTTCGTCGCTGAAGAAACGCGCAATGCCg      >  1:3182067/1‑71 (MQ=255)
                                                              gTTTCCACTACAACTGTGTCACACGTGATCAACAAAaca                                     >  1:3201786/1‑39 (MQ=255)
                                                              gTTTCCACTACAACTGTGTCACACGTGATCAACAAAaca                                     >  1:643253/1‑39 (MQ=255)
                                                              gTTTCCACTACAACTGTGTCACACGTGATCAACAAAACACGTTTCGTCGCTGAAGAAACGCGCAACGCCgt     >  1:427456/1‑71 (MQ=255)
                                                                 tCCACTACAACTGTGTCACACGTGATCAACAAAACACGTTTCGTCGCTgaagaa                   <  1:3176352/54‑1 (MQ=255)
                                                                 tCCACTACAACTGTGTCACACGTGATCAACAAAACACGTTTCGTCGCTgaagaa                   <  1:532597/54‑1 (MQ=255)
                                                                 tCCACTACAACTGTGTCACACGTGATCAACAAAACACGTTTCGTCGCTgaagaa                   <  1:1835937/54‑1 (MQ=255)
                                                                 tCCACTACAACTGTGTCACACGTGATCAACAAAACACGTTTCGTCGCTGAAGAAACGCGCAACGCCGtgtg  >  1:634196/1‑71 (MQ=255)
                                                                   cACTACAACTGTGTCACACGTGATCAACAAAACACGttt                                >  1:3239948/1‑39 (MQ=255)
                                                                   |                                                                    
GAAGTTAGGGTCTGGAGTGAAATGGAATGGCAACAATAAAAGATGTAGCGAAACGAGCAAACGTTTCCACTACAACTGTGTCACACGTGATCAACAAAACACGTTTCGTCGCTGAAGAAACGCGCAACGCCGTGTG  >  W3110S.gb/1739532‑1739667

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: