Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,307,310 C→T 77.8% V24V (GTC→GTT eco → ecotin, a serine protease inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,307,3100CT77.8% 48.0 / 12.1 27V24V (GTC→GTTecoecotin, a serine protease inhibitor
Reads supporting (aligned to +/- strand):  ref base C (4/2);  new base T (7/14);  total (11/16)
Fisher's exact test for biased strand distribution p-value = 1.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

GACCATTCTACCTGCAGTATTGTTTGCCGCTTTCGCTACCACTTCCGCCTGGGCGGCAGAAAGCGTCCAGCCACTGGAAAAAATCGCGCCTTATCCACAAGCTGAAAAAGGGATGAAGCGTCAGGTGATTCAG  >  W3110S.gb/2307244‑2307376
                                                                  |                                                                  
gACCATTCTACCTGCAGTATTGTTTGCCGCTTTCGCTACCACTTCCGCCTGGGCGGCAGAAAGCGTTCAGc                                                                <  1:1382767/71‑1 (MQ=255)
  ccATTCTACCTGCAGTATTGTTTGCCGCTTTCGCTACCACTTCCGCCTGGGCGGCAGAAAGCGTTCAGCCa                                                              <  1:1206034/71‑1 (MQ=255)
               aGTATTGTTTGCCGCTTTCGCTACCACTTCCGCCTGGGCGGCAGAAAGCGTTCAGCCACTGGAAAAAAt                                                   <  1:2077152/69‑1 (MQ=255)
                   ttGTTTGCCGCTTTCGCTACCACTTCCGCCTGGGCGGCAGAAAGCGTTCAGCCACTGGAAAAAATCGCGcc                                             >  1:1546882/1‑71 (MQ=255)
                   ttGTTTGCCGCTTTCGCTACCACTTCCGCCTGGGCGGCAGAAAGCGTTCAGCCACTGGAAAAAATCGCGcc                                             >  1:733888/1‑71 (MQ=255)
                   ttGTTTGCCGCTTTCGCTACCACTTCCGCCTGGGCGGCAGAAAGCGTCCAGCCACTGGaaaaa                                                     <  1:2098556/63‑1 (MQ=255)
                     gTTTGCCGCTTTCGCTACCACTTCCGCCTGGGCGGCAGAAAGCGTTCAGCCACTGGAAAAAATCGCGCCt                                            >  1:12417/1‑70 (MQ=255)
                      tttGCCGCTTTCGCTACCACTTCCGCCTTGGCGGCAGAAAGCGTTCAGCCACTGGAAAAAATCGCGCCtt                                           <  1:1295096/70‑1 (MQ=255)
                      tttGCCGCTTTCGCTACCACTTCCGCCTGGGCGGCAGAAAGCGTTCAGCCACTGGaa                                                        <  1:317805/57‑1 (MQ=255)
                      tttGCCGCTTTCGCTACCACTTCCGCCTGGGCGGCAGAAAGCGTTCAGCCACTGGAAAAAATCGCGCCtt                                           <  1:1542188/70‑1 (MQ=255)
                      tttGCCGCTTTCGCTACCACTTCCGCCTGGGCGGCAGAAAGCGTTCAGCCACTGGAAAAAATCGCGCCtt                                           <  1:2329667/70‑1 (MQ=255)
                        tGCCGCTTTCGCTACCACTTCCGCCTGGGCGGCAGAAAGCGTTCAGCCACTGGAAAAAATCGCGCCTTATc                                        >  1:1680457/1‑71 (MQ=255)
                             cTTTCGCTACCACTTCCGCCTGGGCGGCAGAAAGCGTTCAGCCACTGGAAAAAATCGCGCCTTATCCACaa                                   <  1:495299/71‑1 (MQ=255)
                                   cTACCACTTCCGCCTGGGCGGCAGAAAGCGTTCAGCCACTGGAAAAAATCGCGCCTTATCCACAAGCTGaa                             <  1:2469568/71‑1 (MQ=255)
                                      ccACTTCCGCCTGGGCGGCAGAAAGCGTTCAGCCACTGGAAAAAATcgcg                                               <  1:2391455/50‑1 (MQ=255)
                                                 tGGGCGGCAGAAAGCGTCCAGCCACTGGAAAAAATCGCGCCTTATCCACAAGCTGAAAAAGGGATGAAGCg               >  1:1500515/1‑71 (MQ=255)
                                                   ggcggcAGAAAGCGTTCAGCCACTGGAAAAAATCGCGCCTTATCCACAAGCTGAAAAAGGGATGAAGCGTc             <  1:1474730/71‑1 (MQ=255)
                                                       gcAGAAAGCGTCCAGCCACTGGAAAAAATCGCGCCtt                                           <  1:1960251/37‑1 (MQ=255)
                                                         aGAAAGCGTCCAGCCACTGGAAAAAATCGCGCCTTATCCACAAGCTGAAAAAGGGATGAAGCGTCAGGTGa       >  1:1905181/1‑71 (MQ=255)
                                                           aaaGCGTTCAGCCACTGGAAAAAATCGCGCCTTATCCACAAGCTGAAAAAgg                        <  1:2172572/52‑1 (MQ=255)
                                                           aaaGCGTTCAGCCACTGGAAAAAATCGCGCCTTATCCACAAGCTGAAAAAgg                        <  1:1007251/52‑1 (MQ=255)
                                                           aaaGCGTTCAGCCACTGGAAAAAATCGCGCCTTATCCACAAGCTGAAAAAgg                        <  1:2485057/52‑1 (MQ=255)
                                                             aGCGTTCAGCCACTGGAAAAAATCGCGCCTTATCCACAAGCTGaaaaa                          >  1:285529/1‑48 (MQ=255)
                                                             aGCGTTCAGCCACTGGAAAAAATCGCGCCTTATCCACAAGCTGAAAAAGGGATGAAGCGTCAg           >  1:188977/1‑63 (MQ=255)
                                                             aGCGTTCAGCCACTGGAAAAAATCGCGCCTTATCCACAAGCTGAAAAAGGGATGAAGCGTCAGGTGAtt     >  1:347738/1‑69 (MQ=255)
                                                             aGCGTCCAGCCACTGGAAAAAATCGCGCCTTATCCACAAGCTGAAAAAGGGATGAAGCGTCAGGTGAtt     >  1:1470553/1‑69 (MQ=255)
                                                              gCGTCCAGCCACTGGAAAAAATCGCGCCTTATCCACAAGCTGAAAAAGGGATGAAGCGTCAGGTGATTCAg  >  1:818700/1‑71 (MQ=255)
                                                                  |                                                                  
GACCATTCTACCTGCAGTATTGTTTGCCGCTTTCGCTACCACTTCCGCCTGGGCGGCAGAAAGCGTCCAGCCACTGGAAAAAATCGCGCCTTATCCACAAGCTGAAAAAGGGATGAAGCGTCAGGTGATTCAG  >  W3110S.gb/2307244‑2307376

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: