Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,759,6290GA50.0% 8.6 / 70.3 58A78V (GCG→GTG) selBselenocysteinyl‑tRNA‑specific translation factor
Reads supporting (aligned to +/- strand):  ref base G (14/15);  new base A (14/15);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TCATCAACACGCGCTTCGTCCACGCGATCGGCTTTGGTCAGCGCCACTGTCAGCATCGGGTTACCGGTCAGCTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCGCGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGGCATAGCCGAGATCGATGGTCATGCCGCGC  >  NC_000913/3759491‑3759763
                                                                                                                                          |                                                                                                                                      
tcatcaACACGCGCTTCGTCCACGCGATCGGCTTTGGTCAGCGCCACTGTCAGCATCGGGTTACCGGTCAGCTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAgcg                                                                                                                                        >  2:70084/1‑139 (MQ=255)
 catcaACACGCGCTTCGTCCACGCGATCGGCTTTGGTCAGCGCCACTGTCAGCATCGGGTTACCGGTCAGCTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAgcgc                                                                                                                                       >  2:326490/1‑139 (MQ=255)
   tcaACACGCGCTTCGTCCACGCGATCGGCTTTGGTCAGCGCCACTGTCAGCATCGGGTTACCGGTCAGCTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCACGt                                                                                                                                     >  2:89427/1‑139 (MQ=255)
    caACACGCGCTTCGTCCACGCGATCGGCTTTGGTCAGCGCCACTGTCAGCATCGGGTTACCGGTCAGCTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCGCGTg                                                                                                                                    >  1:710416/1‑139 (MQ=255)
           cgcTTCGTCCACGCGATCGGCTTTGGTCAGCGCCACTGTCAGCATCGGGTTACCGGTCAGCTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCACGTGATCGATa                                                                                                                             <  2:1132951/139‑1 (MQ=255)
            gcTTCGTCCACGCGATCGGCTTTGGTCAGCGCCACTGTCAGCATCGGGTTACCGGTCAGCTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCGCGTGATCGATac                                                                                                                            >  1:327858/1‑139 (MQ=255)
                 gTCCACGCGATCGGCTTTGGTCAGCGCCACTGTCAGCATCGGGTTACCGGTCAGCTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCGCGTGATCGATAccacca                                                                                                                       <  2:327858/139‑1 (MQ=255)
                    cACGCGATCGGCTTTGGTCAGCGCCACTGTCAGCATCGGGTTACCGGTCAGCTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCGCGTGATCGATACCACCAAcg                                                                                                                    <  2:1638281/139‑1 (MQ=255)
                           tCGGCTTTGGTCAGCGCCACTGTCAGCATCGGGTTACCGGTCAGCTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCGCGTGATCGATACCACCAACGCCCGCCa                                                                                                             <  2:991566/139‑1 (MQ=255)
                            cGGCTTTGGTCAGCGCCACTGTCAGCATCGGGTTACCGGTCAGCTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACTAACAGCACGTGATCGATACCACCAACGCCCGCCAg                                                                                                            <  1:1437559/139‑1 (MQ=255)
                            cGGCTTTGGTCAGCGCCACTGTCAGCATCGGGTTACCGGTCAGCTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCGCGTGATCGATACCACCAACGCCCGCCAg                                                                                                            >  2:442480/1‑139 (MQ=255)
                            cGGCTTTGGTCAGCGCCACTGTCAGCATCGGGTTACCGGTCAGCTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCACGTGATCGATACCACCAAcgcc                                                                                                                  >  2:394061/1‑133 (MQ=255)
                            cGGCTTTGGTCAGCGCCACTGTCAGCATCGGGTTACCGGTCAGCTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCACGTGATCGATACCACCAAcgcc                                                                                                                  <  1:394061/133‑1 (MQ=255)
                            cGGCTTTGGTCAGCGCCACTGTCAGCATCGGGTTACCGGTCAGCTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCACGTGATCGATACCACCAACGCCCGCCa                                                                                                             <  1:1411311/138‑1 (MQ=255)
                            cGGCTTTGGTCAGCGCCACTGTCAGCATCGGGTTACCGGTCAGCTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCACGTGATCGATACCACCAACGCCCGCCa                                                                                                             >  2:1411311/1‑138 (MQ=255)
                                 ttGGTCAGCGCCACTGTCAGCATCGGGTTACCGGTCAGCTGCTAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCACGTGATCGATACCACCAACGCCCGCCAGCATGt                                                                                                       >  1:1082590/1‑139 (MQ=255)
                                   ggTCAGCGCCACTGTCAGCATCGGGTTACCGGTCAGCTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGcacg                                                                                                                                      >  1:1534404/1‑106 (MQ=255)
                                     tCAGCGCCACTGTCAGCATCGGGTTACCGGTCAGCTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCACGTGATCGATACCACCAACGCCCGCCAGCATGTTGga                                                                                                   >  1:1651034/1‑139 (MQ=255)
                                         cgcCACTGTCAGCATCGGGTTACCGGTCAGCTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCGCGTGATCGATACCACCAACGCCCGCCAGCATGTTGgaaaga                                                                                               >  1:1347024/1‑139 (MQ=255)
                                         cgcCACTGTCAGCATCGGGTTACCGGTCAGCTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCACGTGATCGATACCACCAACGCCCGCCAGCATGTTGgaaaga                                                                                               <  2:1651034/139‑1 (MQ=255)
                                                          ggTTACCGGTCAGCTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCACGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCggg                                                                              <  2:330025/139‑1 (MQ=255)
                                                          ggTTACCGGTCAGCTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCACGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCggg                                                                              >  2:602109/1‑139 (MQ=255)
                                                             tACCGGTCAGCTGCAAGATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCGCGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAAc                                                                           <  1:1391880/139‑1 (MQ=255)
                                                                  gTCAGCTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCACGTGATc                                                                                                                                 <  1:1186022/80‑1 (MQ=255)
                                                                  gTCAGCTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCACGTGATc                                                                                                                                 >  2:1186022/1‑80 (MQ=255)
                                                                     aGCTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCGCGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATaa                                                                   >  2:622345/1‑139 (MQ=255)
                                                                       cTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCGCGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATaaaa                                                                 >  1:821320/1‑139 (MQ=255)
                                                                       cTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCGCGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATaaaa                                                                 >  1:187992/1‑139 (MQ=255)
                                                                        tGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGACACCACCAACAGCACGTGGTCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAAc                                                                <  2:811072/139‑1 (MQ=255)
                                                                            aaaTCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCGCGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAgg                                                            <  2:1512702/139‑1 (MQ=255)
                                                                                  ccAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCGCGTGATCGATACCACCAACGCCCGCCAGc                                                                                                           >  1:791024/1‑86 (MQ=255)
                                                                                  ccAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCGCGTGATCGATACCACCAACGCCCGCCAGc                                                                                                           <  2:791024/86‑1 (MQ=255)
                                                                                     gATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCGCGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCAcgc                                                       <  2:313382/135‑1 (MQ=255)
                                                                                     gATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCGCGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCAcgc                                                       >  1:313382/1‑135 (MQ=255)
                                                                                       tGCTCGCGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCGCGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCAt                                                 <  1:1593116/139‑1 (MQ=255)
                                                                                        gCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCACGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATc                                                <  1:76009/139‑1 (MQ=255)
                                                                                                gTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCACGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCg                                        >  1:7502/1‑139 (MQ=255)
                                                                                                 tCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCACGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCgg                                       >  2:451537/1‑139 (MQ=255)
                                                                                                 tCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCACGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCgg                                       <  1:85793/139‑1 (MQ=255)
                                                                                                   tgtgCCATCACGCCGTCATCGCACGCCACCACCAACAGCACGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGcc                                     >  1:456661/1‑139 (MQ=255)
                                                                                                     tgCCATCACGCCGTCATCGCACGCCACCACCAACAGCGCGTGATCGATACCACCAACGCCCGCCAGCATGt                                                                                                       >  1:55543/1‑71 (MQ=255)
                                                                                                     tgCCATCACGCCGTCATCGCACGCCACCACCAACAGCGCGTGATCGATACCACCAACGCCCGCCAGCATGt                                                                                                       <  2:55543/71‑1 (MQ=255)
                                                                                                     tgCCATCACGCCGTCATCGCACGCCACCACCAACAGCGCGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAg                                   >  1:329115/1‑139 (MQ=255)
                                                                                                     tgCCATCACGCCGTCATCGCACGCCACCACCAACAGCACGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAg                                   <  1:451537/139‑1 (MQ=255)
                                                                                                          tCACGCCGTCATCGCACTCCACCACCAACAGCGCGTGATCGATACCACCAACGCCCGCCAGCATGCTGGAAAGAAACTTTTCATGACCGGGAACGTCTATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGGc                              <  1:1017167/139‑1 (MQ=255)
                                                                                                          tCACGCCGTCATCGCACGCCACCACCAACAGCGCGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCg                                        <  2:418129/129‑1 (MQ=255)
                                                                                                          tCACGCCGTCATCGCACGCCACCACCAACAGCGCGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCg                                        >  1:418129/1‑129 (MQ=255)
                                                                                                          tCACGCCGTCATCGCACGCCACCACCAACAGCGCGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGGc                              >  2:940841/1‑139 (MQ=255)
                                                                                                          tCACGCCGTCATCGCACGCCACCACCAACAGCACGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGGc                              >  1:649199/1‑139 (MQ=255)
                                                                                                          tCACGCCGTCATCGCACGCCACCACCAACAGCACGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGGc                              >  2:819508/1‑139 (MQ=255)
                                                                                                               ccGTCATCGCACGCCACCACCAACAGCACGTGATCGATACCACCAACGCCCGCCAGCa                                                                                                          <  1:711584/58‑1 (MQ=255)
                                                                                                               ccGTCATCGCACGCCACCACCAACAGCACGTGATCGATACCACCAACGCCCGCCAGCa                                                                                                          >  2:711584/1‑58 (MQ=255)
                                                                                                                cGTCATCGCACGCCACCACCAACAGCACGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGGCATAGcc                        <  2:456661/139‑1 (MQ=255)
                                                                                                                 gTCATCGCACGCCACCACCAACAGCGCGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGGCATAGCCg                       <  2:187992/139‑1 (MQ=255)
                                                                                                                 gTCATCGCACGCCACCACCAACAGCACGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCACGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGGCATAGCCg                       <  2:1082590/139‑1 (MQ=255)
                                                                                                                      cGCACGCCACCACCAACAGCACGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGGCATAGCCg                       >  2:602419/1‑134 (MQ=255)
                                                                                                                      cGCACGCCACCACCAACAGCACGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGGCATAGCCg                       <  1:602419/134‑1 (MQ=255)
                                                                                                                        cACGCCACCACCAACAGCGCGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGGCATAGCCGAGATCGa                <  2:1347024/139‑1 (MQ=255)
                                                                                                                            ccaccaccaACAGCGCGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGGCATAGCCGAGATCGATGGt            >  1:212891/1‑139 (MQ=255)
                                                                                                                                   caACAGCGCGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGGCATAGCCGAGATCGATGGTCATGCcg     <  1:515393/139‑1 (MQ=255)
                                                                                                                                      cAGCGCGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGGCATAGCCGAGATCGATGGTCATGCCgcgc  <  1:70084/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                      
TCATCAACACGCGCTTCGTCCACGCGATCGGCTTTGGTCAGCGCCACTGTCAGCATCGGGTTACCGGTCAGCTGCAAAATCGCCAGATGCTCACGGGTCTGTGCCATCACGCCGTCATCGCACGCCACCACCAACAGCGCGTGATCGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGGCATAGCCGAGATCGATGGTCATGCCGCGC  >  NC_000913/3759491‑3759763

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