Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,987,689 C→G 5.4% R161P (CGC→CCC)  PP_2608 ← shikimate dehydrogenase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,987,6890CG 5.4% 889.0 / 28.7 261R161P (CGC→CCC) PP_2608shikimate dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base C (131/116);  new base G (10/4);  total (141/120)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 277 total aligned reads displayed.

GGGCGGCGGACAACGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGAGGAA  >  NC_002947/2987540‑2987837
                                                                                                                                                     |                                                                                                                                                    
gggCGGCGGACAACGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCAtgc                                                                                                                                                      <  2:6235772/150‑1 (MQ=255)
gggCGGCGGACAACGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCAtgc                                                                                                                                                      >  2:2621290/1‑150 (MQ=255)
gggCGGCGGACAACGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCAtgc                                                                                                                                                      >  2:6087251/1‑150 (MQ=255)
gggCGGCGGACAACGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCAtgc                                                                                                                                                      <  1:3706709/150‑1 (MQ=255)
       ggACAACGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGt                                                                                                                                               <  1:6309248/150‑1 (MQ=255)
       ggACAACGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGt                                                                                                                                               <  1:4987804/150‑1 (MQ=255)
         aCAACGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACGCATGGGTGCGGtat                                                                                                                                             <  1:4847194/150‑3 (MQ=255)
         aCAACGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACGCATGGGTGCGGtat                                                                                                                                             <  2:4810173/150‑3 (MQ=255)
          cAACGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCt                                                                                                                                            <  1:5123272/150‑1 (MQ=255)
            aCGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTgg                                                                                                                                          <  1:6471156/150‑1 (MQ=255)
            aCGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTgg                                                                                                                                          <  1:4311993/150‑1 (MQ=255)
            aCGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCAACCATGCGTGCGGTGCTgg                                                                                                                                          <  1:3308689/150‑1 (MQ=255)
             cGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGTCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTggg                                                                                                                                         <  1:351436/150‑1 (MQ=255)
             cGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTggg                                                                                                                                         <  2:6326389/150‑1 (MQ=255)
             cGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTggg                                                                                                                                         <  2:5566376/150‑1 (MQ=255)
             cGGCAGCTCTGCGCGGGTGCCCAGGCCGACTGGGGACGCGTTGGCCCCCAGGTCGACGTCTTCCAGCCCACGGAACTGTGGGCTGCCCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCAGACCGCTGCCATGCGTGCGGTGCTggg                                                                                                                                         <  2:1988049/150‑1 (MQ=255)
               gCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGt                                                                                                                                       <  2:3282539/150‑1 (MQ=255)
                    tctGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCag                                                                                                                                  >  1:5708356/1‑150 (MQ=255)
                    tctGCACGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCag                                                                                                                                  >  2:4649495/1‑150 (MQ=255)
                     ctGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCaga                                                                                                                                 <  1:543041/150‑1 (MQ=255)
                     ctGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCaga                                                                                                                                 <  1:4665137/150‑1 (MQ=255)
                     ctGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCaga                                                                                                                                 <  1:4603986/150‑1 (MQ=255)
                     ctGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTCGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCaga                                                                                                                                 <  1:4153514/150‑1 (MQ=255)
                       gcgcGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCAGGCGTGCGGTGCTGGGGTCGCAGAgg                                                                                                                               <  2:1243588/150‑1 (MQ=255)
                            ggTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTAt                                                                                                                          <  1:281485/150‑1 (MQ=255)
                                   aTGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCg                                                                                                                   >  2:836092/1‑150 (MQ=255)
                                   aTGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCg                                                                                                                   >  2:5584524/1‑150 (MQ=255)
                                       cGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATcc                                                                                                               >  1:597105/1‑150 (MQ=255)
                                       cGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATcc                                                                                                               <  1:1619470/150‑1 (MQ=255)
                                        gACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATccc                                                                                                              >  2:740793/1‑150 (MQ=255)
                                        gACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATccc                                                                                                              <  2:5635863/150‑1 (MQ=255)
                                          cTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTg                                                                                                            >  1:3462303/1‑150 (MQ=255)
                                           tGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTTCTGACCGTAAGCCCGGGGTAGCCGTTGCCTAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGc                                                                                                           <  2:305034/150‑1 (MQ=255)
                                           tGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACGGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGc                                                                                                           >  2:3529923/1‑150 (MQ=255)
                                              ggACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGACGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTc                                                                                                        >  1:1841831/1‑150 (MQ=255)
                                              ggACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTc                                                                                                        >  1:2737544/1‑150 (MQ=255)
                                                      tGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTg                                                                                                <  1:4188297/150‑1 (MQ=255)
                                                        gCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCa                                                                                              >  1:564082/1‑150 (MQ=255)
                                                         ccaccaGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCAt                                                                                             <  2:3460084/150‑1 (MQ=255)
                                                          caccaGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATa                                                                                            <  2:4167177/150‑1 (MQ=255)
                                                           accaGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAg                                                                                           <  1:536237/150‑1 (MQ=255)
                                                           accaGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAg                                                                                           >  2:6186034/1‑150 (MQ=255)
                                                           accaGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAg                                                                                           <  1:2670502/150‑1 (MQ=255)
                                                                gTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAAt                                                                                      >  2:4093048/1‑150 (MQ=255)
                                                                gTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTAAGCCAGTGCATAGGCAAt                                                                                      <  1:6873714/150‑1 (MQ=255)
                                                                 tCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATg                                                                                     >  1:3072903/1‑150 (MQ=255)
                                                                   gAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGc                                                                                   >  1:6795426/1‑150 (MQ=255)
                                                                   gAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGc                                                                                   <  1:2828555/150‑1 (MQ=255)
                                                                   gAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGc                                                                                   >  1:2138201/1‑150 (MQ=255)
                                                                     aGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCAc                                                                                 <  1:1352049/150‑1 (MQ=255)
                                                                     aGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCAc                                                                                 <  2:4332985/150‑1 (MQ=255)
                                                                      gTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACt                                                                                <  1:6087251/150‑1 (MQ=255)
                                                                         ttCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCCATGGCACTGcc                                                                             >  2:5467426/1‑150 (MQ=255)
                                                                         ttCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGcc                                                                             >  1:6220996/1‑150 (MQ=255)
                                                                         ttCCAGCCCACTGAACTGTGTGCTGACCGAAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGTCGATCCCTGCCTCAGCAAGTGCATAGGCAATGGCACTGcc                                                                             >  1:1532595/1‑150 (MQ=255)
                                                                          tCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCa                                                                            >  1:6723439/1‑150 (MQ=255)
                                                                           tcAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCaa                                                                           >  1:4749743/2‑150 (MQ=255)
                                                                           ccAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCaa                                                                           <  2:1141015/150‑1 (MQ=255)
                                                                           ccAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCaa                                                                           <  2:2204835/150‑1 (MQ=255)
                                                                              gCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAacac                                                                        >  2:568500/1‑150 (MQ=255)
                                                                               cccACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACAcc                                                                       >  1:1031545/1‑150 (MQ=255)
                                                                               cccACTGAACAGAGTGCAGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACAcc                                                                       >  1:1459707/1‑150 (MQ=255)
                                                                                  aCTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACAccgcc                                                                    >  2:6618254/1‑150 (MQ=255)
                                                                                       aCTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAAc                                                               <  2:6516662/150‑1 (MQ=255)
                                                                                         tgtgtgCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCg                                                             >  2:6525199/1‑150 (MQ=255)
                                                                                         tgtgtgCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCg                                                             <  2:1587285/150‑1 (MQ=255)
                                                                                         tgtgtgCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCg                                                             >  2:4907079/1‑150 (MQ=255)
                                                                                          gtgtgCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGa                                                            >  2:593392/1‑150 (MQ=255)
                                                                                             tgCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCa                                                         >  1:4749262/1‑150 (MQ=255)
                                                                                                 gACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAg                                                     >  1:6465259/1‑150 (MQ=255)
                                                                                                          cccGGGGAAGCCGTTGCCGAGCAGTTCGCATCCCGCACCCATGCGTGCGGTGCTGGGGTCGGAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGGCGCAACCGATCACCCGTGGCCGttt                                            <  2:3959227/150‑1 (MQ=255)
                                                                                                          cccGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGttt                                            <  2:1361013/150‑1 (MQ=255)
                                                                                                          cccGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGttt                                            >  2:2906474/1‑150 (MQ=255)
                                                                                                                     cGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCt                                 >  1:3269875/1‑150 (MQ=255)
                                                                                                                          ccGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAg                            >  2:6549765/1‑150 (MQ=255)
                                                                                                                           cGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGc                           >  2:1951057/1‑150 (MQ=255)
                                                                                                                           cGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGc                           >  2:5046405/1‑150 (MQ=255)
                                                                                                                           cGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGc                           >  1:6162859/1‑150 (MQ=255)
                                                                                                                            gAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATATGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGcc                          >  1:1770186/1‑150 (MQ=255)
                                                                                                                            gAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGcc                          >  2:2725958/1‑150 (MQ=255)
                                                                                                                             aGCAGTTCGCATACCGCACCCATGCTTGCGGTGCTGGGGTCGCAGAGGGTAATGATGTCGATCCCTGCCTCAGCCAGTGCATAGCCAATGGCACTGCAAACACCTCCGCAACCGATCACCAGTGCCCGTTTGCCCGCAGGCTCGAAGCCa                         >  1:794349/1‑150 (MQ=255)
                                                                                                                             aGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCa                         <  1:2259098/150‑1 (MQ=255)
                                                                                                                                  ttCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCtt                    >  2:1859672/1‑150 (MQ=255)
                                                                                                                                  ttCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCtt                    <  1:2395554/150‑1 (MQ=255)
                                                                                                                                     cagcaCCGCACGCATGGGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTtgtg                 >  1:1037842/5‑150 (MQ=255)
                                                                                                                                     cagcaCCGCACGCATGGGTGCGGTACTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATg                                                                                     <  2:4494554/78‑1 (MQ=38)
                                                                                                                                     cagcaCCGCACGCATGGGTGCGGTACTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATg                                                                                     >  1:4494554/5‑82 (MQ=38)
                                                                                                                                      cATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTgc                >  1:6562/1‑150 (MQ=255)
                                                                                                                                      cATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTgc                <  2:2903845/150‑1 (MQ=255)
                                                                                                                                         aCCGCACGCATGGGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCa             >  1:2595438/1‑150 (MQ=255)
                                                                                                                                         aCCGCACGCATGGGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCa             >  1:2866846/1‑150 (MQ=255)
                                                                                                                                         aCCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCa             >  2:5717383/1‑150 (MQ=255)
                                                                                                                                          ccGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAg            <  2:5856730/150‑1 (MQ=255)
                                                                                                                                             cACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAgcgc         >  2:281426/1‑150 (MQ=255)
                                                                                                                                                ccATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCTCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGa      >  1:250958/1‑150 (MQ=255)
                                                                                                                                                ccATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGa      <  1:854817/150‑1 (MQ=255)
                                                                                                                                                ccATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGACCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGa      <  1:1983523/150‑1 (MQ=255)
                                                                                                                                                   tgcgtgcgGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGAGGa   >  1:3832090/1‑150 (MQ=255)
                                                                                                                                                    gcgtgcgGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACTAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGAGGaa  <  2:3684952/150‑1 (MQ=255)
                                                                                                                                                    gcgtgcgGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGAGGaa  <  2:4773420/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGGCGGCGGACAACGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGAGGAA  >  NC_002947/2987540‑2987837

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: