Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,981,135 T→G 8.2% intergenic (+31/‑269) 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,1350TG 8.2% 461.9 / 27.3 159intergenic (+31/‑269)rpsA/ihfB30S ribosomal protein S1/integration host factor subunit beta
Reads supporting (aligned to +/- strand):  ref base T (79/67);  new base G (10/3);  total (89/70)
Fisher's exact test for biased strand distribution p-value = 1.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.39e-01

Only 100 of 161 total aligned reads displayed.

ATCGAAAGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCT  >  NC_002947/1980987‑1981284
                                                                                                                                                    |                                                                                                                                                       
aTCGAAAGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTc                                                                                                                                                        <  2:8114/150‑1 (MQ=255)
aTCGAAAGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTc                                                                                                                                                        <  1:4774560/150‑1 (MQ=255)
  cGAAAGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCgg                                                                                                                                                      <  1:3160553/150‑1 (MQ=255)
  cGAAAGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCgg                                                                                                                                                      <  1:2198161/150‑1 (MQ=255)
   gAAAGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCggg                                                                                                                                                     <  2:1137099/150‑1 (MQ=255)
     aaGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTc                                                                                                                                                   <  2:1551366/150‑1 (MQ=255)
       gacgacGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGc                                                                                                                                                 <  2:1247567/150‑1 (MQ=255)
       gacgacGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGc                                                                                                                                                 <  2:4458430/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)
            cgCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCGAGGGTCGCCC‑‑ttt                                                                                                                                          <  2:4257773/150‑1 (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)
                        cgAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCcgc                                                                                                                              >  1:3340425/1‑150 (MQ=255)
                         gAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCcgcg                                                                                                                             >  1:3618043/1‑150 (MQ=255)
                            gCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGAcc                                                                                                                          >  1:4249461/1‑150 (MQ=255)
                             cTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑ttttttt                                                                                                                                      <  2:1943260/137‑1 (MQ=255)
                               aTCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAgg                                                                                                                       >  2:4012636/1‑150 (MQ=255)
                                tCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAggt                                                                                                                      >  2:484338/1‑150 (MQ=255)
                                        cTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTgcggc                                                                                                              <  1:1277825/150‑1 (MQ=255)
                                        cTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTgcggc                                                                                                              <  1:4376185/150‑1 (MQ=255)
                                        cTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTgcggc                                                                                                              <  2:4854083/150‑1 (MQ=255)
                                         tGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTgcggcg                                                                                                             >  1:3687829/1‑150 (MQ=255)
                                         tGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTgcggcg                                                                                                             >  1:4839294/1‑150 (MQ=255)
                                           aaaaaCGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGa                                                                                                           <  2:1623112/150‑1 (MQ=255)
                                                  cTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCAct                                                                                                    <  2:2969936/150‑1 (MQ=255)
                                                   tCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCActg                                                                                                   <  2:4388844/150‑1 (MQ=255)
                                                   tCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCActg                                                                                                   >  1:4559939/1‑150 (MQ=255)
                                                     cGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCActgct                                                                                                 >  1:1967925/1‑150 (MQ=255)
                                                      ggAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCActgctg                                                                                                <  2:3955384/150‑1 (MQ=255)
                                                       gAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTgg                                                                                               >  2:3620398/1‑150 (MQ=255)
                                                       gAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTgg                                                                                               >  2:1613085/1‑150 (MQ=255)
                                                         aGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGaa                                                                                             <  1:4781498/150‑1 (MQ=255)
                                                            tgctgcCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACGCTCGGGTCGCCC‑‑TTTTTTTTGCTCCTCGACCAGGTGGTGCGGCGGTTGGCACTGCTGGAacac                                                                                          <  2:202185/150‑1 (MQ=255)
                                                              ctgcCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGc                                                                                        >  2:2054658/1‑150 (MQ=255)
                                                              ctgcCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGc                                                                                        <  2:515970/150‑1 (MQ=255)
                                                               tgcCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGcc                                                                                       <  2:3270164/150‑1 (MQ=255)
                                                                gcCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCa                                                                                      <  1:914437/150‑1 (MQ=255)
                                                                        caccaTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCGAGGGTCGcc                                                                                                                                                >  2:4558058/1‑86 (MQ=255)
                                                                           caTGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTg                                                                           <  1:4707250/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)
                                                                                  gcgctgcTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGc                                                                    <  2:4581085/150‑1 (MQ=255)
                                                                                    gctgctGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTTTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGccc                                                                  <  1:151418/150‑1 (MQ=255)
                                                                                    gctgctGCGCGAAGCAATGGCCAACCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCCTGAAGCTGGCgtgttcc                                                                  <  1:872058/150‑7 (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)
                                                                                                 gCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCgg                                                     <  1:448079/150‑1 (MQ=255)
                                                                                                   aaTGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGt                                                   >  2:1766631/1‑150 (MQ=255)
                                                                                                       gCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAc                                               >  1:4629924/1‑150 (MQ=255)
                                                                                                        ccAAGCAGAACTGAGTTCTGTTTGATCGGTAAAATGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAcc                                              >  1:5131609/1‑150 (MQ=255)
                                                                                                        ccAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAcc                                              >  1:1065672/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)
                                                                                                               gAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGGTGTAGGAGCGGGTTCACCCGCGAAg                                       >  2:5177492/1‑150 (MQ=255)
                                                                                                                aaCTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGc                                      >  1:2954/1‑150 (MQ=255)
                                                                                                                  cTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAg                                    <  1:4523496/150‑1 (MQ=255)
                                                                                                                  cTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAg                                    >  1:3237831/1‑150 (MQ=255)
                                                                                                                  cTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGAAGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAg                                    >  1:3801271/1‑150 (MQ=255)
                                                                                                                   tGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAgg                                   >  2:3798978/1‑150 (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)
                                                                                                                            tttGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTgg                          >  1:3377265/1‑150 (MQ=255)
                                                                                                                             tttaTCGGTAAAAAGGGCGACCCTCGGCTCCCCC‑‑ATTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGc                         <  1:5293308/147‑1 (MQ=255)
                                                                                                                                   ggTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAccc                                             <  2:3578379/124‑1 (MQ=255)
                                                                                                                                   ggTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCt                   >  2:5394148/1‑150 (MQ=255)
                                                                                                                                     tAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTaa                 <  1:2206808/150‑1 (MQ=255)
                                                                                                                                       aaaaGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGc               <  2:1558872/150‑1 (MQ=255)
                                                                                                                                        aaaGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCg              >  2:3645654/1‑150 (MQ=255)
                                                                                                                                        aaaGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCg              >  1:1607773/1‑150 (MQ=255)
                                                                                                                                        aaaGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCg              >  1:3339214/1‑150 (MQ=255)
                                                                                                                                        aaaGGGCGACCCGAGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCg              >  1:2019452/1‑150 (MQ=255)
                                                                                                                                         aaGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCgg             >  2:2617876/1‑150 (MQ=255)
                                                                                                                                         aaGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCgg             >  2:2336841/1‑150 (MQ=255)
                                                                                                                                         aaGGGCGACCCGAGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCgg             >  1:5434033/1‑150 (MQ=255)
                                                                                                                                         aaGGGCGACCCGAGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGGTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTAGCTAAGCgg             >  2:566074/1‑150 (MQ=255)
                                                                                                                                          aGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGt            >  2:3547238/1‑150 (MQ=255)
                                                                                                                                          aGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGt            >  2:3546870/1‑150 (MQ=255)
                                                                                                                                          aGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGt            <  2:1967925/150‑1 (MQ=255)
                                                                                                                                          aGGGCGACCCGAGGGTCGCCC‑‑TTTGTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGt            >  2:2411843/1‑150 (MQ=255)
                                                                                                                                           gggCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCAc                                               >  2:3257409/1‑114 (MQ=255)
                                                                                                                                            ggCGACCCGAGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAg          >  2:3342729/1‑150 (MQ=255)
                                                                                                                                                 cccTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCAt     >  2:491088/1‑150 (MQ=255)
                                                                                                                                                    tCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCt  >  2:701798/1‑150 (MQ=255)
                                                                                                                                                    tCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCt  <  2:1292191/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                       
ATCGAAAGACGACGCTGAAGAGCGCGAAGCTATCCAGAGCCTGAAAAACGCTCCGGAAGCTGCTGCCGACACCACCATGGCTGCGCTGCTGCGCGAAGCAATGGCCAAGCAGAACTGAGTTCTGTTTGATCGGTAAAAAGGGCGACCCTCGGGTCGCCC‑‑TTTTTTTTGCTCCGCGACCAGGTGGTGCGGCGGATGGCACTGCTGGAACACGCCAGCTTCAGGCTGGCGAAGCCCGTTGTAGGAGCGGGTTCACCCGCGAAGCAGGCGCCACTGGCCTATCTAAGCGGTAGACCATTCT  >  NC_002947/1980987‑1981284

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: