Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 827,265 C→A 6.6% A151A (GCC→GCA ycaC‑I → putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947827,2650CA 6.6% 1012.0 / 41.3 303A151A (GCC→GCAycaC‑Iputative hydrolase
Reads supporting (aligned to +/- strand):  ref base C (131/152);  new base A (9/11);  total (140/163)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 303 total aligned reads displayed.

CCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGC  >  NC_002947/827117‑827412
                                                                                                                                                    |                                                                                                                                                   
ccACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCt                                                                                                                                                    <  1:7506292/150‑1 (MQ=255)
   ccGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGc                                                                                                                                                 <  1:7060200/150‑1 (MQ=255)
     ggCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGa                                                                                                                                               <  1:9047429/150‑1 (MQ=255)
      gCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGAc                                                                                                                                              >  1:1069821/1‑150 (MQ=255)
       cAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACg                                                                                                                                             >  2:5399984/1‑150 (MQ=255)
        aagaagCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGc                                                                                                                                            >  1:2481102/1‑150 (MQ=255)
             gcagcTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGa                                                                                                                                       <  1:5985348/150‑1 (MQ=255)
             gcagcTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCACTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGa                                                                                                                                       <  2:8238540/150‑1 (MQ=255)
              cagcTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAg                                                                                                                                      <  2:6997487/150‑1 (MQ=255)
                gctgctGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGc                                                                                                                                    >  2:7796608/1‑150 (MQ=255)
                  tgctgATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCAACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTg                                                                                                                                  >  1:555666/1‑150 (MQ=255)
                   gctgATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGc                                                                                                                                 >  1:6932556/1‑150 (MQ=255)
                       aTTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCggg                                                                                                                             <  2:177493/150‑1 (MQ=255)
                         tGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGgc                                                                                                                           <  2:8953670/150‑1 (MQ=255)
                         tGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGgc                                                                                                                           <  2:3848959/150‑1 (MQ=255)
                            tGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCa                                                                                                                        >  1:4877366/1‑150 (MQ=255)
                            tGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCa                                                                                                                        >  2:509135/1‑150 (MQ=255)
                                 tgGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACtga                                                                                                                   <  1:828683/150‑1 (MQ=255)
                                 tgGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACtga                                                                                                                   <  2:1406723/150‑1 (MQ=255)
                                   gTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACtgatg                                                                                                                 >  1:352210/1‑150 (MQ=255)
                                    tCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACtgatga                                                                                                                >  2:6619211/1‑150 (MQ=255)
                                    tCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACtgatga                                                                                                                <  1:8565227/150‑1 (MQ=255)
                                    tCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACtgatga                                                                                                                >  2:2847359/1‑150 (MQ=255)
                                     cACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGAc                                                                                                               >  1:5464231/1‑150 (MQ=255)
                                        tGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGAccag                                                                                                            >  2:4346383/1‑148 (MQ=255)
                                        tGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTg                                                                                                            >  2:1062572/1‑150 (MQ=255)
                                           aGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGtt                                                                                                         <  1:4623547/150‑1 (MQ=255)
                                            gTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTc                                                                                                        <  2:393728/150‑1 (MQ=255)
                                             tATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCg                                                                                                       >  2:6233631/1‑150 (MQ=255)
                                               tGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGc                                                                                                     >  2:5899864/1‑150 (MQ=255)
                                                   tGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTgg                                                                                                 >  1:8045943/1‑150 (MQ=255)
                                                   tGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTgg                                                                                                 <  1:4501953/150‑1 (MQ=255)
                                                        ttCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGc                                                                                            >  1:580407/1‑150 (MQ=255)
                                                        ttCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGc                                                                                            <  2:645397/150‑1 (MQ=255)
                                                              gCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTg                                                                                      <  1:5178735/150‑1 (MQ=255)
                                                               cccTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGc                                                                                     >  2:7140844/1‑150 (MQ=255)
                                                                  tGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGc                                                                                                                                                 <  1:7650394/87‑1 (MQ=255)
                                                                  ggTCGGCATGGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCAcc                                                                                  <  2:7273680/149‑1 (MQ=255)
                                                                      ggCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGa                                                                              >  2:847848/1‑150 (MQ=255)
                                                                       gCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGAc                                                                             <  2:4747334/150‑1 (MQ=255)
                                                                        cATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACt                                                                            >  2:7470736/1‑150 (MQ=255)
                                                                         aTTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTg                                                                           <  1:6041256/150‑1 (MQ=255)
                                                                         aTTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTg                                                                           <  1:6633908/150‑1 (MQ=255)
                                                                          ttGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTgg                                                                          >  2:2344203/1‑150 (MQ=255)
                                                                              aggaggGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCa                                                                      <  1:7796599/150‑1 (MQ=255)
                                                                                gaggGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGCTGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAAc                                                                    <  2:7746912/150‑1 (MQ=255)
                                                                                gaggGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAAc                                                                    >  2:2818871/1‑150 (MQ=255)
                                                                                   ggCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGAc                                                                 <  2:8768241/150‑1 (MQ=255)
                                                                                    gtttCGATGTCTTTGTCGTCACAGATGCGTCGGGAACCTTCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGata                                                                >  1:8697904/3‑148 (MQ=2)
                                                                                      ttCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATc                                                              >  2:548276/1‑150 (MQ=255)
                                                                                         gAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGaa                                                           >  1:5044128/1‑150 (MQ=255)
                                                                                         gAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGaa                                                           <  1:3984872/150‑1 (MQ=255)
                                                                                          aGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAg                                                          >  2:8804892/1‑150 (MQ=255)
                                                                                           ggTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAgg                                                         >  1:3923459/1‑150 (MQ=255)
                                                                                            gtctttgtcgtcACAGATGCGTCGGGAACCTTCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAggt                                                        <  2:6220442/144‑2 (MQ=2)
                                                                                               tttgtcgtcACAGATGCGTCGGGAACCTTCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTg                                                     >  2:356605/4‑150 (MQ=2)
                                                                                               ttCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTg                                                     <  1:4365452/150‑1 (MQ=255)
                                                                                                ttgtcgtcACAGATGCGTCGGGAACCTTCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTgg                                                    <  1:5675025/148‑1 (MQ=2)
                                                                                                ttgtcgtcACAGATGCGTCGGGAACCTTCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTgg                                                    <  2:6730346/148‑1 (MQ=2)
                                                                                                tcgtcgTCACCAACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTgg                                                    <  2:555666/150‑1 (MQ=255)
                                                                                                 cgtcgtCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTggg                                                   <  1:2215980/150‑1 (MQ=255)
                                                                                                     gtcACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTAcc                                               <  1:1720886/150‑1 (MQ=255)
                                                                                                     gtcACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTAcc                                               <  1:2515085/150‑1 (MQ=255)
                                                                                                     gtcACAGATGCGTCGGGAACCTTCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTGGGTAcg                                               <  2:7579640/150‑2 (MQ=2)
                                                                                                        aCCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTg                                            >  1:4069318/1‑150 (MQ=255)
                                                                                                           gACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTtc                                         <  2:6868432/150‑1 (MQ=255)
                                                                                                              gCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTcc                                      >  2:5136298/1‑150 (MQ=255)
                                                                                                              gCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTcc                                      >  2:2215754/1‑150 (MQ=255)
                                                                                                                ctctGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATc                                                              <  1:6273047/124‑1 (MQ=255)
                                                                                                                ctctGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCaa                                    >  1:841288/1‑150 (MQ=255)
                                                                                                                  ctGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAAcc                                  <  1:2436816/150‑1 (MQ=255)
                                                                                                                  ctGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAAcc                                  <  1:6151346/150‑1 (MQ=255)
                                                                                                                    ggCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACtgat                                                                                                                  >  1:3772210/1‑68 (MQ=255)
                                                                                                                       atgTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCt                                             <  1:8434193/132‑1 (MQ=255)
                                                                                                                       aCGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATc                             >  2:8569147/1‑150 (MQ=255)
                                                                                                                          ttCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCg                          <  2:4900353/150‑1 (MQ=255)
                                                                                                                          ttCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCg                          <  2:3288749/150‑1 (MQ=255)
                                                                                                                          ttCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCg                          >  2:3097729/1‑150 (MQ=255)
                                                                                                                           tCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTGGGTACGCTGTTTGCAAACCACATTCCgg                         >  2:2655530/1‑150 (MQ=2)
                                                                                                                           tCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTGGGTACGCTGTTTGCAAACCACATTCCgg                         >  1:4679173/1‑150 (MQ=2)
                                                                                                                            cAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGa                        <  1:2614674/150‑1 (MQ=255)
                                                                                                                             aatgAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTGGGTACGCTGTTTGCAAACCACATTCCGGat                       <  1:8812378/147‑2 (MQ=2)
                                                                                                                             aaCGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGAc                       >  2:2452994/1‑150 (MQ=255)
                                                                                                                              aCGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACt                      >  2:4251018/1‑150 (MQ=255)
                                                                                                                               cGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTa                     <  1:6973837/150‑1 (MQ=255)
                                                                                                                                     gACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCaa               <  2:8926503/150‑1 (MQ=255)
                                                                                                                                     gACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGACCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCaa               >  1:2909263/1‑150 (MQ=255)
                                                                                                                                        ccGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCt            >  1:7146419/1‑150 (MQ=255)
                                                                                                                                          gTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCtga          >  1:132083/1‑150 (MQ=255)
                                                                                                                                           tGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATTGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCtgat         <  2:4575296/150‑1 (MQ=255)
                                                                                                                                           tGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCtgat         <  1:5221314/150‑1 (MQ=255)
                                                                                                                                             aCTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCtgatga       >  2:8368563/1‑150 (MQ=255)
                                                                                                                                             aCTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCtgatga       >  2:3758454/1‑150 (MQ=255)
                                                                                                                                              cTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGAc      <  1:6320219/150‑1 (MQ=255)
                                                                                                                                              cTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGAc      <  1:5881400/150‑1 (MQ=255)
                                                                                                                                               tCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGTCCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGTcc     >  2:5816032/1‑150 (MQ=255)
                                                                                                                                                cGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCa    <  1:205273/150‑1 (MQ=255)
                                                                                                                                                cGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCa    <  2:1157643/150‑1 (MQ=255)
                                                                                                                                                cGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCa    <  1:636514/150‑1 (MQ=255)
                                                                                                                                                  gCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGc  >  1:7961883/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCACCGGCAAGAAGCAGCTGCTGATTGCTGGTGTGGTCACTGAAGTATGCGTGGCCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGC  >  NC_002947/827117‑827412

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: