Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,430,396	A→G	9.5%	intergenic (‑78/+83)	PP_4769 ← / ← PP_4770	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,430,396	0	A	G	9.5%	138.5 / 5.5	53	intergenic (‑78/+83)	PP_4769/PP_4770	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/26);  new base G (3/2);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.51e-01

GGCAGCCGTGGGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGA‑TCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCGGCGTA  >  NC_002947/5430249‑5430544                                                                                                                                                     |                                                                                                                                                       ggCAGCCGTGGGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGA‑TCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAt                                                                                                                                                       >  7:252116/1‑149 (MQ=255)          ggtccTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCTCGA‑TCTTTTTGACCTGGGAACCTCATTTTTTGCCTGGAGTGCCGTTGGGCTGCTGCTCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGt                                                                                                                                               <  4:248543/144‑1 (MQ=255)                  cgcACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGA‑TCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTgg                                                                                                                                     <  1:87880/149‑1 (MQ=255)                  cgaaCGGTGTAGTAGAGGTTCTGCAAATCGGCGAACAGCGCGA‑TCTTTTTCCACTGGGCACCTCATTTTTGGCCTGGAGTGCC‑GGTGGC‑GCTTTGCGCCCCTTCGCGGGCACGCCCGCTCCCAAGGGTGTACGCGCG‑ATACCTGTGGGa                                                                                                                                   <  7:48617/146‑1 (MQ=9)                        gtgtAGTAGAGGTTCTGCACATCGGCGAACAGCGCGA‑TCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCgg                                                                                                                               <  7:256116/149‑1 (MQ=255)                           tagtagAGGTTCTGCACATCGGCGAACAGCGCGA‑TCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGc                                                                                                                             <  8:39489/148‑1 (MQ=255)                             gtagAGGTTCTGCACATCGGCGAACAGCGCGA‑TCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACGGGTGTACGCGCG‑ATACCTGTGGGAg                                                                                                                                  <  4:105916/141‑1 (MQ=37)                             gtagAGGTTCTGCACATCGGCGAACAGCGCGA‑TCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACGGGTGTACGCGCG‑ATACCTGTGGGAg                                                                                                                                  >  3:105916/1‑141 (MQ=37)                             gtagAGGTTCTGCACATCGGCGAACAGCGCGA‑TCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTg                                                                                                                          >  3:266625/1‑149 (MQ=255)                              tagAGGTTCTGCACATCGGCGAACAGCGCGATTCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTg                                                                                                                          >  1:46479/1‑149 (MQ=255)                              tagAGGTTCTGCACATCGGCGAACAGCGCGA‑TCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGc                                                                                                                         <  1:139047/149‑1 (MQ=255)                                        gCACATCGGCGAACAGCGCGA‑TCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAggg                                                                                                              >  5:117929/1‑149 (MQ=255)                                         cacaTCGGCGAACAGCGCGA‑TCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAgggg                                                                                                             <  3:351767/149‑1 (MQ=255)                                         cacaTCGGCGAACAGCGCGA‑TCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAg                                                                                                                <  5:388607/146‑1 (MQ=255)                                         cacaTCGGCGAACAGCGCGA‑TCTTTTTCACCTGGGAACCTCATTTTTGGCCGGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAggg                                                                                                              <  4:222475/148‑1 (MQ=255)                                          acaTCGGCGAACAGCGCGA‑TCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGgc                                                                                                            <  3:156162/149‑1 (MQ=255)                                           caTCGGCGAACAGCGCGA‑TCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAg                                                                                                                >  1:110533/1‑144 (MQ=255)                                           caTCGGCGAACAGCGCGA‑TCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAg                                                                                                                <  2:110533/144‑1 (MQ=255)                                           caTCGGCGAACAGCGCGA‑TCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCCTCGCGGGCACCCCCCCGCCCCCCGGGAT‑AGCGGGGGACCCCGGGGGCGGGGGGGTGCCC‑GCGCAgggggg                                                                                                           >  6:207669/1‑147 (MQ=255)                                               ggCGAACAGCGCGA‑TCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCaaa                                                                                                       <  2:88287/149‑1 (MQ=255)                                                gCGAACAGCGCGA‑TCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAg                                                                                                      >  7:54752/1‑149 (MQ=255)                                                             attttttcCCCCGGGGACCCCCATTTTGCGCCTGGATGGCCGGGGCGCGGCGCCCCCGACCTCCGCGGACACCCCCGCTCCACCCCGTCT‑GGCGGGTACCCCGGTGGGTGGGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGaaa                                                                                          <  4:198764/146‑1 (MQ=255)                                                                  tttCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGTTAGTGGGCGTGCCCTGCGAAGGGGCGCAAAGAGCCCCCTGAAATGc                                                                                       <  1:41680/148‑1 (MQ=255)                                                                  tttCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGc                                                                                       >  3:207303/1‑147 (MQ=255)                                                                  tttCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCCTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAAt                                                                                         <  5:207669/145‑1 (MQ=255)                                                                   ttCTCCTGGGAAGCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCca                                                                                     <  4:356660/148‑1 (MQ=255)                                                                        cTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTa                                                                               >  1:96488/1‑149 (MQ=255)                                                                           ggAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCTCAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCG‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTAt                                                                              <  3:34054/147‑1 (MQ=255)                                                                             aaCCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCa                                                                          <  4:266625/149‑1 (MQ=255)                                                                                cTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGCAGGGGCGCAACGCGCCCCCGGAAATGCCAGCAGTATGCCAGAc                                                                       >  8:154969/1‑149 (MQ=255)                                                                                cTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGAc                                                                       <  5:113980/149‑1 (MQ=255)                                                                                cTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACACCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGAc                                                                       >  3:135008/1‑149 (MQ=255)                                                                                cTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCGGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGAc                                                                       >  6:370457/1‑149 (MQ=255)                                                                                 tCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGAcc                                                                      >  5:163332/1‑149 (MQ=255)                                                                                 tCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCCCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGAc                                                                       >  1:115476/1‑148 (MQ=255)                                                                                    tttttGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGt                                                                   <  2:238120/149‑1 (MQ=255)                                                                                        tgtccTGGAGTGCCGGGGGGCGGCTGCGCAGACCTTCGCGGACACGGCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAg                                                                <  6:230030/145‑1 (MQ=255)                                                                                         ggCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGa                                                               <  4:207303/148‑1 (MQ=255)                                                                                           ccTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACg                                                             >  5:367045/1‑148 (MQ=255)                                                                                                     cGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAAtcgtc                                                  <  8:252116/149‑1 (MQ=255)                                                                                                     cGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAAtcgtc                                                  <  3:79088/149‑1 (MQ=255)                                                                                                                         cTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCga                                            <  4:146850/135‑1 (MQ=255)                                                                                                                         cTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCga                                            >  3:146850/1‑135 (MQ=255)                                                                                                                              cgGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGGATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCg                         >  5:190469/1‑149 (MQ=255)                                                                                                                                 acacGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCTGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTg                      >  8:103695/1‑149 (MQ=255)                                                                                                                                  cacGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGc                     <  8:54752/149‑1 (MQ=255)                                                                                                                                   acGCCCGCTCCCACGGGTGTACGCGCG‑ATACCTGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGc                     >  3:370577/1‑148 (MQ=37)                                                                                                                                       ccGCTCCCACGGGTGTACGCGCG‑ATACCTGTGGGAGCGGGCGTGCCC‑GCGACGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCaaaa                >  7:248244/1‑149 (MQ=37)                                                                                                                                             ccACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCg          >  5:73562/1‑149 (MQ=255)                                                                                                                                               acaGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTc        >  7:110340/1‑149 (MQ=255)                                                                                                                                                 aGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCgg      >  2:161321/1‑149 (MQ=255)                                                                                                                                                  ggTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCGGc     <  6:367045/149‑1 (MQ=255)                                                                                                                                                  ggTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCGGc     <  2:115476/149‑1 (MQ=255)                                                                                                                                                     aT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCGGCGTa  >  2:244171/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                       GGCAGCCGTGGGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGA‑TCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCC‑GCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCGGCGTA  >  NC_002947/5430249‑5430544

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑