Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,367,710 A→G 13.3% G272G (GGT→GGC ftsH ← integral membrane ATP‑dependent zinc metallopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,367,7100AG13.3% 800.9 / 88.8 301G272G (GGT→GGCftsHintegral membrane ATP‑dependent zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base A (145/115);  new base G (14/26);  total (159/142)
Fisher's exact test for biased strand distribution p-value = 1.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

Only 100 of 304 total aligned reads displayed.

CCTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACC  >  NC_002947/5367562‑5367859
                                                                                                                                                    |                                                                                                                                                     
ccTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGAc                                                                                                                                                      >  1:3960281/1‑150 (MQ=255)
 cTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGAcc                                                                                                                                                     >  2:2115987/1‑150 (MQ=255)
        tCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGccc                                                                                                                                              <  1:1990567/150‑1 (MQ=255)
          gAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCAt                                                                                                                                            >  2:4116104/1‑150 (MQ=255)
           acgcgggcggGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGCTGGCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCTACTGTTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATg                                                                                                                                           <  1:5103347/148‑1 (MQ=255)
               ggCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCgg                                                                                                                                       <  2:683863/150‑1 (MQ=255)
                gCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGACATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGgc                                                                                                                                      <  2:766382/150‑1 (MQ=255)
                  cGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGgcgc                                                                                                                                    <  2:3147619/150‑1 (MQ=255)
                    gggCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCa                                                                                                                                  <  2:55762/150‑1 (MQ=255)
                     ggCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCAc                                                                                                                                 <  2:8121385/150‑1 (MQ=255)
                        gcagcaGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGAt                                                                                                                              >  1:1098466/1‑150 (MQ=255)
                          agcagcGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATgg                                                                                                                            <  2:2600264/150‑1 (MQ=255)
                          agcagcGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATgg                                                                                                                            >  1:7766877/1‑150 (MQ=255)
                           gcagcGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGc                                                                                                                           >  2:5535828/1‑150 (MQ=255)
                                        gAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTc                                                                                                              >  2:5806615/1‑150 (MQ=255)
                                              aTCCGGACGGTTGGTGACGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTtc                                                                                                        <  2:659853/150‑1 (MQ=255)
                                               tCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTtcg                                                                                                       >  2:614893/1‑150 (MQ=255)
                                                  ggACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGTGACCAACGGCGTCGATTtcgtcg                                                                                                    <  2:8633603/150‑1 (MQ=255)
                                                   gACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGa                                                                                                   <  1:1177047/150‑1 (MQ=255)
                                                         tggtggCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAga                                                                                             >  2:7216889/1‑150 (MQ=255)
                                                          ggtggCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAgat                                                                                            <  2:1017322/150‑1 (MQ=255)
                                                           gtggCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAgatg                                                                                           <  1:5868898/150‑1 (MQ=255)
                                                              gcggcgATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGatgatg                                                                                        >  2:3809524/1‑150 (MQ=255)
                                                                    atgacgatgaTGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACggg                                                                                  <  1:1443595/150‑1 (MQ=255)
                                                                                ccATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTg                                                                      <  1:8100356/150‑1 (MQ=255)
                                                                                ccATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTg                                                                      <  2:449769/150‑1 (MQ=255)
                                                                                ccATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTg                                                                      <  2:7547091/150‑1 (MQ=255)
                                                                                ccATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTg                                                                      <  2:4838822/150‑1 (MQ=255)
                                                                                 catcatTCATCTCGAGGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTgg                                                                     <  2:4319362/150‑1 (MQ=255)
                                                                                 catcatTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTgg                                                                     >  1:6965751/1‑150 (MQ=255)
                                                                                 catcatTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTgg                                                                     >  2:2512623/1‑150 (MQ=255)
                                                                                  atcatTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGc                                                                    >  2:4062971/1‑150 (MQ=255)
                                                                                  atcatTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGc                                                                    <  1:4960848/150‑1 (MQ=255)
                                                                                  atcatTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGc                                                                    <  1:8997115/150‑1 (MQ=255)
                                                                                   tcattcatCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGcc                                                                   >  1:5497659/1‑150 (MQ=255)
                                                                                   tcattcatCTCGAAGCCGTCCAGCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGcc                                                                   >  1:8469302/1‑150 (MQ=255)
                                                                                    cattcatCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCt                                                                  >  1:7956185/1‑150 (MQ=255)
                                                                                         atCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCg                                                             >  2:2935748/1‑150 (MQ=255)
                                                                                          tctcGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGa                                                            >  2:4587353/1‑150 (MQ=255)
                                                                                          tctcGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGa                                                            >  2:3329315/1‑150 (MQ=255)
                                                                                           ctcGCAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCGCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGaa                                                           <  1:5156975/150‑1 (MQ=255)
                                                                                            tcGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTtcgtcg                                                                                                    >  2:7974091/1‑108 (MQ=255)
                                                                                            tcGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAAc                                                          >  1:6702958/1‑150 (MQ=255)
                                                                                            tcGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAAc                                                          <  1:7577698/150‑1 (MQ=255)
                                                                                              gAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACAt                                                        <  2:1308248/150‑1 (MQ=255)
                                                                                              gAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACAt                                                        >  1:7854910/1‑150 (MQ=255)
                                                                                              aaaaccatCCATTTCGACCAGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACAt                                                        <  1:383191/143‑1 (MQ=255)
                                                                                                aaccatCCATTTCGACCAGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGt                                                      >  1:3435216/6‑150 (MQ=255)
                                                                                                  ccatCCATTTCGACCAGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTca                                                    <  1:887679/147‑2 (MQ=255)
                                                                                                  ccatCCATTTCGACCAGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTca                                                    <  1:664039/147‑2 (MQ=255)
                                                                                                   cGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTcgc                                                   <  2:7802991/150‑1 (MQ=255)
                                                                                                    gTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTcgcg                                                  >  2:1860872/1‑150 (MQ=255)
                                                                                                    gTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTcgcg                                                  >  2:8950196/1‑150 (MQ=255)
                                                                                                    atCCATTTCGACCAGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTcacg                                                  >  1:2663815/2‑150 (MQ=255)
                                                                                                      ccATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGca                                                >  2:4022173/1‑150 (MQ=255)
                                                                                                      ccATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGca                                                >  2:3673881/1‑150 (MQ=255)
                                                                                                         tttcgaccaGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGgacgc                                             >  2:1130395/6‑148 (MQ=255)
                                                                                                         tttcgaccaGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGgacgc                                             >  2:2075946/6‑148 (MQ=255)
                                                                                                         tttcgaccaGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGgacgc                                             >  1:9160751/6‑148 (MQ=255)
                                                                                                         tctcCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCacac                                             >  1:5359954/1‑150 (MQ=255)
                                                                                                         tctcCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCacac                                             <  1:5803026/150‑1 (MQ=255)
                                                                                                           tcCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACgg                                           >  1:7894732/1‑150 (MQ=255)
                                                                                                             gaccaGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCt                                         >  2:8261036/2‑150 (MQ=255)
                                                                                                             caccaGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCt                                         <  2:881764/150‑1 (MQ=255)
                                                                                                              accaGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTg                                        >  1:2978429/1‑150 (MQ=255)
                                                                                                              accaGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTg                                        >  2:7447542/1‑150 (MQ=255)
                                                                                                                 agaagcTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCa                                     <  1:2087284/145‑1 (MQ=255)
                                                                                                                  gaagcTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCAc                                    <  1:3207248/146‑1 (MQ=255)
                                                                                                                  gCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCAc                                    <  2:1651888/150‑1 (MQ=255)
                                                                                                                   aagcTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCCCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTCTTCGAACATGTCACGTACGCGGCTTGCAcc                                   <  2:4925940/147‑1 (MQ=255)
                                                                                                                        ggTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAacac                              <  1:7374781/150‑1 (MQ=255)
                                                                                                                          ttGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACAccc                            >  2:6716230/1‑150 (MQ=255)
                                                                                                                           tGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCa                           <  1:877317/150‑1 (MQ=255)
                                                                                                                            gAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATg                                                       >  1:6229187/1‑121 (MQ=255)
                                                                                                                            gAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATg                                                       <  2:6229187/121‑1 (MQ=255)
                                                                                                                            gAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCAc                          >  1:7271662/1‑150 (MQ=255)
                                                                                                                                     cTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTc                 >  2:234839/1‑150 (MQ=255)
                                                                                                                                      tCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTcc                <  1:2364449/150‑1 (MQ=255)
                                                                                                                                      tCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTcc                >  1:3054722/1‑150 (MQ=255)
                                                                                                                                       cACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCa               >  1:2776188/1‑150 (MQ=255)
                                                                                                                                       cACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCa               <  2:1812137/150‑1 (MQ=255)
                                                                                                                                       cACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCa               >  2:3518805/1‑150 (MQ=255)
                                                                                                                                        aCGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGacac                              <  1:3734631/134‑1 (MQ=255)
                                                                                                                                        aCGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCAc              >  2:2434607/1‑150 (MQ=255)
                                                                                                                                          gCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGa            <  2:5160762/150‑1 (MQ=255)
                                                                                                                                           ttcgtcgTTGCCGCCGCCTAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACaaa           <  1:8234012/149‑1 (MQ=255)
                                                                                                                                           ttcgtcgTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACaaa           <  2:6185917/149‑1 (MQ=255)
                                                                                                                                           cTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGaa           <  1:489589/150‑1 (MQ=255)
                                                                                                                                            tcgtcgTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGACACCGACACCCAc                          >  1:4104222/1‑134 (MQ=255)
                                                                                                                                              gtcgtGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTc        <  1:1165232/150‑1 (MQ=255)
                                                                                                                                              gtcgtGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTc        <  1:4035321/150‑1 (MQ=255)
                                                                                                                                               tcgtGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACAccc                            >  2:6458305/1‑129 (MQ=255)
                                                                                                                                               tcgtGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTcc       >  2:2587391/1‑150 (MQ=255)
                                                                                                                                               tcgtGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTcc       >  1:1479414/1‑150 (MQ=255)
                                                                                                                                                cgtGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCg      <  1:6890251/150‑1 (MQ=255)
                                                                                                                                                 gtGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGa     >  1:3619231/1‑150 (MQ=255)
                                                                                                                                                  tGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGaa    >  1:4477550/1‑150 (MQ=255)
                                                                                                                                                   gACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAAc   >  1:215886/1‑150 (MQ=255)
                                                                                                                                                    aCCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAAcc  >  1:5840461/1‑150 (MQ=255)
                                                                                                                                                    aCCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAAcc  <  2:8585494/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCTGGCGGTCGAAGCGGCCGGGGCGCAGCAGCGCAGGGTCGAGCACATCCGGACGGTTGGTGGCGGCGATGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACC  >  NC_002947/5367562‑5367859

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: