Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,787,636	T→C	11.3%	intergenic (+83/‑42)	ahpC → / → ahpF	peroxiredoxin/alkylhydroperoxide reductase small subunit/alkyl hydroperoxide reductase subunit F

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,787,636	0	T	C	11.3%	135.3 / 7.1	54	intergenic (+83/‑42)	ahpC/ahpF	peroxiredoxin/alkylhydroperoxide reductase small subunit/alkyl hydroperoxide reductase subunit F
Reads supporting (aligned to +/- strand):  ref base T (23/24);  new base C (2/4);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGCGCGAAGTCCCGC  >  NC_002947/2787488‑2787779                                                                                                                                                      |                                                                                                                                                  ccGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCt                                                                                                                                                   <  6:6265/149‑1 (MQ=255)  cGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCAGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTg                                                                                                                                                  >  1:207026/1‑149 (MQ=255)    gCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTggg                                                                                                                                                >  2:236898/1‑149 (MQ=255)    gCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCAGGCGAGGTCGCCCggg                                                                                                                                                <  6:219122/149‑1 (MQ=255)        aGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCAGGCGAGGTCGCCCGGGCGtt                                                                                                                                            <  4:17776/149‑1 (MQ=255)          tGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGtttt                                                                                                                                          >  2:262693/1‑149 (MQ=255)            gAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCTAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGttttt                                                                                                                                         >  1:78818/1‑149 (MQ=255)                   ggCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGa                                                                                                                                 <  3:218406/149‑1 (MQ=255)                    gCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAg                                                                                                                                >  1:340659/1‑149 (MQ=255)                                tGGCGCCGTCGCTGGACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGaaaaa                                                                                                                     <  3:355564/147‑1 (MQ=255)                                              gACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAg                                                                                                                    <  4:88578/134‑1 (MQ=255)                                              gACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAg                                                                                                                    >  3:88578/1‑134 (MQ=255)                                              gACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGt                                                                                                     >  7:33352/1‑149 (MQ=255)                                               aCCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCATCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGtt                                                                                                    <  8:75374/149‑1 (MQ=255)                                                 cTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTg                                                                                                   >  4:46511/1‑148 (MQ=255)                                                         c‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACg                                                                                          >  2:232170/1‑149 (MQ=255)                                                           aaGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGTCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGc                                                                                         <  2:207026/149‑1 (MQ=255)                                                                cTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTaaa                                                                                    >  8:236394/1‑149 (MQ=255)                                                                   aGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCg                                                                                 <  7:236394/149‑1 (MQ=255)                                                                    ggTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCg                                                                                 >  5:195785/1‑148 (MQ=255)                                                                      tCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCaa                                                                              >  5:354078/1‑149 (MQ=255)                                                                          gCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑g                                                                          >  4:3971/1‑148 (MQ=255)                                                                            ccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑Gaa                                                                        >  1:88120/1‑148 (MQ=255)                                                                                  cgcgGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTAc                                                                 <  8:33352/149‑1 (MQ=255)                                                                                   gcgggcggTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACt                                                                >  7:178545/1‑149 (MQ=255)                                                                                    cgggcggTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACtt                                                               <  2:78818/149‑1 (MQ=255)                                                                                     gggcggTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTg                                                              <  3:46511/149‑1 (MQ=255)                                                                                           tAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCgg                                                        <  6:354078/149‑1 (MQ=255)                                                                                           tAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCgg                                                        <  1:262693/149‑1 (MQ=255)                                                                                             agcagcCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGt                                                      >  1:185797/1‑149 (MQ=255)                                                                                                     ccAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGc                                              <  1:23384/149‑1 (MQ=255)                                                                                                          gCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATc                                         <  3:3971/149‑1 (MQ=255)                                                                                                            tGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCga                                       <  6:195785/149‑1 (MQ=255)                                                                                                             gAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATc                                         >  3:23393/1‑146 (MQ=255)                                                                                                              aaGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTCGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCg                                        <  1:236898/147‑1 (MQ=255)                                                                                                                   cctaccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGtt                                <  1:232170/149‑1 (MQ=255)                                                                                                                    ctaccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGc                                                          <  1:170678/122‑1 (MQ=255)                                                                                                                    ctaccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGc                                                          >  2:170678/1‑122 (MQ=255)                                                                                                                      accGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCGCGTATGTTGGACGCCACGCTTAATTCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCt                             >  7:157947/1‑149 (MQ=255)                                                                                                                        cGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATc                                   <  7:319678/141‑1 (MQ=255)                                                                                                                        cGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTAGGAGCGGGTCACCCAGCCGATCGAGATc                                   >  8:319678/1‑141 (MQ=255)                                                                                                                          ccccGACAAACCCCCCCCCGACCTCCCCAGCGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTccc                         <  2:140897/149‑1 (MQ=255)                                                                                                                           cccGTCAACCGGCCGGTCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCt                        <  8:178545/149‑1 (MQ=255)                                                                                                                            gcGTCAAACGGCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTAGTTGGAGCGGGTCAGCCAGCCGATCGAGATTGTTGCTTCCCtg                       <  2:340659/148‑2 (MQ=255)                                                                                                                               aaaaaCGCCCAGGCGAGGTCGCCCGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAg                                               <  3:263136/120‑1 (MQ=39)                                                                                                                               aaaaaCGCCCAGGCGAGGTCGCCCGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAg                                               >  4:263136/3‑122 (MQ=39)                                                                                                                                 aaaCGCCCAGGCGAGGTCGCCCGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTcgacga                  >  7:262269/1‑149 (MQ=255)                                                                                                                                     gCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGgc              <  2:88120/149‑1 (MQ=255)                                                                                                                                          ggCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGCGCGaa         <  4:23393/149‑1 (MQ=255)                                                                                                                                              aCCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGCGCGaa         >  7:82965/1‑145 (MQ=255)                                                                                                                                              aCCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGCGCGaa         <  8:82965/145‑1 (MQ=255)                                                                                                                                                 tCGCCTGGGCGTTTTTGTATCCGAGTTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGCGCGAAGTCCcg   >  1:27127/1‑149 (MQ=255)                                                                                                                                                 tCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGCGCGAAGTCCcgc  >  1:243884/1‑149 (MQ=255)                                                                                                                                                 tCGCCCGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGCGCGAAGTCCcgc  <  3:100465/149‑1 (MQ=255)                                                                                                                                                   gCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGCGCGAAGTCCcgc  <  2:27127/147‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  CCGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAG‑TTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGCGCGAAGTCCCGC  >  NC_002947/2787488‑2787779

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