Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,875,612 G→A 64.1% D126N (GAT→AAT)  iap → aminopeptidase in alkaline phosphatase isozyme conversion

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,875,6120GA64.1% 33.4 / 38.4 39D126N (GAT→AAT) iapaminopeptidase in alkaline phosphatase isozyme conversion
Reads supporting (aligned to +/- strand):  ref base G (10/4);  new base A (14/11);  total (24/15)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

ATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGG  >  W3110S.gb/2875547‑2875676
                                                                 |                                                                
aTGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCt                                                             <  1:1983182/71‑1 (MQ=255)
  gAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGa                                                                 >  1:2329086/1‑65 (MQ=255)
        aaaGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGa                                                     >  1:86363/1‑71 (MQ=255)
        aaaGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGa                                                     >  1:748839/1‑71 (MQ=255)
           gcgcCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGc                                                  >  1:1093075/1‑71 (MQ=255)
           gcgcCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGc                                                  >  1:663274/1‑71 (MQ=255)
            cgccgcAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGcc                                                 <  1:2492148/71‑1 (MQ=255)
                     agATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcggc                                        >  1:911470/1‑71 (MQ=255)
                     agATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcggc                                        <  1:2622448/71‑1 (MQ=255)
                     agATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTcggc                                        >  1:2224822/1‑71 (MQ=255)
                       atcatcATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCt                                                             <  1:434368/48‑1 (MQ=255)
                        tcatcaTTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCggg                                     >  1:1637913/1‑71 (MQ=255)
                         catcatTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTg                                                            >  1:2266001/1‑47 (MQ=255)
                         catcatTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGa                                                           >  1:2348832/1‑48 (MQ=255)
                         catcatTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGa                                                     >  1:869225/1‑54 (MQ=255)
                         catcatTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGc                                    >  1:1142343/1‑71 (MQ=255)
                         catcatTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGc                                    >  1:702868/1‑71 (MQ=255)
                         catcatTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGc                                    >  1:2685941/1‑71 (MQ=255)
                           tcatTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCt                                   <  1:2765176/70‑1 (MQ=255)
                            catTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGa                                 >  1:2303889/1‑71 (MQ=255)
                            catTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGa                                 >  1:78161/1‑71 (MQ=255)
                              ttATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCt                                                             <  1:2419900/41‑1 (MQ=255)
                               tATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGt                              >  1:1465360/1‑71 (MQ=255)
                                aTGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGt                              >  1:302760/1‑70 (MQ=255)
                                   gcgcATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTAc                           <  1:1151853/70‑1 (MQ=255)
                                   gcgcATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTAc                           >  1:102655/1‑70 (MQ=255)
                                          tGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGTCAATCTCGGCGGGCTGACGTTACAAGGAATg                   >  1:122573/1‑71 (MQ=255)
                                            gATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATgg                  <  1:1996213/70‑1 (MQ=255)
                                            gATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATgg                  <  1:413503/70‑1 (MQ=255)
                                            gATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATg                   <  1:1436843/69‑1 (MQ=255)
                                            gATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATg                   <  1:821487/69‑1 (MQ=255)
                                                cTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAgg                       <  1:1558599/61‑1 (MQ=255)
                                                cTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatgat             <  1:593767/71‑1 (MQ=255)
                                                cTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatgat             <  1:2397945/71‑1 (MQ=255)
                                                  tACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAgg                       <  1:1953803/59‑1 (MQ=255)
                                                  tACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATaa           >  1:2220413/1‑71 (MQ=255)
                                                   aCGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAc          >  1:1183603/1‑71 (MQ=255)
                                                         cGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcgccg      >  1:964874/1‑69 (MQ=255)
                                                           cTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAgg  >  1:226660/1‑71 (MQ=255)
                                                                 |                                                                
ATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGG  >  W3110S.gb/2875547‑2875676

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: