Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,733,486 C→A 19.0% L896M (CTG→ATG)  lhr → predicted ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,733,4860CA19.0% 188.0 / 21.1 63L896M (CTG→ATG) lhrpredicted ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base C (36/15);  new base A (5/7);  total (41/22)
Fisher's exact test for biased strand distribution p-value = 9.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

GGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCGGGCATCCGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATCCGGGGGAAT  >  W3110S.gb/1733417‑1733537
                                                                     |                                                   
ggCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCGGGCATCCGTACTGTCGc                                                     <  1:1390510/70‑1 (MQ=255)
                 tATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCGGGCATCCGTACTGTCGCTGGACAGCGAGTTACTgc                                   >  1:5031359/1‑71 (MQ=255)
                        cAGAGCGACGCCCCGCTGGCAGAGCGCCGGGCATCCGTACTGTCGCTGGACAGCGAGtt                                        <  1:1145424/59‑1 (MQ=255)
                        cAGAGCGACGCCCCGCTGGCAGAGCGCCGGGCATCCGTACTGTCGCTGGACAGCGAGtt                                        <  1:1923669/59‑1 (MQ=255)
                        cAGAGCGACGCCCCGCTGGCAGAGCGCCGGGCATCCGTACTGTCGCTGGACAGCGAGtt                                        <  1:17470/59‑1 (MQ=255)
                             cGACGCCCCGCTGGCAGAGCGCCGGGCATCCGTACTGTCGCTGGACAGCGAg                                          >  1:2504846/1‑52 (MQ=255)
                             cGACGCCCCGCTGGCAGAGCGCCGGGCATCCGTACTGTCGCTGGACAGCGAg                                          >  1:2397504/1‑52 (MQ=255)
                             cGACGCCCCGCTGGCAGAGCGCCGGGCATCCGTACTGTCGCTGGACAGCGAg                                          >  1:2368734/1‑52 (MQ=255)
                             cGACGCCCCGCTGGCAGAGCGCCGGGCATCCGTACTGTCGCTGGACAGCGAg                                          >  1:3527677/1‑52 (MQ=255)
                             cGACGCCCCGCTGGCAGAGCGCCGGGCATCCGTACTGTCGCTGGACAGCGAg                                          >  1:3647491/1‑52 (MQ=255)
                             cGACGCCCCGCTGGCAGAGCGCCGGGCATCCGTACTGTCGCTGGACAGCGAg                                          >  1:4012104/1‑52 (MQ=255)
                             cGACGCCCCGCTGGCAGAGCGCCGGGCATCCGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCg                       >  1:833278/1‑71 (MQ=255)
                             cGACGCCCCGCTGGCAGAGCGCCGGGCATCCGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCg                       <  1:2970679/71‑1 (MQ=255)
                                 gCCCCGCTGGCAGAGCGCCGGGCATCCGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCg                       >  1:3924016/1‑67 (MQ=255)
                                 gCCCCGCTGGCAGAGCGCCGGGCATCCGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCg                       >  1:1007548/1‑67 (MQ=255)
                                 gCCCCGCTGGCAGAGCGCCGGGCATCCGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCg                       >  1:1396037/1‑67 (MQ=255)
                                             gagCGCCGGGCATGCGTACTGTCGATGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATCCg         <  1:1716819/69‑1 (MQ=255)
                                             gagCGCCGGGCATCCGTACTGTCGATGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATCCg         <  1:4104822/69‑1 (MQ=255)
                                             gagCGCCGGGCATCCGTACTGTCGATGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATCCg         <  1:5504564/69‑1 (MQ=255)
                                             gagCGCCGGGCATCCGTACTGTCGATGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATCCg         <  1:3175065/69‑1 (MQ=255)
                                             gagCGCCGGGCATCCGTACTGTCGATGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATCCg         <  1:420452/69‑1 (MQ=255)
                                             gagCGCCGGGCATCCGTACTGTCGATGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATCCg         <  1:5259002/69‑1 (MQ=255)
                                             gagCGCCGGGCATCCGTACTGTCGATGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATCCg         <  1:2486662/69‑1 (MQ=255)
                                               gcgcCGGGCATCTGTACTGTCGATGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATCCggg       >  1:3910613/1‑69 (MQ=255)
                                               gcgcCGGGCATCCGTACTGTCGATGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATCCggg       >  1:3124840/1‑69 (MQ=255)
                                               gcgcCGGGCATCCGTACTGTCGATGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATCCggg       >  1:390166/1‑69 (MQ=255)
                                               gcgcCGGGCATCCGTACTGTCGATGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATCCggg       >  1:2242094/1‑69 (MQ=255)
                                               gcgcCGGGCATCCGTACTGTCGATGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATCCgg        >  1:621435/1‑68 (MQ=255)
                                                 gcCGGGCATCCGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATCCGGGGGaa   >  1:2738063/1‑71 (MQ=255)
                                                 gcCGGGCATCCGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATCCGGGGGaa   >  1:3920991/1‑71 (MQ=255)
                                                 gcCGGGCATCCGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATCCGGGGGaa   >  1:1978025/1‑71 (MQ=255)
                                                 gcCGGGCATCCGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATCCGGGGGaa   >  1:107785/1‑71 (MQ=255)
                                                        aTCCGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCg              <  1:1398678/53‑1 (MQ=255)
                                                        aTCCGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCg              <  1:4724928/53‑1 (MQ=255)
                                                        aTCCGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCg              <  1:747128/53‑1 (MQ=255)
                                                        aTCCGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCg              <  1:612708/53‑1 (MQ=255)
                                                        aTCCGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCg              <  1:2657894/53‑1 (MQ=255)
                                                        aTCCGTACTGTCGCT‑‑ACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCg              <  1:5165020/51‑1 (MQ=255)
                                                         tCCGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATCCGGGGGAAt  <  1:477507/64‑1 (MQ=255)
                                                         tCCGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATCCGGGGGAAt  <  1:5253284/64‑1 (MQ=255)
                                                         tCCGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATCCGGGGGAAt  <  1:430690/64‑1 (MQ=255)
                                                         tCCGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATCCGGGGGAAt  <  1:2562141/64‑1 (MQ=255)
                                                         tCCGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATCCGGGGGAAt  <  1:3288485/64‑1 (MQ=255)
                                                          ccGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCgg                      >  1:4535961/1‑43 (MQ=255)
                                                          ccGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGAt            >  1:1664774/1‑53 (MQ=255)
                                                          ccGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGAt            >  1:4448957/1‑53 (MQ=255)
                                                          ccGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGAt            >  1:4549293/1‑53 (MQ=255)
                                                          ccGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGAt            >  1:5284266/1‑53 (MQ=255)
                                                          ccGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGAt            >  1:811504/1‑53 (MQ=255)
                                                          ccGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGAt            >  1:1940905/1‑53 (MQ=255)
                                                          ccGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGAt            >  1:5011306/1‑53 (MQ=255)
                                                          ccGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGAt            >  1:2015654/1‑53 (MQ=255)
                                                          ccGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGAt            >  1:57432/1‑53 (MQ=255)
                                                           cGTACTGTGGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATc           >  1:1818581/1‑53 (MQ=255)
                                                           cGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATc           >  1:3236465/1‑53 (MQ=255)
                                                           cGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATc           >  1:4698154/1‑53 (MQ=255)
                                                           cGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATc           >  1:3232880/1‑53 (MQ=255)
                                                           cGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATc           >  1:3226249/1‑53 (MQ=255)
                                                           cGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATc           >  1:296663/1‑53 (MQ=255)
                                                           cGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATc           >  1:2480067/1‑53 (MQ=255)
                                                           cGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATc           >  1:2346073/1‑53 (MQ=255)
                                                           cGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATc           >  1:1842812/1‑53 (MQ=255)
                                                           cGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATc           >  1:1485733/1‑53 (MQ=255)
                                                               cTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATCCgg        >  1:911689/1‑52 (MQ=255)
                                                                     |                                                   
GGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCGGGCATCCGTACTGTCGCTGGACAGCGAGTTACTGCGCAATCTACTCGGACAGGTCGATCCGGGGGAAT  >  W3110S.gb/1733417‑1733537

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: