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,4710AG94.0% 100.1 / ‑2.3 33F9S (TTT→TCT) pepPproline aminopeptidase P II
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (21/10);  total (22/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTT  >  W3110S.gb/3053404‑3053522
                                                                   |                                                   
tGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAGAct                                                  <  1:1110858/71‑1 (MQ=255)
         aaTCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCggg                                         <  1:1362232/71‑1 (MQ=255)
          aTCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGa                                        <  1:2406307/71‑1 (MQ=255)
            cAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGata                                      >  1:1556747/1‑71 (MQ=255)
                              ttGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACActct                    <  1:954600/71‑1 (MQ=255)
                                gCATTTGCTCCACCAGGGCCTGACGGCGACGCTTAGActct                                                >  1:1305455/1‑41 (MQ=38)
                                   tttGCTCCACCAGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCt                 <  1:1353311/69‑1 (MQ=255)
                                     tGCTCCACCAGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGa                                        >  1:1577386/1‑44 (MQ=255)
                                     tGCTCCACCAGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGa                                        >  1:421874/1‑44 (MQ=255)
                                     tGCTCCACCAGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGa                                        >  1:2226925/1‑44 (MQ=255)
                                     tGCTCCACCAGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATATctcactc                              >  1:644590/1‑54 (MQ=255)
                                     tGCTCCACCAGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTa               >  1:1599380/1‑69 (MQ=255)
                                     tGCTCCACCAGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTa               >  1:2032020/1‑69 (MQ=255)
                                     tGCTCCACCAGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTAc              >  1:536685/1‑70 (MQ=255)
                                     tGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACg             >  1:124984/1‑71 (MQ=255)
                                      gCTCCACCAGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGatat                                     >  1:277048/1‑46 (MQ=255)
                                         ccaccaGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGtttt         >  1:161283/1‑71 (MQ=255)
                                          caccaGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGttttt        <  1:1918854/71‑1 (MQ=255)
                                            ccaGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTg      >  1:2024871/1‑71 (MQ=255)
                                             caGGGCCTGACGGCGACGCTGAGACTCTTGCCggg                                         <  1:296390/35‑1 (MQ=255)
                                             caGGGCCTGACGGCGACGCTGAGACTCTTGCCggg                                         <  1:1601643/35‑1 (MQ=255)
                                             caGGGCCTGACGGCGACGCTGAGACTCTTGCCggg                                         <  1:1430878/35‑1 (MQ=255)
                                             aaGGGCCTGACGGCGACGCTGAGACTCTTGCCggg                                         <  1:2410689/34‑1 (MQ=38)
                                              aGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGtt    >  1:2103514/1‑71 (MQ=255)
                                              aGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGtt    >  1:2100508/1‑71 (MQ=255)
                                              aGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGtt    >  1:1645823/1‑71 (MQ=255)
                                               gggCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGtt           <  1:1059019/63‑1 (MQ=255)
                                                  ccTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGtttt  >  1:2510170/1‑69 (MQ=255)
                                                  ccTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGtttt  >  1:948830/1‑69 (MQ=255)
                                                  ccTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGtttt  >  1:2015646/1‑69 (MQ=255)
                                                      aCGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATaa                          >  1:367795/1‑41 (MQ=255)
                                                      aCGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATaa                          >  1:251128/1‑41 (MQ=255)
                                                      aCGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATaa                          >  1:2110601/1‑41 (MQ=255)
                                                                   |                                                   
TGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTT  >  W3110S.gb/3053404‑3053522

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: