Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,030,802 G→T 6.6% W367L (TGG→TTG) ‡ PP_4432 → Xaa‑Pro aminopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,030,8020GT 6.6% 971.1 / 37.3 291W367L (TGG→TTG) ‡PP_4432Xaa‑Pro aminopeptidase
Reads supporting (aligned to +/- strand):  ref base G (136/135);  new base T (12/8);  total (148/143)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 313 total aligned reads displayed.

CACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTG  >  NC_002947/5030654‑5030951
                                                                                                                                                    |                                                                                                                                                     
cacaAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTgg                                                                                                                                                      >  1:7855702/1‑150 (MQ=255)
cacaAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTgg                                                                                                                                                      >  2:3415740/1‑150 (MQ=255)
cacaAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTgg                                                                                                                                                      <  2:2265974/150‑1 (MQ=255)
cacaAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTgg                                                                                                                                                      <  1:3661054/150‑1 (MQ=255)
cacaAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTgg                                                                                                                                                      >  1:8204182/1‑150 (MQ=255)
 acaAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTggg                                                                                                                                                     <  1:481269/150‑1 (MQ=255)
 acaAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTggg                                                                                                                                                     >  2:5885918/1‑150 (MQ=255)
      gACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTc                                                                                                                                                <  1:4696655/150‑1 (MQ=255)
        cAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCga                                                                                                                                              <  2:5315480/150‑1 (MQ=255)
        cAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCga                                                                                                                                              >  1:8758094/1‑150 (MQ=255)
        cAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCga                                                                                                                                              <  2:4000735/150‑1 (MQ=255)
        cAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCga                                                                                                                                              >  2:8205609/1‑150 (MQ=255)
        cAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCga                                                                                                                                              <  1:8168495/150‑1 (MQ=255)
          gCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGCATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCgaga                                                                                                                                            <  1:3170292/150‑1 (MQ=255)
              ttGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTAc                                                                                                                                        <  1:8342972/150‑1 (MQ=255)
                 cGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCga                                                                                                                                     <  1:8920729/150‑1 (MQ=255)
                 cGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCga                                                                                                                                     <  1:5205271/150‑1 (MQ=255)
                   gTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCgaga                                                                                                                                   <  1:1158012/150‑1 (MQ=255)
                   gTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGCATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCgaga                                                                                                                                   <  1:8622103/150‑1 (MQ=255)
                      aCCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCa                                                                                                                                <  2:1954841/150‑1 (MQ=255)
                        caaaTCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATc                                                                                                                              >  1:8073940/3‑150 (MQ=255)
                         cAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATcc                                                                                                                             >  2:4400639/1‑150 (MQ=255)
                          aaTCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCt                                                                                                                            >  2:2967253/1‑150 (MQ=255)
                             cGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGAt                                                                                                                         <  1:2344498/150‑1 (MQ=255)
                              ggTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATc                                                                                                                        >  2:253506/1‑150 (MQ=255)
                              ggTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATc                                                                                                                        >  1:7671189/1‑150 (MQ=255)
                                cctgtccTATCCGCCGGACTGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCAc                                                                                                                      >  1:7004957/7‑150 (MQ=255)
                                  tCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTg                                                                                                                    <  1:5134286/150‑1 (MQ=255)
                                  tCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTg                                                                                                                    <  1:7425607/150‑1 (MQ=255)
                                    aGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGAt                                                                                                                  >  1:6417252/1‑150 (MQ=255)
                                       tATCCGCCGGACTGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGAAAca                                                                                                               <  1:299921/150‑2 (MQ=255)
                                         tCCGCCGGACTGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGAAACAgg                                                                                                             >  1:2537756/1‑150 (MQ=255)
                                         tCCGCCGGACTGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGAAACAgg                                                                                                             >  1:6695905/1‑150 (MQ=255)
                                            gccggactGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGAAACAGGCGt                                                                                                          >  1:8011057/8‑150 (MQ=255)
                                            accaGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGt                                                                                                          >  2:805329/1‑150 (MQ=255)
                                            accaGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGt                                                                                                          >  2:8529455/1‑150 (MQ=255)
                                            accaGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGt                                                                                                          >  1:6754781/1‑150 (MQ=255)
                                              cggactGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGAAACAGGCGTCg                                                                                                        >  1:1566088/6‑150 (MQ=255)
                                              caGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCg                                                                                                        <  2:1341508/150‑1 (MQ=255)
                                               aGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCga                                                                                                       <  2:5737116/150‑1 (MQ=255)
                                               aGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCga                                                                                                       <  1:8436342/150‑1 (MQ=255)
                                                 actGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGAAACAGGCGTCgaat                                                                                                     <  2:4149420/148‑3 (MQ=255)
                                                  ctGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGAAACAGGCGTCgaatg                                                                                                    <  1:3942309/149‑4 (MQ=255)
                                                   tGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGAAACAGGCGTCgaatgc                                                                                                   <  1:1656029/150‑5 (MQ=255)
                                                       gCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCt                                                                                               >  1:7213440/1‑150 (MQ=255)
                                                       gCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCt                                                                                               >  1:548101/1‑150 (MQ=255)
                                                           gcgcACCCTGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTAGGACCTTCCCCTTCATGCCGGGCCGATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTc                                                                                           <  2:1155778/150‑1 (MQ=255)
                                                             gcACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGa                                                                                         <  2:5588792/150‑1 (MQ=255)
                                                             gcACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGa                                                                                         <  1:4318263/150‑1 (MQ=255)
                                                              cACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGAc                                                                                        <  2:8854860/150‑1 (MQ=255)
                                                              cACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGAc                                                                                        <  2:7332081/150‑1 (MQ=255)
                                                                 cATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCGAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGActcc                                                                                     >  2:4676379/1‑150 (MQ=255)
                                                                       ccTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCaa                                                                               <  2:6407774/150‑1 (MQ=255)
                                                                        cTGCGCCCGAGCGACAACAGCGTGCTAAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAg                                                                              >  2:6943189/1‑150 (MQ=255)
                                                                              ccGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTc                                                                        <  2:5094827/150‑1 (MQ=255)
                                                                              ccGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTc                                                                        <  2:6876233/150‑1 (MQ=255)
                                                                              ccGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTc                                                                        >  2:3782811/1‑150 (MQ=255)
                                                                              ccGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTc                                                                        <  2:2352279/150‑1 (MQ=255)
                                                                              ccGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTc                                                                        <  2:6057332/150‑1 (MQ=255)
                                                                                  gCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCa                                                                    >  2:6684392/1‑150 (MQ=255)
                                                                                   cgcctacaGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATTGATGGATGATTGGGGGCTCGAGATTTCCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCaa                                                                   <  1:2821426/145‑1 (MQ=255)
                                                                                   cGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCaa                                                                   <  1:7370602/150‑1 (MQ=255)
                                                                                   cGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCaa                                                                   <  1:3797533/150‑1 (MQ=255)
                                                                                   cGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTAGGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCaa                                                                   <  2:9081540/150‑1 (MQ=255)
                                                                                    gACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAg                                                                  <  2:2186988/150‑1 (MQ=255)
                                                                                           gCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGAc                                                           <  2:7145524/150‑1 (MQ=255)
                                                                                             gTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGa                                                         >  2:1823751/1‑150 (MQ=255)
                                                                                                 tCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAg                                                     >  2:3321255/1‑150 (MQ=255)
                                                                                                   aaaCCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGtt                                                   <  2:5004261/150‑1 (MQ=255)
                                                                                                   aaaCCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGtt                                                   >  1:6899836/1‑150 (MQ=255)
                                                                                                   aaaCCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGtt                                                   >  2:8634492/1‑150 (MQ=255)
                                                                                                    aaCCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGttt                                                  <  2:6145766/150‑1 (MQ=255)
                                                                                                    aaCCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAACTTTTCGTCAGTCCTTGACGAAAAGctt                                                  <  1:2944172/150‑1 (MQ=255)
                                                                                                         ggTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTTGCAAACTTTTCGTCAGTCCTTGACGAAAAGCTTGCGAg                                             <  2:1092556/150‑1 (MQ=38)
                                                                                                          gTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGt                                            >  2:5799428/1‑150 (MQ=255)
                                                                                                          gTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGt                                            >  2:2735569/1‑150 (MQ=255)
                                                                                                          gTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCAATGATACTGGGGTCTAGACGATCTGCTCGACTCAACGCAAGAATTTCGTCAAGGACTGAAGAAAAGTTAGAAAGt                                            >  1:1966295/1‑150 (MQ=255)
                                                                                                           tATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTTGCAAACTTTTCGTCAGTCCTTGACGAAAAGCTTGCGAgga                                           <  2:6557436/150‑3 (MQ=37)
                                                                                                           aaTGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTc                                           <  1:8933819/149‑1 (MQ=255)
                                                                                                               aCCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTg                                       <  2:1719442/150‑1 (MQ=255)
                                                                                                               aCCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTg                                       <  1:957914/150‑1 (MQ=255)
                                                                                                               aCCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTg                                       <  1:7093033/150‑1 (MQ=255)
                                                                                                                 cTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGa                                     >  1:6602247/1‑150 (MQ=255)
                                                                                                                  ttCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAg                                    <  1:2713963/150‑1 (MQ=255)
                                                                                                                    ccACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGtt                                  >  2:232283/1‑150 (MQ=255)
                                                                                                                          cATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTg                            <  1:7138954/150‑1 (MQ=255)
                                                                                                                              ccTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCAc                        >  1:5333301/1‑150 (MQ=255)
                                                                                                                                tGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACgg                      >  1:1761254/1‑150 (MQ=255)
                                                                                                                                tGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACgg                      >  1:6361459/1‑150 (MQ=255)
                                                                                                                                    cTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTAt                  >  2:6190136/1‑150 (MQ=255)
                                                                                                                                    cTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTAt                  >  2:438390/1‑150 (MQ=255)
                                                                                                                                         gatggatgATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCa             <  1:3356342/150‑1 (MQ=255)
                                                                                                                                         gatggatgATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCa             >  2:4503142/1‑150 (MQ=255)
                                                                                                                                           tggatgATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGc           >  1:4486136/1‑150 (MQ=255)
                                                                                                                                           tggatgATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGc           <  2:3069486/150‑1 (MQ=255)
                                                                                                                                             gatgatTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTCTCTTCAGCGa         >  2:7941339/1‑150 (MQ=255)
                                                                                                                                              atgatTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGAc        >  1:5873094/1‑150 (MQ=255)
                                                                                                                                               tgatTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACa       >  1:578938/1‑150 (MQ=255)
                                                                                                                                                   tGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAAtgt   <  1:6252423/150‑1 (MQ=255)
                                                                                                                                                    gggggCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAAtgtg  <  2:6754781/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTG  >  NC_002947/5030654‑5030951

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: