Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,248,207 A→G 10.0% intergenic (+31/+60) ibpA → / ← PP_1983 small heat shock protein IbpA/PAS/PAC/GAF sensor‑containing diguanylate cyclase/phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,248,2070AG10.0% 461.5 / 29.4 160intergenic (+31/+60)ibpA/PP_1983small heat shock protein IbpA/PAS/PAC/GAF sensor‑containing diguanylate cyclase/phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (71/73);  new base G (10/6);  total (81/79)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

Only 100 of 170 total aligned reads displayed.

AAGTGAAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGATGGCA  >  NC_002947/2248058‑2248356
                                                                                                                                                     |                                                                                                                                                     
aaGTGAAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAATATTGGAACCTCCa                                                                                                                                                       >  2:4284771/1‑150 (MQ=255)
aaGTGAAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCa                                                                                                                                                       >  1:4639072/1‑150 (MQ=255)
   tGAAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAATAGGGCACCTCCagag                                                                                                                                                    >  2:4164281/1‑150 (MQ=255)
   tGAAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCagag                                                                                                                                                    >  2:1416889/1‑150 (MQ=255)
   tGAAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCagag                                                                                                                                                    >  1:4450172/1‑150 (MQ=255)
    gAAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAgg                                                                                                                                                   <  1:1621104/150‑1 (MQ=255)
     aaGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGt                                                                                                                                                  <  1:5330774/150‑1 (MQ=255)
       ggCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGc                                                                                                                                                >  1:2250926/1‑150 (MQ=255)
       ggCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGc                                                                                                                                                >  2:1641391/1‑150 (MQ=255)
        gccgccGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGcc                                                                                                                                               <  2:3722448/150‑1 (MQ=255)
        gccgccGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGcc                                                                                                                                               <  2:1790314/150‑1 (MQ=255)
           gccgGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCtt                                                                                                                                            <  1:4683773/150‑1 (MQ=255)
           gccgGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCtt                                                                                                                                            <  2:35847/150‑1 (MQ=255)
               gCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTc                                                                                                                                        <  2:466977/150‑1 (MQ=255)
                 cTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTGGTTCGt                                                                                                                                      >  2:4594289/1‑150 (MQ=255)
                    gCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCt                                                                                                                                   >  2:1647072/1‑150 (MQ=255)
                     ccAATGGCCTGCTCCGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAAGGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCtt                                                                                                                                  <  2:3570892/150‑1 (MQ=255)
                         tGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCa                                                                                                                              <  1:5006491/150‑1 (MQ=255)
                            ccTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGc                                                                                                                           <  1:1555139/150‑1 (MQ=255)
                                cTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTc                                                                                                                                        <  1:556633/133‑1 (MQ=255)
                                cTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTc                                                                                                                                        >  2:556633/1‑133 (MQ=255)
                                cTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTc                                                                                                                       >  1:102454/1‑150 (MQ=255)
                                   aGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCcg                                                                                                                    >  1:3888702/1‑150 (MQ=255)
                                    gCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCcgc                                                                                                                   <  2:2728854/150‑1 (MQ=255)
                                     cATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGCCCt                                                                                                                                             <  2:1324397/123‑1 (MQ=255)
                                     cATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGCCCt                                                                                                                                             >  1:1324397/1‑123 (MQ=255)
                                      aTCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGt                                                                                                                 >  1:4338075/1‑150 (MQ=255)
                                        cGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAg                                                                                                               <  2:3211123/150‑1 (MQ=255)
                                        cGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAg                                                                                                               <  2:3888702/150‑1 (MQ=255)
                                           ccTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTg                                                                                                            >  1:229739/1‑150 (MQ=255)
                                             tgctgcGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCGATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTgcg                                                                                                          >  2:2668487/1‑150 (MQ=255)
                                                  cgcATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGcc                                                                                                     <  1:3137518/150‑1 (MQ=255)
                                                     aTCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTc                                                                                                  <  2:3909390/150‑1 (MQ=255)
                                                      tCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCa                                                                                                 <  2:2250926/150‑1 (MQ=255)
                                                         tGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGCCCtttt                                                                                                                                          <  1:4466969/106‑1 (MQ=255)
                                                         tGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGc                                                                                              <  2:2376646/150‑1 (MQ=255)
                                                           ccTGAAGAACCCAAGGCCAAGCGCATTCCGATCAAGGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGa                                                                                            <  1:4857802/150‑1 (MQ=255)
                                                            cTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGAt                                                                                           <  2:4954553/150‑1 (MQ=255)
                                                             tGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATc                                                                                          <  1:890077/150‑1 (MQ=255)
                                                              gaagaaGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCa                                                                                         <  2:4450172/150‑1 (MQ=255)
                                                              gaagaaGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCa                                                                                         >  2:1246353/1‑150 (MQ=255)
                                                                agaagCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGt                                                                                       >  1:1048746/1‑150 (MQ=255)
                                                                 gaagCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTa                                                                                      <  2:3432429/150‑1 (MQ=255)
                                                                     ccaaggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCt                                                                                  <  2:386082/150‑1 (MQ=255)
                                                                      caaggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTg                                                                                 >  2:518453/1‑150 (MQ=255)
                                                                           ccaagCGCATTCCGATCAATGGTGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGcc                                                                                                                                               <  2:3915036/83‑1 (MQ=39)
                                                                           ccaagCGCATTCCGATCAATGGTGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGcc                                                                                                                                               >  1:3915036/1‑83 (MQ=39)
                                                                           ccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAATGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGgc                                                                            >  1:1231603/1‑150 (MQ=255)
                                                                             aagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGgcgc                                                                          >  2:2629489/1‑150 (MQ=255)
                                                                             aagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGgcgc                                                                          >  2:4730865/1‑150 (MQ=255)
                                                                               gcgcATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCt                                                                        >  1:1976790/1‑150 (MQ=255)
                                                                                cgcATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTc                                                                       <  2:916075/150‑1 (MQ=255)
                                                                                cgcATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTc                                                                       <  2:637816/150‑1 (MQ=255)
                                                                                    ttCAGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGAGGCCGACATGCGATCAGTAATCTGACGGCGCCTCCAGc                                                                   >  1:3665548/1‑150 (MQ=255)
                                                                                      ccGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCaa                                                                 <  2:2818489/150‑1 (MQ=255)
                                                                                        gATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACa                                                               >  2:1951444/1‑150 (MQ=255)
                                                                                                gCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGa                                                       <  2:4867915/150‑1 (MQ=255)
                                                                                                  gAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGGAGGTGCGCGGCCg                                                                                                    >  2:1622371/1‑103 (MQ=255)
                                                                                                  gAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGCAGGTGCGCGGCCg                                                                                                    <  1:1622371/103‑1 (MQ=255)
                                                                                                   aaaaaCCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGcc                                                         >  1:3544379/1‑145 (MQ=255)
                                                                                                   aaaaaCCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGcc                                                         <  2:3544379/145‑1 (MQ=255)
                                                                                                     aaaCCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTc                                                  >  1:4015773/1‑150 (MQ=255)
                                                                                                        ccAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcgg                                               <  1:2629489/150‑1 (MQ=255)
                                                                                                         cAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGTCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcggt                                              >  2:3063521/1‑150 (MQ=255)
                                                                                                         cAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCa                                                                                         >  1:2922209/1‑107 (MQ=255)
                                                                                                         cAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCa                                                                                         >  1:2922238/1‑107 (MQ=255)
                                                                                                         cAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCa                                                                                         <  2:2922238/107‑1 (MQ=255)
                                                                                                          aGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcggta                                             >  1:5186375/1‑150 (MQ=255)
                                                                                                          aGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcggta                                             >  1:3338365/1‑150 (MQ=255)
                                                                                                             gcTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcggtaacg                                          <  1:493346/150‑1 (MQ=255)
                                                                                                              cTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcggtaacgg                                         <  2:3396485/150‑1 (MQ=255)
                                                                                                               tGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcggtaacggt                                        >  1:2174071/1‑150 (MQ=255)
                                                                                                                    tGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATgg                                   >  1:650805/1‑150 (MQ=255)
                                                                                                                     gATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGc                                  >  1:4168566/1‑150 (MQ=255)
                                                                                                                     gATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGc                                  <  2:1976790/150‑1 (MQ=255)
                                                                                                                             tGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCtgtg                          <  2:102454/150‑1 (MQ=255)
                                                                                                                              gCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGc                         >  1:1205155/1‑150 (MQ=255)
                                                                                                                               cAACGAAGAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCa                        >  2:5028396/1‑150 (MQ=255)
                                                                                                                               cAACGAAAAAAAGGGCCCCTCCAGAGGTGCCCTTTTGCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCGGTGCa                        <  2:4863788/150‑1 (MQ=255)
                                                                                                                               cAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTACTGGCCGGCTGTGCa                        >  1:5165479/1‑150 (MQ=255)
                                                                                                                               cAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCa                        <  1:1647533/150‑1 (MQ=255)
                                                                                                                               cAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCa                        >  2:4558255/1‑150 (MQ=255)
                                                                                                                                aaCGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAg                       >  2:2841235/1‑150 (MQ=255)
                                                                                                                                aaCGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAg                       <  1:1641391/150‑1 (MQ=255)
                                                                                                                                 aCGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcgg                                               >  2:5007220/1‑125 (MQ=255)
                                                                                                                                 aCGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGa                      >  2:1836115/1‑150 (MQ=255)
                                                                                                                                   gAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATa                    >  1:628478/1‑150 (MQ=255)
                                                                                                                                    aaaaaaaGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATaa                   >  2:345143/1‑150 (MQ=255)
                                                                                                                                       aaaaGGGCACCTCTGGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAAccc                >  1:3390/1‑150 (MQ=255)
                                                                                                                                         aaGGGCACCTCTGGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTg              >  1:3866462/1‑150 (MQ=255)
                                                                                                                                          aGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTgg             >  2:2287550/1‑150 (MQ=255)
                                                                                                                                          aGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTgg             <  1:4284771/150‑1 (MQ=255)
                                                                                                                                            ggCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTa           >  2:2450108/1‑150 (MQ=255)
                                                                                                                                               aCCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAg        >  2:5319060/1‑150 (MQ=255)
                                                                                                                                                 cTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTCACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGAt      >  1:3782346/1‑150 (MQ=255)
                                                                                                                                                 cTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGAt      >  2:3046728/1‑150 (MQ=255)
                                                                                                                                                  tCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGATg     >  1:3334634/1‑150 (MQ=255)
                                                                                                                                                  tCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGATg     >  1:3280525/1‑150 (MQ=255)
                                                                                                                                                   ccAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGATgg    >  1:1878810/1‑150 (MQ=255)
                                                                                                                                                     agagGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGATGGCa  >  2:1289184/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGTGAAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGATGGCA  >  NC_002947/2248058‑2248356

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