Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,794,292 C→A 7.8% C45* (TGC→TGA) ‡ PP_2447 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,794,2920CA 7.8% 155.7 / 3.1 51C45* (TGC→TGA) ‡PP_2447hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/22);  new base A (2/2);  total (27/24)
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.49e-01

CCGCTGTCGGACCTATGATCGGTGGCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATGGCGCTGGCGGCAGTGCTGTCGCAGCAGGGGGCGCCTATTACGCCGGGC  >  NC_002947/2794144‑2794439
                                                                                                                                                     |                                                                                                                                                   
ccgcTGTCGGACCTATGATCGGTGGCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTgc                                                                                                                                                     <  6:301683/149‑1 (MQ=255)
  gcTGTCGGACCTATGATCGGTGGCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTgcgc                                                                                                                                                   <  1:59399/149‑1 (MQ=255)
  gcTGTCGGACCTATGATCGGTGGCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTgcgc                                                                                                                                                   >  2:297227/1‑149 (MQ=255)
   cTGTCGGACCTATGATTGGTGGCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGcc                                                                                                                                                  >  2:44111/1‑149 (MQ=255)
      tCGGACCTATGATCGGTGGCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGgc                                                                                                                                               >  5:360883/1‑149 (MQ=255)
      tCGGACCTATGATCGGTGGCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGgc                                                                                                                                               >  3:102129/1‑149 (MQ=255)
              aTGATCGGTGGCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGcc                                                                                                                                       >  3:117890/1‑149 (MQ=255)
              aTGATCGGTGGCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGcc                                                                                                                                       >  2:253680/1‑149 (MQ=255)
              aTGATCGGTGGCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCGCAGCCGgccaggg                                                                                                                                          >  4:259481/1‑141 (MQ=255)
              aTGATCGGTGGCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCGCAGCCGgccaggg                                                                                                                                          <  3:259481/146‑6 (MQ=255)
                        gCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCGCAGCCGgcca                                                                                                                                             >  8:143258/1‑131 (MQ=255)
                        gCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCGCAGCCGgcca                                                                                                                                             <  7:143258/133‑3 (MQ=255)
                          cAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTc                                                                                                                             >  3:80270/1‑147 (MQ=255)
                           agagCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAg                                                                                                                          >  8:92809/1‑149 (MQ=255)
                            gagCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAgg                                                                                                                         <  2:145400/149‑1 (MQ=255)
                             agCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAgg                                                                                                                         <  1:256777/149‑1 (MQ=255)
                              gCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGc                                                                                                                       >  8:177631/1‑149 (MQ=255)
                                cGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCt                                                                                                                     <  5:370311/149‑1 (MQ=255)
                                  ggCGCTGCCTGCCTTGCGTTGGATCAATCGTGTGTCCATGGATCAAAATTCCTTCAAGGGTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGc                                                                                                                   <  7:285307/149‑1 (MQ=255)
                                  ggCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGc                                                                                                                   >  8:28943/1‑149 (MQ=255)
                                      ctgcctgcCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCa                                                                                                               <  5:173273/149‑1 (MQ=255)
                                         cctgcctgGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGc                                                                                                             <  1:253680/148‑1 (MQ=255)
                                                  cGTTTGTACATGGGATTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCGACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGGAACCTCAAGGGCCTGCCGCAGCGAGTCGAGAt                                                                                                   <  5:109485/149‑1 (MQ=255)
                                                    ttGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATca                                                                                                 <  1:297227/149‑1 (MQ=255)
                                                        aaCAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATcagcag                                                                                             >  4:18057/1‑149 (MQ=255)
                                                         aCAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCagcagc                                                                                            <  4:80270/149‑1 (MQ=255)
                                                           aaGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTTCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGt                                                                                          >  3:232183/1‑149 (MQ=255)
                                                                       aTGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCg                                                                              <  4:102129/149‑1 (MQ=255)
                                                                         ggAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTg                                                                            >  7:282046/1‑149 (MQ=255)
                                                                          gAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCTGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTAGCGTg                                                                            <  6:360883/148‑1 (MQ=255)
                                                                                     cTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAAcc                                                                >  1:52358/1‑149 (MQ=255)
                                                                                     cTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAAc                                                                 <  1:173610/148‑1 (MQ=255)
                                                                                      ttCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAAc                                                                       >  8:93932/1‑141 (MQ=255)
                                                                                      ttCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAAc                                                                       <  7:93932/141‑1 (MQ=255)
                                                                                      ttCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACca                                                               >  1:333324/1‑149 (MQ=255)
                                                                                          aGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGc                                                           >  2:311981/1‑149 (MQ=255)
                                                                                           ggCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCg                                                          >  2:217646/1‑149 (MQ=255)
                                                                                              tGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCgggg                                                       >  5:70158/1‑149 (MQ=255)
                                                                                               gCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGc                                                      <  3:18057/149‑1 (MQ=255)
                                                                                                  ccccGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCAt                                                   <  4:117890/149‑1 (MQ=255)
                                                                                                    ccGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATgg                                                 >  2:107257/1‑149 (MQ=255)
                                                                                                    ccGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATgg                                                 >  8:86084/1‑149 (MQ=255)
                                                                                                     cggcggCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATGgc                                                <  7:92809/149‑1 (MQ=255)
                                                                                                              cTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATGGCGCTggcgg                                        <  1:311981/148‑1 (MQ=255)
                                                                                                                  ctgcctgcctgcTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATGGCGCTGGCGGCAGTg                                   <  7:177631/149‑1 (MQ=255)
                                                                                                                      ctgcctgcTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATGGCGCTGGCGGCAg                                     >  4:45193/1‑143 (MQ=255)
                                                                                                                      ctgcctgcTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATGGCGCTGGCGGCAg                                     <  3:45193/143‑1 (MQ=255)
                                                                                                                                gctgCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATGGCGCTGGCGGCAGTGCTGTCGCAGCAggg                     >  2:300939/1‑149 (MQ=255)
                                                                                                                                   gCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATGGCGCTGGCGGCAGTGCTGTCGCAGCAgg                      >  3:262207/1‑145 (MQ=255)
                                                                                                                                   gCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATGGCGCTGGCGGCAGTGCTGTCGCAGCAgg                      <  4:262207/145‑1 (MQ=255)
                                                                                                                                     ccACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTTCCGCAGCGCGGCGAGCTCAGCAGCGTGCCGGTCTACCGGGGCAACGCCCACCCCCGCGGGGCCAAGGCCCTGGCGGCCGTGCGTTCGCAGCAGGGGgcg                 >  8:145904/1‑149 (MQ=255)
                                                                                                                                      cACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATGGCGCTGGCGGCAGTGCTGTCGCAGCAGGGGgcgc                <  3:243138/149‑1 (MQ=255)
                                                                                                                                       aCCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATGGCGCTGGCGGCAGTGCTGTCGCAGCAGGGGGCGcc               >  5:218182/1‑149 (MQ=255)
                                                                                                                                          cTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATGGCGCTGGCGGCAGTGCTGTCGCAGCAGGGGGCGCCTAt            <  4:232183/149‑1 (MQ=255)
                                                                                                                                                    gcgcCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATGGCGCTGGCGGCAGTGCTGTCGCAGCAGGGGGCGCCTATTACGCCGGGc  >  4:52756/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCGCTGTCGGACCTATGATCGGTGGCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATGGCGCTGGCGGCAGTGCTGTCGCAGCAGGGGGCGCCTATTACGCCGGGC  >  NC_002947/2794144‑2794439

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: