Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 178,902 C→T 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,9020CT88.9% 171.4 / 12.1 63L84L (TTG→TTAmtn5'‑methylthioadenosine/S‑adenosylhomocysteine nucleosidase
Reads supporting (aligned to +/- strand):  ref base C (4/3);  new base T (28/28);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

CGGACAGCCTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCAGCGACT  >  NC_000913/178806‑178998
                                                                                                |                                                                                                
cGGACAGCCTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAAc                                                                                               <  1:1072660/100‑1 (MQ=255)
 gaaCAGCCTAGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACg                                                                                              <  1:649843/98‑1 (MQ=255)
  gACAGCCTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGt                                                                                             >  1:1564291/1‑100 (MQ=255)
   aCAGCCTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGtt                                                                                            <  1:1629344/100‑1 (MQ=255)
    cAGCCTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTg                                                                                           <  2:1443534/100‑1 (MQ=255)
      gCCTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGt                                                                                         >  2:34894/1‑100 (MQ=255)
         tGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGcc                                                                                      <  1:1928867/100‑1 (MQ=255)
          ggTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCa                                                                                     >  2:1473958/1‑100 (MQ=255)
           gTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAg                                                                                    >  2:51693/1‑100 (MQ=255)
            tAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAgg                                                                                   <  1:1534831/100‑1 (MQ=255)
               cTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCa                                                                                >  1:852612/1‑100 (MQ=255)
               cTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCa                                                                                >  1:1066045/1‑100 (MQ=255)
               cTGACCGTATTCATAACCAAATGCCGGGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGcc                                                                                 >  2:1187754/1‑99 (MQ=255)
                tGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAAAGTTGGTGCCAGGCCa                                                                                >  1:1557610/1‑99 (MQ=255)
                 gACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCAcc                                                                              <  1:1585345/100‑1 (MQ=255)
                     gTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCa                                                                          <  2:225339/100‑1 (MQ=255)
                          cATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACCACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAAcc                                                                     >  1:809354/1‑100 (MQ=255)
                          cATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAAcc                                                                     >  1:1779039/1‑100 (MQ=255)
                          cATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAAcc                                                                     >  2:1239944/1‑100 (MQ=255)
                          cATAACCAAATGCAGTGACATCCGCGTCGTGATAACGTACTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAAcc                                                                     >  1:1200574/1‑100 (MQ=255)
                            tAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCgg                                                                   >  1:1874024/1‑100 (MQ=255)
                            tAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCgg                                                                   <  2:1213585/100‑1 (MQ=255)
                             aaCCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGgt                                                                  <  2:1017007/100‑1 (MQ=255)
                                cAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGtt                                                               >  2:1256235/1‑100 (MQ=255)
                                cAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGtt                                                               <  1:1270645/100‑1 (MQ=255)
                                cAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGtt                                                               <  2:852613/100‑1 (MQ=255)
                                cAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGtt                                                               <  2:622758/100‑1 (MQ=255)
                                cAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGtt                                                               >  2:866218/1‑100 (MQ=255)
                                cAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGtt                                                               >  1:264260/1‑100 (MQ=255)
                                 aaaTGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTta                                                              <  2:81399/100‑1 (MQ=255)
                                       cGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATc                                                        <  2:1405611/100‑1 (MQ=255)
                                       cGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATc                                                        >  2:1683195/1‑100 (MQ=255)
                                        gTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATca                                                       <  2:1378607/100‑1 (MQ=255)
                                          gACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATcaca                                                     <  2:716272/100‑1 (MQ=255)
                                                  cgtcgtGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCtt                                             >  2:631007/1‑100 (MQ=255)
                                                   gtcgtGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCtt                                             >  2:728165/1‑99 (MQ=255)
                                                    tcgtGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGc                                           <  1:866217/100‑1 (MQ=255)
                                                      gtGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAg                                         >  2:531347/1‑100 (MQ=255)
                                                      gtGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAg                                         >  2:657879/1‑100 (MQ=255)
                                                         aTAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAgtgt                                      >  2:727569/1‑100 (MQ=255)
                                                           aaCGTGCTGCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACAGCTGGCTTGCAGTGtt                                     <  1:518489/99‑1 (MQ=255)
                                                            aCGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTcc                                   <  2:1874036/100‑1 (MQ=255)
                                                            aCGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTcc                                   <  1:946682/100‑1 (MQ=255)
                                                             cGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCa                                  >  1:1197594/1‑100 (MQ=255)
                                                             cGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCa                                  >  2:516938/1‑100 (MQ=255)
                                                             cGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCa                                  >  2:186000/1‑100 (MQ=255)
                                                             cGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCa                                  >  1:1836818/1‑100 (MQ=255)
                                                                gCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACAGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACa                               >  1:737450/1‑100 (MQ=255)
                                                                 cTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAg                              <  1:51693/100‑1 (MQ=255)
                                                                  ttCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGc                             <  2:1631908/100‑1 (MQ=255)
                                                                      tCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAg                         >  2:301882/1‑100 (MQ=255)
                                                                         gagaCAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCGGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTgg                      <  1:1464794/100‑1 (MQ=255)
                                                                         gagaCAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTgg                      <  2:1557621/100‑1 (MQ=255)
                                                                           gaCAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCa                    >  1:918663/1‑100 (MQ=255)
                                                                           gaCAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCa                    <  1:1233667/100‑1 (MQ=255)
                                                                              aaCGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCAccc                 <  1:285925/100‑1 (MQ=255)
                                                                                       gCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCa        >  2:1738012/1‑100 (MQ=255)
                                                                                       gCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCa        >  2:138682/1‑100 (MQ=255)
                                                                                       gCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCa        >  2:877551/1‑100 (MQ=255)
                                                                                        cccACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCAg       <  1:1473950/100‑1 (MQ=255)
                                                                                          cACTTTTACCGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTACTAATCACATCTGGCTTGCACTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCAGCg     <  1:1928689/100‑1 (MQ=255)
                                                                                          cACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCAGCg     <  1:1256231/100‑1 (MQ=255)
                                                                                          cACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAACTGGCACCCAGCGCCGCAGCg     <  1:1187752/100‑1 (MQ=255)
                                                                                             ttttAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCAGCGACt  >  2:1426334/1‑100 (MQ=255)
                                                                                             ttttAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCAGCGACt  >  1:79472/1‑100 (MQ=255)
                                                                                                |                                                                                                
CGGACAGCCTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCAGCGACT  >  NC_000913/178806‑178998

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