Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,356,405 G→C 7.0% G94A (GGC→GCC)  PP_1181 → two‑component system response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,356,4050GC 7.0% 237.7 / 5.8 76G94A (GGC→GCC) PP_1181two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base G (36/33);  new base C (2/3);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.40e-01

ATCATCGACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATCGA  >  NC_002947/1356257‑1356546
                                                                                                                                                    |                                                                                                                                             
atcatcGACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAggg                                                                                                                                               <  1:246719/149‑1 (MQ=255)
 tcatcGACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGc                                                                                                                                              <  5:73748/149‑1 (MQ=255)
  catcGACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGcc                                                                                                                                             >  6:52147/1‑149 (MQ=255)
    tcGACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTg                                                                                                                                           >  3:328282/1‑149 (MQ=255)
    tcGACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTg                                                                                                                                           <  5:52147/149‑1 (MQ=255)
      gACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTgg                                                                                                                                          <  2:270581/148‑1 (MQ=255)
       aCAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGAc                                                                                                                                        <  6:366594/149‑1 (MQ=255)
        cagccagcCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGCCCTGGAcc                                                                                                                                       <  6:105273/149‑1 (MQ=255)
         agccagccCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCt                                                                                                                                      <  6:344159/149‑1 (MQ=255)
         agccagccCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGCCCTGGACCt                                                                                                                                      <  1:24997/149‑1 (MQ=255)
          gccagccCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGCCCTGGACCTg                                                                                                                                     <  5:315396/149‑1 (MQ=255)
           ccagccCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTgg                                                                                                                                    >  6:210508/1‑149 (MQ=255)
           ccagccCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTgg                                                                                                                                    <  5:276684/149‑1 (MQ=255)
           ccagccCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTgg                                                                                                                                    >  1:208980/1‑149 (MQ=255)
            cagccCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGc                                                                                                                                              >  3:234404/1‑138 (MQ=255)
            cagccCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGc                                                                                                                                              <  4:234404/138‑1 (MQ=255)
                        gTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCgacgac                                                                                                                         <  4:328282/147‑1 (MQ=255)
                        gTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACt                                                                                                                        <  2:122451/148‑1 (MQ=255)
                        gTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTa                                                                                                                       >  2:78560/1‑149 (MQ=255)
                         tGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCTGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTAc                                                                                                                      <  4:3490/149‑1 (MQ=255)
                              cTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGgc                                                                                                                                  <  3:123879/132‑1 (MQ=255)
                              cTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGgc                                                                                                                                  >  4:123879/1‑132 (MQ=255)
                                     tGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAAccc                                                                                                          >  7:67138/1‑149 (MQ=255)
                                         gCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCcgacga                                                                                                                          >  8:306574/1‑129 (MQ=255)
                                         gCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCcgacga                                                                                                                          <  7:306574/129‑1 (MQ=255)
                                         gCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTAc                                                                                                                      >  3:138241/1‑133 (MQ=255)
                                         gCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTAc                                                                                                                      <  4:138241/133‑1 (MQ=255)
                                          cTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTgc                                                                                                      <  8:67138/148‑1 (MQ=255)
                                           tCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGCCTACGTCTGCAAACCCGTgcgc                                                                                                    >  5:121649/1‑149 (MQ=255)
                                                ggCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACg                                                                                               >  4:67468/1‑149 (MQ=255)
                                                 tcGAGGACGGCCTGATCATCTGCCGCCGGGTGCGCAGCCAGTACCCGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGc                                                                                              <  6:121649/148‑1 (MQ=255)
                                                     ggACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTtgc                                                                                          >  4:276294/1‑149 (MQ=255)
                                                           ccTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCcgac                                                                                                                            <  3:234656/109‑1 (MQ=255)
                                                           ccTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCcgac                                                                                                                            >  4:234656/1‑109 (MQ=255)
                                                            cTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTgg                                                                                    >  7:321741/1‑148 (MQ=255)
                                                               aGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCAc                                                                                 <  5:210508/148‑1 (MQ=255)
                                                                   tCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGTCCAGGCCCTggac                                                                                                                                  >  6:236275/1‑93 (MQ=255)
                                                                   tCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGTCCAGGCCCTggac                                                                                                                                  <  5:236275/95‑3 (MQ=255)
                                                                   tCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGAc                                                                                                                                        <  6:107069/89‑1 (MQ=255)
                                                                   tCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGAc                                                                                                                                        >  5:107069/1‑89 (MQ=255)
                                                                   tCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTtg                                                                                           <  3:333620/134‑1 (MQ=255)
                                                                   tCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTtg                                                                                           >  4:333620/1‑134 (MQ=255)
                                                                   tCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATc                                                                            >  5:7750/1‑149 (MQ=255)
                                                                       ccgccgGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAgg                                                                        >  8:185483/1‑149 (MQ=255)
                                                                            gggTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGAGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAAACCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCct                                                                    >  1:263206/1‑148 (MQ=255)
                                                                                gcgcAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCtgctgc                                                               >  3:258894/1‑149 (MQ=255)
                                                                                 cgcAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCtgctgc                                                               <  8:321741/148‑1 (MQ=255)
                                                                                  gcAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCTCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTgcgc                                                             <  2:208980/149‑1 (MQ=255)
                                                                                  gcAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTgcgc                                                             >  2:41096/1‑149 (MQ=255)
                                                                                  gcAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTgcgc                                                             >  3:239349/1‑149 (MQ=255)
                                                                                      ccAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCa                                                         >  2:129131/1‑149 (MQ=255)
                                                                                          tACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCg                                                      >  3:342229/1‑148 (MQ=255)
                                                                                             cAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAAc                                                  >  2:177098/1‑149 (MQ=255)
                                                                                              aGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAAcc                                                 <  6:7750/149‑1 (MQ=255)
                                                                                                gggCCGATCCTCATGCTCACGGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGt                                               <  3:276294/149‑1 (MQ=255)
                                                                                                gggCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGt                                               <  4:239349/149‑1 (MQ=255)
                                                                                                gggCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCATGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGt                                               >  2:121267/1‑149 (MQ=255)
                                                                                                 ggCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTg                                              <  2:53638/149‑1 (MQ=255)
                                                                                                  gCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTg                                              >  6:253608/1‑148 (MQ=255)
                                                                                                     gATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACa                                          <  5:253608/149‑1 (MQ=255)
                                                                                                      aTCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAg                                         >  3:10333/1‑149 (MQ=255)
                                                                                                          tCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAg                                     >  4:39679/1‑149 (MQ=255)
                                                                                                             ctctctCCGCGCGCTGTGCTTAACTGGACCTGGTCTAGGGCCTGGACGTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAgcgc                                  <  6:309743/147‑1 (MQ=255)
                                                                                                                  aCCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGa                             >  1:207504/1‑149 (MQ=255)
                                                                                                                    cgcgcgCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACTTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGGAGCCAGGCCCTCCTGCGCCGCAGCGAAACCGTGGACCGCAAgcgc                                  >  5:309743/1‑142 (MQ=255)
                                                                                                                    cgcgcgCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGAcc                           >  4:19641/1‑149 (MQ=255)
                                                                                                                      cgcgcAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTg                         <  1:41096/149‑1 (MQ=255)
                                                                                                                           aGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCtt                    <  4:342229/149‑1 (MQ=255)
                                                                                                                             tgatgaACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGTTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCg                  <  3:67468/149‑1 (MQ=255)
                                                                                                                               atgaACTGGACCAGGTCCAGGGCCTGGAGCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGt                >  1:212892/1‑149 (MQ=255)
                                                                                                                                tgaACTGGCCCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCCGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTg               >  4:51553/1‑149 (MQ=255)
                                                                                                                                   aCTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCTGTgcgc            <  4:258894/149‑1 (MQ=255)
                                                                                                                                      cgcccAGGTCCAGGCCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTgg                                                                                    >  1:212244/4‑74 (MQ=255)
                                                                                                                                        aCCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCg        >  4:244039/1‑148 (MQ=255)
                                                                                                                                            ggTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATCg   <  3:39679/149‑1 (MQ=255)
                                                                                                                                             gTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATCGa  <  1:177098/149‑1 (MQ=255)
                                                                                                                                             gTCCAGGCCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATCGa  >  8:339381/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
ATCATCGACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATCGA  >  NC_002947/1356257‑1356546

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