Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,055,004 T→G 5.4% L15R (CTG→CGG)  oprB‑III → carbohydrate‑selective porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,055,0040TG 5.4% 850.5 / 21.7 263L15R (CTG→CGG) oprB‑IIIcarbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base T (176/71);  new base G (9/7);  total (185/78)
Fisher's exact test for biased strand distribution p-value = 2.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 295 total aligned reads displayed.

CCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGCCCACGGTG  >  NC_002947/4054855‑4055153
                                                                                                                                                     |                                                                                                                                                     
  ggaaaaCCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCGTGCCCGCGGATGCCCCTCAAAGAAAACCCGAGATTCTGTAATGCCATCCGCCAGTCGCATAACCCCCAACCGGCTGctggccctgg                                                                                                                                                     >  1:5559922/3‑150 (MQ=38)
   tAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggccctggc                                                                                                                                                    <  1:4268791/150‑1 (MQ=255)
       aCCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCGTGCCCGCGGATGCCCCGCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCCCCCGGCTGCTGGCCCGGGCCAGc                                                                                                                                                >  2:4776273/1‑150 (MQ=38)
        ccGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCAGGGCCAGca                                                                                                                                               <  1:6466426/150‑1 (MQ=255)
        ccGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCGACAAGGCCCCTGGAGGAGCGGGCTTGCCCGCGAAAGCCCCTCAAAGACCACCCGAGATTCTGTAATGCCCAACGCCATTCGCATCACCCCCACCCGGCGGCTGGCCCGGGCCCGcc                                                                                                                                               >  2:4707787/1‑149 (MQ=37)
         cgcgTCGCCTGCTTCGCGGGCAAGCCCGCACCTACACGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGcac                                                                                                                                              <  1:3562749/150‑1 (MQ=255)
                ctggCTTCGCGGGCACGCCCGCGCAAACACGGCCCGTTTAGGCGGGGGGTTGCCCGCGAATGCCCCTCAAACAAACCCCGATATTATGTAAGGCCAGCCGAAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGcc                                                                                                                                       <  1:113702/147‑1 (MQ=255)
                ccagCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGcc                                                                                                                                       <  2:2615028/147‑1 (MQ=255)
                     ttCGCGGGCAAGCCCGCTCCTACAATGCCCCTGTAGGAGCGGGCTTGGCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATTCCAACCGCAATTCGCATAAACCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGcccattc                                                                                                                                  >  2:4591785/1‑146 (MQ=38)
                         cgGGAACGCCCGCCCCTCCCCGGCCCCTTTTGGAGCGGTCTGGCCCGCGAATGCCCCCCAAAGACAACCCGAGATTCTTTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGc                                                                                                                              <  1:1165302/150‑1 (MQ=255)
                          gggCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcg                                                                                                                             >  1:4832344/1‑150 (MQ=255)
                           ggCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCTGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCCGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcga                                                                                                                            >  1:979677/1‑150 (MQ=255)
                           ggCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcga                                                                                                                            <  2:1639131/150‑1 (MQ=255)
                            gCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTCATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcgac                                                                                                                           >  1:5061206/1‑150 (MQ=255)
                              aaGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcgacga                                                                                                                         >  2:4007917/1‑150 (MQ=255)
                                  ccGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtctc                                                                                                                     >  1:3784664/1‑150 (MQ=255)
                                  ccGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtctc                                                                                                                     <  2:2249757/150‑1 (MQ=255)
                                   cGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCa                                                                                                                    >  2:5638227/1‑150 (MQ=255)
                                    gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCACTTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCAt                                                                                                                   >  1:492863/1‑150 (MQ=255)
                                    gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCAt                                                                                                                   >  2:3845582/1‑150 (MQ=255)
                                    gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCAt                                                                                                                   <  2:4123227/150‑1 (MQ=255)
                                    gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCAt                                                                                                                   >  2:769213/1‑150 (MQ=255)
                                    gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCAt                                                                                                                   >  2:494467/1‑150 (MQ=255)
                                        cTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAcc                                                                                                               >  2:6560492/1‑150 (MQ=255)
                                         tACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCCGCACCACTGCCCTTGCCAGCGACGATCTCATGAccc                                                                                                              >  2:4280197/1‑150 (MQ=255)
                                         tACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAccc                                                                                                              >  1:6344966/1‑150 (MQ=255)
                                         tACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAccc                                                                                                              >  2:3481783/1‑150 (MQ=255)
                                          aCAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCg                                                                                                             >  1:6024142/1‑150 (MQ=255)
                                            atgcccccTGTTGGCGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgca                                                                                                           <  1:806054/146‑1 (MQ=255)
                                            aaGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgca                                                                                                           >  2:5522993/1‑150 (MQ=255)
                                            aaGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgca                                                                                                           >  1:4983175/1‑150 (MQ=255)
                                             aGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAgcagcagg                                                                                                                                          <  1:3166116/118‑3 (MQ=255)
                                               gCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgcagca                                                                                                        >  2:5355641/1‑150 (MQ=255)
                                               gCCCCTGTAGGAGCGGGCTAGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCAACACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgcagca                                                                                                        >  2:4378876/1‑150 (MQ=255)
                                                 cccTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAACGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCAcc                                                                                                      <  2:972209/150‑1 (MQ=255)
                                                 cccTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCAcc                                                                                                      >  1:1419797/1‑150 (MQ=255)
                                                    tGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTg                                                                                                   >  2:2400651/1‑150 (MQ=255)
                                                        ggAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCg                                                                                               <  2:4446537/150‑1 (MQ=255)
                                                           gCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTg                                                                                            >  1:3429964/1‑150 (MQ=255)
                                                                 ttGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTggg                                                                                      >  1:1149900/1‑150 (MQ=255)
                                                                     ccGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGggcgg                                                                                  >  2:6367177/1‑150 (MQ=255)
                                                                     ccGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGggcgg                                                                                  >  2:1115269/1‑150 (MQ=255)
                                                                          caTGCCCCTCAAAGACAACCCGAGAGTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCGGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTgc                                                                             <  2:5110607/149‑1 (MQ=255)
                                                                             gCCCCTCAAAGCCAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGcc                                                                          <  1:2139212/150‑1 (MQ=255)
                                                                               cccTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCAt                                                                        >  1:5348057/1‑150 (MQ=255)
                                                                               cccTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCAt                                                                        >  2:2960634/1‑150 (MQ=255)
                                                                               cccTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCAt                                                                        >  1:5393384/1‑150 (MQ=255)
                                                                                  tCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCaa                                                                     >  1:2404694/1‑150 (MQ=255)
                                                                                  tCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCaa                                                                     >  2:3626744/1‑150 (MQ=255)
                                                                                  tCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCaa                                                                     >  1:3494996/1‑150 (MQ=255)
                                                                                    aaaGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACt                                                                   >  1:6097255/1‑150 (MQ=255)
                                                                                      aGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCg                                                                 <  1:1655156/150‑1 (MQ=255)
                                                                                       gACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTTCGCCATCAACTCga                                                                >  2:3411415/1‑150 (MQ=255)
                                                                                        aCAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCgaa                                                               <  2:1474080/150‑1 (MQ=255)
                                                                                          aaCCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCgaaga                                                             >  2:1280531/1‑150 (MQ=255)
                                                                                          aaCCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCgaaga                                                             >  1:3257959/1‑150 (MQ=255)
                                                                                            cccGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGaagaag                                                           >  2:5759600/1‑150 (MQ=255)
                                                                                            cccGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGaagaag                                                           >  2:5592206/1‑150 (MQ=255)
                                                                                            cccGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGaagaag                                                           >  2:2684500/1‑150 (MQ=255)
                                                                                             ccGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAagaaga                                                          >  1:6406684/1‑150 (MQ=255)
                                                                                             ccGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAagaaga                                                          >  1:6286522/1‑150 (MQ=255)
                                                                                             ccGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAagaaga                                                          >  2:5504340/1‑150 (MQ=255)
                                                                                             ccGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAagaaga                                                          >  1:1577062/1‑150 (MQ=255)
                                                                                              cGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGAc                                                         >  1:6250067/1‑150 (MQ=255)
                                                                                              cGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGAc                                                         >  1:4557115/1‑150 (MQ=255)
                                                                                               gagaTTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACg                                                        >  1:1910676/1‑150 (MQ=255)
                                                                                                  aTTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCg                                                     >  1:5740266/1‑150 (MQ=255)
                                                                                                   ttCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGt                                                    >  1:5958951/1‑150 (MQ=255)
                                                                                                      tGTAATGCCATCCGCAATTCGCATCACCCCCCCCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTGGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCaa                                                 <  2:1762629/150‑1 (MQ=255)
                                                                                                        tAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGa                                               >  1:2781789/1‑150 (MQ=255)
                                                                                                        tAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGa                                               >  1:3874282/1‑150 (MQ=255)
                                                                                                           tGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCa                                            >  1:2682675/1‑150 (MQ=255)
                                                                                                           tGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCa                                            >  1:6581285/1‑150 (MQ=255)
                                                                                                             ccATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCAcc                                          <  2:3895053/150‑1 (MQ=255)
                                                                                                                   gCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGAc                                    >  1:1263244/1‑150 (MQ=255)
                                                                                                                   gCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGAc                                    >  1:6686466/1‑150 (MQ=255)
                                                                                                                   gCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGAc                                    >  1:3658567/1‑150 (MQ=255)
                                                                                                                     aaTTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTa                                  <  2:3215989/150‑1 (MQ=255)
                                                                                                                            aTCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGt                           <  1:2112033/150‑1 (MQ=255)
                                                                                                                            aTCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGt                           <  2:2440603/150‑1 (MQ=255)
                                                                                                                                 ccccACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAAc                      >  2:3126260/1‑150 (MQ=255)
                                                                                                                                   ccACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGAGTCAAGATCACCGGCGACTACAGCGGTGAAACCg                    >  2:6683235/1‑150 (MQ=255)
                                                                                                                                    cACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGc                   <  2:2324213/150‑1 (MQ=255)
                                                                                                                                    cACCCTGCGGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGc                   <  1:4776273/150‑1 (MQ=255)
                                                                                                                                    cACCCGGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGc                   <  1:1983953/150‑1 (MQ=255)
                                                                                                                                       ccTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTa                >  1:1743102/1‑150 (MQ=255)
                                                                                                                                       agtgctgctgGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTa                >  1:227273/3‑150 (MQ=255)
                                                                                                                                        gtggtgctgGCCAGGGCCAGCAGCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTAc               >  2:6234546/5‑150 (MQ=255)
                                                                                                                                        ctgctgctgGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTAc               >  1:1819081/1‑150 (MQ=255)
                                                                                                                                         tggtgctgGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACa              >  1:1356688/4‑150 (MQ=255)
                                                                                                                                         tggtgctgGCCAGGGCCAGCAGCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGCCTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGACACCGCCTACa              >  2:115024/4‑150 (MQ=255)
                                                                                                                                         tggtgctgGCCAGGGCCAGCAGCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTAGAAGAATACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACa              >  1:4651195/4‑150 (MQ=255)
                                                                                                                                            tgctgGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACg           >  2:3949278/1‑150 (MQ=255)
                                                                                                                                               tggccagggccAGCAGCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGccc        >  2:6128507/1‑150 (MQ=255)
                                                                                                                                               tggccagggccAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGccc        >  2:2061970/1‑150 (MQ=255)
                                                                                                                                                  ccctggccAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGc                   <  1:4278510/136‑1 (MQ=255)
                                                                                                                                                  ccctggccAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGCCCACg     >  2:6777675/1‑150 (MQ=255)
                                                                                                                                                   cctggccAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGCCCACgg    >  2:6736507/1‑150 (MQ=255)
                                                                                                                                                     tggccAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGCCCACGGTg  >  2:5872278/1‑150 (MQ=255)
                                                                                                                                                     tggccAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGCCCACGGTg  >  1:4094196/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGCCCACGGTG  >  NC_002947/4054855‑4055153

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