Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,981,136 C→A 8.6% intergenic (+32/‑268) 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,1360CA 8.6% 502.0 / 24.2 158intergenic (+32/‑268)rpsA/ihfB30S ribosomal protein S1/integration host factor subunit beta
Reads supporting (aligned to +/- strand):  ref base C (78/66);  new base A (11/3);  total (89/69)
Fisher's exact test for biased strand distribution p-value = 9.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 166 total aligned reads displayed.

ATCGAAAGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTT  >  NC_002947/1980987‑1981285
                                                                                                                                                     |                                                                                                                                                       
aTCGAAAGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTc                                                                                                                                                         <  2:8114/150‑1 (MQ=255)
aTCGAAAGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTc                                                                                                                                                         <  1:4774560/150‑1 (MQ=255)
  cGAAAGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCgg                                                                                                                                                       <  1:2198161/150‑1 (MQ=255)
   gAAAGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCggg                                                                                                                                                      <  2:1137099/150‑1 (MQ=255)
    aaaGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGt                                                                                                                                                     >  1:5027860/1‑150 (MQ=255)
     aaGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTc                                                                                                                                                    <  2:1551357/150‑1 (MQ=255)
     aaGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTc                                                                                                                                                    <  1:334134/150‑1 (MQ=255)
        acgacgCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGcc                                                                                                                                                 <  1:399202/150‑1 (MQ=255)
        acgacgCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGcc                                                                                                                                                 <  1:399104/150‑1 (MQ=255)
        acgacgATGAAGAGCGAGAAGATATCAAGAGCCTGAAAAACGCTACGGAAGCTGCTGCCGAAACCAACATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGTCTACCCTCGCGTCGcc                                                                                                                                                 >  1:4378353/1‑150 (MQ=255)
          gacgCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑t                                                                                                                                             >  2:5102783/1‑150 (MQ=255)
             gCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑tttt                                                                                                                                          >  2:1195802/1‑150 (MQ=255)
             gCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCGAGGGTCGCCC‑‑tttt                                                                                                                                          <  2:923505/150‑1 (MQ=255)
                gAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGGTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTAGCCC‑‑ttttttt                                                                                                                                       >  2:3797260/1‑150 (MQ=255)
                 aaGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑tttttttt                                                                                                                                      <  1:2699061/150‑1 (MQ=255)
                  agagCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTg                                                                                                                                     >  2:2704256/1‑150 (MQ=255)
                  agagCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTg                                                                                                                                     >  1:1603517/1‑150 (MQ=255)
                   gagCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGc                                                                                                                                    >  2:4368324/1‑150 (MQ=255)
                        cgAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCcgc                                                                                                                               >  1:3340425/1‑150 (MQ=255)
                         gAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCcgcg                                                                                                                              >  1:3618043/1‑150 (MQ=255)
                            gCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGAcc                                                                                                                           >  1:4249461/1‑150 (MQ=255)
                               aTCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAgg                                                                                                                        >  2:4012636/1‑150 (MQ=255)
                                        cTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTgcggc                                                                                                               <  1:1277825/150‑1 (MQ=255)
                                        cTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTgcggc                                                                                                               <  2:4854083/150‑1 (MQ=255)
                                        cTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTgcggc                                                                                                               <  1:4376288/150‑1 (MQ=255)
                                           aaaaaCGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGa                                                                                                            <  2:1623112/150‑1 (MQ=255)
                                           aaaaaCGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGa                                                                                                            <  1:4735434/150‑1 (MQ=255)
                                                   tCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCActg                                                                                                    <  2:4388844/150‑1 (MQ=255)
                                                   tCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCActg                                                                                                    >  1:4559939/1‑150 (MQ=255)
                                                    ccGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCActgc                                                                                                   >  1:888128/1‑150 (MQ=255)
                                                     cGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCActgct                                                                                                  >  1:1967925/1‑150 (MQ=255)
                                                      ggAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCActgctg                                                                                                 <  2:3955384/150‑1 (MQ=255)
                                                            tgctgcCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACGCTCGGGTCGCCC‑‑TTTTTTTTGCTCCTCGACCAGGTGGTGCGGCGGTTGGCACTGCTGGAacac                                                                                           <  2:202185/150‑1 (MQ=255)
                                                              ctgcCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGc                                                                                         <  2:515970/150‑1 (MQ=255)
                                                              ctgcCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGc                                                                                         >  2:2054658/1‑150 (MQ=255)
                                                               tgcCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGcc                                                                                        <  2:3270059/150‑1 (MQ=255)
                                                                       ccaccaTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCGTCAg                                                                                >  1:4704466/1‑150 (MQ=255)
                                                                        caccaTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCGAGGGTCGcc                                                                                                                                                 <  1:4558058/86‑1 (MQ=255)
                                                                           caTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTg                                                                            <  2:4137652/150‑1 (MQ=255)
                                                                           caTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTCTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCCTGAAGCTg                                                                            <  1:544534/150‑1 (MQ=255)
                                                                            aTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTgg                                                                           >  1:1233057/1‑150 (MQ=255)
                                                                             tGGCTGCGCTGCTGCTCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGTTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCCTGAAGCTGGc                                                                          <  1:4096833/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)
                                                                                    gctgctGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGccc                                                                   <  2:4801938/150‑1 (MQ=255)
                                                                                       gctgcGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGtt                                                                >  1:3724811/1‑150 (MQ=255)
                                                                                       gctgcGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGtt                                                                <  1:3322089/150‑1 (MQ=255)
                                                                                       gctgcGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGtg                                                                >  2:1523013/1‑149 (MQ=255)
                                                                                       gctgcGCGAAGCAATGGCCAAGCAGAACTGAGGTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGtt                                                                <  1:4056358/150‑1 (MQ=255)
                                                                                        ctgcGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCCTTTTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGt                                                                 >  1:1445491/1‑150 (MQ=255)
                                                                                          gcgcGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTa                                                             <  1:5102783/150‑1 (MQ=255)
                                                                                            gcgAAGCAATGGCCAAGCAGACCTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTGTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGGTGTAgg                                                           >  2:3015552/1‑150 (MQ=255)
                                                                                            gcgAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAgg                                                           >  2:873389/1‑150 (MQ=255)
                                                                                              gAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAg                                                         >  1:1525489/1‑150 (MQ=255)
                                                                                              gAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAg                                                         >  1:2057133/1‑150 (MQ=255)
                                                                                                 gCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCgg                                                      <  1:448079/150‑1 (MQ=255)
                                                                                                  cAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAATGGCGACCCTCGGGTCGCCA‑‑TTTTTTTTGCTCAGAGACCAGGAGGTGCGGAGGATTGCACTGCTGGAACACGCAAGATTCAGGCTGGCGAAGCCCGTTGTAGGAGCggg                                                     >  2:3207891/1‑150 (MQ=255)
                                                                                                   aaTGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGt                                                    >  2:1766631/1‑150 (MQ=255)
                                                                                                       gCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAc                                                <  1:2336399/150‑1 (MQ=255)
                                                                                                        ccAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAcc                                               >  1:1065672/1‑150 (MQ=255)
                                                                                                         cAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAccc                                              <  1:1774915/150‑1 (MQ=255)
                                                                                                           aGCAGACCTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCcgc                                            <  1:4112803/150‑1 (MQ=255)
                                                                                                            gCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCcgcg                                           >  1:220979/1‑150 (MQ=255)
                                                                                                              aGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGaa                                         >  1:1052612/1‑150 (MQ=255)
                                                                                                                  cTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAg                                     >  1:3237831/1‑150 (MQ=255)
                                                                                                                   tGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAgg                                    >  2:3798978/1‑150 (MQ=255)
                                                                                                                          tGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACt                             >  2:2164513/1‑150 (MQ=255)
                                                                                                                          tGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACt                             >  2:3997994/1‑150 (MQ=255)
                                                                                                                             ttGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGc                          <  2:1603517/150‑1 (MQ=255)
                                                                                                                              tGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGcc                         >  2:3587912/1‑150 (MQ=255)
                                                                                                                                   ggTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAccc                                              >  1:3578379/1‑124 (MQ=255)
                                                                                                                                   ggTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAccc                                              <  2:3578379/124‑1 (MQ=255)
                                                                                                                                   ggTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCt                    >  2:5394148/1‑150 (MQ=255)
                                                                                                                                   ggTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCt                    >  2:1857940/1‑150 (MQ=255)
                                                                                                                                       aaaaGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGGTCCGCGACCAGGTGGTGGGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTTGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGc                <  2:4604968/150‑1 (MQ=255)
                                                                                                                                       aaaaGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGc                <  2:1558872/150‑1 (MQ=255)
                                                                                                                                        aaaGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCg               >  2:3645654/1‑150 (MQ=255)
                                                                                                                                        aaaGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCg               >  1:1607790/1‑150 (MQ=255)
                                                                                                                                        aaaGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCg               >  1:4811178/1‑150 (MQ=255)
                                                                                                                                         aaGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCgg              >  2:2336841/1‑150 (MQ=255)
                                                                                                                                         aaGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCgg              >  2:2617876/1‑150 (MQ=255)
                                                                                                                                         aaGGGCGACCCGAGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCgg              >  1:5434033/1‑150 (MQ=255)
                                                                                                                                         aaGGGCGACCCGAGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGGTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTAGCTAAGCgg              >  2:566074/1‑150 (MQ=255)
                                                                                                                                          aGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGt             <  2:1967925/150‑1 (MQ=255)
                                                                                                                                          aGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGt             >  2:3547238/1‑150 (MQ=255)
                                                                                                                                          aGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCCTGAAGCTGGCg                                                                         >  2:5174606/1‑90 (MQ=255)
                                                                                                                                          aGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCCTGAAGCTGGCg                                                                         <  1:5174606/90‑1 (MQ=255)
                                                                                                                                          aGGGCGACCCGAGGGTCGCCC‑‑TTTGTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGt             >  2:2411843/1‑150 (MQ=255)
                                                                                                                                           gggCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAc                                                <  1:3257409/114‑1 (MQ=255)
                                                                                                                                           gggCGACCCGAGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTa            >  2:2927198/1‑150 (MQ=255)
                                                                                                                                           gggCGACCCGAGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTa            >  1:3059883/1‑150 (MQ=255)
                                                                                                                                            ggCGACCCGAGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAg           >  2:3342729/1‑150 (MQ=255)
                                                                                                                                            ggCGACCCGAGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAg           >  2:3342788/1‑150 (MQ=255)
                                                                                                                                                 cccTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCAt      >  2:491002/1‑150 (MQ=255)
                                                                                                                                                    tCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCt   <  2:1292191/150‑1 (MQ=255)
                                                                                                                                                    tCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCt   >  2:701798/1‑150 (MQ=255)
                                                                                                                                                     cGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCtt  <  2:269718/150‑1 (MQ=255)
                                                                                                                                                     cGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCtt  <  1:3546870/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
ATCGAAAGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCTT  >  NC_002947/1980987‑1981285

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: