Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,932,050 C→G 7.0% G476A (GGC→GCC)  PP_3466 ← ABC transporter permease/ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,932,0500CG 7.0% 182.8 / 3.7 57G476A (GGC→GCC) PP_3466ABC transporter permease/ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (31/22);  new base G (2/2);  total (33/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 = 9.63e-01

GCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGTTCTGCGCCAGGGTGCCACGGAACAACACGATGT  >  NC_002947/3931909‑3932197
                                                                                                                                             |                                                                                                                                                   
gCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGc                                                                                                                                                >  3:151496/1‑147 (MQ=255)
  cgcgATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGGCGAGCGCGGCCAGCagtt                                                                                                                                            <  5:141418/149‑4 (MQ=255)
      aTCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATc                                                                                                                                        >  8:243692/1‑149 (MQ=255)
          cTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAgg                                                                                                                                    <  2:84156/149‑1 (MQ=255)
           ttCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGTCGAGCGCGGCCAGCagtt                                                                                                                                            <  4:200760/140‑4 (MQ=255)
           ttCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGGCGAGCGCGGCCAGCagtt                                                                                                                                            >  3:200760/1‑137 (MQ=255)
           ttCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGc                                                                                                                                   >  7:231717/1‑149 (MQ=255)
               cGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAg                                                                                                                               >  7:6062/1‑149 (MQ=255)
                 cTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCg                                                                                                                             <  4:215446/149‑1 (MQ=255)
                      tGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAgg                                                                                                                        <  5:155479/149‑1 (MQ=255)
                                 aCTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAg                                                                                                             >  2:182611/1‑149 (MQ=255)
                                   tGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGc                                                                                                           >  3:137739/1‑149 (MQ=255)
                                       gAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAgg                                                                                                                                    >  4:9592/1‑120 (MQ=255)
                                       gAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAgg                                                                                                                                    <  3:9592/120‑1 (MQ=255)
                                       gAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTc                                                                                                       >  6:14909/1‑149 (MQ=255)
                                       gAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTc                                                                                                       <  7:243692/149‑1 (MQ=255)
                                           tGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGc                                                                                                   <  8:6062/149‑1 (MQ=255)
                                                tCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCgcgc                                                                                              >  3:188578/1‑149 (MQ=255)
                                                tCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCgcgc                                                                                              >  3:5857/1‑149 (MQ=255)
                                                  gTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCgcgc                                                                                              <  7:32144/147‑1 (MQ=255)
                                                                  gCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTc                                                                            <  1:182611/149‑1 (MQ=255)
                                                                   ccGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCAACCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTc                                                                            >  7:18939/1‑148 (MQ=255)
                                                                             cAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCAc                                                                 >  7:189259/1‑149 (MQ=255)
                                                                                   tGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGTCAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCCCCACCCGTTCCAGTTCATCACGGCTgg                                                           >  6:192979/1‑149 (MQ=255)
                                                                                   tGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCCCCCGTTCCAGTTCATCACGGCTgg                                                           >  7:248268/1‑149 (MQ=255)
                                                                                   tGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTgg                                                           >  6:82366/1‑149 (MQ=255)
                                                                                   tGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTgg                                                           >  1:148086/1‑149 (MQ=255)
                                                                                     cACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGc                                                                                                                           >  2:183193/1‑83 (MQ=255)
                                                                                     cACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGc                                                                                                                           <  1:183193/83‑1 (MQ=255)
                                                                                        gAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCg                                                      >  6:23804/1‑149 (MQ=255)
                                                                                            aaTGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCa                                                  >  8:34384/1‑149 (MQ=255)
                                                                                             aTGGCCATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCAc                                                 <  5:82366/149‑1 (MQ=255)
                                                                                              tGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACg                                                <  4:151496/149‑1 (MQ=255)
                                                                                                 cGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCACAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACGAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCa                                             <  8:189259/149‑1 (MQ=255)
                                                                                                        ttgttTCTGGCCCCCGGCCAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCa                                      <  5:14909/149‑1 (MQ=255)
                                                                                                        ttgttTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGc                                        >  5:57856/1‑147 (MQ=255)
                                                                                                        ttgttTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCa                                      <  5:23804/149‑1 (MQ=255)
                                                                                                         tgttTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAg                                     >  4:56997/1‑149 (MQ=255)
                                                                                                           tttCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAgg                                    <  3:56997/148‑1 (MQ=255)
                                                                                                             tCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGtt                                  <  7:34384/148‑1 (MQ=255)
                                                                                                                 gCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGTTCTgcg                             >  8:70868/1‑149 (MQ=255)
                                                                                                                 gCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGTTCTgcg                             >  3:242869/1‑149 (MQ=255)
                                                                                                                 gCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGTTCTgc                              <  4:188578/148‑1 (MQ=255)
                                                                                                                  cccccGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGTTCTgcgc                            >  4:163463/1‑149 (MQ=255)
                                                                                                                           aaCTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGTTCTGCGCCAGGGTGcc                   >  5:78616/1‑149 (MQ=255)
                                                                                                                           aaCTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGTTCTGCGCCAGGGTGcc                   >  7:48433/1‑149 (MQ=255)
                                                                                                                            cttgctgGGCGAGCGCGGCCAGCAGGCGATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGTTCTGCGCCAGGGTGCCa                  >  1:148675/3‑149 (MQ=255)
                                                                                                                             ctgctgGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGTTCTGCGCCAGGGTGCCAc                 <  4:5857/149‑1 (MQ=255)
                                                                                                                             ctgctgGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGTTCTGCGCCAGGGTGCCAc                 >  2:229482/1‑149 (MQ=255)
                                                                                                                              tgctgGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGTTCTGCGCCAGGGTGCCACg                >  2:12966/1‑149 (MQ=255)
                                                                                                                                ctgGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTTTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGTTCTGCGCCAGGGTGCCACgg               <  2:148086/148‑1 (MQ=255)
                                                                                                                                ctgGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGTTCTGCGCCAGGGTGCCACGGa              <  3:163463/149‑1 (MQ=255)
                                                                                                                                   gCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGTTCTGCGCCAGGGTGCCACGGaa             <  1:12966/147‑1 (MQ=255)
                                                                                                                                      gcgcTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGTTCTGCGCCAGGGTGCCACGGAacaaca        >  7:237907/1‑149 (MQ=255)
                                                                                                                                          tCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCa                                      >  2:216387/1‑115 (MQ=255)
                                                                                                                                          tCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCa                                      <  1:216387/115‑1 (MQ=255)
                                                                                                                                            gCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGTTCTGCGCCAGGGTGCCACGGAACAACACGATGt  >  7:37561/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
GCCGCGATCACTTCACGCTCGGTGGCCTCATCCACTGCCGAAGTGGCCTCGTCCAGCACCAGGATCGCCGGGGCCTGCAACACTGCACGAGCAATGGCGATGCGTTGTTTCTGGCCCCCGGACAACTGCTGGCCGCGCTCGCCCAGCAAGCCATCCAGGCCCAGCGGCAGGCTGTCGACCAGGCTGTCCAGCCGCGCCAGGCGCACCACCCGTTCCAGTTCATCACGGCTGGCCTCGGGCACGCCATAAGCCAGGTTCTGCGCCAGGGTGCCACGGAACAACACGATGT  >  NC_002947/3931909‑3932197

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: