Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,406,214 G→A 6.9% P132S (CCG→TCG)  PP_4750 ← amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,406,2140GA 6.9% 171.0 / 2.3 58P132S (CCG→TCG) PP_4750amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (24/30);  new base A (2/2);  total (26/32)
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.13e-01

ACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCGAACTGGCCGGCGGGAATGGCCGCGAACGCACCGCGGTACACCTCGGACATATGCGCCGCCGTGATCAGCCCCAGGCCAACCACCGCTG  >  NC_002947/5406066‑5406360
                                                                                                                                                    |                                                                                                                                                   
aCGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCgg                                                                                                                                                     >  4:96209/1‑149 (MQ=255)
      aGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAacac                                                                                                                                               <  4:62952/149‑1 (MQ=255)
       ggCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAAcacc                                                                                                                                              <  5:103782/149‑1 (MQ=255)
              aaaGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATc                                                                                                                                         <  5:47833/147‑1 (MQ=255)
                       cTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGtt                                                                                                                              <  3:295323/149‑1 (MQ=255)
                       cTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGCTGGGATTGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGtt                                                                                                                              >  5:310130/1‑149 (MQ=255)
                              aCCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGTTACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGgcgc                                                                                                                      <  3:249305/149‑1 (MQ=255)
                               ccTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGcc                                                                                                                     <  4:311288/149‑1 (MQ=255)
                                         cTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAgg                                                                                                               <  3:63716/145‑1 (MQ=255)
                                         cTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAgg                                                                                                               >  4:63716/1‑145 (MQ=255)
                                         cTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTaa                                                                                                           <  8:325119/149‑1 (MQ=255)
                                          tGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAg                                                                                                          >  6:122771/1‑149 (MQ=255)
                                          tGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAg                                                                                                          >  3:315548/1‑149 (MQ=255)
                                            aaaGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGc                                                                                                         <  4:247980/148‑1 (MQ=255)
                                             aaGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGcac                                                                                                       <  3:219887/149‑1 (MQ=255)
                                              aGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCAGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGcac                                                                                                       <  3:307188/148‑1 (MQ=255)
                                               gCTGATTTCACCAACGCCAATCGTCGAGGCGACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGcaca                                                                                                      >  1:30670/1‑148 (MQ=255)
                                                 tGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAg                                                                                                   >  3:165448/1‑149 (MQ=255)
                                                      tCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAggcgg                                                                                               <  4:315548/148‑1 (MQ=255)
                                                        aCCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCtt                                                                                            >  2:281790/1‑149 (MQ=255)
                                                         ccAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTc                                                                                           >  5:289874/1‑149 (MQ=255)
                                                          taaCGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCg                                                                                          <  8:254647/148‑1 (MQ=255)
                                                          cAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTc                                                                                           >  2:130676/1‑148 (MQ=255)
                                                          cAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCg                                                                                          <  8:280349/149‑1 (MQ=255)
                                                              gCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGTCTTCGAACt                                                                                      <  3:96209/149‑1 (MQ=255)
                                                               ccAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCGAACTg                                                                                     >  2:167359/1‑149 (MQ=255)
                                                                 aaTCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCGAACTGGc                                                                                   >  2:9205/1‑149 (MQ=255)
                                                                      tcgAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCCTTAATAGCGCCGACAGGTTAAGCACATAGGCGGGTTGGAACTGGCcggc                                                                               >  1:115778/1‑149 (MQ=255)
                                                                        gAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCGAACTGGCcggcgg                                                                             <  7:315083/148‑1 (MQ=255)
                                                                              aCGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCGAACTGGCCGGCGGGAATGGc                                                                      >  6:241816/1‑149 (MQ=255)
                                                                              aCGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACACACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCg                                                                                          >  1:207878/1‑129 (MQ=255)
                                                                              aCGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACACACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCg                                                                                          <  2:207878/129‑1 (MQ=255)
                                                                                  ccGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCGAACTGGCCGGCGGGAATGGCcgc                                                                   >  1:201588/1‑148 (MQ=255)
                                                                                       tCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCGAACTGGCCGGCGGGAATGGCCGCGAACGc                                                             <  1:281790/149‑1 (MQ=255)
                                                                                         cTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCGAACTGGCCGGCGGGAATGGCCGCGAACGCAc                                                           >  2:257534/1‑149 (MQ=255)
                                                                                             aGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCGAACTGGCCGGCGGGAATGGc                                                                      >  2:131576/1‑134 (MQ=255)
                                                                                             aGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCGAACTGGCCGGCGGGAATGGc                                                                      <  1:131576/134‑1 (MQ=255)
                                                                                                  gCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCGAACTGGCCGGCGGGAATGGCCGCGAACGCACCGCGGTaca                                                  >  5:131192/1‑149 (MQ=255)
                                                                                                    cAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGcaca                                                                                                      <  3:17565/95‑1 (MQ=255)
                                                                                                    cAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTCCGGCAACACCACATCGAAAAACCCTT‑GCGGCGCCGACAGCTTAAGcaca                                                                                                      >  4:17565/1‑95 (MQ=255)
                                                                                                     aaTGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGc                                                                                                         >  5:155405/1‑91 (MQ=255)
                                                                                                     aaTGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGc                                                                                                         <  6:155405/91‑1 (MQ=255)
                                                                                                      aTGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCGAACTGGCCGGCGGGAATGGCCGCGAACGCACCGCGGTACACCTc                                              <  2:201588/149‑1 (MQ=255)
                                                                                                            tAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCGAACTGGCCGGCGGGAATGGCCGCGAACGCACCGCGGTACACCTCGGACat                                        <  1:167359/149‑1 (MQ=255)
                                                                                                            tAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGACGACAGGTTAAGCACATAGGCGGCTTCGAACTGGCCGGCGGGAATGGCCGCGAACGCACCGCGGTACACCTCGGACat                                        <  5:241816/149‑1 (MQ=255)
                                                                                                              ggTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCGAACTGGCCGGCGGGAATGGCCGCGAACGCACCGCGGTACACCTCGGACatat                                      <  4:165448/149‑1 (MQ=255)
                                                                                                               gTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCGAACTGGCCGGCGGGAATGGCCGCGAACGCACCGCGGTACACCCCGGACATATg                                     >  8:247459/1‑149 (MQ=255)
                                                                                                                tGGCTGAGGTGGGAATGACGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCGAACTGGCCGGCGGGAATGGCCGCGAACGCACCGCGGTACACCTCGGACatat                                      <  5:33281/147‑1 (MQ=255)
                                                                                                                tGGCTGAGGTGGGAATGACGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCGAACTGGCCGGCGGGAATGGCCGCGAACGCACCGCGGTACACCTCGGACatat                                      >  6:33281/1‑147 (MQ=255)
                                                                                                                                  cGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCGAACTGGCCGGCGGGAATGGCCGCGAACGCACCGCGGTACACCTCGGACATATGCGCCGCCGTGATCAGcccc                  <  1:130676/149‑1 (MQ=255)
                                                                                                                                    aTACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCGAACTGGCCGGCGGGAATGGCCGCGAACGCACCGCGGTACACCCCGGACATATGCGCCGCCGTGATCAGCCCCAg                <  7:247459/149‑1 (MQ=255)
                                                                                                                                     tACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCGAACTGGCCGGCGGGAATGGCCGCGAACGCACCGCGGTACACCTCGGACATATGCGCCACCGTGATCAGCCCCAgg               <  6:289874/149‑1 (MQ=255)
                                                                                                                                          aCCAATTGCGACAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCGAACTGGCCGGCGGGAATGGCCGCGAACGCACCGCGGTACACCTCGGACATATGCGCCGCCGTGATCAGCCCCAgg               <  7:38826/144‑1 (MQ=255)
                                                                                                                                          aCCAATTGCGACAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCGAACTGGCCGGCGGGAATGGCCGCGAACGCACCGCGGTACACCTCGGACATATGCGCCGCCGTGATCAGCCCCAgg               <  7:43946/144‑1 (MQ=255)
                                                                                                                                          aCCAATTGCGACAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCGAACTGGCCGGCGGGAATGGCCGCGAACGCACCGCGGTACACCTCGGACATATGCGCCGCCGTGATCAGCCCCAgg               >  8:38826/1‑144 (MQ=255)
                                                                                                                                          aCCAATTGCGACAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCGAACTGGCCGGCGGGAATGGCCGCGAACGCACCGCGGTACACCTCGGACATATGCGCCGCCGTGATCAGCCCCAgg               >  8:43946/1‑144 (MQ=255)
                                                                                                                                           ccAATTGCTGCAACACCACATCGAAAAACAGTT‑GCGGCGCCGACAGGTTAAGAACATAGGCGGCTTCGAACTGGCCGGCGGGAATGGCCGCGAACGCACCGCGGTACACCTCGGACATATGCGCCGCCGTGATCAGCCCCAGGCCAac          >  5:35242/1‑148 (MQ=255)
                                                                                                                                                 gCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCGAACTGGCCGGCGGGAATGGCCGCGAACGCACCGCGGTACACCTCGGACATATGCGCCGCCGTGATCAGCCCCAGGCCAACCCCCGCt   >  3:175691/1‑149 (MQ=255)
                                                                                                                                                  cGGCAACACCACATCGAATAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCGAACTTGCCGGCGGGAATGGCCGCGAACGCACCGCGGTACACCTCGGACATATGCGCCGCCGTGATCAGCCCCATGCCAA‑CACCGCTg  <  6:310130/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACGCTCAGGCCCTGAAAGGTTTGCTGTGAAACCTGGTAAGCCTGAAAGCTGATTTCACCAACGCCAATCGTCGAGGCCACGGCCGAGTCCTTGAGCAGGCCAATGGCGTAGGTGGCTGAGGTGGGAATGGCGATACGCACCAATTGCGGCAACACCACATCGAAAAACCGTT‑GCGGCGCCGACAGGTTAAGCACATAGGCGGCTTCGAACTGGCCGGCGGGAATGGCCGCGAACGCACCGCGGTACACCTCGGACATATGCGCCGCCGTGATCAGCCCCAGGCCAACCACCGCTG  >  NC_002947/5406066‑5406360

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: ‑