Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 858,325 C→G 8.7% A168P (GCA→CCA)  phr ← deoxyribodipyrimidine photolyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947858,3250CG 8.7% 439.1 / 24.1 138A168P (GCA→CCA) phrdeoxyribodipyrimidine photolyase
Reads supporting (aligned to +/- strand):  ref base C (64/62);  new base G (4/8);  total (68/70)
Fisher's exact test for biased strand distribution p-value = 3.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 142 total aligned reads displayed.

CAGGAAGCGTGACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTGCGCCTGGATGGCCGACTTTT  >  NC_002947/858176‑858474
                                                                                                                                                     |                                                                                                                                                     
cAGGAAGCGTGACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGCTTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTGCACCCGCCGTGc                                                                                                                                                       >  2:3873005/1‑150 (MQ=255)
  ggAAGCGTGACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCa                                                                                                                                                     >  2:4851207/1‑150 (MQ=255)
   gAAGCGTGACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAg                                                                                                                                                    <  1:2612083/150‑1 (MQ=255)
    aaGCGTGACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAgg                                                                                                                                                   <  2:1552019/150‑1 (MQ=255)
        gTGACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTTGTCGCGCAATGCGGGCGTTGTTTTTCCGAACCCCTCTACATGTTGCGTGATTGGGTCGCTGCCGATGTTCATCGGGGCCTGGCGCTTCACCCGCCCTGGCActgcgg                                                                                                                                               <  1:4640759/150‑1 (MQ=255)
        gTGACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCCTGGCAcggcgg                                                                                                                                               <  1:2298670/150‑1 (MQ=255)
         tGACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCCGGGGTAGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCAACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCCAGCCAcggcggg                                                                                                                                              >  1:3886505/1‑150 (MQ=255)
         tGACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCCTGGCAcggcggg                                                                                                                                              <  2:1633842/150‑1 (MQ=255)
             aGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAgg                                                                                                                                          >  1:1139559/1‑150 (MQ=255)
             aGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAgg                                                                                                                                          <  1:2785811/150‑1 (MQ=255)
             aGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCCTGGCAcggcgggtgaa                                                                                                                                          <  1:1570638/150‑5 (MQ=255)
               gCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCt                                                                                                                                        >  2:609315/1‑150 (MQ=255)
                   gACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCgg                                                                                                                                    <  1:4055900/150‑1 (MQ=255)
                    aCTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGt                                                                                                                                   >  1:3923516/1‑150 (MQ=255)
                    aCTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGt                                                                                                                                   <  1:3629193/150‑1 (MQ=255)
                    aCTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGt                                                                                                                                   <  1:905707/150‑1 (MQ=255)
                          ccTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAgggcggg                                                                                                                                              >  2:2882375/1‑133 (MQ=255)
                          ccTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAgggcggg                                                                                                                                              <  1:2882375/133‑1 (MQ=255)
                           cTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATg                                                                                                                            >  2:2190141/1‑150 (MQ=255)
                           cTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATg                                                                                                                            <  1:2472690/150‑1 (MQ=255)
                            tCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGc                                                                                                                           >  1:3404609/1‑150 (MQ=255)
                                cTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCaa                                                                                                                       >  1:3669021/1‑150 (MQ=255)
                                cTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCaa                                                                                                                       >  1:1088933/1‑150 (MQ=255)
                                  tCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGCCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAgg                                                                                                                     <  1:429733/150‑1 (MQ=255)
                                           ggCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGtcttc                                                                                                            >  2:416476/1‑150 (MQ=255)
                                           ggCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGtcttc                                                                                                            <  2:4871576/150‑1 (MQ=255)
                                             ccAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTcttctt                                                                                                          <  2:2126214/150‑1 (MQ=255)
                                             ccAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTcttctt                                                                                                          <  2:2812598/150‑1 (MQ=255)
                                             ccAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTcttctt                                                                                                          <  2:3056904/150‑1 (MQ=255)
                                              cAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTg                                                                                                         >  2:1388333/1‑150 (MQ=255)
                                               aGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGa                                                                                                        <  2:2576016/150‑1 (MQ=255)
                                                gTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGaa                                                                                                       >  1:575913/1‑150 (MQ=255)
                                                 tGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAAc                                                                                                      <  2:3490928/150‑1 (MQ=255)
                                                  ggTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACt                                                                                                     <  1:964432/150‑1 (MQ=255)
                                                   gTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTg                                                                                                    <  1:2608297/150‑1 (MQ=255)
                                                     cgcgCACTGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGc                                                                                                  <  2:4444160/150‑1 (MQ=255)
                                                     cgcgCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGc                                                                                                  <  2:2267412/150‑1 (MQ=255)
                                                      gcgcAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCt                                                                                                 <  2:2131847/150‑1 (MQ=255)
                                                      gcgcAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCt                                                                                                 >  2:5084013/1‑150 (MQ=255)
                                                      gcgcAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCt                                                                                                 <  2:2132072/150‑1 (MQ=255)
                                                          aaTGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAac                                                                                             >  2:2549821/1‑150 (MQ=255)
                                                             gCGGGCGTGGTTTTTCCGAACCCCTCCACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACAcc                                                                                          <  1:343169/150‑1 (MQ=255)
                                                                  cGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGaa                                                                                     <  1:2459208/150‑1 (MQ=255)
                                                                  cGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGaa                                                                                     <  1:3638979/150‑1 (MQ=255)
                                                                      gTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTaa                                                                                 <  1:4636969/150‑1 (MQ=255)
                                                                       tttttCCGAACCCCTCTACATGTTGAGGGATAGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAAAACATGGAAGTAAt                                                                                >  2:3985274/1‑150 (MQ=255)
                                                                        ttttCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATc                                                                               <  1:1722322/150‑1 (MQ=255)
                                                                         tttCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGTTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCg                                                                              <  1:2205537/150‑1 (MQ=255)
                                                                         tttCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCg                                                                              >  1:1826305/1‑150 (MQ=255)
                                                                             cGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCAc                                                                          >  2:907835/1‑150 (MQ=255)
                                                                              gAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACt                                                                         >  2:4454559/1‑150 (MQ=255)
                                                                                aCCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGc                                                                       <  1:4953617/150‑1 (MQ=255)
                                                                                aCCCCACTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCTCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCACAAGCCAGGTCTTCTTGAACTCGCTGAACACCTGGAAGTAATCGCCACTGc                                                                       <  2:3899600/150‑1 (MQ=255)
                                                                                  cccTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCga                                                                     >  2:4434587/1‑150 (MQ=255)
                                                                                  cccTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCga                                                                     >  1:4436394/1‑150 (MQ=255)
                                                                                  cccTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCga                                                                     >  1:1534207/1‑150 (MQ=255)
                                                                                    ctctACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGTGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCgagt                                                                   >  2:1637963/1‑150 (MQ=255)
                                                                                    ctctACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCgagt                                                                   >  2:298976/1‑150 (MQ=255)
                                                                                    ctctACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCgagt                                                                   >  2:298927/1‑150 (MQ=255)
                                                                                      ctACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCgagtga                                                                 <  2:1224406/150‑1 (MQ=255)
                                                                                      ctACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCgagtga                                                                 <  2:1224311/150‑1 (MQ=255)
                                                                                         cATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTa                                                              <  2:4694505/150‑1 (MQ=255)
                                                                                             ttGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGt                                                          >  2:3011390/1‑150 (MQ=255)
                                                                                                 gggATTGGGTCGCTGCCGATGGTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTTCGGTGccc                                                      <  2:3844621/150‑1 (MQ=255)
                                                                                                   gATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCcgg                                                    >  2:662057/1‑150 (MQ=255)
                                                                                                   gATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCcgg                                                    >  1:3947502/1‑150 (MQ=255)
                                                                                                   gATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCcgg                                                    >  2:745710/1‑150 (MQ=255)
                                                                                                    aTTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCcgga                                                   >  1:2104509/1‑150 (MQ=255)
                                                                                                     ttGGGTCGCTGCCGATGTTCAGCGTGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCcggac                                                  <  1:1637963/150‑1 (MQ=255)
                                                                                                     ttGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCcggac                                                  >  1:2881110/1‑150 (MQ=255)
                                                                                                        ggTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCcggacgga                                               <  1:2348803/150‑1 (MQ=255)
                                                                                                        ggTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGAGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTCGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCcggacgga                                               >  2:4439800/1‑150 (MQ=255)
                                                                                                         gTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGgaa                                              <  1:2839009/150‑1 (MQ=255)
                                                                                                         gTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGgaa                                              <  2:309712/150‑1 (MQ=255)
                                                                                                         gTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGgaa                                              <  1:2627451/150‑1 (MQ=255)
                                                                                                         gTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGgaa                                              <  1:1033770/150‑1 (MQ=255)
                                                                                                         gTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGgaa                                              >  2:3967665/1‑150 (MQ=255)
                                                                                                                 cGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCt                                      >  1:3743/1‑150 (MQ=255)
                                                                                                                  gATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTg                                     >  2:4789394/1‑150 (MQ=255)
                                                                                                                       tCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGa                                <  2:1595830/150‑1 (MQ=255)
                                                                                                                            ggggCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGACCTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTgg                           <  2:1972530/150‑1 (MQ=255)
                                                                                                                            ggggCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTgg                           <  2:376551/150‑1 (MQ=255)
                                                                                                                            ggggCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTgg                           >  2:3719176/1‑150 (MQ=255)
                                                                                                                                  tGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAGGAACAGCTGGTCGAGGTGGCTGTgc                     <  2:3923516/150‑1 (MQ=255)
                                                                                                                                  tGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTgc                     >  1:4921339/1‑150 (MQ=255)
                                                                                                                                      gcTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTGCGCCt                 >  1:1872363/1‑150 (MQ=255)
                                                                                                                                      gcTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTGCGCCt                 >  1:4449957/1‑150 (MQ=255)
                                                                                                                                       cTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCCGATGCTCAAGGCCGGTCTTCTTGAACTGGCTGAACACCTGGACGTAATCGCCCCTGCGAGTGAGTCCGGTGCCCGGACGGAAGACCAGCTGGTCGAGGTGGCTGTGCGCCTg                >  2:5006309/1‑150 (MQ=255)
                                                                                                                                         tCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTGCGCCTGGa              >  1:3413319/1‑150 (MQ=255)
                                                                                                                                          tgcccGCCCTGGCACGGCGGGTAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGt                             <  1:3823974/132‑1 (MQ=255)
                                                                                                                                          tgcccGCCCTGGCACGGCGGGTAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGt                             <  1:3824053/132‑1 (MQ=255)
                                                                                                                                          tgcccGCCCTGGCACGGCGGGTAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGt                             >  2:3823974/3‑134 (MQ=255)
                                                                                                                                          tgcccGCCCTGGCACGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTGCGCCTGGAt             >  2:808630/3‑150 (MQ=255)
                                                                                                                                           aCCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTGCGCCTGGATg            >  2:2656504/1‑150 (MQ=255)
                                                                                                                                           aCCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGGCGAGGTGGCTGTGCGCCTGGATg            >  2:3398637/1‑150 (MQ=255)
                                                                                                                                            cccGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTGCGCCTGGATgg           >  2:1785200/1‑150 (MQ=255)
                                                                                                                                              cgccgTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCTCTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTGCGCCTGGATGGcc         <  2:3886505/150‑1 (MQ=255)
                                                                                                                                              cgccgTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTGCGCCTGGATGGcc         <  1:1358363/150‑1 (MQ=255)
                                                                                                                                                     ccAGGGCGTGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTGCGCCTGGATGGACGACtatt  >  1:4089234/1‑150 (MQ=255)
                                                                                                                                                     ccAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTGCGCCTGGATGGCCGACtttt  <  1:1388333/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGGAAGCGTGACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTGCGCCTGGATGGCCGACTTTT  >  NC_002947/858176‑858474

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: