Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,162,146 C→G A162G (GCA→GGA)  btuB → vitamin B12/cobalamin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,162,1460CG100.0% 2095.1 / NA 507A162G (GCA→GGA) btuBvitamin B12/cobalamin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (256/251);  total (256/251)

Only 100 of 554 total aligned reads displayed.

ATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCAT  >  NC_000913/4162049‑4162243
                                                                                                 |                                                                                                 
   cGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTg                                                                                             <  1:4233975/101‑1 (MQ=255)
      gggCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTgnnn                                                                                          >  2:10015908/1‑98 (MQ=255)
        gCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGaa                                                                                        >  1:6096648/1‑101 (MQ=255)
          cgcgCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCTCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGc                                                                                      <  2:3360362/101‑1 (MQ=255)
            cgcTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCaa                                                                                    >  1:10865751/1‑101 (MQ=255)
            cgcTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCaa                                                                                    <  1:6457049/101‑1 (MQ=255)
            cgcTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGGGGGGAAGCaa                                                                                    >  2:789082/1‑101 (MQ=255)
             gcTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCaaa                                                                                   >  1:1359341/1‑100 (MQ=255)
             gcTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAAt                                                                                   >  1:2536441/1‑101 (MQ=255)
              cTCCGCTGTTTATGGTTCCGATGCGATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATa                                                                                  >  1:9127512/1‑101 (MQ=255)
              cTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATa                                                                                  >  1:7432368/1‑101 (MQ=255)
               tCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGACCCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAg                                                                                 <  1:3248770/101‑1 (MQ=255)
               tCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAg                                                                                 <  1:11828030/101‑1 (MQ=255)
               tCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAg                                                                                 >  2:6631843/1‑101 (MQ=255)
               tCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAg                                                                                 <  1:8129974/101‑1 (MQ=255)
                nnncTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGt                                                                                <  2:7585539/98‑1 (MQ=255)
                ccGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGt                                                                                <  1:5370315/101‑1 (MQ=255)
                 cGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGtt                                                                               >  1:1205107/1‑101 (MQ=255)
                   cTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGAGGGGAAGCAATAGTTAt                                                                             >  1:7789686/1‑101 (MQ=255)
                    tGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAAGGGAAATTTCAGGAGGGAGGGGAAGCAATAGTTATc                                                                            >  2:6052901/1‑101 (MQ=255)
                    tGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATc                                                                            >  1:13351860/1‑101 (MQ=255)
                    tGTTTATGGTTCCGATGCAATAGGCGGGGTGGGGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGGGGGGAAGCAATAGTTAta                                                                            >  2:3834057/1‑100 (MQ=255)
                     gTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCa                                                                           >  1:8485841/1‑101 (MQ=255)
                     gTTTATGGTTCCGATGCAATAGGCGGGGTGGGGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGGGGGGAAGCAATAGTTATCa                                                                           >  2:2245937/1‑101 (MQ=255)
                      tttATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAg                                                                          >  1:7205006/1‑101 (MQ=255)
                       ttATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGa                                                                         >  1:3470112/1‑101 (MQ=255)
                       ttATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGa                                                                         >  1:1553652/1‑101 (MQ=255)
                        tATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGaa                                                                        >  2:1159404/1‑101 (MQ=255)
                         aTGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGCGGGTGGGGAAGCAATAGTTATCAGAAc                                                                       <  1:13085198/101‑1 (MQ=255)
                          tGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACt                                                                      <  1:14564014/101‑1 (MQ=255)
                           ggTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAGCCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTa                                                                     <  1:10554777/101‑1 (MQ=255)
                            gTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTat                                                                    >  1:12227534/1‑101 (MQ=255)
                               nnnaTGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatgat                                                                 <  2:10225134/98‑1 (MQ=255)
                                cGATGCAATAGGCGGGGGGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGCGGGGAAGCAATAGTTATCAGAACTaggatg                                                                >  2:2559415/1‑101 (MQ=255)
                                 gATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatgannn                                                               >  2:8955484/1‑98 (MQ=255)
                                  aTGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtc                                                              >  2:2738621/1‑101 (MQ=255)
                                    gCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtctc                                                            >  2:7087989/1‑101 (MQ=255)
                                     nnntAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTctct                                                           <  2:11846464/98‑1 (MQ=255)
                                      nnnaGGCGGGGTGGTGAATATCATCACGACGCGCGATGACCCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTa                                                          <  2:8209587/98‑1 (MQ=255)
                                      aaTAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtctnnn                                                          >  2:14160587/1‑98 (MQ=255)
                                      aaTAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTa                                                          >  1:281620/1‑101 (MQ=255)
                                       aTAGGCGGGGTTGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTAc                                                         <  2:2877853/101‑1 (MQ=255)
                                       aTAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTAc                                                         <  1:14135694/101‑1 (MQ=255)
                                        tAGGCGGGGTGGGGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGGGGGGAAGCAATAGTTATCAGAACTATGATGTctctnnn                                                        >  2:14448075/1‑98 (MQ=255)
                                          ggCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTGATCAGAACTATGATGTCTCTACgcg                                                      >  2:4582204/1‑100 (MQ=255)
                                            nnnggtggtGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcagt                                                    <  2:9161567/98‑2 (MQ=255)
                                             ggggTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcagca                                                   <  2:3210989/101‑1 (MQ=255)
                                               nnnggtgAATATCATCACGACGCGCGATGAACCCGGAACGGCAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcagcacc                                                 <  2:13099563/98‑3 (MQ=255)
                                               ggtggtGAATATCGTCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcagcnnn                                                 >  2:13858976/1‑98 (MQ=255)
                                                gtggtgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaaca                                                <  1:2876711/101‑1 (MQ=255)
                                                  ggtgAATATCATCCCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCaacaac                                              <  1:4702183/101‑1 (MQ=255)
                                                  ggtgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCaacaac                                              <  1:11162298/101‑1 (MQ=255)
                                                    tgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTg                                            <  1:2192823/101‑1 (MQ=255)
                                                    tgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTg                                            >  1:5246614/1‑101 (MQ=255)
                                                     gAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCaacaacnnn                                           >  2:12492947/1‑98 (MQ=255)
                                                     gAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTgg                                           <  1:7047380/101‑1 (MQ=255)
                                                     gAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTgg                                           >  1:394215/1‑101 (MQ=255)
                                                      aaTATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaa                                                  >  1:10134900/1‑93 (MQ=255)
                                                       atatCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGGGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTgggg                                         >  1:13806817/1‑101 (MQ=255)
                                                        tatCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTggnnn                                        >  2:7470941/1‑98 (MQ=255)
                                                        nnncatcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTggggg                                        <  2:10988167/98‑1 (MQ=255)
                                                          tcatcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTggggnnn                                      >  2:13525448/1‑98 (MQ=255)
                                                          tcatcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGa                                       <  1:5026580/100‑1 (MQ=255)
                                                           catcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTgggggnnn                                     >  2:10327637/1‑98 (MQ=255)
                                                              nnnggcgcgcgATGAACGCGGAACGGAAATTTCAGGAGCGTCGGGAAGCAATATTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGa                                  <  2:10549704/96‑1 (MQ=255)
                                                              caCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTAACAGAACTATGATGTCTCTACGCAGCAACAACTGCGGGATannn                                  >  2:7973082/1‑98 (MQ=255)
                                                               nnnacgCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGac                                 <  2:11231625/98‑1 (MQ=255)
                                                               acgacgCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGac                                 <  2:1923239/101‑1 (MQ=255)
                                                                 gacgCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGacac                               >  1:12531001/1‑101 (MQ=255)
                                                                  nnncgcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACg                              <  2:14848090/98‑1 (MQ=255)
                                                                    gcgcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACggg                            <  2:4438612/101‑1 (MQ=255)
                                                                     cgcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGt                           >  2:5118521/1‑101 (MQ=255)
                                                                     cgcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGGCTCTACGCAGCAACAACTGGGGGATAAGACACGGGt                           >  1:3721443/1‑101 (MQ=255)
                                                                      gcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTa                          >  1:8787722/1‑101 (MQ=255)
                                                                        gATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAAc                        >  2:6110303/1‑101 (MQ=255)
                                                                         nnnaaCCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTaagg                       <  2:12227534/98‑3 (MQ=255)
                                                                         aTGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACACCTGGGGGATAAGACACGGGTAACg                       <  1:7908703/101‑1 (MQ=255)
                                                                         aTGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACg                       <  1:5584796/101‑1 (MQ=255)
                                                                          tGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGc                      <  1:11632904/101‑1 (MQ=255)
                                                                           nnncccGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCt                     <  2:9609111/98‑1 (MQ=255)
                                                                           gAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCt                     >  1:2676790/1‑101 (MQ=255)
                                                                             aCCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGt                   <  2:4131000/101‑1 (MQ=255)
                                                                              nnnggaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGtt                  <  2:11084859/98‑1 (MQ=255)
                                                                              cccGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGtt                  >  1:930539/1‑101 (MQ=255)
                                                                               ccGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTg                 >  1:10333366/1‑101 (MQ=255)
                                                                               ccGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTg                 >  1:9710755/1‑101 (MQ=255)
                                                                                cggaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGtnnn                >  2:8595395/1‑98 (MQ=255)
                                                                                 ggaacggaaATTTCAGGAGGGTGGGGAAGCAATAGGTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATATGACACGGGTAACGCTGTTggg               >  2:5705999/1‑101 (MQ=255)
                                                                                  gaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGc              <  1:6727686/101‑1 (MQ=255)
                                                                                  gaacggaaATTTCAGGAGGGTGGGGAAGCAAGAGTTATCAGACCTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGc              <  1:2015710/101‑1 (MQ=255)
                                                                                   aacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCg             <  2:5496821/101‑1 (MQ=255)
                                                                                   aacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGAGATAAGACACGGGTAACGCTGTTggnnn             >  2:12862051/1‑98 (MQ=255)
                                                                                      ggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGAtt          >  2:731003/1‑101 (MQ=255)
                                                                                      ggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGAtt          <  2:6386228/101‑1 (MQ=255)
                                                                                       gaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTa         <  1:584495/101‑1 (MQ=255)
                                                                                       gaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTa         >  1:1674684/1‑101 (MQ=255)
                                                                                        taaTTTCAGGCGGGTGGGGAAGCAATAGTTATCAGAACTATGTTGTCTCTACGCAGCAATAACTGGGGGATAAAACACGGGTACTGCGGTTGAGCGATTAt        >  2:6896862/2‑101 (MQ=255)
                                                                                        aaaTTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTAt        >  1:11380404/1‑101 (MQ=255)
                                                                                            ttCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGccc    <  2:3615763/101‑1 (MQ=255)
                                                                                            ttCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGccc    <  2:5282144/101‑1 (MQ=255)
                                                                                                 |                                                                                                 
ATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCAT  >  NC_000913/4162049‑4162243

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