Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,554,258	G→T	7.7%	S223* (TCG→TAG)	PP_4041 ←	glycoside hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,554,258	0	G	T	7.7%	161.3 / 2.5	52	S223* (TCG→TAG)	PP_4041	glycoside hydrolase
Reads supporting (aligned to +/- strand):  ref base G (23/25);  new base T (2/2);  total (25/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCAGGGTCATGGTCGCGTCGCGCAGCCAGCAGAAGCGGTAGTCCCAGTTGCGCTCGCCGCCCACACGCTCGGGCAGCGACGTGGTGACGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGTCGATGTTGCAGTCAGGCTGCACCGGCAAGTGCGAGGGCTGGTGACGCAGG  >  NC_002947/4554113‑4554398                                                                                                                                                   |                                                                                                                                              cagcagGGTCATGGTCGCGTCGCGCAGCCAGCAGAAGCGGTAGTCCCAGTTGCGCTCGCCGCCCACACGCTCGGGCAGCGACGTGGTGACGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCagcgagc                                                                                                                                            >  1:273952/1‑149 (MQ=255) cagcagGGTCATGGTCGCGTCGCGCAGCCAGCAGAAGCGGTAGTCCCAGTTGCGCTCGCCGCCCACACGCTCGGGCAGCGACGTGGTGACGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCagcgagc                                                                                                                                            >  3:153999/1‑149 (MQ=255)              gtcgcgtcgcGCAGCCAGCAGAAGCGGTAGTCCCAGTTGCGCTCGCCGCCCACACGCTCGGGCAGCGACGTGGTGACGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGc                                                                                                                              >  8:181584/1‑149 (MQ=255)                   gtcgcgCAGCCAGCAGAAGCGGTAGTCCCAGTTGCGCTCGCCGCCCACACGCTCGGGCTGCGACGTGGTGACGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGtat                                                                                                                         <  4:248197/149‑1 (MQ=255)                   gtcgcgCAGCCAGCAGAAGCGGTAGTCCCAGTTGCGCTCGCCGCCCACACGCTCGGGCAGCGACGTGGTGACGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCTCC‑AGGGCGGta                                                                                                                          >  7:139908/1‑148 (MQ=255)                   gtcgcgCAGCCAGCAGAAGCGGTAGTCCCAGTTGCGCTCGCCGCCCACACGCTCGGGCAGCGACGTGGTGACGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGtat                                                                                                                         <  8:66731/149‑1 (MQ=255)                      gcgcAGCCAGCAGAAGCGGTAGTCCCAGTTGCGCTCGCCGCCCACACGCTCGGGCAGCGACGTGGTGACGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGCGCCGCTCGCTGCTCAccctcaag                                                                                                                              >  4:173585/1‑136 (MQ=255)                      gcgcAGCCAGCAGAAGCGGTAGTCCCAGTTGCGCTCGCCGCCCACACGCTCGGGCAGCGACGTGGTGACGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGCGCCGCTCGCTGCTCAccctcaag                                                                                                                              <  3:173585/142‑7 (MQ=255)                         cagccagcAGAAGCGGTAGTCCCAGTTGCGCTCGCCGCCCACACGCTCGGGCAGCGACGTGGTGACGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTc                                                                                                                   <  4:210408/149‑1 (MQ=255)                         cagccagcAGAAGCGGTAGTCCCAGTTGCGCTCGCCGCCCACACGCTCGGGCAGCGACGTGGTGACGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTc                                                                                                                   <  4:206458/149‑1 (MQ=255)                             cagcagAAGCGGTAGTCCCAGTTGCGCTCGCCGCCCACACGCTCGGGCAGCGACGTGGTGACGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTAc                                                                                                               >  6:310922/1‑149 (MQ=255)                                cagAAGCGGTAGTCCCAGTTGCGCTCGCCGCCCACACGCTCGGGCAGCGACGTGGTGACGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTg                                                                                                            >  5:17891/1‑149 (MQ=255)                                    aGCGGTAGTCCCAGTTGCGCTCGCCGCCCACACGCTCGGGCAGCGACGTGGTGACGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCggg                                                                                                        <  7:181584/149‑1 (MQ=255)                                              ccAGTTGCGCTCGCCGCCCACACGCTCGGGCAGCGACGTGGTGACGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCgg                                                                                              >  3:101834/1‑149 (MQ=255)                                              ccAGTTGCGCTCGCCGCCCACACGCTCGGGCAGCGACGTGGTGACGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCgg                                                                                              >  3:127390/1‑149 (MQ=255)                                              ccAGTTGCGCTCGCCGCCCACACGCTCGGGCAGCGACGTGGTGACGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCgg                                                                                              >  7:87745/1‑149 (MQ=255)                                                aGTTGCGCTCGCCGCCCACACGCTCGGGCAGCGACGTGGTGACGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGc                                                                                             >  4:122482/1‑148 (MQ=255)                                                    gcgcTCGCCGCCCACACGCTCGGGCAGCGACGTGGTGACGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAgg                                                                                        >  2:182139/1‑149 (MQ=255)                                                       gtCGCCGGCCACCCGCTCGGGCAGCGACGTGGTGACCGCAGCGACAATTCC‑GCCGGCGGGCTCGTAGGTCATTGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCTGTCGGCGAAGGCCt                                                                                     <  4:37980/148‑1 (MQ=255)                                                             gccCACACGCTCGGGCAGCGACGTGGTGACGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCAtt                                                                               <  4:252688/149‑1 (MQ=255)                                                                 acacGCTCGGGCAGCGACGTGGTGACGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTc                                                                           <  6:259318/149‑1 (MQ=255)                                                                 acacGCTCGGGCAGCGACGTGGTGACGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTc                                                                           <  4:50063/149‑1 (MQ=255)                                                                           gCAGCGACGTGGTGACGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGttt                                                                 <  2:273952/149‑1 (MQ=255)                                                                                gACGTGGTGACGGCAGCGACAATACCGGCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTAc                                                                                                               >  1:160905/1‑99 (MQ=255)                                                                                   gtggtgACGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCTAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGc                                                         >  6:144247/1‑149 (MQ=255)                                                                                    tggtgACGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGAGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTCTTTGCCAGGGCt                                                        >  2:242289/1‑149 (MQ=255)                                                                                         aCGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGt                                                    <  1:182139/148‑1 (MQ=255)                                                                                         aCGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGt                                                    <  4:101834/148‑1 (MQ=255)                                                                                         aCGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGt                                                    <  3:175612/148‑1 (MQ=255)                                                                                           ggCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGTCGa                                                 >  3:70686/1‑149 (MQ=255)                                                                                           ggCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGTCGa                                                 >  3:198419/1‑149 (MQ=255)                                                                                            gCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGTCGAt                                                <  3:122482/149‑1 (MQ=255)                                                                                                cGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGTCGATGTTg                                            >  8:180181/1‑149 (MQ=255)                                                                                                         c‑gccgGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGTCGATg                                               >  8:279196/1‑137 (MQ=255)                                                                                                         c‑gccgGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGTCGATg                                               <  7:279196/137‑1 (MQ=255)                                                                                                         c‑gccgGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGTCGATGTTGCAGTCAGGc                                   <  5:310922/149‑1 (MQ=255)                                                                                                         c‑gccgGTGGGAGAGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGTCGATGTTGCAGTCAGGc                                   >  2:77784/1‑149 (MQ=255)                                                                                                            ccgGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTc                                                                                                                   >  4:245198/1‑67 (MQ=255)                                                                                                            ccgGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTc                                                                                                                   <  3:245198/67‑1 (MQ=255)                                                                                                             aggTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGTCGATGTTGCAGTCAGGCTGc                                <  4:103731/148‑1 (MQ=255)                                                                                                                  ggAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGTCGATGTTGCAGTCAGGCTGCACCg                            <  1:304423/148‑1 (MQ=255)                                                                                                                        tAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGTCGATGTTGCAGTCAGGCTGCACCGGCAAGTg                     <  1:77784/149‑1 (MQ=255)                                                                                                                         aGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGTCGATGTTGCAGTCAGGCTGCACCGGCAAGTGc                    >  8:67106/1‑149 (MQ=255)                                                                                                                           gTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACCGAGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGTCGATGTTGCAGTCAGGCTGCACCGGCAAGTGc                    >  1:49688/1‑148 (MQ=255)                                                                                                                             cATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGTCGATGTTGCAGTCAGGCTGCACCGGCAAGCGCGAgg                >  6:238229/1‑149 (MQ=255)                                                                                                                              aTGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGTCGATGTTGCAGTCAGGCTGCACCGGCAAGTGCGAggg               >  8:205119/1‑149 (MQ=255)                                                                                                                              aTGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGTCGATGTTGCAGTCAGGCTGCACCGGCAAGTGCGAgg                <  6:17891/148‑1 (MQ=255)                                                                                                                              aTGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGTCGATGTTGCAGTCAGGCTGCACCGGCAAGTGCGAgg                <  4:70686/148‑1 (MQ=255)                                                                                                                                  ccTTGAGGGTGAGCAGCTAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGTCGATGTTGCAGTCAGGCTGCACCGGCAAGTGCGAGGGCTgg           <  5:144247/149‑1 (MQ=255)                                                                                                                                   cTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGTCGATGTTGCAGTCAGGCTGCACCGGCAAGTGCGAGGGCTGGt          <  4:153999/149‑1 (MQ=255)                                                                                                                                    ttGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGTCGATGTTGCAGTCAGGCTGCACCGGCAAGTGCGAg                 >  2:96968/1‑141 (MQ=255)                                                                                                                                    ttGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCAATGCTCGCAAGTGTTTGCCAGGGCTTGGTCGATGTTGCAGTCAGGCTGCACCGGCAAGTGCGAGGGCTGGt          >  1:96609/1‑148 (MQ=255)                                                                                                                                    ttGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGCCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGTCGATGTTGCAGTCAGGCTGCACCGGCAAGTGCGAg                 <  1:96968/141‑1 (MQ=255)                                                                                                                                           tGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGTCGATGTTGCAGTCAGGCTGCACCGGCAAGTGCGAGGGCTGGTGACGCAgg  <  2:49688/149‑1 (MQ=255)                                                                                                                                                agcgagcgGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGTCGATGTTGCAGTCAg                                     <  7:184720/109‑1 (MQ=255)                                                                                                                                                agcgagcgGCGAACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGTCGATGTTGCAGTCAg                                     >  8:184720/1‑109 (MQ=255)                                                                                                                                                   |                                                                                                                                              CAGCAGGGTCATGGTCGCGTCGCGCAGCCAGCAGAAGCGGTAGTCCCAGTTGCGCTCGCCGCCCACACGCTCGGGCAGCGACGTGGTGACGGCAGCGACAATACC‑GCCGGTGGGAGCGTAGGTCATGGCCTTGAGGGTGAGCAGCGAGCGGCGCACC‑AGGGCGGTATACGGTCCTACCTGCGGGCAGCGGTCGGCGAAGGCCTGCCATTGCTCGCAAGTGTTTGCCAGGGCTTGGTCGATGTTGCAGTCAGGCTGCACCGGCAAGTGCGAGGGCTGGTGACGCAGG  >  NC_002947/4554113‑4554398

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