Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,053,471 A→G 21.9% 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,4710AG21.9% 61.0 / 15.4 32F9S (TTT→TCT) pepPproline aminopeptidase P II
Reads supporting (aligned to +/- strand):  ref base A (14/11);  new base G (3/4);  total (17/15)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

CAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTA  >  W3110S.gb/3053406‑3053529
                                                                 |                                                          
cagcaAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAGActct                                                       >  1:1426299/1‑71 (MQ=255)
        aTCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGa                                               <  1:1637929/71‑1 (MQ=255)
        aTCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGa                                               <  1:2041545/71‑1 (MQ=255)
               cggcgCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAGACTCTTGCCgg                                                 <  1:1399241/62‑1 (MQ=255)
                ggcgCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATctcac                                       >  1:846450/1‑71 (MQ=255)
                  cgcTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATctcact                                      >  1:715716/1‑70 (MQ=255)
                        cGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAac                                <  1:238713/70‑1 (MQ=255)
                            ttGCATTTGCTTCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCgg                                                 <  1:927070/49‑1 (MQ=255)
                              gCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACtctc                          >  1:2047541/1‑70 (MQ=255)
                              gCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGcc                                                   <  1:2548243/45‑1 (MQ=255)
                                 tttGCTCCACCAGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATAtct                                          <  1:1143268/51‑1 (MQ=255)
                                 tttGCTCCACCAGGGCCTGACGGCGACGCTGAAAct                                                         <  1:1131857/36‑1 (MQ=255)
                                 tttGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGa                                               >  1:777473/1‑46 (MQ=255)
                                  ttGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTAc                     >  1:842593/1‑71 (MQ=255)
                                   tGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGa                                               >  1:975040/1‑44 (MQ=255)
                                   tGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGa                                               >  1:447436/1‑44 (MQ=255)
                                   tGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACg                    >  1:1878286/1‑71 (MQ=255)
                                          ccaGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGtttttt              >  1:2212908/1‑70 (MQ=255)
                                          ccaGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCAt                                   <  1:2175644/49‑1 (MQ=255)
                                           caGGGCCTGACGGCGACGCTGAGACTCTTGCCGGGATAtct                                          <  1:1627529/41‑1 (MQ=255)
                                                ccTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGttt                 >  1:1084801/1‑61 (MQ=255)
                                                ccTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGttt                 >  1:332748/1‑61 (MQ=255)
                                                 cTGACGGCGACGCTGAGACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGtttt         <  1:521581/68‑1 (MQ=255)
                                                    aCGGCGACGCTGAAACTCTTGCCGGGATATctcac                                       >  1:676847/1‑35 (MQ=255)
                                                    aCGGCGACGCTGAAACTCTTGCCGGGATATctcac                                       >  1:2525078/1‑35 (MQ=255)
                                                    aCGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGt                   >  1:2170688/1‑55 (MQ=255)
                                                    aCGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGt            >  1:928121/1‑62 (MQ=255)
                                                    aCGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAgtg    <  1:1383562/70‑1 (MQ=255)
                                                     cGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTa  <  1:1722605/71‑1 (MQ=255)
                                                      ggCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAac                                >  1:837399/1‑40 (MQ=255)
                                                              tGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTAc                     <  1:2015519/43‑1 (MQ=255)
                                                               gAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCtt                       <  1:302328/40‑1 (MQ=255)
                                                                 |                                                          
CAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTA  >  W3110S.gb/3053406‑3053529

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: