Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,189,339 A→G 100% L160S (TTA→TCA)  ycfD ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,189,3390AG95.7% 69.4 / ‑3.8 23L160S (TTA→TCA) ycfDconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (15/7);  total (16/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AAGAATATCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTAACAGATCCGGATGTGGGCAGTGCTGTTTCATTTGCAGCTTTTCGCCCACTCGCCAGCGACGACGTCC  >  W3110S.gb/1189281‑1189406
                                                          |                                                                   
aaGAATATCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTGACAGATCCGGAt                                                         >  1:1845622/1‑71 (MQ=255)
           cAGGCTCCAGCTCTTCATCGATGATGGCTACGAACGGATCGACCTGTGACAGATCCGGATGTGGGCAGTGc                                              <  1:1389224/71‑1 (MQ=255)
                tCCAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTGACAGATCCGGATGTGGGc                                                   <  1:701288/61‑1 (MQ=255)
                tCCAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTGACAGATCCGGATGTGGGCAGTGCTGttt                                         >  1:2489062/1‑71 (MQ=255)
                 ccAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTGACAGATCCGGa                                                          >  1:946681/1‑53 (MQ=255)
                 ccAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTGACAGATCCGGa                                                          >  1:2385823/1‑53 (MQ=255)
                     ctctTCATCGATGATGGCTTCGAACGGATCGACCTGTGACAGATCCGGATGTGGGCAGTGCTGTTTc                                        <  1:2185872/67‑1 (MQ=255)
                          cATCGATGATGGCTTCGAACGGATCGACCTGTGACAGATCCGGATGTGGGCAGTGCTGTTTCATTTGCAGc                               >  1:2273910/1‑71 (MQ=255)
                          cATCGATGATGGCTTCGAACGGATCGACCTGTGACAGATCCGGATGTGGGCAGTGCTGTTTCACTTGCAGc                               >  1:2554057/1‑71 (MQ=255)
                            tCGATGATGGCTTCGAACGGATCGACCTGTGACAGATCCGGa                                                          <  1:2317929/42‑1 (MQ=255)
                              gatgatGGCTTCGAACGGATCGCCCTGTGACAGATCCGGATGTGGGCAGTGCTGTTTCATTTGCAGCtttt                           >  1:2662597/1‑71 (MQ=255)
                              gatgatGGCTTCGAACGGATCGACCTGTGACAGATCCGGATGTGGGCAGTGCTGTTTCATTTGCAGCttt                            >  1:714349/1‑70 (MQ=255)
                               atgatgGCTTCGAACGGATCGACCTGTGACAGATCCGGa                                                          <  1:1649049/39‑1 (MQ=255)
                                  atgGCTTCGAACGGATCGACCTGTGACAGATCCGGAtgt                                                       >  1:1119540/1‑39 (MQ=255)
                                  atgGCTTCGAACGGATCGACCTGTGACAGATCCGGATGTGGGCAGTGCTGTTTCATTTGCAGCTTTTCGc                        >  1:1369685/1‑70 (MQ=255)
                                    ggCTTCGAACGGATCGACCTGTGACAGATCCGGATGTGGGCAGTGCTGTTTCAt                                      >  1:629813/1‑54 (MQ=255)
                                           aaCGGATCGACCTGTGACAGATCCGGATGTGGGCAGTGCTGTTTCATTTGCAGCTTTTc                          <  1:2775070/59‑1 (MQ=255)
                                              ggATCGACCTGTGACAGATCCGGATGTGGGCAGTGCTGTTTCAt                                      >  1:522393/1‑44 (MQ=255)
                                              ggATCGACCTGTGACAGATCCGGATGTGGGCAGTGCTGTTTCAt                                      >  1:62011/1‑44 (MQ=255)
                                              ggATCGACCTGTGACAGATCCGGATGTGGGCAGTGCTGTTTCAt                                      >  1:979737/1‑44 (MQ=255)
                                                  cGACCTGTGACAGATCCGGATGTGGGCAGTGCTGTTTCATTTGCAGCTCTTCGCCCACTCGCCAGcgacga       >  1:2139846/1‑71 (MQ=255)
                                                      cTGTGACAGATCCGGATGTGGGCAGTGCTGTTTCATTTGCAGCTTTTCGCCCACTCGCCAGc            <  1:291168/62‑1 (MQ=255)
                                                       tGTAACAGATCCGGATGTGGGCAGTGCTGTTTCATTTGCAGCTTTTCGCCCACTCGCCAGCGACGACGTcc  >  1:548149/1‑71 (MQ=255)
                                                          |                                                                   
AAGAATATCACCAGGCTCCAGCTCTTCATCGATGATGGCTTCGAACGGATCGACCTGTAACAGATCCGGATGTGGGCAGTGCTGTTTCATTTGCAGCTTTTCGCCCACTCGCCAGCGACGACGTCC  >  W3110S.gb/1189281‑1189406

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: