Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,875,612 G→A 28.6% 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,6120GA28.6% 78.6 / 33.8 49D126N (GAT→AAT) iapaminopeptidase in alkaline phosphatase isozyme conversion
Reads supporting (aligned to +/- strand):  ref base G (17/18);  new base A (6/8);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 7.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTA  >  W3110S.gb/2875544‑2875680
                                                                    |                                                                    
cTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCg                                                                      <  1:2943336/69‑1 (MQ=255)
cTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCg                                                                      <  1:529477/69‑1 (MQ=255)
cTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCg                                                                      <  1:476449/69‑1 (MQ=255)
  cATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAAt                                                                    >  1:558267/1‑69 (MQ=255)
  cATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAAt                                                                    >  1:178671/1‑69 (MQ=255)
   aTGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGAt                                                                    <  1:1518897/68‑1 (MQ=255)
     gAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCt                                                                 <  1:1690623/69‑1 (MQ=255)
     gAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCt                                                                 <  1:3179284/69‑1 (MQ=255)
     gAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCt                                                                 <  1:1917766/69‑1 (MQ=255)
     gAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCt                                                                 <  1:2256365/69‑1 (MQ=255)
     gAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCt                                                                 <  1:2664622/69‑1 (MQ=255)
      aaGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTg                                                                >  1:1880077/1‑69 (MQ=255)
              gcgcCGCAGCAGATCATCATTATGGCGCATCTTGATACTTACGCCCCGCTGAGCGa                                                                     <  1:3337308/56‑1 (MQ=255)
              gcgcCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGa                                                                     <  1:3617497/56‑1 (MQ=255)
              gcgcCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGa                                                                     <  1:3073045/56‑1 (MQ=255)
                gccgcAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGc                                                      <  1:2523575/69‑1 (MQ=255)
                gccgcAGCAGATCATCATTATGGCGAATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGc                                                      <  1:3595151/69‑1 (MQ=255)
                          atcatcATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTcggc                                            <  1:3155090/69‑1 (MQ=255)
                           tcatcaTTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTcggcg                                           >  1:363657/1‑69 (MQ=255)
                            catcatTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCt                                                                 <  1:1164571/46‑1 (MQ=255)
                             atcatTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCt                                                                 >  1:2703660/1‑45 (MQ=255)
                               catTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAAtctc                                               >  1:1808850/1‑61 (MQ=255)
                               catTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCt                                       >  1:3415699/1‑69 (MQ=255)
                               catTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCt                                       <  1:514859/69‑1 (MQ=255)
                               catTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCt                                       <  1:725746/69‑1 (MQ=255)
                               catTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCt                                       >  1:2011926/1‑69 (MQ=255)
                                atTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTg                                      >  1:3008096/1‑69 (MQ=255)
                                atTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTg                                      >  1:2833572/1‑69 (MQ=255)
                                atTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTg                                      >  1:75249/1‑69 (MQ=255)
                                atTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTg                                      >  1:667176/1‑69 (MQ=255)
                                atTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTg                                      >  1:1738/1‑69 (MQ=255)
                                atTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTg                                      >  1:1635746/1‑69 (MQ=255)
                                   aTGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGAc                                    >  1:617040/1‑68 (MQ=255)
                                   aTGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGAc                                    >  1:3360291/1‑68 (MQ=255)
                                           tCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAgg                           <  1:2900013/69‑1 (MQ=255)
                                           tCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAgg                           >  1:3303735/1‑69 (MQ=255)
                                           tCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAgg                           <  1:3450820/69‑1 (MQ=255)
                                           tCTGGATACTTACGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAgg                           <  1:36648/69‑1 (MQ=255)
                                             tGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGaa                         >  1:243102/1‑69 (MQ=255)
                                             tGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGaa                         >  1:1219580/1‑69 (MQ=255)
                                                 tACTTACGCCCCGCTGAGCGATGCTGACGCCGATTCCAATCTCGGCGGGCTGACGTTACAAGGAATg                       <  1:58640/67‑1 (MQ=255)
                                                     tACGCCGCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatgat                 <  1:2643932/69‑1 (MQ=255)
                                                     tACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatgat                 <  1:366503/69‑1 (MQ=255)
                                                     tACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatgat                 <  1:868447/69‑1 (MQ=255)
                                                       cGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGTCTGACGTTACAAGGAATGgat                    >  1:1435277/1‑64 (MQ=255)
                                                       cGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgat                    <  1:165721/64‑1 (MQ=255)
                                                       cGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgat                    <  1:2828305/64‑1 (MQ=255)
                                                       cGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgat                    >  1:3460001/1‑64 (MQ=255)
                                                       cGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgat                    >  1:3670687/1‑64 (MQ=255)
                                                       cGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgat                    <  1:2824026/64‑1 (MQ=255)
                                                       cGCCCCGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgat                    >  1:2834733/1‑64 (MQ=255)
                                                            cGCTGAGCAATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAAcgccg          <  1:937004/69‑1 (MQ=255)
                                                                    gATGCTGACGCCGATGCCAATCTCGGCGTGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTa  <  1:377312/69‑1 (MQ=255)
                                                                    gATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatga                  <  1:540027/53‑1 (MQ=255)
                                                                    gATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatga                  <  1:1747703/53‑1 (MQ=255)
                                                                    gATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGgatga                  <  1:2803764/53‑1 (MQ=255)
                                                                    gATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTa  <  1:1043418/69‑1 (MQ=255)
                                                                    gATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTa  <  1:3637460/69‑1 (MQ=255)
                                                                    gATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTa  <  1:327702/69‑1 (MQ=255)
                                                                    |                                                                    
CTCATGAAGGCAAAGCGCCGCAGCAGATCATCATTATGGCGCATCTGGATACTTACGCCCCGCTGAGCGATGCTGACGCCGATGCCAATCTCGGCGGGCTGACGTTACAAGGAATGGATGATAACGCCGCAGGTTTA  >  W3110S.gb/2875544‑2875680

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: