Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,537,405 +GC intergenic (+61/‑394) b3407 → / → feoA hypothetical protein/ferrous iron transport protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,537,4041.C100.0% 743.2 / NA 210intergenic (+60/‑395)b3407/feoAhypothetical protein/ferrous iron transport protein A
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (100/110);  total (100/110)
*NC_0009133,537,4042.G100.0% 743.0 / NA 210intergenic (+60/‑395)b3407/feoAhypothetical protein/ferrous iron transport protein A
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (100/110);  total (100/110)

Only 100 of 223 total aligned reads displayed.

GCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCA‑‑CCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATA  >  NC_000913/3537312‑3537487
                                                                                               ||                                                                                   
 cgcTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCA‑‑c                                                                                    <  2:447363/93‑1 (MQ=255)
 cgcTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCA‑‑c                                                                                    <  1:5286850/93‑1 (MQ=255)
   cTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCA‑‑cgc                                                                                  <  2:5028082/93‑3 (MQ=255)
           ggcAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTg                                                                            <  1:3354424/93‑1 (MQ=38)
             cAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGc                                                                          <  1:5620444/93‑1 (MQ=37)
             cAATGGGAAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGc                                                                          <  2:1429567/93‑1 (MQ=37)
              aaTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCa                                                                         >  1:5170770/1‑93 (MQ=37)
               aTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCaa                                                                        >  2:3476575/1‑93 (MQ=37)
               aTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCaa                                                                        >  2:3180887/1‑93 (MQ=37)
               aTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCaa                                                                        >  2:5660918/1‑93 (MQ=37)
               aTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCaa                                                                        <  2:4345682/93‑1 (MQ=37)
               aTGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCaa                                                                        >  1:682374/1‑93 (MQ=37)
                tGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCaac                                                                       >  2:5309415/1‑93 (MQ=37)
                tGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCaac                                                                       >  1:4347092/1‑93 (MQ=37)
                 gggCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAAATTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGGaacg                                                                      >  1:459355/1‑93 (MQ=37)
                  ggCAAAAAACGTTAAACGCCCGTACTGGCCTACGTTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCaacga                                                                     <  1:2616029/93‑1 (MQ=37)
                  ggCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCaacga                                                                     >  1:5077125/1‑93 (MQ=37)
                  ggCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCaacga                                                                     >  2:24101/1‑93 (MQ=37)
                  ggCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCaacga                                                                     >  2:3441875/1‑93 (MQ=37)
                   gCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCaacgaa                                                                    <  1:798121/93‑1 (MQ=38)
                   gCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCaacgaa                                                                    <  1:3281990/93‑1 (MQ=37)
                   gCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCaacgaa                                                                    <  2:2927979/93‑1 (MQ=38)
                    cAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCaacgaac                                                                   >  1:4922903/1‑93 (MQ=37)
                    cAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCaacgaac                                                                   <  2:354690/93‑1 (MQ=37)
                     aaaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCaacgaacg                                                                  >  2:3618386/1‑93 (MQ=37)
                      aaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGc                                                                 >  2:5860271/1‑93 (MQ=37)
                      aaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGc                                                                 >  1:5716726/1‑93 (MQ=37)
                       aacaCGTTAAACGCGCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGcc                                                                <  2:4330066/90‑1 (MQ=25)
                       aaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGcc                                                                <  1:4750676/93‑1 (MQ=37)
                        aaaCGTTAACCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCt                                                               <  1:1850463/93‑1 (MQ=37)
                        aaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCt                                                               >  2:1026243/1‑93 (MQ=37)
                        aaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCt                                                               >  2:3354975/1‑93 (MQ=37)
                        aaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCt                                                               <  1:1119326/93‑1 (MQ=37)
                        aaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCt                                                               <  1:4594771/93‑1 (MQ=37)
                          aCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTg                                                             >  2:1251271/1‑93 (MQ=37)
                            gtaaaaCGCCCTTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGcc                                                           <  1:1446408/90‑1 (MQ=25)
                             ttAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCt                                                          <  1:5240140/93‑1 (MQ=37)
                              tAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTg                                                         <  1:3668540/93‑1 (MQ=37)
                               attcGCCCGTACTGGCCTACGGTTCGTATTTGCACGATATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTGGCCCGa                                                        >  1:5493614/4‑93 (MQ=21)
                                aaCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGAc                                                       >  1:3148921/1‑93 (MQ=38)
                                 aCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACa                                                      <  1:2563187/93‑1 (MQ=37)
                                 aCGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTCGCCTGACa                                                      >  1:4827028/1‑93 (MQ=37)
                                   ccccGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTCACATa                                                    <  1:2092264/92‑1 (MQ=37)
                                    cccGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATaa                                                   <  2:6027164/93‑1 (MQ=37)
                                     ccGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATaaa                                                  >  2:4852685/1‑93 (MQ=38)
                                     ccGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATaaa                                                  <  2:5528975/93‑1 (MQ=37)
                                     ccGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATaaa                                                  >  2:750530/1‑93 (MQ=38)
                                     ccGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATaaa                                                  >  1:2352926/1‑93 (MQ=38)
                                      cGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATaaaa                                                 <  2:1826088/93‑1 (MQ=37)
                                       gTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAg                                                <  1:1626188/93‑1 (MQ=37)
                                        tACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGt                                               <  2:4585664/93‑1 (MQ=37)
                                         aCTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTg                                              >  2:3998498/1‑93 (MQ=37)
                                         aCTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTg                                              <  2:831931/93‑1 (MQ=37)
                                          cTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGc                                             <  2:1063538/93‑1 (MQ=37)
                                          cTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGc                                             >  1:4411616/1‑93 (MQ=37)
                                           tGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGcc                                            >  2:4175126/1‑93 (MQ=37)
                                            ggCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCg                                           <  1:169777/93‑1 (MQ=37)
                                            ggCCTACGGTTCGAATTTGCACGAAATCATAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCg                                           >  2:394628/1‑93 (MQ=38)
                                             gCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCgg                                          >  1:5852215/1‑93 (MQ=38)
                                             gCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCgg                                          <  1:1887480/93‑1 (MQ=37)
                                             gCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCgg                                          >  2:5137737/1‑93 (MQ=38)
                                                tACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGgaga                                       >  1:2997306/1‑93 (MQ=38)
                                                 aCGGTTTGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGaa                                      <  2:1011971/93‑1 (MQ=37)
                                                 aCGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGaa                                      <  2:5783659/93‑1 (MQ=37)
                                                   ggTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAata                                    <  2:4130753/93‑1 (MQ=37)
                                                    gTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAatat                                   <  1:268056/93‑1 (MQ=37)
                                                     ttCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATAtc                                  <  2:892032/93‑1 (MQ=37)
                                                     ttCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATAtc                                  <  1:3892657/93‑1 (MQ=37)
                                                      tCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATAtct                                 >  2:3752002/1‑93 (MQ=38)
                                                       cGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATAtctc                                <  2:1373785/93‑1 (MQ=37)
                                                        gAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTcc                               <  2:5438628/93‑1 (MQ=37)
                                                        gAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTcc                               <  2:676215/93‑1 (MQ=37)
                                                        gAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTcc                               <  2:2095278/93‑1 (MQ=37)
                                                         aaTTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCg                              <  2:369657/93‑1 (MQ=37)
                                                             tGCAGGAAAGGGTAGGGCAGATAGGGCCT‑‑CCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAt                          <  1:1232609/93‑1 (MQ=255)
                                                             tGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAt                          <  1:2411557/93‑1 (MQ=38)
                                                              gCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAtt                         >  2:4886693/1‑93 (MQ=38)
                                                                aCGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAtttt                       <  2:5746260/93‑1 (MQ=38)
                                                                 cGAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCAttttt                      <  2:1762415/93‑1 (MQ=38)
                                                                  gAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTa                     <  1:2212357/93‑1 (MQ=38)
                                                                  gAAATCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTa                     <  1:2798749/93‑1 (MQ=38)
                                                                   aaaTCGTAGGACAGATAAGGCGTCGACG‑‑CCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTAt                    <  1:4494802/93‑1 (MQ=37)
                                                                     aTCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTc                  <  2:301591/93‑1 (MQ=38)
                                                                      tCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTcc                 >  1:605439/1‑93 (MQ=38)
                                                                      tCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTcc                 >  1:3975061/1‑93 (MQ=38)
                                                                      tCGTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTcc                 >  1:3083943/1‑93 (MQ=38)
                                                                        gTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTTTTctcc               >  1:1197968/1‑90 (MQ=37)
                                                                        gTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCcac               <  2:2848195/93‑1 (MQ=38)
                                                                        gTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCcac               >  2:5495642/1‑93 (MQ=38)
                                                                        gTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCcac               <  1:2474206/93‑1 (MQ=38)
                                                                        gTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCcac               <  1:2307559/93‑1 (MQ=38)
                                                                        gTAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCcac               <  1:1138669/93‑1 (MQ=38)
                                                                         tAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCcaca              >  1:413429/1‑93 (MQ=38)
                                                                         tAGGGCAGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCcaca              <  1:3451684/93‑1 (MQ=38)
                                                                             gaaGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCa          <  1:3584488/91‑1 (MQ=38)
                                                                               aGATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCaaa        <  1:5675702/93‑1 (MQ=38)
                                                                                gATAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAAc       >  1:5912889/1‑93 (MQ=39)
                                                                                 aTAAGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACt      >  1:4652218/1‑93 (MQ=39)
                                                                                   aaGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCa    >  1:4662936/1‑93 (MQ=39)
                                                                                    aGGCGT‑‑TCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCat   <  2:1992872/93‑1 (MQ=39)
                                                                                               ||                                                                                   
GCGCTGGCGGCGGCAATGGGCAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGT‑‑TCA‑‑CCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATATCTCCGGCATTTTTATTCCACAGCCAAACTCATA  >  NC_000913/3537312‑3537487

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: