Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 178,902 C→T 100% L84L (TTG→TTA mtn ← 5'‑methylthioadenosine/S‑adenosylhomocysteine nucleosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913178,9020CT98.7% 254.9 / ‑4.6 79L84L (TTG→TTAmtn5'‑methylthioadenosine/S‑adenosylhomocysteine nucleosidase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (38/39);  total (39/40)
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.67e-01

CCGGACAGCCTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCAGCGAC  >  NC_000913/178805‑178997
                                                                                                 |                                                                                               
ccGGACAGCCTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTaa                                                                                               >  2:1615874/1‑100 (MQ=255)
 cGGACAGCCTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAAc                                                                                              <  1:1028394/100‑1 (MQ=255)
  ggACAGCCTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACg                                                                                             >  2:183902/1‑100 (MQ=255)
     cAGCCTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCTTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGt                                                                                            >  2:696486/1‑98 (MQ=255)
        ccTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTg                                                                                       <  1:382675/100‑1 (MQ=255)
        ccTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTg                                                                                       >  2:2049999/1‑100 (MQ=255)
        ccTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTg                                                                                       <  2:1002495/100‑1 (MQ=255)
         cTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGc                                                                                      >  2:2242934/1‑100 (MQ=255)
           ggTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCa                                                                                    >  2:2089112/1‑100 (MQ=255)
            gTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAg                                                                                   <  1:1151098/100‑1 (MQ=255)
             tAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGTCAgg                                                                                  <  1:183904/100‑1 (MQ=255)
             tAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAgg                                                                                  <  1:1677538/100‑1 (MQ=255)
               aCTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTCTAACGTTGGTGCCAGGcc                                                                                >  2:2053576/1‑100 (MQ=255)
                    ccGTATTCATAACCAAATGCCGTGACATCCGCGTCGTAATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGcccccggcagga                                                                           >  2:1369665/1‑90 (MQ=255)
                     cGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGc                                                                          <  1:1346259/100‑1 (MQ=255)
                      gTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCa                                                                         <  1:352455/100‑1 (MQ=255)
                       tATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAg                                                                        <  1:2227731/100‑1 (MQ=255)
                         ttCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGaa                                                                      <  2:151428/100‑1 (MQ=255)
                          tCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAAc                                                                     <  1:448013/100‑1 (MQ=255)
                            aTAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAAcc                                                                    >  2:1827527/1‑99 (MQ=255)
                            aTAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCg                                                                   <  1:2180581/100‑1 (MQ=255)
                              aaCCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGgt                                                                 <  1:534591/100‑1 (MQ=255)
                               aCCAAATGCCGTGACATCCGCGTCGTGATCACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTACCGTTGGTGCCAGGCCACCGGCAGAACCGgtg                                                                >  2:1940983/1‑100 (MQ=255)
                                 cAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGtt                                                              <  1:2409797/100‑1 (MQ=255)
                                  aaaTGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTta                                                             >  1:1357015/1‑100 (MQ=255)
                                  aaaTGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTta                                                             >  1:1356386/1‑100 (MQ=255)
                                      gCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTtaata                                                          <  2:1169837/99‑1 (MQ=255)
                                       ccGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTaataat                                                        <  1:2303778/100‑1 (MQ=255)
                                       ccGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTaataat                                                        >  2:901801/1‑100 (MQ=255)
                                        cGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATc                                                       <  1:961185/100‑1 (MQ=255)
                                        cGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATc                                                       >  1:432878/1‑100 (MQ=255)
                                         gttaCATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATca                                                      <  2:2104982/97‑1 (MQ=255)
                                         gTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATca                                                      <  1:738178/100‑1 (MQ=255)
                                         gTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATca                                                      <  1:1132879/100‑1 (MQ=255)
                                         gTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATca                                                      >  1:1046629/1‑100 (MQ=255)
                                          tGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTCCGTGCCAGGCCACCTGCAGAACCGGTGTTACTAATcac                                                     >  2:972516/1‑100 (MQ=255)
                                           gACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATcaca                                                    <  1:219901/100‑1 (MQ=255)
                                             cATCCGCGTCGTGATAACGTGCTTCGGCCGAGACAACGATCTCTCCCACTTTTCACGTTGGTGCCAGGCCACCGGCTGACCCGGTGTTAATAATcac                                                     >  2:888351/1‑97 (MQ=255)
                                                ccGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTgg                                               >  1:2039697/1‑100 (MQ=255)
                                                ccGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTgg                                               <  1:921383/100‑1 (MQ=255)
                                                ccGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTgg                                               >  2:827292/1‑100 (MQ=255)
                                                     tcgtGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGc                                          <  2:2039683/100‑1 (MQ=255)
                                                       gtGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAg                                        >  1:1113941/1‑100 (MQ=255)
                                                       gtGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAg                                        <  2:742744/100‑1 (MQ=255)
                                                          aTAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAgtgt                                     >  1:1652849/1‑100 (MQ=255)
                                                            aaCGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTc                                   <  2:124974/100‑1 (MQ=255)
                                                             aCGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTcc                                  >  1:975240/1‑100 (MQ=255)
                                                              cGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCa                                 <  1:719886/100‑1 (MQ=255)
                                                               gTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCaa                                <  1:1197589/100‑1 (MQ=255)
                                                               gTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCa                                 <  2:1545155/99‑1 (MQ=255)
                                                                tGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTACCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAAc                               <  1:836863/100‑1 (MQ=255)
                                                                 tcTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGACCAACGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGCTCCAACa                              >  2:1716560/2‑100 (MQ=255)
                                                                 ccgatctCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACa                              <  2:539760/94‑1 (MQ=255)
                                                                  cTTTGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATGACATCTGGCTTGCAGTGTTCCAACAg                             <  1:1558046/100‑1 (MQ=255)
                                                                  cTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAg                             >  2:1974164/1‑100 (MQ=255)
                                                                   ttCGTCCGAGACAACGATATCGCCCAGTTTGAACGTTGGTGCCAGGCCAGCGGCAGAACCGGTGTTGATAATCGCATCTGGGTTGCAGAGTTCAAACAg                             >  2:1291642/1‑99 (MQ=255)
                                                                   ttCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGc                            >  1:1566835/1‑100 (MQ=255)
                                                                        ccGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACaagatcgg                       >  1:539767/1‑93 (MQ=255)
                                                                          gagaCAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTgg                     <  2:563866/100‑1 (MQ=255)
                                                                          gagaCAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTgg                     <  1:1789956/100‑1 (MQ=255)
                                                                          gagaCAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTgg                     <  2:1199080/100‑1 (MQ=255)
                                                                            gaCAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCTTCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCa                   >  2:1769596/1‑100 (MQ=255)
                                                                              cAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGTCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTg                      >  2:1069694/1‑95 (MQ=255)
                                                                              cAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCAcc                 >  1:854770/1‑100 (MQ=255)
                                                                                aCGATAGCGCCCACTTTTAACGTGGGTGCCAGGCCACCGGCAGAACCGCTGTTAATAATTGCCTCTGGCTTTCAGTGTTCTCACAGCaagctcg                     >  2:2161096/1‑89 (MQ=255)
                                                                                  gATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAgc             <  1:2164198/100‑1 (MQ=255)
                                                                                     atCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGcgcc          >  1:1297321/1‑100 (MQ=255)
                                                                                       cGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAACCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGcgccgc        >  2:1130267/1‑100 (MQ=255)
                                                                                        tcccACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCa       >  2:893082/2‑100 (MQ=255)
                                                                                        gCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGcgccgta       <  1:1735717/100‑3 (MQ=255)
                                                                                        gCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCa       >  2:892889/1‑100 (MQ=255)
                                                                                        gCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCa       >  1:2504759/1‑100 (MQ=255)
                                                                                         cccACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGTTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCAg      >  2:2341470/1‑100 (MQ=255)
                                                                                         cccACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCAg      <  1:1841116/100‑1 (MQ=255)
                                                                                         cccACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCAg      <  1:334239/100‑1 (MQ=255)
                                                                                         cccACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCAg      >  1:2101654/1‑100 (MQ=255)
                                                                                          ccACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCAGc     <  2:752405/100‑1 (MQ=255)
                                                                                          ccACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCAGc     <  1:1977343/100‑1 (MQ=255)
                                                                                             cTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCAGCGAc  >  2:1970772/1‑100 (MQ=255)
                                                                                              ttttAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCAGCGAc  <  1:1096889/99‑1 (MQ=255)
                                                                                                 |                                                                                               
CCGGACAGCCTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCAGCGAC  >  NC_000913/178805‑178997

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: