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% 876.9 / NA 205A162G (GCA→GGA) btuBvitamin B12/cobalamin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (98/107);  total (98/107)

Only 100 of 228 total aligned reads displayed.

TATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCAT  >  NC_000913/4162048‑4162243
                                                                                                  |                                                                                                 
tatCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGGTGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTcaggag                                                                                                 >  2:3019107/1‑101 (MQ=255)
tatCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTcaggag                                                                                                 <  1:2882991/101‑1 (MQ=255)
 atCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGGGAATATCATCACGACGCGCGGTGAACCCGGAACGGAAAATTCTGGAgg                                                                                                >  2:1307789/1‑101 (MQ=255)
  tCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAggg                                                                                               <  1:5007008/101‑1 (MQ=255)
   ccGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAgggg                                                                                              >  2:1767630/1‑100 (MQ=255)
    cGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAggggg                                                                                             >  2:274175/1‑99 (MQ=255)
    cGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTg                                                                                             >  1:4197837/1‑101 (MQ=255)
    cGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTg                                                                                             >  2:3988575/1‑101 (MQ=255)
     gTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCAAGACGGGCGATGAAGCCGGAACGGAAATTTCAGgagggagg                                                                                            >  2:4200979/1‑101 (MQ=255)
     gTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGGGGTGAATTTCATCACGACGCCCGATGAACCCGGAACGGAAATTTCAGGGGGGTgg                                                                                            >  2:3547566/1‑101 (MQ=255)
      tGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTggg                                                                                           <  2:65926/101‑1 (MQ=255)
       gngCNNCGCTCCGCNGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATAATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGggggggc                                                                                          >  2:1111118/1‑100 (MQ=255)
        ggCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGa                                                                                         <  2:27172/101‑1 (MQ=255)
        ggCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGa                                                                                         <  2:2523369/101‑1 (MQ=255)
         gCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGaa                                                                                        >  2:519066/1‑101 (MQ=255)
          ccGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAg                                                                                       <  2:2807936/101‑1 (MQ=255)
          ccGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAg                                                                                       >  1:3020593/1‑101 (MQ=255)
          ccGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAg                                                                                       <  1:3628153/101‑1 (MQ=255)
           cgcgcgCCGCTGTTTATAGTTCCGATGCAATAGGCGGGGTGGTGCATATCATCACGACGCGAGATGAACCCGGAACGGAAATTTCAGGAGGGTCGGGaatc                                                                                      <  2:535769/101‑3 (MQ=255)
            gcgcTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGGATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCa                                                                                     >  2:6956342/1‑101 (MQ=255)
              gcTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAAt                                                                                   >  1:1953622/1‑101 (MQ=255)
              gcTCCGCTGTTTATGGTTCCGATGCAATAGGCGGAGGGGGGAATATCATGACGACGCCCGATTAACCCGGAACGGAGATTTCAGGAGGGTGGGGAAGaaaa                                                                                   >  2:2168064/1‑100 (MQ=255)
               cTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATa                                                                                  <  1:1416699/101‑1 (MQ=255)
                 ccGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGt                                                                                <  1:730538/101‑1 (MQ=255)
                 ccGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGt                                                                                >  1:6374572/1‑101 (MQ=255)
                 ccGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGGGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTAAGGAGGGCGGGGAAGCAATAGt                                                                                >  1:912852/1‑101 (MQ=255)
                 ccGCTGGTTATGGTACCGATTCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAAGATGGTGGGAAAGCAATAGt                                                                                <  2:722602/101‑1 (MQ=255)
                  cGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGtt                                                                               <  2:5831986/101‑1 (MQ=255)
                    cTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTAt                                                                             >  1:1202973/1‑101 (MQ=255)
                    cTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTAt                                                                             <  2:4222586/101‑1 (MQ=255)
                     tGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATc                                                                            <  2:32371/101‑1 (MQ=255)
                      gTTTATGGGTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCa                                                                           >  1:5786491/1‑101 (MQ=255)
                       tttATGGTTCCGATGCAATAGGCGGGGGGGTGAATATCATCACGACGCGCGATGAACCCGGACCGGAAATTTTAGGAGGGCGGGGAAGCAATAGTTATCAg                                                                          >  2:1158330/1‑101 (MQ=255)
                        ttATGGTTCCGATGCAATAGGCGGGGTGGTTAATATCATCACGACTCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGa                                                                         <  1:4690885/101‑1 (MQ=255)
                        ttATGGTTCCGATGCAATAGGCGGGGTGGGGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGa                                                                         >  2:1918612/1‑101 (MQ=255)
                         tctGGTTCGGATGCAATAGGCGGGGTTGTGAATATCATCTCGACGCGCGATGAACCCGCAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGaa                                                                        <  1:1431422/99‑1 (MQ=255)
                         tATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGaa                                                                        >  1:3829303/1‑101 (MQ=255)
                           tGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACt                                                                      >  1:3347178/1‑101 (MQ=255)
                             gTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTat                                                                    >  1:3370703/1‑101 (MQ=255)
                                ccGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatgat                                                                 <  1:1027181/101‑1 (MQ=255)
                                 cGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTatgatg                                                                <  1:6872620/101‑1 (MQ=255)
                                  gATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGt                                                               >  2:1569318/1‑101 (MQ=255)
                                  gATGCAAGAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGCACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGt                                                               <  1:4420938/101‑1 (MQ=255)
                                   aTGCAATAGGCGGGGTGTTGAATATCATCACGACGCGCGATGAGCCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtc                                                              <  2:4317197/101‑1 (MQ=255)
                                    tGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtct                                                             <  1:1139874/101‑1 (MQ=255)
                                    tGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGAAGtct                                                             <  2:2451082/101‑1 (MQ=255)
                                    tGCAATAGGCGGGGCGGGGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGtcc                                                             >  2:3892140/1‑100 (MQ=255)
                                      cAATAGGCGGGGTGGTGAATATCATCACGACGCGAGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTctct                                                           <  1:667139/101‑1 (MQ=255)
                                      cAATAGGCGGGGGGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTctct                                                           >  1:2268037/1‑101 (MQ=255)
                                       aaTAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAGCGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTa                                                          >  1:4874843/1‑101 (MQ=255)
                                       aaTAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTa                                                          >  2:5098439/1‑101 (MQ=255)
                                       aaTAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTa                                                          <  1:1675593/101‑1 (MQ=255)
                                        aTAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTAc                                                         <  2:5781742/101‑1 (MQ=255)
                                        aTAGGCGGGGTGGTAAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCGAc                                                         <  2:4069008/101‑1 (MQ=255)
                                          aGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgc                                                       >  2:226993/1‑101 (MQ=255)
                                             cGGGGCGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGCGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcagc                                                    >  1:6200089/1‑101 (MQ=255)
                                              ggggTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcagca                                                   >  1:5280124/1‑101 (MQ=255)
                                               gggTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaa                                                  <  1:2684218/101‑1 (MQ=255)
                                                ggtggtGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaac                                                 >  1:4647045/1‑101 (MQ=255)
                                                ggtggtGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaac                                                 >  2:4121529/1‑101 (MQ=255)
                                                 gtggtgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaaca                                                <  1:1553026/101‑1 (MQ=255)
                                                   ggtgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCCGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCaacaac                                              >  1:6517477/1‑101 (MQ=255)
                                                   ggtgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCaacaac                                              >  2:644644/1‑101 (MQ=255)
                                                    gtgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACt                                             <  2:5043804/101‑1 (MQ=255)
                                                     tgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTg                                            <  1:3811441/101‑1 (MQ=255)
                                                     ggAATATCATCACGACGCGCGATGACCCCGGAACGGAAATTTCAGGAGCGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCCACGCAGCAACAACTg                                            <  2:4377624/100‑1 (MQ=255)
                                                      gAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTgg                                           <  1:4732258/101‑1 (MQ=255)
                                                        atatCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTgggg                                         >  1:2557285/1‑101 (MQ=255)
                                                         tatCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTggggg                                        <  1:2135215/101‑1 (MQ=255)
                                                          atcatcACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGa                                       <  1:281869/101‑1 (MQ=255)
                                                           tcatcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGaa                                      >  2:678074/1‑100 (MQ=255)
                                                           tcatcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGAt                                      >  1:3268361/1‑101 (MQ=255)
                                                           tcatcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGAt                                      >  1:3670030/1‑101 (MQ=255)
                                                            catcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATa                                     <  1:3498144/101‑1 (MQ=255)
                                                             atcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATaa                                    >  2:3819230/1‑101 (MQ=255)
                                                              tcaCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAg                                   <  2:1477207/101‑1 (MQ=255)
                                                               caCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGa                                  <  1:1918612/101‑1 (MQ=255)
                                                                acgacgCGCGATGAACCCGGAACGGAAATTTCAGGAGGGGGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGac                                 >  1:3766682/1‑101 (MQ=255)
                                                                  gacgCGCGATGAACCCGGAACGTAAATTTCAGGAGTGTGGGGAAGCAATAGTTATCAGAACTAGGATGTCTCTACGCAGCAACAACTGGGGGATAAGacac                               <  2:3115994/101‑1 (MQ=255)
                                                                  gacgCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGacac                               >  1:5425211/1‑101 (MQ=255)
                                                                   gcgcgcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACg                              <  2:5242798/100‑1 (MQ=255)
                                                                   acgCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACg                              <  2:3444914/101‑1 (MQ=255)
                                                                   acgCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACg                              <  2:4829095/101‑1 (MQ=255)
                                                                     gcgcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACggg                            >  2:1439097/1‑101 (MQ=255)
                                                                     gcgcgATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACggg                            >  2:6484836/1‑101 (MQ=255)
                                                                       gcgATGAACCCGGAACGGAAATTTCAGGTGGGGGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTa                          >  1:1937244/1‑101 (MQ=255)
                                                                          aTGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTaa                         >  2:5533934/1‑99 (MQ=255)
                                                                           tGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGc                      >  2:2216328/1‑101 (MQ=255)
                                                                              aCCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGt                   <  1:1111118/101‑1 (MQ=255)
                                                                               cccGGAACGGAAATTTCAGGAGGGGGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGtt                  >  1:6055761/1‑101 (MQ=255)
                                                                                ccGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTg                 >  1:5550749/1‑101 (MQ=255)
                                                                                acggaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTg                 <  2:2198857/100‑1 (MQ=255)
                                                                                  ggaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTggg               <  1:5098439/101‑1 (MQ=255)
                                                                                   gaacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGc              >  1:4944576/1‑101 (MQ=255)
                                                                                    aacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCg             >  2:3105303/1‑101 (MQ=255)
                                                                                    aacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCg             >  1:6977863/1‑101 (MQ=255)
                                                                                        gaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTa         >  2:501533/1‑101 (MQ=255)
                                                                                        gaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTa         <  2:6327721/101‑1 (MQ=255)
                                                                                         aaaTTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTAt        <  1:871534/101‑1 (MQ=255)
                                                                                              tCAGGAGGGTGGGAAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCa   <  1:3964954/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
TATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCAT  >  NC_000913/4162048‑4162243

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