Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,766,473 A→G D248G (GAT→GGT)  b2632 → putative GTP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,766,4730AG99.0% 244.6 / NA 95D248G (GAT→GGT) b2632putative GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (40/54);  minor base C (1/0);  total (41/54)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TACTGCAACATCCCTTTCGGACATGGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGATACCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGT  >  NC_000913/2766404‑2766542
                                                                     |                                                                     
tACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCAt                                                                 >  2:4008873/1‑76 (MQ=255)
tACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCAt                                                                 >  2:1268590/1‑76 (MQ=255)
 aCTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATa                                                                <  1:3037021/76‑1 (MQ=255)
  cTGCAACATCCCTTTCGAGCAGCGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATAc                                                               <  2:709420/76‑1 (MQ=255)
  cTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATAc                                                               <  1:115097/76‑1 (MQ=255)
  cTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATAc                                                               >  2:2472989/1‑76 (MQ=255)
   tGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGAACATTTCGGTGAAACCGTAGGTGCCATAcc                                                              >  2:2660946/1‑75 (MQ=255)
    gCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTc                                                             >  2:3612892/1‑76 (MQ=255)
     cAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCg                                                            <  1:2361065/76‑1 (MQ=255)
      caCATCCCTTTCGAACAACTGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGa                                                           <  2:4027875/75‑1 (MQ=255)
      aaCATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGa                                                           <  2:1085241/76‑1 (MQ=255)
       aCATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGAt                                                          <  2:3339371/76‑1 (MQ=255)
       aCATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGAt                                                          <  2:2527547/76‑1 (MQ=255)
       aCATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGAt                                                          <  1:3957390/76‑1 (MQ=255)
        cATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATa                                                         <  2:3727075/76‑1 (MQ=255)
         aTCCCTTTCGACCAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATAc                                                        <  2:2691947/76‑1 (MQ=255)
         aTCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATAc                                                        <  1:3202396/76‑1 (MQ=255)
          tCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACg                                                       <  2:1960571/76‑1 (MQ=255)
            ccTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACgga                                                     >  2:1763192/1‑75 (MQ=255)
             cTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTa                                                    <  1:2493435/76‑1 (MQ=255)
                         ggTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTccc                                        <  1:3846797/76‑1 (MQ=255)
                         ggTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTccc                                        <  2:3597408/76‑1 (MQ=255)
                         ggTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTccc                                        >  1:2374563/1‑76 (MQ=255)
                         ggTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTccc                                        >  2:2229011/1‑76 (MQ=255)
                         ggTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTccc                                        >  2:951300/1‑76 (MQ=255)
                          ttAGCCCGGGAGCAGGCACGTGACGATTTCGGTGACACCGTAGGTTCGATACTCGATACTGTAAGCACCTTTcccc                                       >  1:1696229/2‑76 (MQ=255)
                          gTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTcccc                                       >  1:1666668/1‑76 (MQ=255)
                          gTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTcccc                                       <  2:4520238/76‑1 (MQ=255)
                           tAGCCCGGGCGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATGCTCGATCCGGTAAGCACCTTTCCCCt                                      <  2:1236306/76‑1 (MQ=255)
                           tAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTCAGCACCTTTCCCCt                                      >  2:3334695/1‑76 (MQ=255)
                            aGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAGCCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCtt                                     <  1:410439/76‑1 (MQ=255)
                            aGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTATGCACCTTTCCCCtt                                     <  2:1716474/76‑1 (MQ=255)
                            aGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCtt                                     <  1:3121443/76‑1 (MQ=255)
                            aGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCtt                                     <  2:160211/76‑1 (MQ=255)
                            aGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACAGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCtt                                     >  1:3369602/1‑76 (MQ=255)
                             gCCCGGTAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCCCCTTTCCCCTTa                                    <  1:3034632/76‑1 (MQ=255)
                             gCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTa                                    >  1:3800015/1‑76 (MQ=255)
                               ccGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAg                                                  >  2:4519304/1‑60 (MQ=255)
                               ccGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAg                                                  <  1:4519304/60‑1 (MQ=255)
                               ccGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTAtt                                  >  2:2225406/1‑76 (MQ=255)
                               ccGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTAtt                                  <  2:2721754/76‑1 (MQ=255)
                               ccGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGCGCCATACTCGATACGGTAAGCACCTTTCCCCTTAtc                                  >  2:3959884/1‑75 (MQ=255)
                                 gggAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTcc                                >  1:2733852/1‑76 (MQ=255)
                                 gggAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTcc                                >  1:3996105/1‑76 (MQ=255)
                                  ggAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTccc                               <  2:1943917/76‑1 (MQ=255)
                                  ggAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTccc                               <  1:4039616/76‑1 (MQ=255)
                                    aGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGc                             >  2:3065232/1‑76 (MQ=255)
                                    aGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGc                             >  1:3704531/1‑76 (MQ=255)
                                    aGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGc                             >  1:457018/1‑76 (MQ=255)
                                     gCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGcc                            <  1:2781154/76‑1 (MQ=255)
                                     gCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGcc                            >  1:3205665/1‑76 (MQ=255)
                                     gCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGcc                            <  2:1947266/76‑1 (MQ=255)
                                     gCAGGCACGTGACGATTTCGGTGAAACAGCAGCCCCCATACCCGTGACGGCAAGCACCTTGCCCCTCATCcccccc                            >  2:3390780/1‑76 (MQ=255)
                                      cAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGccc                           >  2:1895930/1‑76 (MQ=255)
                                       aGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGcccc                          >  1:418017/1‑76 (MQ=255)
                                       aGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGcccc                          <  1:882607/76‑1 (MQ=255)
                                       aGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGcccc                          >  2:1898960/1‑76 (MQ=255)
                                       aGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGcccc                          <  1:3449198/76‑1 (MQ=255)
                                        ggCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCg                         <  1:4022283/76‑1 (MQ=255)
                                        ggCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCg                         <  1:2310609/76‑1 (MQ=255)
                                         gCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCgg                        <  1:528057/76‑1 (MQ=255)
                                           ccGTGACGATTCCGGTGCAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTCTTCCCGCCCCGGTg                      <  1:582230/75‑1 (MQ=255)
                                           agttGACGATTTCGGCGCGTCCGTAGGTTCCCTACTCGGTACGGTAACCTCCTTTCCCCTTATTCCCGCCCCCGTg                      <  1:2409557/73‑1 (MQ=255)
                                            cGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGc                     <  1:4481295/76‑1 (MQ=255)
                                            cGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGc                     <  1:3640951/76‑1 (MQ=39)
                                            cGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGc                     <  1:4335613/76‑1 (MQ=39)
                                            cGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGCTGc                     >  1:723092/1‑76 (MQ=255)
                                            cGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGGAAGCACCTTTCCCCTTATTCCCGCCCCGGTGc                     <  2:714909/76‑1 (MQ=38)
                                             gTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCg                    >  1:4441985/1‑76 (MQ=38)
                                             gTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCg                    >  1:2858359/1‑76 (MQ=38)
                                              tGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCgg                   <  2:3636495/76‑1 (MQ=38)
                                              tGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCgg                   >  1:1151299/1‑76 (MQ=38)
                                              ggACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCgg                   <  2:2700958/75‑1 (MQ=38)
                                               gACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGa                  >  2:1516985/1‑76 (MQ=38)
                                               gACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGa                  <  1:3771041/76‑1 (MQ=38)
                                                  gATTTCGGTGAAGTCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGGCCa               <  1:3379290/76‑1 (MQ=37)
                                                  gATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCa               <  1:2398316/76‑1 (MQ=38)
                                                  gATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCa               >  1:4481794/1‑76 (MQ=38)
                                                   aTTTCGGTGAAACCGTAGGTGCGATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCat              <  1:259387/76‑1 (MQ=38)
                                                    tttCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCatt             <  1:2330394/76‑1 (MQ=39)
                                                    tttCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGAACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCatt             >  2:4729488/1‑76 (MQ=38)
                                                     ttCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCatta            >  1:4211106/1‑76 (MQ=39)
                                                      tCGGTGAAACCTTGGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCattat           <  2:4644652/76‑1 (MQ=38)
                                                      tCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCattat           <  1:944953/76‑1 (MQ=255)
                                                      tCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCattat           <  1:973733/76‑1 (MQ=255)
                                                       tggTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCattatt          <  2:3369062/75‑1 (MQ=38)
                                                       cGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCattatt          >  2:3684880/1‑76 (MQ=39)
                                                        ggTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTc         >  2:525688/1‑76 (MQ=38)
                                                          tGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAg       <  2:3889317/76‑1 (MQ=255)
                                                           gAAACCGGAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAgg      <  2:2134572/76‑1 (MQ=39)
                                                            aaaCCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGc     >  2:1429158/1‑76 (MQ=38)
                                                             aaCCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGcc    >  2:4002111/1‑76 (MQ=38)
                                                              aCCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCg   >  1:2975093/1‑76 (MQ=38)
                                                               ccGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGt  >  1:3676304/1‑76 (MQ=38)
                                                               ccGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGt  <  2:837168/76‑1 (MQ=255)
                                                                     |                                                                     
TACTGCAACATCCCTTTCGGACATGGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGATACCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGT  >  NC_000913/2766404‑2766542

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