Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,695,020 C→G 7.3% A329G (GCC→GGC)  wspC → biofilm formation methyltransferase WspC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,695,0200CG 7.3% 168.3 / 3.6 55A329G (GCC→GGC) wspCbiofilm formation methyltransferase WspC
Reads supporting (aligned to +/- strand):  ref base C (25/26);  new base G (2/2);  total (27/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 = 3.20e-01

CGTTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCA  >  NC_002947/1694876‑1695164
                                                                                                                                                  |                                                                                                                                                 
cgtTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCgg                                                                                                                                               <  8:139791/149‑1 (MQ=255)
cgtTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCgg                                                                                                                                               <  6:116713/149‑1 (MQ=255)
  ttCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCt                                                                                                                                             >  5:205929/1‑149 (MQ=255)
                 gCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCACCGGCATTGGCCAGCCGGGCAATGCCggc                                                                                                                                <  4:189552/149‑2 (MQ=38)
                 gCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGc                                                                                                                              >  6:174137/1‑149 (MQ=255)
                  cACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGTCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑ACGGCATTGGCCAGCCGGGCAATGCCggccag                                                                                                                             <  1:273260/149‑5 (MQ=255)
                                     cAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTa                                                                                                            <  4:111179/146‑1 (MQ=255)
                                     cAGTGTGACGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTa                                                                                                            >  3:111179/1‑146 (MQ=255)
                                        tgtgCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAAt                                                                                                      <  2:226521/149‑1 (MQ=255)
                                          tgCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGc                                                                                                    <  5:134447/149‑1 (MQ=255)
                                            ccGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGcc                                                                                                   >  1:209400/1‑148 (MQ=255)
                                                 gCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGc                                                                                              >  4:199568/1‑148 (MQ=255)
                                                   cTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGc                                                                                                          >  2:218542/1‑134 (MQ=255)
                                                   cTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGc                                                                                                          <  1:218542/134‑1 (MQ=255)
                                                   cTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCt                                                                                             <  5:278410/147‑1 (MQ=255)
                                                   cTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCt                                                                                             >  6:278410/1‑147 (MQ=255)
                                                   cTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTAc                                                                                           >  5:25147/1‑149 (MQ=255)
                                                   cTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTAc                                                                                           <  3:70095/149‑1 (MQ=255)
                                                      gTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCCAGGCCCGTAGCGAATGCCAACGCTAc                                                                                           <  1:14636/147‑1 (MQ=255)
                                                      gTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCt                                                                                         >  8:310934/1‑148 (MQ=255)
                                                       ttCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGa                                                                                       >  5:210946/1‑149 (MQ=255)
                                                           cggcgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTCAGCCa                                                                                   <  7:310934/149‑1 (MQ=255)
                                                            ggcgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAAt                                                                                                      <  8:58145/129‑1 (MQ=255)
                                                            ggcgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAAt                                                                                                      >  7:58145/1‑129 (MQ=255)
                                                                    ctgcCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACCCTACCTGAGCCAGTACccgcc                                                                          >  8:149299/1‑149 (MQ=255)
                                                                               ggccaggccaGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCa                                                                                                                    >  6:38834/1‑96 (MQ=255)
                                                                               ggccaggccaGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCa                                                                                                                    <  5:38834/96‑1 (MQ=255)
                                                                                gccaggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgt                                                              <  6:25147/149‑1 (MQ=255)
                                                                                gccaggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgt                                                              <  6:164757/149‑1 (MQ=255)
                                                                                gccaggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgt                                                              >  2:2979/1‑149 (MQ=255)
                                                                                 ccaggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgtg                                                             <  3:199568/149‑1 (MQ=255)
                                                                                  caggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgtgt                                                            <  2:209400/149‑1 (MQ=255)
                                                                                       cagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGtactac                                                       >  1:337222/1‑149 (MQ=255)
                                                                                          cgcgAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActgg                                                    <  6:210946/149‑1 (MQ=255)
                                                                                             gAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActggctg                                                 >  7:123598/1‑149 (MQ=255)
                                                                                               gCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTggg                                               >  8:297403/1‑149 (MQ=255)
                                                                                                 ccGCGAGGACAGCCTCGTGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCATCGAGCAGGGCCGTAGCTAATGCCAACGCTACCTGAGCCAGTACCGGCCGTCGGCACAGGTGTACTACTGGCTGGGtt                                             <  1:146461/149‑1 (MQ=255)
                                                                                                 ccGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGtt                                             >  3:149924/1‑149 (MQ=255)
                                                                                                     gAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCt                                         >  8:134527/1‑149 (MQ=255)
                                                                                                       gCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCa                                       <  1:28510/149‑1 (MQ=255)
                                                                                                       gCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCa                                       <  5:75482/149‑1 (MQ=255)
                                                                                                               aCGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATAcc                               >  3:136798/1‑149 (MQ=255)
                                                                                                                  tgggggCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGaa                            <  4:149924/148‑1 (MQ=255)
                                                                                                                     cggCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGc                         >  2:38072/2‑149 (MQ=255)
                                                                                                                         cAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATg                     >  8:237912/1‑149 (MQ=255)
                                                                                                                         cAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATg                     <  6:8783/149‑1 (MQ=255)
                                                                                                                              aGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCgcag                <  7:149299/149‑1 (MQ=255)
                                                                                                                               gTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCgcag                >  6:64279/1‑148 (MQ=255)
                                                                                                                                   cTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAgg            <  7:297403/148‑1 (MQ=255)
                                                                                                                                      g‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCt        >  8:223263/1‑149 (MQ=255)
                                                                                                                                          cGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGAt                      <  5:236827/133‑1 (MQ=255)
                                                                                                                                          cGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGAt                      >  6:236827/1‑133 (MQ=255)
                                                                                                                                            gCATTGGCCAGCCGGGCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGc    >  2:224633/1‑149 (MQ=255)
                                                                                                                                            gCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGc    <  7:223263/149‑1 (MQ=255)
                                                                                                                                              aTTGGCCAGCCGGGCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCa  >  1:137765/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CGTTCACGACGTTGCCCGCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCA  >  NC_002947/1694876‑1695164

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: