Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 961,720 A→T 9.1% L232Q (CTG→CAG)  selD ← selenide, water dikinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947961,7200AT 9.1% 142.1 / 5.5 55L232Q (CTG→CAG) selDselenide, water dikinase
Reads supporting (aligned to +/- strand):  ref base A (22/27);  new base T (3/2);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01

CCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGTGGTGAGG  >  NC_002947/961573‑961866
                                                                                                                                                    |                                                                                                                                                  
ccGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCCCGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGCCCAGGTGGCCGAGCAg                                                                                                                                                   >  4:163111/1‑149 (MQ=255)
ccGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAg                                                                                                                                                   >  4:85308/1‑149 (MQ=255)
ccGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCCACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAg                                                                                                                                                   <  7:141304/149‑1 (MQ=255)
 cGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCTGCTCGGCCACCTGGTCGAGCTgg                                                                                                                                                  <  2:71098/149‑1 (MQ=255)
  gTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTGACCTTCGGCCAGCTCGACCAGGTGGCCGAGCAgg                                                                                                                                                  <  1:32464/149‑1 (MQ=255)
  gTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGGGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGc                                                                                                                                                 <  1:167820/149‑1 (MQ=255)
     gCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACAGTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCg                                                                                                                                               <  3:16088/148‑1 (MQ=255)
      cTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTACGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCg                                                                                                                                               >  6:53101/1‑147 (MQ=255)
      cTGTAGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGAGTAGTTCAGGTGCGCTGTGAGGCAGCTG‑CCTTCGGCCAGCTCGACAAGGTGGCCGAGCAGGCCg                                                                                                                                               >  8:115006/1‑147 (MQ=255)
            tAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGt                                                                                                                                       <  1:25084/149‑1 (MQ=255)
                  cgcAGGGTTCCACCGGGGATGCAGCCTACGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTc                                                                                                                                 <  5:53101/149‑1 (MQ=255)
                   gcAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCg                                                                                                                                >  4:39543/1‑149 (MQ=255)
                         ttCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATg                                                                                                                          <  6:42126/149‑1 (MQ=255)
                         ttCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATg                                                                                                                          >  7:22238/1‑149 (MQ=255)
                         ttCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATg                                                                                                                          >  3:20499/1‑149 (MQ=255)
                            cACCGGTGATGCAGCCTTCGGCCATGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGcc                                                                                                                       <  7:115006/149‑1 (MQ=255)
                             aCCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGCCCAGGTGGCCGCGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCt                                                                                                                      >  4:28037/1‑149 (MQ=255)
                                     tGCAGCCATCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCACGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGAcc                                                                                                                 >  4:19618/1‑146 (MQ=255)
                                        aGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTc                                                                                                            >  1:52860/1‑148 (MQ=255)
                                        aGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTc                                                                                                            >  6:258071/1‑148 (MQ=255)
                                                  cAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCg                                                                                                  >  8:57725/1‑148 (MQ=255)
                                                    tgTAATGTGCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAg                                                                                               <  1:257185/148‑1 (MQ=255)
                                                    ggTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAg                                                                                               >  3:214272/1‑149 (MQ=255)
                                                                cagtAGGAAGGCTTGGCACGGCGGCGGAGTTAAGGTGCGCGGTGAGGCAGCTG‑CCTTCGGCATGATAGACGATGTGGCCGAGCAGGCCGAAGCCTGTGAATTCGGTCATTGCATTGACCCCGTCCAGTTTGCCGAAGCGGCTTCCGGGc                                                                                   <  4:20499/146‑1 (MQ=255)
                                                                cacTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGCCGTCGGTCATGGCCTTTCCCCCGTCCACTTTTCCGAAGCGGCTGCCGGGc                                                                                   >  5:9892/1‑149 (MQ=255)
                                                                 acTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCg                                                                                  <  5:63452/149‑1 (MQ=255)
                                                                      gcaggcGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCCGTTTGCCGAAGCGGCTGCCGGGCGTGtt                                                                              >  1:252531/1‑148 (MQ=255)
                                                                       caggcGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGa                                                                            <  5:258071/149‑1 (MQ=255)
                                                                         ggcGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAgg                                                                          <  7:57725/149‑1 (MQ=255)
                                                                             tGGCACGGCTGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGc                                                                      <  8:22238/149‑1 (MQ=255)
                                                                             tGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGc                                                                      <  3:19618/149‑1 (MQ=255)
                                                                                  cggcggCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATc                                                                 <  4:214272/149‑1 (MQ=255)
                                                                                  cggcggCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCACGACCACGTGGCCGCGCCCCCCGAAGACCGGGGCCTAGCGGAGGGCCGAGGCCCCCTCCCCGTTGCCGACACGGCTGCCGGGCCGGTGGAGGGTGgaccag                                                                 >  2:124958/1‑146 (MQ=255)
                                                                                         gTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTcgc                                                          <  3:85308/149‑1 (MQ=255)
                                                                                          tAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCa                                                                                                          <  6:209645/100‑1 (MQ=255)
                                                                                          tAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCa                                                                                                          >  5:209645/1‑100 (MQ=255)
                                                                                            gTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAg                                                       <  3:39543/149‑1 (MQ=255)
                                                                                                  gggcgcGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGtt                                                                              <  5:254317/118‑1 (MQ=255)
                                                                                                  gggcgcGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGtt                                                                              >  6:254317/3‑120 (MQ=255)
                                                                                                       cgGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTgg                                            <  6:112733/149‑1 (MQ=255)
                                                                                                       cgGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTgg                                            <  4:255487/149‑1 (MQ=255)
                                                                                                            tggCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGGCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTg                                       <  1:124958/148‑1 (MQ=255)
                                                                                                                cGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGcaca                                   <  2:52860/149‑1 (MQ=255)
                                                                                                                            gCCTGCTCGGCCACCTGGTCGAGCTGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCAGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGc                        >  7:75211/1‑149 (MQ=255)
                                                                                                                             cctgCTCGGCCACCTGGTCGAGCTGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGcc                       >  2:266128/4‑149 (MQ=255)
                                                                                                                             ccAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGcc                       >  1:144884/1‑149 (MQ=255)
                                                                                                                             ccAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGcc                       <  3:163111/149‑1 (MQ=255)
                                                                                                                               tgCTCGGCCACCTGGTCGAGCTGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCctt                     >  8:6455/2‑149 (MQ=255)
                                                                                                                                         ggTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCgg             >  6:117885/1‑147 (MQ=255)
                                                                                                                                         ggTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCgg             >  6:113692/1‑147 (MQ=255)
                                                                                                                                              ccGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGG‑CCTGCACACGCAGCTTGGCCTTCTTTTCGGCAgtggt      >  3:63369/1‑148 (MQ=255)
                                                                                                                                                 aGCTGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGTGGTGAg   <  7:6455/149‑1 (MQ=255)
                                                                                                                                                 aGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTgg                                            <  4:158551/108‑1 (MQ=255)
                                                                                                                                                 aGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTgg                                            >  3:158551/1‑108 (MQ=255)
                                                                                                                                                  gCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGTGGTGAgg  <  3:28037/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGCGGTGAGGCCGCTG‑CCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAGTGGTGAGG  >  NC_002947/961573‑961866

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