Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 328,688 T→C 6.8% I188V (ATC→GTC)  PP_0270 ← sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947328,6880TC 6.8% 169.7 / 2.3 59I188V (ATC→GTC) PP_0270sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (26/29);  new base C (2/2);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

CCGCCTGGCGCACTTCCTCCGGGCCATCCAGCGCCAGGGGAGGCTGCTCGAGGTTGCTGCCCAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCAACGGTCAGCCAGTGGCCATCGGCCAGCGGCACCGCGAACTGGATATGCGCCTGGGCCGT  >  NC_002947/328541‑328832
                                                                                                                                                   |                                                                                                                                                 
ccgccTGGCGCACTTCCTCCGGGCCATCCAGCGCCAGGGGAGGCTGCTCGAGGTTGCTGCCCAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATg                                                                                                                                                  >  5:134387/1‑149 (MQ=255)
ccgccTGGCGCACTTCCTCCGGGCCATCCAGCGCCAGGGGAGGCTGCTCGAGGTTGCTGCCCAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATg                                                                                                                                                  >  2:306663/1‑149 (MQ=255)
 cgccTGGCGCACTTCCTCCGGGCCATCCAGCGCCAGGGGAGGCTGCTCGAGGTTGCTGCCCAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATgc                                                                                                                                                 <  4:139540/149‑1 (MQ=255)
             ttCCTCCGGGCCATCCAGCGCCAGGGGAGGCTGCTCGAGGTTGCTGCCCAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATaca                                                                                                                                     >  2:343446/1‑149 (MQ=255)
               cctccGGGCCATCCAGCGCCAGGGGAGGCTGCTCGAGGTTGCTGCCCAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATacac                                                                                                                                    <  3:132934/148‑1 (MQ=255)
               cctccGGGCCATCCAGCGCCAGGGGAGGCTGCTCGAGGTTGCTGCCCAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATaca                                                                                                                                     <  6:31017/147‑1 (MQ=255)
               cctccGGGCCATCCAGCGCCAGGGGAGGCTGCTCGAGGTTGCTGCCCAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGACGCGCAGCGCATACAcc                                                                                                                                   >  8:20282/1‑149 (MQ=255)
               cctccGGGCCATCCAGCGCCAGGGGAGGCTGCTCGAGGTTGCTGCCCAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGACGCGCAGCGCATACAcc                                                                                                                                   >  8:14704/1‑149 (MQ=255)
                   cGGGCCATCCAGCGCCAGGGGAGGCTGCTCGAGGTTGCTGCCCAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCgct                                                                                                                              <  6:2786/149‑2 (MQ=255)
                      gCCATCCAGCGCCAGGGGAGGCTGCTCGAGGTTGCTGCCCAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCgcagca                                                                                                                           <  1:343446/149‑1 (MQ=255)
                      gCCATCCAGCGCCAGGGGAGGCTGCTCGAGGTTGCTGCCCAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCgcag                                                                                                                             >  6:169572/1‑147 (MQ=255)
                          tCCAGCGCCAGGGGAGGCTGCTCGAGGTTGCTGCCCAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCAt                                                                                                                                        <  7:265051/133‑1 (MQ=255)
                          tCCAGCGCCAGGGGAGGCTGCTCGAGGTTGCTGCCCAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCAt                                                                                                                                        >  8:265051/1‑133 (MQ=255)
                            cAGCGCCAGGGGAGGCTGCTCGAGGTTGCTGCCCAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAAt                                                                                                                     >  4:73399/1‑149 (MQ=255)
                                  cAGGGGAGGCTGCTCGAGGTTGCTGCCCAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGa                                                                                                               <  1:375693/149‑1 (MQ=255)
                                  cAGGGGAGGCTGCTCGAGGTTGCTGCCCAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGa                                                                                                               <  4:57470/149‑1 (MQ=255)
                                     gggAGGCTGCTCGAGGTTGCTGCCCAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGACGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCa                                                                                                            <  7:14704/149‑1 (MQ=255)
                                     gggAGGCTGCTCGAGGTTGCTGCCCAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGACGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCa                                                                                                            <  7:20282/149‑1 (MQ=255)
                                        aGGCTGCTCGAGGTTGCTGCCCAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCa                                                                                                         <  6:117374/149‑1 (MQ=255)
                                                        ctgcCCAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACTGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTcgc                                                                                          <  1:227604/148‑1 (MQ=255)
                                                         tgcCCAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTaa                                                                                                                      >  6:99028/1‑119 (MQ=255)
                                                         tgcCCAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTaa                                                                                                                      <  5:99028/119‑1 (MQ=255)
                                                         tgcCCAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTcgccg                                                                                        >  1:190328/1‑149 (MQ=255)
                                                          gcCCAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTcgccgc                                                                                       >  3:83253/1‑149 (MQ=255)
                                                            ccAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTc                                                                                     >  7:42342/1‑149 (MQ=255)
                                                            ccAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTc                                                                                     >  6:44494/1‑149 (MQ=255)
                                                             cAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCg                                                                                                                                 <  6:304082/105‑1 (MQ=255)
                                                             cAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCg                                                                                                                                 >  5:304082/1‑105 (MQ=255)
                                                                     gTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACAACAGCAACACCGCGACGATGCGCAGCGCATACACCCGCCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTg                                                                             <  1:285061/149‑1 (MQ=255)
                                                                      tGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTg                                                                             >  3:355805/1‑147 (MQ=255)
                                                                          gCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTc                                                                       <  4:83253/149‑1 (MQ=255)
                                                                                gCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGcc                                                                 >  8:149026/1‑149 (MQ=255)
                                                                                           cAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCAACGGt                                                        >  8:19147/1‑147 (MQ=255)
                                                                                                ggCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCAACGGTCAGCCAg                                                 >  7:5116/1‑149 (MQ=255)
                                                                                                   gcAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCAACGGTCAGCCAGTgg                                              <  4:6273/149‑1 (MQ=255)
                                                                                                   gcAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACAACGCGACGATGCGAATCGCATACACCCGCAGCAGGTAATCGATGATAAGCACACAGGCCGACTCGCAGATCCAGCCATGCAGTTCCTGGCCCGCAACGGTCAGCCAGTgg                                              >  3:307465/1‑149 (MQ=255)
                                                                                                         acaGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCAACGGTCAGCCAGTGGCCATCg                                        <  7:149026/149‑1 (MQ=255)
                                                                                                          caGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCAACGGTCAGCCAGTGGCCATCgg                                       >  5:32825/1‑149 (MQ=255)
                                                                                                             gcgcACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCAACGGTCAGCCAGTGGCCATCGGCCa                                    <  1:306663/149‑1 (MQ=255)
                                                                                                               gcACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCAACGGTCAGCCAGTGGCCATCGGCCAGc                                  <  5:44494/149‑1 (MQ=255)
                                                                                                               gcACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCAACGGTCAGCCAGTGGCCATCGGCCAGc                                  >  3:28390/1‑149 (MQ=255)
                                                                                                                   ggCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCAACGGTCAGCCAGTGGCCATCGGCCAGCGGCa                              >  5:151026/1‑149 (MQ=255)
                                                                                                                   ggCAACCAGGCACACAAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCACCGGTCAGCCAGTGGCCATCGGCCAGCGGCa                              >  5:355327/1‑149 (MQ=255)
                                                                                                                     cAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCAACGGTCAGCCAGTGGCCATCg                                        >  5:154088/1‑137 (MQ=255)
                                                                                                                     cAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCAACGGTCAGCCAGTGGCCATCg                                        <  6:154088/137‑1 (MQ=255)
                                                                                                                       aCCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCAACGGTCAGCCAGTGGCCATCGGCCAGCGGCACcgc                          <  7:19147/149‑1 (MQ=255)
                                                                                                                        ccAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCAACGGTCAGCCAGTGGCCATCGGCCAGCGGCACcgcg                         <  6:151026/149‑1 (MQ=255)
                                                                                                                        ccAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCAACGGTCAGCCAGAGGCCATCGGCCAGCGGCACcgc                          >  7:87440/1‑148 (MQ=255)
                                                                                                                           ggCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCAACGGTCAGCCAGTGGCCATCGGCCAGCGGCACCGCGAAc                      >  1:252146/1‑149 (MQ=255)
                                                                                                                             cacaCCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCAACGGTCAGCCAGTGGCCATCGGCCAGCGGCACCGCGAACTg                    <  2:190328/149‑1 (MQ=255)
                                                                                                                              acacCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCAACGGTCAGCCAGTGGCCATCGGCCAGCGGCACCGCGAACTgg                   <  4:355805/149‑1 (MQ=255)
                                                                                                                              acacCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCAACGGTCAGCCAGTGGCCATCGGCCAGCGGCACCGCGAACTgg                   <  3:73399/149‑1 (MQ=255)
                                                                                                                               cacCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCAACGGTCAGCCAGTGGCCATCGGCCAGCGGCACCGCGAACTGGa                  <  5:169572/149‑1 (MQ=255)
                                                                                                                                    gCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCAACGGTCAGCCAGTGGCCATCGGCCAGCGGCACCGCGAACTGGATATgc             <  2:252146/149‑1 (MQ=255)
                                                                                                                                         acCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCAACGGTCAGCCAGTGGCCATCGGCCAGCGGCACCGCGAACTGGATATGCGCCTg        >  2:22871/1‑149 (MQ=255)
                                                                                                                                          ccGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCAACGGTCAGc                                                    >  7:242407/1‑104 (MQ=255)
                                                                                                                                          ccGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCAACGGTCAGc                                                    <  8:242407/104‑1 (MQ=255)
                                                                                                                                          ccGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCAACGGTCAGCCAGTGGCCATCGGCCAGCGGCACCGCGAACTGGATATGCGCCTgg       <  6:134387/149‑1 (MQ=255)
                                                                                                                                               acgaTGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCAACGGTCAGCCAGTGGCCATCGGCCAGCGGCACCGCGAACTGGATATGCGCCTGGGCCGt  <  2:155194/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCGCCTGGCGCACTTCCTCCGGGCCATCCAGCGCCAGGGGAGGCTGCTCGAGGTTGCTGCCCAGGCCGCGTGCGGCATCGGCCAGGCGCCGCAACGGGCGCAGGCACAGGCGCACGGCAACCAGGCACACCAGCAACACCGCGACGATGCGCAGCGCATACACCCGCAGCAGGTAATCGCTGATCAGCACCCAGGCCGACTCGCCGCTCCAGCCCTGCAGTTCCTGGCCCGCAACGGTCAGCCAGTGGCCATCGGCCAGCGGCACCGCGAACTGGATATGCGCCTGGGCCGT  >  NC_002947/328541‑328832

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