Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,053,471 A→G 100% F9S (TTT→TCT)  pepP ← proline aminopeptidase P II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,053,4710AG86.7% 36.4 / ‑0.4 15F9S (TTT→TCT) pepPproline aminopeptidase P II
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (10/3);  total (10/5)
Fisher's exact test for biased strand distribution p-value = 9.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAG  >  W3110S.gb/3053426‑3053530
                                             |                                                           
cTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCAt                                    <  1:1731816/71‑1 (MQ=255)
cTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCAt                                    <  1:1140919/71‑1 (MQ=255)
cTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCAt                                    <  1:361749/71‑1 (MQ=255)
               tGCTCCACCAGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACg                     >  1:986777/1‑71 (MQ=255)
                        aGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGtt            >  1:1382072/1‑71 (MQ=255)
                        aGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGtt            >  1:347930/1‑71 (MQ=255)
                          ggCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGtttt          >  1:1483168/1‑71 (MQ=255)
                          ggCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGtttt          >  1:627392/1‑71 (MQ=255)
                          ggCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGtttt          >  1:693193/1‑71 (MQ=255)
                           gCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATa                                   <  1:775561/45‑1 (MQ=255)
                                aCGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGttt           >  1:548481/1‑64 (MQ=255)
                                aCGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAgtg     >  1:215111/1‑70 (MQ=255)
                                aCGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAgtg     >  1:377377/1‑70 (MQ=255)
                                aCGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAgtg     >  1:92831/1‑70 (MQ=255)
                                  ggCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTag  <  1:612307/71‑1 (MQ=255)
                                             |                                                           
CTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAG  >  W3110S.gb/3053426‑3053530

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: ‑