Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 1,595,970 (C)8→9 coding (3831/3978 nt) ydeK ← predicted lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,595,9621.C75.0% 26.3 / 4.2 12coding (3839/3978 nt)ydeKpredicted lipoprotein
Reads supporting (aligned to +/- strand):  ref base . (3/0);  new base C (4/5);  total (7/5)
Fisher's exact test for biased strand distribution p-value = 2.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

ACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAG  >  W3110S.gb/1595895‑1596020
                                                                    |                                                          
acCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccc                                                         >  1:191870/1‑71 (MQ=255)
acCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccc                                                         >  1:628946/1‑71 (MQ=255)
acCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccc                                                         >  1:2812498/1‑71 (MQ=255)
acCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccc                                                         >  1:2903769/1‑71 (MQ=255)
acCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccc                                                         >  1:343821/1‑71 (MQ=255)
       ttATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccccccccc                                                   <  1:355519/70‑2 (MQ=255)
                      tgtTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCcgtcg                                       <  1:1920788/68‑1 (MQ=255)
                      tgtTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCcgtcg                                       <  1:2871427/68‑1 (MQ=255)
                      tgtTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCcgtcg                                       <  1:3680605/68‑1 (MQ=255)
                      tgtTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTccc                                                         >  1:3907654/1‑49 (MQ=255)
                          gACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGa                               >  1:1076009/1‑71 (MQ=255)
                          gACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGa                               >  1:2214830/1‑71 (MQ=255)
                          gACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGa                               >  1:502219/1‑71 (MQ=255)
                           aCTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCcgtcg                                       <  1:765294/63‑1 (MQ=255)
                                   tatTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGa                               <  1:2981732/63‑1 (MQ=255)
                                             agagTCGGAACCGAAAGCAATTTCCCCCCCCCCCGATATCcgtc                                        >  1:1919148/1‑44 (MQ=255)
                                             agagTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgttg             >  1:3559048/1‑71 (MQ=255)
                                             agagTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCgttg             >  1:2546221/1‑71 (MQ=255)
                                                         gAAAGCAATTTCCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAg  >  1:1304769/1‑70 (MQ=255)
                                                                    |                                                          
ACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAG  >  W3110S.gb/1595895‑1596020

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