Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 426,471 G→A 100% R112H (CGT→CAT)  tgt → tRNA‑guanine transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913426,4710GA96.4% 261.0 / ‑0.0 83R112H (CGT→CAT) tgttRNA‑guanine transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (39/41);  total (40/43)
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.44e-01

ATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTTGATG  >  NC_000913/426374‑426566
                                                                                                 |                                                                                               
aTGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATa                                                                                               <  2:2048299/100‑1 (MQ=255)
 tGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATaa                                                                                              >  1:813949/1‑100 (MQ=255)
 tGCAGTGCAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACCGGGCGTGCACTTCCATa                                                                                               >  1:534198/1‑99 (MQ=255)
    aGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAAccc                                                                                           <  1:603590/100‑1 (MQ=255)
            gggCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACg                                                                                   >  2:460532/1‑100 (MQ=255)
             ggCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGGGCACTTCCATAACCCGATCAACgg                                                                                  >  2:1844314/1‑100 (MQ=255)
              gCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGc                                                                                 <  2:588851/100‑1 (MQ=255)
              gCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGc                                                                                 <  2:1461623/100‑1 (MQ=255)
              gCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGc                                                                                 <  1:2177270/100‑1 (MQ=255)
              gCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGc                                                                                 <  1:1832434/100‑1 (MQ=255)
               ccGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCg                                                                                <  1:40630/100‑1 (MQ=255)
               ccGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCg                                                                                >  2:1487717/1‑100 (MQ=255)
               ccGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCg                                                                                >  2:2061418/1‑100 (MQ=255)
                     cTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGa                                                                          >  2:118040/1‑100 (MQ=255)
                         ccGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCATCGGCGATCCGAtttt                                                                      >  2:2559622/1‑100 (MQ=255)
                         ccGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGAtttt                                                                      <  2:1390890/100‑1 (MQ=255)
                           gACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTcc                                                                    <  1:907048/100‑1 (MQ=255)
                            aCTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCt                                                                   <  1:172566/100‑1 (MQ=255)
                             cTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTc                                                                  <  2:582356/100‑1 (MQ=255)
                                cggcggCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGAt                                                               <  2:179113/100‑1 (MQ=255)
                                 ggcggcTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATc                                                              >  2:1116103/1‑100 (MQ=255)
                                 ggcggcTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATcc                                                             >  1:2057723/1‑100 (MQ=255)
                                  gcggcTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATcc                                                             <  1:226213/100‑1 (MQ=255)
                                    ggcTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTg                                                           >  2:1359378/1‑100 (MQ=255)
                                    ggcTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTg                                                           >  1:152220/1‑100 (MQ=255)
                                     gcTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGa                                                          <  2:2409643/100‑1 (MQ=255)
                                      cTTCCAGGTTTTCAGCCTTTGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATcc                                                             <  1:2155458/96‑1 (MQ=255)
                                       ttCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGaaa                                                        <  1:39450/100‑1 (MQ=255)
                                        tCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGaaa                                                        >  2:195865/1‑99 (MQ=255)
                                         ccAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGaaaaa                                                      >  2:1727008/1‑100 (MQ=255)
                                          cAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAAt                                                     >  1:1054831/1‑100 (MQ=255)
                                          cAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAAt                                                     >  2:220600/1‑100 (MQ=255)
                                           aGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATc                                                    >  1:1502432/1‑100 (MQ=255)
                                              tCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAAt                                                 <  1:2404893/100‑1 (MQ=255)
                                              tCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAAt                                                 <  2:313892/100‑1 (MQ=255)
                                               cTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATg                                                <  2:933636/100‑1 (MQ=255)
                                               cTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATg                                                <  2:156416/100‑1 (MQ=255)
                                                 tCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGga                                              <  2:391807/100‑1 (MQ=255)
                                                 tCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGga                                              >  2:1899025/1‑100 (MQ=255)
                                                  cAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGgag                                             <  2:617593/100‑1 (MQ=255)
                                                   aGCCTTGGCGATATTCGTAAAATCTCCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGgaga                                            <  1:847184/100‑1 (MQ=255)
                                                    gCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGAt                                           >  1:1837598/1‑100 (MQ=255)
                                                      cTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTc                                         <  2:363436/100‑1 (MQ=255)
                                                       ttGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCa                                        <  1:490902/100‑1 (MQ=255)
                                                        tGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAg                                       >  2:498775/1‑100 (MQ=255)
                                                          gCGATATTCGTAAAATCACCGAACAGGCCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTa                                     <  1:874854/100‑1 (MQ=255)
                                                           cGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTAc                                    >  2:466961/1‑100 (MQ=255)
                                                               atTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATc                                <  1:1530922/100‑1 (MQ=255)
                                                                 tagTAAAATCACCCAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCtt                              <  2:1635192/98‑1 (MQ=255)
                                                                 tCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCtt                              >  1:77366/1‑100 (MQ=255)
                                                                  cGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTg                             <  2:1502441/100‑1 (MQ=255)
                                                                  cGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTg                             >  2:29897/1‑100 (MQ=255)
                                                                   gTAAAATCGCCGAGCGGGGCGTGTACTTCCATAGCCCGAGCAACGTCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACTATCTTgg                            <  1:483164/100‑1 (MQ=255)
                                                                   gTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTgg                            >  1:898277/1‑100 (MQ=255)
                                                                   gTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTgg                            >  2:54910/1‑100 (MQ=255)
                                                                   gTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATGCAGTACGATCTTgg                            <  2:125151/100‑1 (MQ=255)
                                                                      aaaTCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGACCTTGGtt                          >  1:359389/1‑99 (MQ=255)
                                                                       aaTCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCg                        >  1:136214/1‑100 (MQ=255)
                                                                        agcACCCAACCTGGCGTGCACTTCCATAACCCGATCAACGTCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTTGTTCgg                       <  2:2188718/98‑1 (MQ=255)
                                                                        aTCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCgg                       <  1:1186220/100‑1 (MQ=255)
                                                                        aTCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCgg                       <  2:952400/100‑1 (MQ=255)
                                                                          cACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGa                      <  2:479862/99‑1 (MQ=255)
                                                                           aCCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCATTGGAGATTCAGTACGATCTTGGTTCGGata                    >  1:1207443/1‑100 (MQ=255)
                                                                           aCCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTCCGATCTTGGTTCGGata                    >  2:214804/1‑100 (MQ=255)
                                                                           aCCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCAGata                    >  2:1866518/1‑100 (MQ=255)
                                                                             cGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATc                  >  2:2329857/1‑100 (MQ=255)
                                                                              gAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCg                 >  2:926209/1‑100 (MQ=255)
                                                                              gAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCg                 >  1:165064/1‑100 (MQ=255)
                                                                              gAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCg                 >  1:425851/1‑100 (MQ=255)
                                                                               aaCAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGt                >  2:2329763/1‑100 (MQ=255)
                                                                               aaCAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGt                <  1:286200/100‑1 (MQ=255)
                                                                               aaCAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGt                >  2:785949/1‑100 (MQ=255)
                                                                               aaCAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGt                <  2:1136706/100‑1 (MQ=255)
                                                                                aCAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTc               <  2:156526/100‑1 (MQ=255)
                                                                                 cAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTACTCGATACTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTcn              <  1:743760/100‑2 (MQ=255)
                                                                                    ggCGTGCACTTCAATAACCCGATCAACGGAGATCCGATTTTCCTCGATACTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGa           >  1:1170966/1‑100 (MQ=255)
                                                                                     gCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGAt          <  2:814333/100‑1 (MQ=255)
                                                                                      cGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATc         >  1:1736566/1‑100 (MQ=255)
                                                                                       gTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCt        <  2:1881468/100‑1 (MQ=255)
                                                                                         gCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGACCCTGAAAAATCAATGGAGATTCAGTCCGATCTTGGTTCGGATATCGTCATGATCttt      <  2:874597/100‑1 (MQ=255)
                                                                                          cACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtg     >  2:1201188/1‑100 (MQ=255)
                                                                                            cTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtgat   >  1:916236/1‑100 (MQ=255)
                                                                                             ttCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtgatg  <  1:1039774/100‑1 (MQ=255)
                                                                                                 |                                                                                               
ATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTTGATG  >  NC_000913/426374‑426566

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