Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 739,200 T→C 7.9% P72P (CCA→CCG) ‡ pilD ← leader peptidase/N‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947739,2000TC 7.9% 214.6 / 6.8 76P72P (CCA→CCG) ‡pilDleader peptidase/N‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (37/33);  new base C (3/3);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

CCACTACCAGCGACAGCAGCGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCCCACCAGCA  >  NC_002947/739052‑739348
                                                                                                                                                    |                                                                                                                                                    
ccACTACCAGCGACAGCAGCGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAgtgt                                                                                                                                                      <  5:227366/149‑1 (MQ=255)
    tACCAGCGACAGCAGCGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTggg                                                                                                                                                   >  4:154092/1‑148 (MQ=255)
      ccAGCGACAGCAGCGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCAc                                                                                                                                                >  4:228228/1‑149 (MQ=255)
      ccAGCGACAGCAGCGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCAc                                                                                                                                                >  5:4966/1‑149 (MQ=255)
      ccAGCGACAGCAGCGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCCCACACTGCGCGCAc                                                                                                                                                >  1:244338/1‑149 (MQ=255)
        aGCGACAGCAGCGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCg                                                                                                                                              >  3:157558/1‑149 (MQ=255)
         gCGACAGCAGCGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCgg                                                                                                                                             >  8:202727/1‑149 (MQ=255)
         gCGACAGCAGCGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCgg                                                                                                                                             >  5:151611/1‑149 (MQ=255)
            aCAGCAGCGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAg                                                                                                                                          <  5:235891/149‑1 (MQ=255)
             cagcagCGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAgg                                                                                                                                         >  3:225512/1‑149 (MQ=255)
             cagcagCGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCCCACACTGCGCGCACCggattc                                                                                                                                         <  2:244338/149‑5 (MQ=255)
                 agcGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCggg                                                                                                                                     >  3:244559/1‑149 (MQ=255)
                   cgcgcTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGggca                                                                                                                                   <  3:228228/149‑1 (MQ=255)
                     cgcTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGggcagg                                                                                                                                 <  6:252427/149‑1 (MQ=255)
                          gCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGcac                                                                                                                              <  7:202727/147‑1 (MQ=255)
                          gCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAg                                                                                                                            >  4:194123/1‑149 (MQ=255)
                          gCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAg                                                                                                                            <  6:4966/149‑1 (MQ=255)
                             aCTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAgg                                                                                                                                         <  5:203347/133‑1 (MQ=255)
                             aCTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAgg                                                                                                                                         >  6:203347/1‑133 (MQ=255)
                                 ccACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCCCACACTGCGCGCACCgg                                                                                                                                             <  5:29491/125‑1 (MQ=255)
                                 ccACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCCCACACTGCGCGCACCgg                                                                                                                                             >  6:29491/1‑125 (MQ=255)
                                  cACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAAc                                                                                                                    >  7:48307/1‑149 (MQ=255)
                                       ccGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAg                                                                                                                            <  4:31433/136‑1 (MQ=255)
                                       ccGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAg                                                                                                                            >  3:31433/1‑136 (MQ=255)
                                        cGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGt                                                                                                              <  8:221575/149‑1 (MQ=255)
                                              aCGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGtgctgc                                                                                                         <  3:154092/148‑1 (MQ=255)
                                                  aGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGt                                                                                                    >  5:175763/1‑149 (MQ=255)
                                                   ggCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTc                                                                                                   <  8:48307/149‑1 (MQ=255)
                                                   ggCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTc                                                                                                   <  3:244832/149‑1 (MQ=255)
                                                    gCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGggca                                                                                                                                   <  4:26331/116‑1 (MQ=255)
                                                    gCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGggca                                                                                                                                   >  3:23011/1‑116 (MQ=255)
                                                    gCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGggca                                                                                                                                   <  4:23011/116‑1 (MQ=255)
                                                    gCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGggca                                                                                                                                   >  3:26331/1‑116 (MQ=255)
                                                           cGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGcc                                                                                            <  3:194123/148‑1 (MQ=255)
                                                             gTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCg                                                                                                  <  3:19237/140‑1 (MQ=255)
                                                             gTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCg                                                                                                  >  4:19237/1‑140 (MQ=255)
                                                             gTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGcccc                                                                                          >  3:61815/1‑148 (MQ=255)
                                                                       gACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTg                                                                               >  5:58180/1‑149 (MQ=255)
                                                                            aCAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAGTACCTCTTGCGCCt                                                                          >  3:135483/1‑149 (MQ=255)
                                                                             cAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTc                                                                         <  4:244559/149‑1 (MQ=255)
                                                                              aGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCa                                                                        >  3:230229/1‑149 (MQ=255)
                                                                                cgcccgccCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACg                                                                      >  6:196960/1‑149 (MQ=255)
                                                                                     gcccAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGcc                                                                 <  3:172958/149‑1 (MQ=255)
                                                                                          aGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACTAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGc                                                            <  7:145827/149‑1 (MQ=255)
                                                                                           ggCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGtgc                                                                                                            <  1:207475/100‑1 (MQ=255)
                                                                                           ggCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGtgc                                                                                                            >  2:207475/1‑100 (MQ=255)
                                                                                           ggCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGcc                                                           <  5:196960/149‑1 (MQ=255)
                                                                                                 gTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGtg                                                                                                             >  5:195882/1‑93 (MQ=255)
                                                                                                 gTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGtg                                                                                                             <  6:195882/93‑1 (MQ=255)
                                                                                                   aGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCa                                                                                               <  8:128189/105‑1 (MQ=255)
                                                                                                   aGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCa                                                                                               >  7:128189/1‑105 (MQ=255)
                                                                                                    gCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCa                                                  <  4:157558/149‑1 (MQ=255)
                                                                                                     cTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCAt                                                 >  6:256222/1‑149 (MQ=255)
                                                                                                      tgatgaCCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCaa                                                                                        <  4:27668/109‑1 (MQ=255)
                                                                                                      tgatgaCCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCaa                                                                                        >  3:27668/1‑109 (MQ=255)
                                                                                                      tgatgaCCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCAt                                                 >  1:90537/1‑148 (MQ=255)
                                                                                                           aCCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCgg                                           <  6:151611/149‑1 (MQ=255)
                                                                                                             cGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCgg                                           >  5:148159/1‑147 (MQ=255)
                                                                                                                    ttCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGAt                                   <  6:148159/148‑1 (MQ=255)
                                                                                                                    ttCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATa                                  >  3:234261/1‑149 (MQ=255)
                                                                                                                    agcatccACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATa                                  <  8:7601/144‑1 (MQ=255)
                                                                                                                       tcCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGa                                    >  6:86701/1‑144 (MQ=255)
                                                                                                                       tcCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGa                                    <  5:86701/144‑1 (MQ=255)
                                                                                                                         ccACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACctct                                                                                 <  4:202212/97‑1 (MQ=255)
                                                                                                                         ccACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACctct                                                                                 >  3:202212/1‑97 (MQ=255)
                                                                                                                         ccACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATacac                               >  6:208771/1‑147 (MQ=255)
                                                                                                                         ccACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATacac                               >  8:107203/1‑147 (MQ=255)
                                                                                                                         ccACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCa                             >  8:103719/1‑149 (MQ=255)
                                                                                                                         ccACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCa                             <  4:234261/149‑1 (MQ=255)
                                                                                                                           acgcacgAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGc                           <  6:58180/149‑1 (MQ=255)
                                                                                                                                 gAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGtt                      >  6:104012/1‑148 (MQ=255)
                                                                                                                                 gAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGtt                      >  5:88657/1‑148 (MQ=255)
                                                                                                                                   aTCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGa                    >  1:140222/1‑148 (MQ=255)
                                                                                                                                    cccGGTGCCCACACTGCGCGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTg                     <  7:239416/145‑1 (MQ=255)
                                                                                                                                    cccGGTGCCCACACTGCGCGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTg                     >  8:239416/2‑146 (MQ=255)
                                                                                                                                          gcgcgcAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTg            <  7:103719/149‑1 (MQ=255)
                                                                                                                                          gcgcgcAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTg            <  7:107203/149‑1 (MQ=255)
                                                                                                                                                    tgGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCCCACcagca  <  4:225512/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCACTACCAGCGACAGCAGCGCGCTGGCCACTTCCACTACCGGGTAACGCAGGCTGATGCGGTTTTTGCAGGACGAACAGCGCCCGCCCAAGGCCAGGTAGCTGATGACCGGAATGTTCTCCCACGCACGAATCCGGTGCGCGCAGTGTGGGCACCGGGAGGCGGGCAGGCACAGGTCGAAACGCGCGTGCTGCGTGGTCGGCAGCCCCAATACCTCTTGCGCCTCACGCTGCCATTGCCGCTCCAGCATGATCGGCAGGCGATACACCAGCACGTTGATAAAGCTGCCCACCAGCA  >  NC_002947/739052‑739348

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: