Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,053,087 T→A 29.7% Y137F (TAT→TTT)  pepP ← proline aminopeptidase P II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,053,0870TA29.7% 70.4 / 35.9 47Y137F (TAT→TTT) pepPproline aminopeptidase P II
Reads supporting (aligned to +/- strand):  ref base T (14/19);  new base A (8/6);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 5.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TGCGGTGAGATTTTGCCGCGAACCTTTACGCAGTTTTTCCAGCGCACTGTTCACGATTACATCAGCATATGCATATTCGCCCTGGGCATGGTAAACCACATCCAGGCCGTTAAGTAGTTGATAAAGTTGCTGATTGAT  >  W3110S.gb/3053020‑3053157
                                                                   |                                                                      
tgcggtgAGATTTTGCCGCGAACCTTTACGCAGTTTTTCCAGCGCACTGTTCACGATTACATCAgcatatg                                                                     <  1:1885360/71‑1 (MQ=255)
  cggtgAGATTTTGCCGCGTACCTTTACGCAGTTTTTCCAGCGCACTGTTCACGATTACATCAgcatatgca                                                                   <  1:1948994/71‑1 (MQ=255)
  cggtgAGATTTTGCCGCGAACCTTTACGCAGTTTTTCCAGCGCACTGTTCACGATTACATCAgcatatgca                                                                   >  1:3094693/1‑71 (MQ=255)
  cggtgAGATTTTGCCGCGAACCTTTACGCAGTTTTTCCAGCGCACTGTTCACGATTACATCAgcatatgca                                                                   >  1:1548510/1‑71 (MQ=255)
  cggtgAGATTTTGCCGCGAACCTTTACGCAGTTTTTCCAGCGCACTGTTCACGATTACATCAgcatatgca                                                                   >  1:1622733/1‑71 (MQ=255)
  cggtgAGATTTTGCCGCGAACCTTTACGCAGTTTTTCCAGCGCACTGTTCACGATTACATCAgcatatgca                                                                   >  1:3752664/1‑71 (MQ=255)
           tttGCCGCGAACCTTTACGCAGTTTTTCCAGCGCACTGTTCACGATTACATCAgcata                                                                       >  1:1355928/1‑58 (MQ=255)
           tttGCCGCGAACCTTTACGCAGTTTTTCCAGCGCACTGTTCACGATTACATCAGCAAATGCATATTCGccc                                                          >  1:475115/1‑71 (MQ=255)
           tttGCCGCGAACCTTTACGCAGTTTTTCCAGCGCACTGTTCACGATTACATCAGCAAATGCATATTCGccc                                                          >  1:1926650/1‑71 (MQ=255)
            ttGACGCGAACCTTTACGCAGTTTTTCCAGCGCACTGTTCACGATTACATCAGCATATGCATATTAGCCCt                                                         >  1:2512147/1‑71 (MQ=255)
                  cgAACCTTTACGCAGTTTTTCCAGCGCACTGTTCACGATTACATCAgcatatg                                                                     <  1:3565007/53‑1 (MQ=255)
                      ccTTTACGCAGTTTTTCCAGCGCACTGTTCACGATTACATCAGCATATGCATATTCGCCCTGGGCATGGt                                                >  1:3418881/1‑70 (MQ=255)
                      ccTTTACGCAGTTTTTCCAGCGCACTGTTCACGATTACATCAGCATATGCATATTCGCCCTGGGCATGGTa                                               <  1:3147513/71‑1 (MQ=255)
                      ccTTTACGCAGTTTTTCCAGCGCACTGTTCACGATTACATCAGCATATGCATATTCGCCCTGGGCATGGTa                                               <  1:3234897/71‑1 (MQ=255)
                        tttACGCAGTTTTTCCAGCGCACTGTTCACGATTACATCAGCATATGCATATTCGCCCTggg                                                      <  1:2259841/62‑1 (MQ=255)
                          tACGCAGTTTTTCCAGCGCACTGTTCACGATTACATCAGCATATGCATATTCGCCCTGGGCATGGTaaa                                             >  1:2698254/1‑69 (MQ=255)
                               aGTTTTTCCAGCGCACTGTTCACGATTACATCAGCATATGCATATTCGCCCTGGGCATGGTAAACCACATc                                      <  1:455364/71‑1 (MQ=255)
                               aGTTTTTCCAGCGCACTGTTCACGATTACATCAGCAAATGCATATTCGCCCTg                                                        >  1:2047341/1‑53 (MQ=255)
                                 tttttCCAGCGCACTGTTCACGATTACATCAgcat                                                                        <  1:320915/35‑1 (MQ=255)
                                 tttttCCAGCGCACTGTTCACGATTACATCAGCATATGCATATTCGCCCTGGGCATGGTAAACCACATcc                                     <  1:1629379/70‑1 (MQ=255)
                                 tttttCCAGCGCACTGTTCACGATTACATCAGCATATGCATATTCGCCCTGGGCATGGTAAACCACATcc                                     <  1:354897/70‑1 (MQ=255)
                                 tttttCCAGCGCACTGTTCACGATTACATCAGCATATGCATATTCGCCCTGGGCATGGTAAACCACATcc                                     <  1:3841753/70‑1 (MQ=255)
                                 tttttCCAGCGCACTGTTCACGATTACATCAGCATATGCATATTCGCCCTGGGCATGGTAAACCACATcc                                     <  1:502776/70‑1 (MQ=255)
                                 tttatCCAGCGCACTGTTCACGATTACATCAGCATATGCATATTCGCCCTGGGCATGGTAAACCACATcc                                     <  1:2558724/70‑1 (MQ=255)
                                   tttCCAGCGCACTGTTCACGATTACATCAGCATATGCATATTCGCCCTGGGCATGGTAAACCACATc                                      >  1:1394703/1‑67 (MQ=255)
                                   tttCCAGCGCACTGTTCACGATTACATCAGCATATGCATATTCGCCCTGGGCATGGTAAACCACATc                                      >  1:2810377/1‑67 (MQ=255)
                                   tttCCAGCGCACTGTTCACGATTACATCAGCATATGCATATTCGCCCTGGGCATGGTAAACCACATCCAgg                                  >  1:3085803/1‑71 (MQ=255)
                                   tttCCAGCGCACTGTTCACGATTACATCAGCAAATGCATATTCGCCCTGGGCATGGTAAACCACATCCAgg                                  >  1:3329839/1‑71 (MQ=255)
                                        aGCGCACTGTTCACGATTACATCAGCATATGCATATTCGCCCTGGGCATGGTAAACCACATCCAGGCCGtt                             <  1:3500674/71‑1 (MQ=255)
                                        aGCGCACTGTTCACGATTACATCAGCATATGCATATTCGCCCTGGGCATGGTAAACCACATCCAGGCCGtt                             <  1:350992/71‑1 (MQ=255)
                                         gcgcACTGTTCACGATTACATCAGCATATGCATATTCGCCCTGGGCATGGTAAACca                                          <  1:2759713/57‑1 (MQ=255)
                                         gcgcACTGTTCACGATTACATCAGCATATGCATATTCGCCCTGGGCATGGTAAACca                                          <  1:3857668/57‑1 (MQ=255)
                                         gcgcACTGTTCACGATTACATCAGCATATGCATATTCGCCCTGGGCATGGTAAACca                                          <  1:2079452/57‑1 (MQ=255)
                                          cgcACTGTTCACGATTACATCAGCAAATGCATATTCGCCCTGGGCATGGTAAACCACATCCAGGCCGtt                             >  1:255326/1‑69 (MQ=255)
                                           gcACTGTTCACGATTACATCAGCATATGCATATTCGCCCTGGGCATGGTAAACCACATCCAGGCCGTTAag                          <  1:3785778/71‑1 (MQ=255)
                                           gcACTGTTCACGATTACATCAGCATATGCATATTCGCCCTGGGCATGGTAAACCACATCCAGGCCGTTAag                          <  1:959456/71‑1 (MQ=255)
                                             aCTGTTCACGATTACATCAGCATATGCATATTCGCCCTGGGCATGGTAAACCACATCCAGGCCGTTAagta                        <  1:3523494/71‑1 (MQ=255)
                                             aCTGTTCACGATTACATCAGCATATGCATATTCGCCCTGGGCATGGTAAACCACATCCAGGCCGTTAagta                        <  1:2075116/71‑1 (MQ=255)
                                              cTGTTCACGATTACATCAGCAAATGCATATTCGCCCTGGGCATGGTaa                                              <  1:2510110/48‑1 (MQ=255)
                                              cTGTTCACGATTACATCAGCAAATGCATATTCGCCCTGGGCATGGTAAACCACATc                                      <  1:3788474/56‑1 (MQ=255)
                                               tGTTCACGATTACATCAGCAAATGCATATTCGCCCTGGGCATGGTAAACCACATc                                      <  1:1979979/55‑1 (MQ=255)
                                                gTTCACGATTACATCAGCATATGCATATTCGCCCTGGGCATGGTAAACCACATCCAGGCCGTTAAGTAGtt                     >  1:2145686/1‑71 (MQ=255)
                                                gTTCACGATTACATCAGCAAATGCATATTCGCCCTGGGCATGGTAAACCACATCCa                                    >  1:2363465/1‑56 (MQ=255)
                                                gTTCACGATTACATCAGCAAATGCATATTCGCCCTGGGCATGGTAAACCACATCCa                                    >  1:2598623/1‑56 (MQ=255)
                                                   cACGATTACATCAGCATATGCATATTCGCCCTGGGCa                                                    >  1:1813313/1‑37 (MQ=255)
                                                   cACGATTACATCAGCATATGCATATTCGCCCTGGGCa                                                    >  1:3744824/1‑37 (MQ=255)
                                                    aCGATTACATCAGCAAATGCATATTCGCCCTGGGCATGGTAAACCACATCCAGGCCGTTAAGTAGTTGATa                 <  1:3171134/71‑1 (MQ=255)
                                                     cGATTACATCAGCATATGCATATTCGCCCTGGGCATGGTAAACCACATCCAGGCCGTTAAGTAGTTGATaa                <  1:335302/71‑1 (MQ=255)
                                                       aTTACATCAGCAAATGCATATTCGCCCTGGGCATGGTAAACCACATCCAGGCCGTTAagtag                       <  1:2624329/62‑1 (MQ=255)
                                                       aTTACATCAGCAAATGCATATTCGCCCTGGGCATGGTAAACCACATCCAGGCCGTTAagtag                       <  1:1726709/62‑1 (MQ=255)
                                                         tACATCAGCATATGCATATTCGCCCTGGGCATGGTAAACCACATCCAGGCCGTTaa                           >  1:3378213/1‑56 (MQ=255)
                                                         tACATCAGCATATGCATATTCGCCCTGGGCATGGTAAACCACATCCAGGCCGTTaa                           >  1:2245411/1‑56 (MQ=255)
                                                         tACATCAGCATATGCATATTCGCCCTGGGCATGGTAAACCACATCCAGGCCGTTaa                           >  1:1247852/1‑56 (MQ=255)
                                                         tACATCAGCAAATGCATATTCGCCCTGGGCATGGTAAACCACATCCAGGCCGTTAAGTAGtt                     >  1:3078735/1‑62 (MQ=255)
                                                             tCAGCATATGCATATTCGCCCTGGGCATGGTAAACCACATCCAGGCCGTTAAGTAGTTGATAAAGTTGCt         >  1:1477368/1‑70 (MQ=255)
                                                             tCAGCATATGCATATTCGCCCTGGGCATGGTAAACCACATCCAGGCCGTTAAGTAGTTGATAAAGTTGCt         >  1:260719/1‑70 (MQ=255)
                                                                 catatgcatatTCGCCCTGGGCATGGTAAACCACATCCAGGCCGTTAAGTAGTTGATAAAGTTGCtgatt     <  1:55340/70‑1 (MQ=255)
                                                                   tatgcatatTCGCCCTGGGCATGGTAAACCACATc                                      <  1:3071427/35‑1 (MQ=255)
                                                                   tatgcatatTCGCCCTGGGCATGGTAAACCACATc                                      <  1:419773/35‑1 (MQ=255)
                                                                   tatgcatatTCGCCCTGGGCATGGTAAACCACATCGAGGCCGTTAAGTAGTTGATAAAGTTGCt         <  1:87917/64‑1 (MQ=255)
                                                                   tatgcatatTCGCCCTGGGCATGGTAAACCACATCCAGGCCGTTAAGTAGTTGATAAAGTTGCtgattgat  >  1:2458983/1‑71 (MQ=255)
                                                                   tatgcatatTCGCCCTGGGCATGGTAAACCACATCCAGGCCGTTAAGTAGTTGATAAAGTTGCtgattgat  >  1:2902148/1‑71 (MQ=255)
                                                                   tatgcatatTCGCCCTGGGCATGGTAAACCACATCCAGGCCGTTAAGTAGTTGATAAAGTTGCtgattgat  >  1:664345/1‑71 (MQ=255)
                                                                   tatgcatatTCGCCCTGGGCATGGTAAACCACATCCAGGCCGTTAAGTAGTTGATAAAGTTGCtgattgat  >  1:3060571/1‑71 (MQ=255)
                                                                   tatgcatatTCGCCCTGGGCATGGTAAACCACATCCAGGCCGTTAAGTAGTTGATAAAGTTGCt         <  1:2225784/64‑1 (MQ=255)
                                                                   tatgcatatTCGCCCTGGGCATGGTAAACCACATCCAGGCCGTTAAGTAGTTGATAAAGTTGCt         <  1:1626926/64‑1 (MQ=255)
                                                                   |                                                                      
TGCGGTGAGATTTTGCCGCGAACCTTTACGCAGTTTTTCCAGCGCACTGTTCACGATTACATCAGCATATGCATATTCGCCCTGGGCATGGTAAACCACATCCAGGCCGTTAAGTAGTTGATAAAGTTGCTGATTGAT  >  W3110S.gb/3053020‑3053157

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: