Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,462,102 A→G 42.9% F53L (TTC→CTC)  sixA ← phosphohistidine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,462,1020AG42.9% 2.7 / 24.0 21F53L (TTC→CTC) sixAphosphohistidine phosphatase
Reads supporting (aligned to +/- strand):  ref base A (4/8);  new base G (4/5);  total (8/13)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

CACTGGAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCAGGCGAGATTC  >  W3110S.gb/2462042‑2462171
                                                            |                                                                     
cACTGGAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAGCGGGCTCAcc                                                             <  1:967031/71‑1 (MQ=255)
cACTGGAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAGCGGGCTCAcc                                                             <  1:1693096/71‑1 (MQ=255)
  cTGGAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAGCGGGCTCACCAg                                                           <  1:387170/71‑1 (MQ=255)
  cTGGAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAg                                                           <  1:658870/71‑1 (MQ=255)
     gAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGaa                                                                       >  1:1313085/1‑56 (MQ=255)
     gAAGGCAGGTTCAGACAATCCACTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGaa                                                                       >  1:412779/1‑56 (MQ=255)
         gCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGt                                                    <  1:362303/71‑1 (MQ=255)
            ggTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCg                                                 <  1:946101/71‑1 (MQ=255)
                cAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAGCGGGCTCACCAGAACGCGTTCGAttt                                             >  1:723144/1‑71 (MQ=255)
                 aGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAGCGGGCTCACCAGAACGCGTTCGAttt                                             <  1:1664803/70‑1 (MQ=255)
                 aGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGAttt                                             <  1:118701/70‑1 (MQ=255)
                               ttcttcCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCtt                              >  1:242359/1‑71 (MQ=255)
                                 cttcCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCAc                                                              >  1:1281247/1‑37 (MQ=255)
                                 cttcCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGAt                                               >  1:1481032/1‑52 (MQ=255)
                                 cttcCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGAt                                               >  1:585989/1‑52 (MQ=255)
                                 cttcCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCttt                             <  1:1203453/70‑1 (MQ=255)
                                     cAGTGTTTGCTCGGCTCGCAGGAGCGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCa                           <  1:113390/68‑1 (MQ=255)
                                      aGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCa                       <  1:1527180/71‑1 (MQ=255)
                                           ttGCTCGGCTCGCAGGAGCGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCg                  >  1:674093/1‑71 (MQ=255)
                                                 ggctcgCAGGAGCGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCAt              >  1:581876/1‑69 (MQ=255)
                                                 ggctcgCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCa            <  1:370102/71‑1 (MQ=255)
                                                      gCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCAGGc         <  1:1545008/69‑1 (MQ=255)
                                                         ggAGCGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCAGGAGAGAt    >  1:783872/1‑71 (MQ=255)
                                                            aCGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGATCGCCATCAGGCGAGATTc  <  1:158188/70‑1 (MQ=255)
                                                            |                                                                     
CACTGGAAGGCAGGTTCAGACAATCCCCTACTTCTTCCAGTGTTTGCTCGGCTCGCAGGAACGGGCTCACCAGAACGCGTTCGATTTCCACTTTTTGACCTTTCAGCCAGTTCGCCATCAGGCGAGATTC  >  W3110S.gb/2462042‑2462171

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: