Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,274,854 G→A 58.3% T108M (ACG→ATG)  ychN ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,274,8540GA58.3% 6.8 / 26.4 24T108M (ACG→ATG) ychNconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (3/7);  new base A (6/8);  total (9/15)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

CCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTG  >  W3110S.gb/1274788‑1274911
                                                                  |                                                         
cccGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTc                                                         <  1:1029312/69‑1 (MQ=255)
cccGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTc                                                         <  1:63376/69‑1 (MQ=255)
cccGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTc                                                         <  1:1255941/69‑1 (MQ=255)
                  aGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACaa                                      <  1:947235/70‑1 (MQ=255)
                  aGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACaa                                      <  1:135970/70‑1 (MQ=255)
                         cccTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAg                                              >  1:637200/1‑55 (MQ=255)
                         cccTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCATCCATTGCGCCAGTTCCACAAGAGTACCg                              <  1:1485711/71‑1 (MQ=255)
                           cTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCATCCATTGCGCCAGTTCCACAAGAGTACAGAt                            >  1:1113009/1‑71 (MQ=255)
                            tCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCATc                                                         <  1:1185799/41‑1 (MQ=255)
                             cAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCATCCATTGCGCCAGTTCCACAAGAGTACCGAttt                          >  1:936238/1‑71 (MQ=255)
                              aaaCATTAAAATGTGAGCACTTTATCGGCTGACAGCATCCATTGCGCCAGTTCCACAAGAGTACCGAtt                           <  1:1089084/69‑1 (MQ=255)
                                aCATTAAAATGTGAGCACTTTATCGGCTGACAGCATCCa                                                       >  1:1125103/1‑39 (MQ=255)
                                aCATTAAAATGTGAGCACTTTATCGGCTGACAGCATCCa                                                       >  1:594824/1‑39 (MQ=255)
                                  aTTAAAATGTGAGCACTTTATCGGCTGACAGCATCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCAc                      >  1:1308797/1‑70 (MQ=255)
                                   ttAAAATGTGAGCACTTTATCGGCTGACAGCATcc                                                        <  1:300699/35‑1 (MQ=255)
                                     aaaaTGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCa                                               >  1:472642/1‑42 (MQ=255)
                                             aGCACTTTATCGGCTGACAGCATCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCatcaatca           <  1:897662/70‑1 (MQ=255)
                                                  tttATCGGCTGACAGCATCCATTGCGCCAGTTCCACAAGAGTAcc                               <  1:1229459/45‑1 (MQ=255)
                                                  tttATCGGCTGACAGCATCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAgg       >  1:137982/1‑69 (MQ=255)
                                                    tATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGAt                            >  1:1284842/1‑46 (MQ=255)
                                                    tATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAg    <  1:908573/70‑1 (MQ=255)
                                                      tCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAgtg  <  1:805195/70‑1 (MQ=255)
                                                      tCGGCTGACAGCATCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAgtg  <  1:791648/70‑1 (MQ=255)
                                                             aCAGCATCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAgg       <  1:1028513/58‑1 (MQ=255)
                                                                  |                                                         
CCCGCATAGTGGGGCGTCAGACGCCCCCTCAAACATTAAAATGTGAGCACTTTATCGGCTGACAGCGTCCATTGCGCCAGTTCCACAAGAGTACCGATTTCCACCCCATCAATCAGAGGAAGTG  >  W3110S.gb/1274788‑1274911

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: