Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,427,357 A→T 19.0% D8V (GAT→GTT)  ydaY → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,427,3570AT19.0% 43.9 / 6.0 21D8V (GAT→GTT) ydaYhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (11/6);  new base T (2/2);  total (13/8)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGT  >  W3110S.gb/1427289‑1427423
                                                                    |                                                                  
cTCAAAATGGGTTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGAt                                                                   >  1:1012634/1‑70 (MQ=255)
cTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATc                                                                  <  1:2409437/71‑1 (MQ=255)
      aTGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAg                                                                >  1:849127/1‑67 (MQ=255)
      aTGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATc                                                            >  1:2385459/1‑71 (MQ=255)
            tCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGTTCAGTATCGCTccc                                                      >  1:2208975/1‑71 (MQ=255)
                ttCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGTTCAGTATCGCTCCCGt                                                    <  1:650096/69‑1 (MQ=255)
                ttCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGt                                                    <  1:1130367/69‑1 (MQ=255)
                ttCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGt                                                    <  1:2408780/69‑1 (MQ=255)
                        aTTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTa                                           <  1:2045595/70‑1 (MQ=255)
                           gTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATcc                                       >  1:1028044/1‑71 (MQ=255)
                            ttCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGTTCAGTATCGCTCCCGTAATGCAtt                                            <  1:1904999/65‑1 (MQ=255)
                              ctgtctgGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTc                                    >  1:1406613/1‑71 (MQ=255)
                              ctgtctgGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTc                                    >  1:1371291/1‑71 (MQ=255)
                                             gATGAGCAGAAGCAGCGATAACGTTCAGTATCGCTCCCGTAATGCATTAATCCGTCGCcaca                              >  1:1073722/1‑62 (MQ=255)
                                               tGAGCAGAAGCAGAGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATg                   <  1:82827/71‑1 (MQ=255)
                                                     gAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGc              >  1:1044899/1‑70 (MQ=255)
                                                          agcGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGa                 >  1:1130345/1‑62 (MQ=255)
                                                          agcGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGa                 >  1:488217/1‑62 (MQ=255)
                                                           gcGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTg       <  1:1364163/71‑1 (MQ=255)
                                                                aaCGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGt  >  1:1823498/1‑71 (MQ=255)
                                                                 aCGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCa     >  1:1465413/1‑67 (MQ=255)
                                                                    |                                                                  
CTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGT  >  W3110S.gb/1427289‑1427423

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: