Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 426,471 G→A 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,4710GA94.9% 243.9 / 2.6 78R112H (CGT→CAT) tgttRNA‑guanine transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (43/31);  total (45/33)
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.88e-01

ATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTTGATGAG  >  NC_000913/426374‑426568
                                                                                                 |                                                                                                 
aTGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATa                                                                                                 <  2:574906/100‑1 (MQ=255)
 tGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTaa                                                                                                >  1:750000/1‑100 (MQ=255)
 tGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATaa                                                                                                >  2:763504/1‑100 (MQ=255)
 tGCAGTGGAAGGGGCCGATCATCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATaa                                                                                                >  1:1545785/1‑100 (MQ=255)
    aGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAAccc                                                                                             <  2:1772739/100‑1 (MQ=255)
         aaGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCa                                                                                        >  1:1290536/1‑100 (MQ=255)
            gggCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATTAACg                                                                                     >  2:551626/1‑100 (MQ=255)
            gggCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACg                                                                                     <  1:1036815/100‑1 (MQ=255)
             ggCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACgg                                                                                    >  2:1010636/1‑100 (MQ=255)
             ggCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACgg                                                                                    >  2:976281/1‑100 (MQ=255)
             ggCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACgg                                                                                    >  1:195188/1‑100 (MQ=255)
             ggCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACgg                                                                                    <  1:507920/100‑1 (MQ=255)
             ggCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACg                                                                                     >  2:483597/1‑99 (MQ=255)
               ccGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCg                                                                                  >  1:581197/1‑100 (MQ=255)
                    tcttccgatcTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCg                                                                             <  1:383/91‑1 (MQ=255)
                      tCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGAt                                                                           >  2:382405/1‑100 (MQ=255)
                       cACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGCTATTCGTAAAATCACCGAACAGGACGTGCACTTCCATAAACCGATCAACGACAATCCGAtt                                                                          >  2:601383/1‑100 (MQ=255)
                       cACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGAtt                                                                          >  1:843583/1‑100 (MQ=255)
                         ccGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGAtttt                                                                        >  1:678652/1‑100 (MQ=255)
                           gACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTcc                                                                      >  2:933242/1‑100 (MQ=255)
                             cTCCGGCGGCTTCCAGGTCTTCAGCCTTGTCGATATTCTTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTc                                                                    <  2:1408782/100‑1 (MQ=255)
                             cTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTc                                                                    <  2:195188/100‑1 (MQ=255)
                             cTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTc                                                                    <  2:216595/100‑1 (MQ=255)
                               ccGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGa                                                                  <  2:1159961/100‑1 (MQ=255)
                               ccGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGagatcggaaga                                                                  >  2:383/1‑90 (MQ=255)
                               ccGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGa                                                                  >  1:375590/1‑100 (MQ=255)
                                cggcggCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGAt                                                                 >  2:262796/1‑100 (MQ=255)
                                 ggcggcTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATc                                                                >  1:918856/1‑100 (MQ=255)
                                     gcTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGa                                                            >  1:1744455/1‑100 (MQ=255)
                                     gcTTCCAGGTCTTCAGCCTTGGCGACCTTCGTAAAATCACCGAACAGGCCGTGCACCTCCATCACCAGCTCAACGGCGATACGATTTTCCTCGATCCtca                                                            <  1:1058595/100‑3 (MQ=255)
                                      cTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGaa                                                           >  1:1272342/1‑100 (MQ=255)
                                       ttCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGaaa                                                          <  2:154543/100‑1 (MQ=255)
                                        tCCAGGGCTTCAGCCTTGGCGATATTCTTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGaaaa                                                         <  2:1595395/100‑1 (MQ=255)
                                         ccAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGaaaaa                                                        >  2:556904/1‑100 (MQ=255)
                                          cAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAAt                                                       <  2:1152695/100‑1 (MQ=255)
                                           aGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATc                                                      <  1:153677/100‑1 (MQ=255)
                                              tCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAAt                                                   <  2:750004/100‑1 (MQ=255)
                                              tCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAAt                                                   <  1:201313/100‑1 (MQ=255)
                                               cTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATg                                                  >  1:427790/1‑100 (MQ=255)
                                                ttCAGTCTTGGCGATTTTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATg                                                  <  2:1022721/99‑1 (MQ=255)
                                                ttCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATgg                                                 <  2:354180/100‑1 (MQ=255)
                                                ttCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATgg                                                 <  2:496085/100‑1 (MQ=255)
                                                 tCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGga                                                <  1:255775/100‑1 (MQ=255)
                                                    gCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGAt                                             >  1:25464/1‑100 (MQ=255)
                                                    gCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGAt                                             >  1:1362974/1‑100 (MQ=255)
                                                    gCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGAt                                             >  1:879332/1‑100 (MQ=255)
                                                        tGGCGATATTCGTAAAATCACCGAACAGGACGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAg                                         <  1:263439/100‑1 (MQ=255)
                                                          gCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTa                                       >  2:637609/1‑100 (MQ=255)
                                                          gCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGACCCTGAAAAATCAATGGAGATTCAGTa                                       <  2:1102805/100‑1 (MQ=255)
                                                            gATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACg                                     >  2:875938/1‑100 (MQ=255)
                                                               atTCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATc                                  <  1:987328/100‑1 (MQ=255)
                                                                ttCGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAAATTCAGTACGATCt                                 >  2:184114/1‑100 (MQ=255)
                                                                 tCGTAAAATCACCGAACAGGGCGTGCACTTCCATACCCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATAGAGATTCAGTACGATCtt                                >  2:1765867/1‑100 (MQ=255)
                                                                  cGTAAAATCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCtt                                >  2:1475939/1‑99 (MQ=255)
                                                                     aaaaTCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGtt                            >  1:1567577/1‑100 (MQ=255)
                                                                        aTCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCgg                         <  1:1492161/100‑1 (MQ=255)
                                                                         tCACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGa                        >  1:1389880/1‑100 (MQ=255)
                                                                          cACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGat                       >  1:690062/1‑100 (MQ=255)
                                                                          cACCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCACTACGATCTTGGTTCGGat                       >  1:689918/1‑100 (MQ=255)
                                                                           aCCGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGata                      >  2:177420/1‑100 (MQ=255)
                                                                            ccGAACAGGGCGTGCACTTCCATAACCCTATTAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGatat                     <  2:680215/100‑1 (MQ=255)
                                                                            ccGAACAGGGCGTGCACTTCCATAACCCGATCAACGCCGATCCGATTTTCCTCGATCCTGAAAAAAAAATATAGATTCAGTAAAATATTGGTTCGGatat                     >  2:1128249/1‑100 (MQ=255)
                                                                             cGAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATc                    >  2:315005/1‑100 (MQ=255)
                                                                              gAACAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCg                   >  2:979120/1‑100 (MQ=255)
                                                                                aCAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTc                 <  1:1817988/100‑1 (MQ=255)
                                                                                 cAGGGCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCa                >  2:825703/1‑100 (MQ=255)
                                                                                   gggCGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATg              >  2:1033862/1‑100 (MQ=255)
                                                                                      cGTGCACTTCCATAACCCGATCAACGGCGATCCGCTTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATCTCGCCATGATc           <  1:953197/100‑1 (MQ=255)
                                                                                      cGTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATc           >  1:891278/1‑100 (MQ=255)
                                                                                       gTGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCt          >  1:6808/1‑100 (MQ=255)
                                                                                        tGCACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCtt         <  1:551622/100‑1 (MQ=255)
                                                                                          cACTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtg       <  1:720139/100‑1 (MQ=255)
                                                                                           aCTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTAGTTCGGATATCGTCATGATCTTtga      >  2:1205003/1‑100 (MQ=255)
                                                                                            cTTCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtgat     <  2:825138/100‑1 (MQ=255)
                                                                                             ttCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtgatg    <  2:131176/100‑1 (MQ=255)
                                                                                             ttCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtgatg    <  1:715790/100‑1 (MQ=255)
                                                                                              tCCATAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTtgatga   <  1:347725/100‑1 (MQ=255)
                                                                                                cGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTTGATGAg  >  1:1052856/1‑99 (MQ=255)
                                                                                                 |                                                                                                 
ATGCAGTGGAAGGGGCCGATCCTCACCGACTCCGGCGGCTTCCAGGTCTTCAGCCTTGGCGATATTCGTAAAATCACCGAACAGGGCGTGCACTTCCGTAACCCGATCAACGGCGATCCGATTTTCCTCGATCCTGAAAAATCAATGGAGATTCAGTACGATCTTGGTTCGGATATCGTCATGATCTTTGATGAG  >  NC_000913/426374‑426568

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: ‑