Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,625,476 T→C L635L (TTA→TTG)  dcp ← dipeptidyl carboxypeptidase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,625,4760TC96.5% 206.7 / ‑2.1 57L635L (TTA→TTG) dcpdipeptidyl carboxypeptidase II
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (24/31);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

TTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTTAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGTAATAACCTGCAGCATATCCGCCACCAAAAATATGGGCGAAATAACTGCTGCGATAGCGTGGCGGTA  >  NC_000913/1625342‑1625609
                                                                                                                                      |                                                                                                                                     
ttCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATc                                                                                                                                   >  1:715925/1‑139 (MQ=255)
ttCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCT‑GAAAGGATCGCCCCCCGAAAACGCAACCCATTTTCCCGCGTCAATcc                                                                                                                                  >  1:1031595/1‑139 (MQ=255)
 tCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATcc                                                                                                                                  >  1:61282/1‑139 (MQ=255)
  cAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCg                                                                                                                                 >  2:176996/1‑139 (MQ=255)
  cAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCg                                                                                                                                 >  2:917435/1‑139 (MQ=255)
     gCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGccc                                                                                                                              >  2:813399/1‑139 (MQ=255)
          gATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTc                                                                                                                         >  1:56914/1‑139 (MQ=255)
                  agcaTCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAAcc                                                                                                                 <  1:562663/139‑1 (MQ=255)
                   gcaTCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCa                                                                                                                <  1:478438/139‑1 (MQ=255)
                   gcaTCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCa                                                                                                                <  1:367227/139‑1 (MQ=255)
                      tCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTTAATCCGcc                                                                                                                               >  2:159440/1‑121 (MQ=255)
                      tCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTTAATCCGcc                                                                                                                               <  1:159440/121‑1 (MQ=255)
                       cGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGa                                                                                                            >  1:483875/1‑139 (MQ=255)
                       cGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGa                                                                                                            <  1:896761/139‑1 (MQ=255)
                                 ttAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACcatcat                                                                                                  <  2:426067/139‑1 (MQ=255)
                                       gCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCa                                                                                            >  1:263782/1‑139 (MQ=255)
                                         cTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCAGCCAAc                                                                                          >  2:86712/1‑139 (MQ=255)
                                                   cATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTcc                                                                                >  1:1081702/1‑139 (MQ=255)
                                                         cGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGAt                                                                          <  2:263782/139‑1 (MQ=255)
                                                          gATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATa                                                                         <  1:162410/139‑1 (MQ=255)
                                                           aTACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGGCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATaa                                                                        <  1:680668/139‑1 (MQ=255)
                                                           aTACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATaa                                                                        <  1:917435/139‑1 (MQ=255)
                                                               aGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGt                                                                    >  1:482947/1‑139 (MQ=255)
                                                                ggCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGta                                                                   <  1:813399/139‑1 (MQ=255)
                                                                  cGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGtaat                                                                 >  1:604453/1‑139 (MQ=255)
                                                                  cGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCAACATCGGCCAACATTTGCGTCCACAGATAAGCGtaat                                                                 <  2:61282/139‑1 (MQ=255)
                                                                   gTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGtaata                                                                <  2:705442/139‑1 (MQ=255)
                                                                   gTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGtaata                                                                <  2:604453/139‑1 (MQ=255)
                                                                    ttCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGtaataa                                                               >  1:322304/1‑139 (MQ=255)
                                                                         gATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGTAATAACCTgc                                                          <  1:596967/139‑1 (MQ=255)
                                                                               cGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGcc                                                                                                                               <  2:982929/64‑1 (MQ=255)
                                                                               cGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGcc                                                                                                                               >  1:982929/1‑64 (MQ=255)
                                                                                    ttACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTcaaca                                                                                                                     <  2:662336/69‑1 (MQ=255)
                                                                                    ttACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTcaaca                                                                                                                     >  1:662336/1‑69 (MQ=255)
                                                                                      aCCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGTAATAACCTGCAGCATATCCGcca                                             <  2:483875/139‑1 (MQ=255)
                                                                                         tctGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGTAATAACCTGCAGCATATCCGccacca                                          >  1:374225/1‑139 (MQ=255)
                                                                                            ggAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGTAATAACCTGCAGCATATCCGCCACCaaaa                                       >  1:46597/1‑139 (MQ=255)
                                                                                                  gATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTg                                                                                                                            >  2:1002044/1‑48 (MQ=255)
                                                                                                  gATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTg                                                                                                                            <  1:1002044/48‑1 (MQ=255)
                                                                                                  gATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGTAATAACCTGCAGCATATCCGCCACCAAAAATATgg                                 <  1:156715/139‑1 (MQ=255)
                                                                                                  gATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGTAATAACCTGCAGCATATCCGCCACCAAAAATATgg                                 >  2:434333/1‑139 (MQ=255)
                                                                                                  gATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGTAATAACCTGCAGCATATCCGCCACCAAAAATATgg                                 >  1:119601/1‑139 (MQ=255)
                                                                                                  gATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGTAATAACCTGCAGCATATCCGCCACCAAAAATATgg                                 >  2:156715/1‑139 (MQ=255)
                                                                                                       ccTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTcaa                                                                                                                       >  2:437728/1‑48 (MQ=255)
                                                                                                       ccTCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTcaa                                                                                                                       <  1:437728/48‑1 (MQ=255)
                                                                                                         tCGCGAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGTAATAACCTGCAGCATATCCGCCACCAAAAATATGGGCGAAAt                          <  2:46597/139‑1 (MQ=255)
                                                                                                            cgAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGTAATAACCTGCAGCATATCCGCCACCAAAAATATGGGCGAAATAAc                       <  2:482947/139‑1 (MQ=255)
                                                                                                            cgAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGTAATAACCTGCAGCATATCCGCCACCAAAAATATGGGCGAAATAAc                       >  1:61048/1‑139 (MQ=255)
                                                                                                            cgAAAACGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGTAATAACCTGCAGCATATCCGCCACCAAAAATATGGGCGAAATAAc                       <  2:715925/139‑1 (MQ=255)
                                                                                                                  cGCAGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGTAATAACCTGCAGCATATCCGCCACCAAAAATATGGGCGAAATAACtgctgc                 <  2:322304/139‑1 (MQ=255)
                                                                                                                     aGCCCATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGTAATAACCTGCAGCATATCCGCCACCAAAAATATGGGCGAAATAACTGCTGCGAt              >  1:413635/1‑139 (MQ=255)
                                                                                                                       cccATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGTAATAACCTGCAGCATATCCGCCACCAAAAATATGGGCGAAATAACTGCTGCGATAg            <  2:61048/139‑1 (MQ=255)
                                                                                                                         cATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGTAATAACCTGCAGCATATCCGCCACCAAAAATATGGGCGAAATAACTGCTGCGATAGCg          <  2:374225/139‑1 (MQ=255)
                                                                                                                         cATTTTCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGTAATAACCTGCAGCATATCCGCCACCAAAAATATGGGCGAAATAACTGCTGCGATAGCg          <  2:747361/139‑1 (MQ=255)
                                                                                                                           ttttCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGTAATAACCTGCAGCATATCCGCCACCAAAAATATGGGCGAAATAACTGCTGCGATAGCGTg        <  2:1081702/139‑1 (MQ=255)
                                                                                                                              tCACGCGTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGTAATAACCTGCAGCATATCCGCCACCAAAAATATGGGCGAAATAACTGCTGCGATAGCGTGGCg     <  2:56914/139‑1 (MQ=255)
                                                                                                                                 cgcgTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGTAATAACCTGCAGCATATCCGCCACCAAAAATATGGGCGAAATAACTGCTGCGATAGCGTGGCGGta  <  2:1031595/139‑1 (MQ=255)
                                                                                                                                 cgcgTCAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGTAATAACCTGCAGCATATCCGCCACCAAAAATATGGGCGAAATAACTGCTGCGATAGCGTGGCGGta  <  2:995993/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                     
TTCAAGCCACGATGTTGCAGCATCGGCATAATCTTAGGTGCCTTACCGCGCCATTGTCGATACAGGCGTTCCAGATCTTCGCTGTTACCTCTGGAAAGGATCGCCTCGCGAAAACGCAGCCCATTTTCACGCGTTAATCCGCCCTGCTCAACAAACCACTGATAACCATCATCGGCCAACATTTGCGTCCACAGATAAGCGTAATAACCTGCAGCATATCCGCCACCAAAAATATGGGCGAAATAACTGCTGCGATAGCGTGGCGGTA  >  NC_000913/1625342‑1625609

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑