Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,981,205 C→G 5.1% intergenic (+101/‑199) rpsA → / → ihfB 30S ribosomal protein S1/integration host factor subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,981,2050CG 5.1% 562.6 / 13.6 168intergenic (+101/‑199)rpsA/ihfB30S ribosomal protein S1/integration host factor subunit beta
Reads supporting (aligned to +/- strand):  ref base C (92/67);  new base G (2/7);  total (94/74)
Fisher's exact test for biased strand distribution p-value = 4.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 168 total aligned reads displayed.

ACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGTTCAAATCCTTCCCGGCGT  >  NC_002947/1981057‑1981351
                                                                                                                                                      |                                                                                                                                                  
accaccATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCa                                                                                                                                                   >  2:3563385/1‑150 (MQ=255)
     caTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTg                                                                                                                                              <  2:4137652/150‑1 (MQ=255)
     caTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTg                                                                                                                                              <  1:4707250/150‑1 (MQ=255)
          ctgcgctgcTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGaa                                                                                                                                         <  2:5196518/150‑1 (MQ=255)
           tgcgctgcTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCCTGAAGCTGGCgtgt                                                                                                                                        <  2:4529274/150‑4 (MQ=255)
            gcgctgcTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGc                                                                                                                                       <  2:4581085/150‑1 (MQ=255)
              gctgctGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTTTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGccc                                                                                                                                     <  1:151418/150‑1 (MQ=255)
              gctgctGCGCGAAGCAATGGCCAACCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCCTGAAGCTGGCgtgttcc                                                                                                                                     <  1:872058/150‑7 (MQ=255)
               ctgctgCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCCTGAAGCTGGCgtgttcca                                                                                                                                    <  1:4085336/150‑8 (MQ=255)
                 gctgcGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGtg                                                                                                                                  >  2:1523013/1‑149 (MQ=255)
                  ctgcGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGt                                                                                                                                   >  1:1445491/1‑150 (MQ=255)
                  ctgcGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTg                                                                                                                                 <  2:392626/150‑1 (MQ=255)
                   tgcGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGt                                                                                                                                <  1:1195802/150‑1 (MQ=255)
                     cgcgAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAg                                                                                                                              <  2:580550/150‑1 (MQ=255)
                      gcgAAGCAATGGCCAAGCAGACCTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTGTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGGTGTAgg                                                                                                                             >  2:3015552/1‑150 (MQ=255)
                        gAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAg                                                                                                                           >  1:1525489/1‑150 (MQ=255)
                        gAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAg                                                                                                                           >  1:2057133/1‑150 (MQ=255)
                                 gCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAc                                                                                                                  >  1:4629924/1‑150 (MQ=255)
                                  ccAAGCAGAACTGAGTTCTGTTTGATCGGTAAAATGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAcc                                                                                                                 >  1:5131609/1‑150 (MQ=255)
                                   cAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAccc                                                                                                                <  1:1774915/150‑1 (MQ=255)
                                    aaGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCcg                                                                                                               <  1:3291809/150‑1 (MQ=255)
                                     aaaaGAACTAAGTTCTGATTGATCGGAAAAAATGTAGACCCTAGGGTCGCCC‑‑ATTTTTTATATACCCGAAAAGGTGGTGCGTCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTATTAGGAGCGGGTTCACCcgc                                                                                                              >  2:4801915/4‑150 (MQ=255)
                                     aGCAGACCTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCcgc                                                                                                              <  1:4112803/150‑1 (MQ=255)
                                      gCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCcgc                                                                                                              >  1:2284195/1‑150 (MQ=255)
                                        aGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGaa                                                                                                           >  1:1052612/1‑150 (MQ=255)
                                          aaCTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGc                                                                                                         >  1:2954/1‑150 (MQ=255)
                                            cTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAg                                                                                                       <  1:4523496/150‑1 (MQ=255)
                                                 ttCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGcc                                                                                                  >  2:2710191/1‑150 (MQ=255)
                                                    tGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACt                                                                                               >  2:2164513/1‑150 (MQ=255)
                                                    tGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACt                                                                                               >  2:3997994/1‑150 (MQ=255)
                                                        tGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGcc                                                                                           >  2:3587912/1‑150 (MQ=255)
                                                             ggTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCt                                                                                      >  2:5394148/1‑150 (MQ=255)
                                                             ggTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCt                                                                                      >  2:1857940/1‑150 (MQ=255)
                                                               tAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTaa                                                                                    <  1:2206808/150‑1 (MQ=255)
                                                                 aaaaGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGc                                                                                  <  2:1558872/150‑1 (MQ=255)
                                                                  aaaGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCg                                                                                 >  1:3339214/1‑150 (MQ=255)
                                                                  aaaGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCg                                                                                 >  1:1607773/1‑150 (MQ=255)
                                                                   aaGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCgg                                                                                >  2:2336841/1‑150 (MQ=255)
                                                                   aaGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCgg                                                                                >  2:2618156/1‑150 (MQ=255)
                                                                   aaGGGCGACCCGAGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCgg                                                                                >  1:5434033/1‑150 (MQ=255)
                                                                    aGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGt                                                                               <  2:1967925/150‑1 (MQ=255)
                                                                    aGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGt                                                                               >  2:3546870/1‑150 (MQ=255)
                                                                    aGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCCTGAAGCTGGCg                                                                                                                                           >  2:5174606/1‑90 (MQ=255)
                                                                    aGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCCTGAAGCTGGCg                                                                                                                                           <  1:5174606/90‑1 (MQ=255)
                                                                     gggCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAc                                                                                                                  <  1:3257409/114‑1 (MQ=255)
                                                                     gggCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAc                                                                                                                  >  2:3257409/1‑114 (MQ=255)
                                                                     gggCGACCCGAGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTa                                                                              >  1:3059883/1‑150 (MQ=255)
                                                                     gggCGACCCGAGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTa                                                                              >  2:2927198/1‑150 (MQ=255)
                                                                      ggCGACCCGAGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAg                                                                             >  2:3342788/1‑150 (MQ=255)
                                                                           cccTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCAt                                                                        >  2:491088/1‑150 (MQ=255)
                                                                               cGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCtt                                                                    <  1:3997994/150‑1 (MQ=255)
                                                                               cGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCtt                                                                    <  1:3546870/150‑1 (MQ=255)
                                                                               cGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCtt                                                                    <  1:3547238/150‑1 (MQ=255)
                                                                                  gTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGaa                                                                 <  1:1613085/150‑1 (MQ=255)
                                                                                  gTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGaa                                                                 <  2:4839294/150‑1 (MQ=255)
                                                                                  gTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGaa                                                                 <  1:4368324/150‑1 (MQ=255)
                                                                                  gTCGCCA‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGaa                                                                 <  1:2927198/150‑1 (MQ=255)
                                                                                   tagCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAAt                                                                <  1:3797260/148‑1 (MQ=255)
                                                                                     gCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTa                                                              <  2:888128/150‑1 (MQ=255)
                                                                                     gCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGGTGTAGGAGCGGGTTCACCCGCGAAGCCGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTa                                                              >  1:3456555/1‑150 (MQ=255)
                                                                                      ccc‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTAc                                                             <  2:3179217/150‑1 (MQ=255)
                                                                                      ccc‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTAc                                                             >  1:5448067/1‑150 (MQ=255)
                                                                                       cc‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACt                                                            >  1:5297768/1‑150 (MQ=255)
                                                                                       cc‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACt                                                            <  1:780886/150‑1 (MQ=255)
                                                                                           ttttttttGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAc                                                                                                                  <  2:648801/94‑1 (MQ=255)
                                                                                           ttttttttGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGACCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTc                                                          >  1:5051864/1‑150 (MQ=255)
                                                                                              tttttGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAAt                                                       <  1:891063/150‑1 (MQ=255)
                                                                                              tttttGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAAt                                                       <  2:3992975/150‑1 (MQ=255)
                                                                                                tttGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCt                                                     <  2:2752153/150‑1 (MQ=255)
                                                                                                 ttGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCtt                                                    >  1:4468355/1‑150 (MQ=255)
                                                                                                     tCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAAt                                                >  1:3951928/1‑150 (MQ=255)
                                                                                                     tCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAAt                                                >  1:1883835/1‑150 (MQ=255)
                                                                                                       cgcgACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCCATCTTGAAttt                                              >  2:5284338/1‑150 (MQ=255)
                                                                                                          gACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAAtttttt                                           <  1:484338/150‑1 (MQ=255)
                                                                                                          gACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAAtttttt                                           <  2:3618043/150‑1 (MQ=255)
                                                                                                          gACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAAtttttt                                           >  1:2538731/1‑150 (MQ=255)
                                                                                                            ccAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTtta                                         >  2:3473828/1‑150 (MQ=255)
                                                                                                             cAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCcg                                                                                                               >  2:2417479/1‑79 (MQ=255)
                                                                                                             cAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCcg                                                                                                               <  1:2417479/79‑1 (MQ=255)
                                                                                                              aGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAAttttttt                                          >  2:3725537/1‑147 (MQ=255)
                                                                                                              aGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAAttttttt                                          <  1:3725537/147‑1 (MQ=255)
                                                                                                                tttgtgCGGCGGATGTCACTGCTGTAACACGCCAGCTTGAGGCTGGAGAATACATTTTTATGATAGTGTTCTAAGGCGAATCAGGGGACACTGGCATATATATGAGGTTGTCCATTATTGAATTAGTTCTATCTTGAATTTTTTTATTaa                                     <  2:4378353/147‑1 (MQ=255)
                                                                                                                gtggtgCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTaa                                     <  1:2710191/150‑1 (MQ=255)
                                                                                                                 tggtgCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTaat                                    >  2:3760310/1‑149 (MQ=255)
                                                                                                                       ggcggATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGa                              <  1:5394148/150‑1 (MQ=255)
                                                                                                                         cggATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATc                            >  1:4850902/1‑150 (MQ=255)
                                                                                                                                 aCTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTCGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGt                    >  2:366849/1‑150 (MQ=255)
                                                                                                                                  ctgctgGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTATACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGtt                   >  2:2460276/1‑150 (MQ=255)
                                                                                                                                      tgGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTATGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTATGTCAAGATCCACCCTGTTCaat               >  2:2260013/1‑149 (MQ=255)
                                                                                                                                      tgGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAATATCCACCCTGTTCaaa               >  2:1407432/1‑150 (MQ=255)
                                                                                                                                      tgGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGTTCaaa               >  1:3630250/1‑150 (MQ=255)
                                                                                                                                      tgGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGTTCaaa               >  2:76552/1‑150 (MQ=255)
                                                                                                                                       ggAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTATTAAGTCAAGATCCACCCTGTTCAAATc             <  1:4012636/150‑1 (MQ=255)
                                                                                                                                         aaCACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGTTCAAATcc            >  2:4956118/1‑150 (MQ=255)
                                                                                                                                          acacGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGTTCAAATCCt           <  2:2284195/150‑1 (MQ=255)
                                                                                                                                           ttcGCCAGCCTGAAGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGTTCAAATCCtt          >  1:1055876/3‑150 (MQ=255)
                                                                                                                                                cAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGTTCAAATCCTTCCCgg     >  1:1274741/1‑150 (MQ=255)
                                                                                                                                                cAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGTTCAAATCCTTCCCgg     <  2:4704466/150‑1 (MQ=255)
                                                                                                                                                    ttCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAAt                                                                >  1:1761800/1‑87 (MQ=255)
                                                                                                                                                    ttCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAAt                                                                <  2:1761800/87‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
ACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTTGAATTACTTCAATCTTGAATTTTTTTATTAAGTCAAGATCCACCCTGTTCAAATCCTTCCCGGCGT  >  NC_002947/1981057‑1981351

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: