Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,194,976 A→G 7.1% Q372R (CAG→CGG)  uxpA → lipoprotein UxpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,194,9760AG 7.1% 158.3 / 2.5 57Q372R (CAG→CGG) uxpAlipoprotein UxpA
Reads supporting (aligned to +/- strand):  ref base A (26/26);  new base G (2/2);  total (29/28)
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.60e-01

GGGCGTTTTTTCCCGGCCTTTGAACAGCAACTGGCGCCTGCCGCGCAGCTGCGTGTCGCGCAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTTGCTGCCGGGCTTGCGTTATGACTACCCGCTGGCAGCGGGGGCGGCCATTACCCGCGAA  >  NC_002947/1194840‑1195121
                                                                                                                                        |                                                                                                                                                  
gggCGTTTTTTCCCGGCCTTTGAACAGCAACTGGCGCCTGCCGCGCAGCTGCGTGTCGCGCAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGgcg                                                                                                                                       >  2:231574/1‑149 (MQ=255)
      tttttCCCGGCCTTTGAACAGCAACTGGCGCCTGCCGCGCAGCTGCGTGTCGCGCAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGc                                                                                                                                                 >  6:346857/1‑134 (MQ=255)
      tttttCCCGGCCTTTGAACAGCAACTGGCGCCTGCCGCGCAGCTGCGTGTCGCGCAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGc                                                                                                                                                 <  5:346857/134‑1 (MQ=255)
       ttttCCCGGCCTTTGAACAGCAACTGGCGCCTGCCGCGCAGCTGCGTGTCGCGCAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGa                                                                                                                                 >  5:21769/1‑148 (MQ=255)
             cGGCCTTTGAACAGCAACTGGCGCCTGCCGCGCAGCTGCGTGTCGCGCAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGt                                                                                                                          <  7:283405/149‑1 (MQ=255)
             cGGCCTTTGAACAGCAACTGGCGCCTGCCGCGCAGCTGCGTGTCGCGCAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGt                                                                                                                          >  3:191172/1‑149 (MQ=255)
              ggCCTTTGAACAGCAACTGGCGCCTGCCGCGCAGCTGCGGGTCGCGCAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTg                                                                                                                         <  8:194195/149‑1 (MQ=255)
               gCCTTTGAACAGCAACTGGCGCCTGCCGCGCAGCTGCGTGTCGCGCAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTgg                                                                                                                        <  8:42697/149‑1 (MQ=255)
                  tttGAACAGCAACTGGCGCCTGCCGCGCAGCTGCGTGTCGCGCAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGAGCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCgg                                                                                                                     >  1:284634/1‑149 (MQ=255)
                       aCAGCAACTGGCGCCTGCCGCGCAGCTGCGTGTCGCGCAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCt                                                                                                                >  7:2231/1‑149 (MQ=255)
                         aGCAACTGGCGCCTGCCGCGCAGCTGCGTGTCGCGCAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTgg                                                                                                              >  2:251603/1‑149 (MQ=255)
                         aGCAACTGGCGCCTGCCGCGCAGCTGCGTGTCGCGCAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTgg                                                                                                              <  8:2231/149‑1 (MQ=255)
                         aGCAACTGGCGCCTGCCGCGCAGCTGCGTGTCGCGCAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTgg                                                                                                              <  1:251603/149‑1 (MQ=255)
                           cAACTGGCGCCTGCCGCGCAGCTGCGTGTCGCGCAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGa                                                                                                            >  2:45026/1‑149 (MQ=255)
                                 gcgcCTGCCGCGCAGCTGCGTGTCGCGCAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCTCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGActgc                                                                                                        <  4:191172/147‑1 (MQ=255)
                                 gcgcCTGCCGCGCAGCTGCGTGTCGCGCAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGActgcc                                                                                                       <  6:21769/148‑1 (MQ=255)
                                    ccTGCCGCGCAGCTGCGTGTCGCGCAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCctgct                                                                                                   >  5:201787/1‑149 (MQ=255)
                                        ccGCGCAGCTGCGTGTCGCGCAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCa                                                                                                >  2:252301/1‑148 (MQ=255)
                                          gcgcAGCTGCGTGTCGCGCAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCCGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATc                                                                                              >  8:282649/1‑148 (MQ=255)
                                           cgcAGCTGCGTGTCGCGCAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCaa                                                                                            <  7:21726/149‑1 (MQ=255)
                                           cgcAGCTGCGTGTCGCGCAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCaa                                                                                            <  7:27001/149‑1 (MQ=255)
                                              aGCTGCGTGTCGCGCAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCagggc                                                                                         >  8:226874/1‑145 (MQ=255)
                                              aGCTGCGTGTCGCGCAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGcc                                                                                         <  8:143713/149‑1 (MQ=255)
                                               gCTGCGTGTCGCGCAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGcc                                                                                         <  1:252301/149‑1 (MQ=255)
                                                   cGTGTCGCGCAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGGTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTg                                                                                      <  1:110647/147‑1 (MQ=255)
                                                          cgcAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCcgacga                                                                              >  4:36082/1‑148 (MQ=255)
                                                                  cAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATAGCAGCAt                                                                     <  5:18942/149‑1 (MQ=255)
                                                                     gcgcAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATagca                                                                         <  5:160499/142‑1 (MQ=255)
                                                                     gcgcAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATagca                                                                         >  6:160499/1‑142 (MQ=255)
                                                                            tGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTtg                                                           >  8:133725/1‑149 (MQ=255)
                                                                                ccAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTtgctgc                                                       >  6:124476/1‑149 (MQ=255)
                                                                                 cAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTTGCTGcc                                                      >  7:86407/1‑149 (MQ=255)
                                                                                 cAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTTGCTGcc                                                      >  5:349010/1‑149 (MQ=255)
                                                                                 cAGCAACGGGCCGACCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCtgctgc                                                                                                 >  5:31552/1‑106 (MQ=255)
                                                                                 cAGCAACGGGCCGACCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCtgctgc                                                                                                 <  6:31552/106‑1 (MQ=255)
                                                                                  aGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTTGCTGCCg                                                     <  6:201787/149‑1 (MQ=255)
                                                                                             ggccACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTTGCTGCCGGGCTTGCGTTa                                          <  3:36082/149‑1 (MQ=255)
                                                                                              gccACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTTGCTGCCGGGCTTGCGTTAt                                         <  5:124476/149‑1 (MQ=255)
                                                                                                cACGCGGCCTGGCTCGACCAGCCGCTGGCCCGTGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTTGCTGCCGGGCTTGCGTTATGa                                       >  6:77191/1‑149 (MQ=255)
                                                                                                 aCGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTACATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTTGCTGCCGGGCTTGCGTTATGAc                                      >  4:320491/1‑149 (MQ=255)
                                                                                                    cgGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTTGCTGCCGGGCTTGCGTTATGACTAc                                   <  1:231574/149‑1 (MQ=255)
                                                                                                           gCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTTGCTGCCGGGCTTGCGTTATGACTACCCGCTgg                            <  8:86407/149‑1 (MQ=255)
                                                                                                             tCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTtg                                                           <  6:58356/116‑1 (MQ=255)
                                                                                                             tCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTtg                                                           >  5:58356/1‑116 (MQ=255)
                                                                                                             tCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTTGCTGCCGGGCTTGCGTTATGACTACCCGCTgg                            >  4:20814/1‑147 (MQ=255)
                                                                                                                   aGCCGCTGGCCCGCGCGCCCCGGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTTGCTGCCGGGCTTg                                               >  4:120766/1‑122 (MQ=255)
                                                                                                                   aGCCGCTGGCCCGCGCGCCCCGGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTTGCTGCCGGGCTTg                                               >  4:115847/1‑122 (MQ=255)
                                                                                                                   aGCCGCTGGCCCGCGCGCCCCGGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTTGCTGCCGGGCTTg                                               <  3:120766/122‑1 (MQ=255)
                                                                                                                   aGCCGCTGGCCCGCGCGCCCCGGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTTGCTGCCGGGCTTg                                               <  3:115847/122‑1 (MQ=255)
                                                                                                                    gCCGCTGGCCCGCGGGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTTGCTGCCGGGCTTGCGTTATGACTACCCGCTGGCAGCggggg                   >  7:300532/1‑149 (MQ=255)
                                                                                                                    gCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGAAACCCGTGGCGGTAGCTGGGACTGCCGGCAGCATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTTGCTGCCGGGCTTGCGTTATGACTACCCGCTGGCAGCggggg                   <  7:226874/149‑1 (MQ=255)
                                                                                                                            cccGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTTGCTGCCGGGCTTGCGTTATGACTACCCGCTGGCAGCGGGGGCGGCCAt            >  6:341572/1‑148 (MQ=255)
                                                                                                                             ccGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTTGCTGCCGGGCTTGCGTTATGACTACCCGCTGGCAGCGGGGGCGGCCATTa          >  2:344270/1‑149 (MQ=255)
                                                                                                                                cgcgcCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTTGCTGCCGGGCTTGCGTTATGACTACCCGCTGGCAGCGGGGGCGGCCATTAccc       <  6:349010/149‑1 (MQ=255)
                                                                                                                                  cgcCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTTGCTGCCGGGCTTGCGTTATGACTACCCGCTGGCAGCGGGGGCGGCCATTACCcgc     <  3:98895/149‑1 (MQ=255)
                                                                                                                                    ccccAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTCGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTTGCTGCCGGGCTTGCGTTATGACTACCCGCTGGCAGCGGGGGCGGCCATTACCCGCGa   <  5:341572/149‑1 (MQ=255)
                                                                                                                                      ccAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTTGCTGCCGGGCTTGCGTTATGACTACCCGCTGGCAGCGGGGGCGGCCATTACCAGCGaa  >  6:267743/1‑148 (MQ=255)
                                                                                                                                        |                                                                                                                                                  
GGGCGTTTTTTCCCGGCCTTTGAACAGCAACTGGCGCCTGCCGCGCAGCTGCGTGTCGCGCAATTGCAGGCGCAACTGCGCCAGCAACGGGCCGGCCACGCGGCCTGGCTCGACCAGCCGCTGGCCCGCGCGCCCCAGGCCTTGTGGCGCCGCGATACCCGTGGCGGTAGCTGGGACTGCCTGCTGCATCAAGCCTTGGCCGACGATAGCAGCATGCCCGTGTTGCTGCCGGGCTTGCGTTATGACTACCCGCTGGCAGCGGGGGCGGCCATTACCCGCGAA  >  NC_002947/1194840‑1195121

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