Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,486,594 G→A 17.0% R123Q (CGG→CAG) ‡ yhdZ → amino acid ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,486,5940GA17.0% 490.2 / 85.2 211R123Q (CGG→CAG) ‡yhdZamino acid ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (82/93);  new base A (21/15);  total (103/108)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 214 total aligned reads displayed.

GATGGCGTGGAACTGACCAACGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGT  >  NC_002947/1486446‑1486743
                                                                                                                                                    |                                                                                                                                                     
gATGGCGTGGAACTGACCAACGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACgg                                                                                                                                                      >  1:337634/1‑150 (MQ=255)
   ggCGTGGAACTGACCAACGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGgcgc                                                                                                                                                   >  2:1675735/1‑150 (MQ=255)
      gTGGAACTGACCAACGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAg                                                                                                                                                >  2:3786457/1‑150 (MQ=255)
         gAACTGACCAACGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGcc                                                                                                                                             <  1:2476173/150‑1 (MQ=255)
          aaCTGACCAACGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCg                                                                                                                                            >  1:1887338/1‑150 (MQ=255)
           aCTGACCAACGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGa                                                                                                                                           <  1:3797845/150‑1 (MQ=255)
                 cAACGACCTCACGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCCCTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAAt                                                                                                                                     >  2:3012476/1‑150 (MQ=255)
                 cAACGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCaa                                                                                                                                                 <  1:2050689/138‑1 (MQ=255)
                 cAACGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCaa                                                                                                                                                 >  2:2050689/1‑138 (MQ=255)
                 cAACGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGACAt                                                                                                                                     >  1:3936871/1‑150 (MQ=255)
                 cAACGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAAt                                                                                                                                     >  1:2810448/1‑150 (MQ=255)
                      aCCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCa                                                                                                                                <  2:653104/150‑1 (MQ=255)
                       ccTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCAt                                                                                                                               <  2:1470408/150‑1 (MQ=255)
                            aGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACt                                                                                                                          >  1:341076/1‑150 (MQ=255)
                             gCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTa                                                                                                                         >  1:3675211/1‑150 (MQ=255)
                                   cGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGa                                                                                                                   <  1:66550/150‑1 (MQ=255)
                                         cATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCgtgt                                                                                                             <  2:3588123/150‑1 (MQ=255)
                                           tCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAAt                                                                                                                                     <  1:1523960/124‑1 (MQ=255)
                                           tCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGc                                                                                                           <  2:13635/150‑1 (MQ=255)
                                           tCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGc                                                                                                           <  1:399316/150‑1 (MQ=255)
                                                  cgAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGACCTTGCCATGCACTACCTGGAGCGTGTGCGGATTcc                                                                                                    >  2:1576671/1‑150 (MQ=255)
                                                  cgAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTcc                                                                                                    <  1:632761/150‑1 (MQ=255)
                                                   gAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCg                                                                                                   <  2:5048061/150‑1 (MQ=255)
                                                     ggTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGa                                                                                                 <  2:1595694/150‑1 (MQ=255)
                                                     ggTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGa                                                                                                 <  2:858217/150‑1 (MQ=255)
                                                        gggCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCa                                                                                              >  2:1999828/1‑150 (MQ=255)
                                                        gggCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCa                                                                                              >  2:2896954/1‑150 (MQ=255)
                                                        gggCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCa                                                                                              <  1:446790/150‑1 (MQ=255)
                                                        gggCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCa                                                                                              <  1:4384404/150‑1 (MQ=255)
                                                         ggCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAg                                                                                             <  1:4577599/150‑1 (MQ=255)
                                                            aTGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGca                                                                                          <  1:812719/150‑1 (MQ=255)
                                                            aTGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGca                                                                                          >  1:3852245/1‑150 (MQ=255)
                                                                 gtTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACaa                                                                                     <  2:5084836/150‑1 (MQ=255)
                                                                   tcgaGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGt                                                                                   >  2:3991609/4‑150 (MQ=255)
                                                                    ccAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTa                                                                                  >  2:4918032/1‑150 (MQ=255)
                                                                     cAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTAc                                                                                 >  2:9679/1‑150 (MQ=255)
                                                                     cAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTAc                                                                                 >  1:4880784/1‑150 (MQ=255)
                                                                     cAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTAc                                                                                 >  1:3881487/1‑150 (MQ=255)
                                                                      aGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTAcc                                                                                >  2:2611644/1‑150 (MQ=255)
                                                                      aGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTAcc                                                                                >  1:576858/1‑150 (MQ=255)
                                                                        cACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCt                                                                              <  1:1320347/150‑1 (MQ=255)
                                                                         aCTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAGGTACCCTg                                                                             <  2:885736/150‑1 (MQ=255)
                                                                           ttCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGc                                                                           >  2:4106954/1‑150 (MQ=255)
                                                                             cAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCa                                                                         <  2:82694/150‑1 (MQ=255)
                                                                             cAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCa                                                                         <  2:4850423/150‑1 (MQ=255)
                                                                             cAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCa                                                                         >  2:4315207/1‑150 (MQ=255)
                                                                              aaCCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAgg                                                                                            <  1:49115/130‑1 (MQ=255)
                                                                                ccTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGtt                                                                      <  1:2607589/150‑1 (MQ=255)
                                                                                ccTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGtt                                                                      <  1:1516416/150‑1 (MQ=255)
                                                                                   gTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTc                                                                   >  1:2301246/1‑150 (MQ=255)
                                                                                    gtCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCt                                                                  <  2:3518678/149‑1 (MQ=255)
                                                                                       ccGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTggc                                                               >  1:3123228/1‑150 (MQ=255)
                                                                                       ccGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTggc                                                               >  2:791836/1‑150 (MQ=255)
                                                                                           aCCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGa                                                                                                 >  2:2012713/1‑112 (MQ=255)
                                                                                           aCCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGa                                                                                                 <  1:2012713/112‑1 (MQ=255)
                                                                                               gAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCcagca                                                       <  1:2764824/150‑1 (MQ=255)
                                                                                                      cTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCgt                                                <  2:2320031/150‑1 (MQ=255)
                                                                                                       tGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGc                                                                                                                 <  1:4450845/84‑1 (MQ=255)
                                                                                                       tGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGc                                                                                                                 >  2:4450845/1‑84 (MQ=255)
                                                                                                       tGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCgtg                                               >  2:4715230/1‑150 (MQ=255)
                                                                                                             aCTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCa                                         <  1:3041159/150‑1 (MQ=255)
                                                                                                               tGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAgg                                                                                                                                         >  2:3534885/1‑52 (MQ=255)
                                                                                                               tGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCAtt                                       >  1:4413739/1‑150 (MQ=255)
                                                                                                               tGCACCCTGGCACCGATTTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGAtt                                       >  2:4601252/1‑150 (MQ=255)
                                                                                                               tGCACCCTGGCACCGATCTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGAtt                                       >  2:1559855/1‑150 (MQ=255)
                                                                                                                 cACCCTGGCACCGATTTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGc                                     <  1:4711153/150‑1 (MQ=255)
                                                                                                                  aCCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGcc                                    >  1:1825382/1‑150 (MQ=255)
                                                                                                                  aCCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGcc                                    >  1:4202704/1‑150 (MQ=255)
                                                                                                                    ccTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCcg                                  >  1:4840464/1‑150 (MQ=255)
                                                                                                                    ccTGGCACCGATTTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCcg                                  >  2:141445/1‑150 (MQ=255)
                                                                                                                     cTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCcgc                                 >  2:2663623/1‑150 (MQ=255)
                                                                                                                      tGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCcgcg                                >  1:1584968/1‑150 (MQ=255)
                                                                                                                      tGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCcgcg                                >  1:3735782/1‑150 (MQ=255)
                                                                                                                          gcCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGtt                            >  1:3021277/1‑150 (MQ=255)
                                                                                                                          gcCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGtt                            >  2:4020437/1‑150 (MQ=255)
                                                                                                                          accGATTTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCcgtgcgct                            >  2:2530342/2‑148 (MQ=255)
                                                                                                                          accGATCTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCcgtgcgct                            >  1:1020636/2‑148 (MQ=255)
                                                                                                                            cGATTTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCcgtgcgcttt                          <  2:1116892/150‑5 (MQ=255)
                                                                                                                            cGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGGGTGCGGATTCCGGAGCAGGCACACACGTTCCCTGGCGAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTtgt                          <  2:3333480/150‑1 (MQ=255)
                                                                                                                             gATCTGGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCcgtgcgctttg                         <  2:475438/150‑6 (MQ=255)
                                                                                                                              atctgGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGc                        <  2:2869687/147‑1 (MQ=255)
                                                                                                                                 tgGGTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATg                     <  2:180405/150‑1 (MQ=255)
                                                                                                                                  gggTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGa                    >  2:1961696/1‑150 (MQ=255)
                                                                                                                                   ggTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGaa                   <  2:4678371/150‑1 (MQ=255)
                                                                                                                                   ggTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGaa                   >  1:4276450/1‑150 (MQ=255)
                                                                                                                                   ggTACGCAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGaa                   <  1:4753915/150‑1 (MQ=255)
                                                                                                                                    gtacgcAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGaac                  <  1:1210195/147‑2 (MQ=255)
                                                                                                                                    gtacgcAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGaac                  <  2:1008994/147‑2 (MQ=255)
                                                                                                                                     tGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGc                 >  1:1147073/1‑150 (MQ=255)
                                                                                                                                        gcAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAAcccga              >  2:4905989/1‑148 (MQ=255)
                                                                                                                                        gcAAGATGCCAAAGAAGAAGGCGGAAGAAATCGCCATGCATTATCTGGCGCGCGTGAAAATTCCGGAACAGGCCAACAAATATCCGGGCCAGCTTTCAGGCGGCCAGCAGCAGCGTGTGGCGATTGCCCGTGCGCTTTGCATGAAcccga              >  1:92478/1‑148 (MQ=255)
                                                                                                                                         cAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCaa             <  2:576858/150‑1 (MQ=255)
                                                                                                                                         cAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCaa             <  1:3877619/150‑1 (MQ=255)
                                                                                                                                          aaGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAg            >  2:455498/1‑150 (MQ=255)
                                                                                                                                             aTGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATc         >  1:4067279/1‑150 (MQ=255)
                                                                                                                                               gCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCAt       <  2:1864286/150‑1 (MQ=255)
                                                                                                                                               gCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCAt       <  2:4405083/150‑1 (MQ=255)
                                                                                                                                                  aCGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCt    >  1:77629/1‑150 (MQ=255)
                                                                                                                                                  aCGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCt    >  1:1260478/1‑150 (MQ=255)
                                                                                                                                                   cGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTg   <  2:3199818/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GATGGCGTGGAACTGACCAACGACCTCAAGCAGATCGAGGCCATCCGCCGCGAGGTGGGCATGGTGTTCCAGCACTTCAACCTGTTCCCGCACCTGAGCATCCTGGAGAACTGCACCCTGGCGCCGATGTGGGTGCGCAAGATGCCACGGCGCAAGGCCGAGGAAATTGCCATGCACTACCTGGAGCGTGTGCGGATTCCGGAGCAGGCACACAAGTACCCTGGCCAGTTGTCTGGCGGCCAGCAGCAGCGTGTGGCCATTGCCCGCGCGTTGTGCATGAAGCCCAAGATCATGCTGT  >  NC_002947/1486446‑1486743

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