Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,759,734 G→A A43V (GCC→GTC)  selB ← selenocysteyl‑tRNA‑specific translation elongation factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,759,7340GA100.0% 1161.4 / NA 283A43V (GCC→GTC) selBselenocysteyl‑tRNA‑specific translation elongation factor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (135/148);  total (135/148)

Only 100 of 283 total aligned reads displayed.

CGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGGCATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGCCTGCAATAAGGTTGTTTTGCCGTG  >  NC_000913/3759636‑3759831
                                                                                                  |                                                                                                 
 gATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACAt                                                                                                <  1:4932488/101‑1 (MQ=255)
 gATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACAt                                                                                                >  1:4932968/1‑101 (MQ=255)
  aTACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATa                                                                                               <  2:4932690/101‑1 (MQ=255)
       accaACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCga                                                                                          >  2:4932973/1‑101 (MQ=255)
        ccaACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAg                                                                                              <  2:4932590/96‑1 (MQ=255)
          aaCGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGAt                                                                                       >  1:4932978/1‑101 (MQ=255)
           aCGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATc                                                                                      >  1:4932982/1‑101 (MQ=255)
             gcccgccAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGa                                                                                    <  1:4932440/101‑1 (MQ=255)
              cccgccAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGAt                                                                                   >  2:4932984/1‑101 (MQ=255)
                 gccAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGt                                                                                >  2:4932989/1‑101 (MQ=255)
                  ccAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTc                                                                               <  2:4932553/101‑1 (MQ=255)
                   cAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCa                                                                              >  1:4932991/1‑101 (MQ=255)
                   cAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCa                                                                              >  2:4932993/1‑101 (MQ=255)
                    aGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGcc                                                                                            <  2:4932514/86‑1 (MQ=255)
                       aTGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGcc                                                                          >  1:4933001/1‑101 (MQ=255)
                       aTGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGcc                                                                          >  1:4932998/1‑101 (MQ=255)
                         gTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCcgc                                                                        >  1:4933008/1‑101 (MQ=255)
                         gTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCcgc                                                                        >  1:4933006/1‑101 (MQ=255)
                         gTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCcgc                                                                        <  2:4932616/101‑1 (MQ=255)
                          ttGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCcgcg                                                                       >  2:4933011/1‑101 (MQ=255)
                           tGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCgcgc                                                                      >  1:4933012/1‑101 (MQ=255)
                             gaaagaaaCTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCtt                                                                    <  1:4932311/101‑1 (MQ=255)
                             gaaagaaaCTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCtt                                                                    <  1:4932622/101‑1 (MQ=255)
                              aaagaaaCTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCttt                                                                   >  2:4933016/1‑101 (MQ=255)
                                agaaaCTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCttttt                                                                 >  2:4933020/1‑101 (MQ=255)
                                agaaaCTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCttttt                                                                 >  1:4933021/1‑101 (MQ=255)
                                 gaaaCTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCtttttt                                                                <  2:4932884/101‑1 (MQ=255)
                                 gaaaCTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCtttttt                                                                <  1:4932566/101‑1 (MQ=255)
                                 gaaaCTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCtttttt                                                                <  1:4932628/101‑1 (MQ=255)
                                  aaaCTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCttttttt                                                               >  2:4933025/1‑101 (MQ=255)
                                   aaCTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTttc                                                              >  1:4933027/1‑101 (MQ=255)
                                    aCTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTttct                                                             <  1:4932607/101‑1 (MQ=255)
                                    aCTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTttct                                                             <  2:4932769/101‑1 (MQ=255)
                                    aCTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTttct                                                             <  2:4932630/101‑1 (MQ=255)
                                     cTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTttctt                                                            <  2:4932783/101‑1 (MQ=255)
                                     cTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTttctt                                                            <  2:4932609/101‑1 (MQ=255)
                                       tttCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTcc                                                          >  1:4933031/1‑101 (MQ=255)
                                        ttCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTc                                                                               >  2:4933036/1‑79 (MQ=255)
                                        ttCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCg                                                         >  1:4933034/1‑101 (MQ=255)
                                          cATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTttctt                                                            >  2:4933040/1‑96 (MQ=255)
                                          cATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGc                                                       <  2:4932708/101‑1 (MQ=255)
                                           aTGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCa                                                      >  1:4933041/1‑101 (MQ=255)
                                           aTGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCa                                                      <  2:4932152/101‑1 (MQ=255)
                                            tGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAg                                                     >  2:4933045/1‑101 (MQ=255)
                                             gACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGa                                                    >  1:4933046/1‑101 (MQ=255)
                                             gACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGa                                                    <  1:4932522/101‑1 (MQ=255)
                                                cGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACgg                                                 >  1:4933050/1‑101 (MQ=255)
                                                 gggAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGt                                                <  1:4932618/101‑1 (MQ=255)
                                                     aCGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGc                                            <  1:4932541/101‑1 (MQ=255)
                                                     aCGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGc                                            <  2:4932601/101‑1 (MQ=255)
                                                      cGTCGATAAAACCAGGCACGCGACCCTCCGGCTGCGGCCAGTAGACATAGCCga                                                                                          >  1:4933056/1‑54 (MQ=39)
                                                      cGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCa                                           >  2:4933055/1‑101 (MQ=255)
                                                       gTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCAt                                          <  2:4932464/101‑1 (MQ=255)
                                                         cGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCAttt                                        >  1:4933062/1‑101 (MQ=255)
                                                         cGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCAttt                                        >  2:4933061/1‑101 (MQ=255)
                                                           aTAAAGCCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTAc                                      <  2:4932697/101‑1 (MQ=255)
                                                           aTAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTAc                                      >  1:4933066/1‑101 (MQ=255)
                                                           aTAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTAc                                      <  2:4932434/101‑1 (MQ=255)
                                                             aaaaCCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCa                                                                              >  2:4933070/1‑59 (MQ=255)
                                                             aaaaCCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTAc                                      >  1:4933069/1‑99 (MQ=255)
                                                             aaaaCCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGc                                    <  1:4932857/101‑1 (MQ=255)
                                                              aaaCCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGcc                                   <  2:4932903/101‑1 (MQ=255)
                                                               aaCCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCAttt                                        <  1:4932774/95‑1 (MQ=255)
                                                               aaCCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCa                                  <  1:4932600/101‑1 (MQ=255)
                                                                aCCAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAg                                 <  2:4932846/101‑1 (MQ=255)
                                                                  cAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCttttt                                                                 <  1:4932393/67‑1 (MQ=255)
                                                                  cAGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTa                               <  1:4932650/101‑1 (MQ=255)
                                                                   aGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTaa                              <  1:4932826/101‑1 (MQ=255)
                                                                   aGGCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTaa                              >  2:4933074/1‑101 (MQ=255)
                                                                     gCACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATc                            <  2:4932803/101‑1 (MQ=255)
                                                                      cACGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCa                                              <  1:4932656/82‑1 (MQ=255)
                                                                       aCGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGc                          <  2:4932728/101‑1 (MQ=255)
                                                                       aCGCGACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGc                          <  2:4932673/101‑1 (MQ=255)
                                                                         gcgACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCa                                  >  2:4933078/1‑91 (MQ=255)
                                                                          cgACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGCCTg                       <  1:4932791/101‑1 (MQ=255)
                                                                          cgACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGCCTg                       <  1:4932681/101‑1 (MQ=255)
                                                                           gACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGCCTGc                      >  2:4933081/1‑101 (MQ=255)
                                                                           gACCATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGCCTGc                      >  1:4933080/1‑101 (MQ=255)
                                                                             ccATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGCCTGCaa                    <  2:4932636/101‑1 (MQ=255)
                                                                              cATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCTCCTGCAAt                   >  1:4933086/1‑101 (MQ=255)
                                                                              cATCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGCCTGCAAt                   >  2:4933087/1‑101 (MQ=255)
                                                                               aTCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGCCTGCAATa                  <  1:4932716/101‑1 (MQ=255)
                                                                                tCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGCCTGCAATaa                 <  1:4932546/101‑1 (MQ=255)
                                                                                tCCGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGCCTGCAATaa                 <  1:4932755/101‑1 (MQ=255)
                                                                                  cGGCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGCCTGCAATAAgg               <  1:4932699/101‑1 (MQ=255)
                                                                                   ggCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGCCTGCAATAAGgt              >  1:4933092/1‑101 (MQ=255)
                                                                                    gCTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGCCTGCAATAAGgtt             <  1:4932669/101‑1 (MQ=255)
                                                                                     cTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCAt                                          >  2:4933093/1‑71 (MQ=255)
                                                                                     cTGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGCCTGCAATAAGgttg            <  1:4932582/101‑1 (MQ=255)
                                                                                      tGCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGCCTGCAATAAGgttgt           <  1:4932975/101‑1 (MQ=255)
                                                                                       gCGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGCCTGCAATAAGgttgtt          <  1:4932555/101‑1 (MQ=255)
                                                                                        cGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGCCTGCAATaa                 <  1:4932722/93‑1 (MQ=255)
                                                                                        cGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGCCTGCAATAAGGTTGttt         <  2:4932661/101‑1 (MQ=255)
                                                                                        cGGCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGCCTGCAATAAGGTTGttt         <  1:4932943/101‑1 (MQ=255)
                                                                                         ggCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGCCTGCAATAAGGTTGtttt        <  2:4932540/101‑1 (MQ=255)
                                                                                         ggCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGCCTGCAATAAGGTTGtttt        >  2:4933096/1‑101 (MQ=255)
                                                                                          gCCAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGCCTGCAATAAGGTTGTTTTg       <  2:4932744/101‑1 (MQ=255)
                                                                                            cAGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGCCTGCAATAAGgttgtt          >  1:4933100/1‑96 (MQ=255)
                                                                                             aGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGCCTGCAATAAGGTTGTTTTg       >  1:4933104/1‑98 (MQ=255)
                                                                                             aGTAGACATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGCCTGCAATAAGGTTGTTTTGCCg    >  2:4933103/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
CGATACCACCAACGCCCGCCAGCATGTTGGAAAGAAACTTTTCATGACCGGGAACGTCGATAAAACCAGGCACGCGACCATCCGGCTGCGGCCAGTAGGCATAGCCGAGATCGATGGTCATGCCGCGCTTTTTTTCTTCCGGCAGACGGTCAGCATTTACGCCAGTAATCGCCTGCAATAAGGTTGTTTTGCCGTG  >  NC_000913/3759636‑3759831

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