Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,971,956	A→C	8.8%	intergenic (‑391/+26)	PP_4379 ← / ← flgL	3‑oxoacyl‑ACP synthase/flagellar hook‑associated protein FlgL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,971,956	0	A	C	8.8%	171.8 / 7.6	57	intergenic (‑391/+26)	PP_4379/flgL	3‑oxoacyl‑ACP synthase/flagellar hook‑associated protein FlgL
Reads supporting (aligned to +/- strand):  ref base A (25/27);  new base C (3/2);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

GCACTGGACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTA  >  NC_002947/4971809‑4972093                                                                                                                                                    |                                                                                                                                           gcACTGGACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCgcag                                                                                                                                           >  6:263851/1‑149 (MQ=255) gcACTGGACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCgcag                                                                                                                                           <  5:263851/149‑1 (MQ=255)       gACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCTGCGCAG‑ccc                                                                                                                                    >  8:71839/1‑149 (MQ=255)         ccccGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑ccc                                                                                                                                    >  7:40181/1‑147 (MQ=255)          cccGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAg                                                                                                                                 >  4:211568/1‑149 (MQ=255)           ccGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCTGCGCAG‑ccctatc                                                                                                                                <  7:71839/149‑5 (MQ=255)           ccGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAg                                                                                                                                 >  3:99610/1‑148 (MQ=255)            cGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCTGCGCAG‑ccctatcg                                                                                                                               <  8:239704/149‑6 (MQ=255)             gAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGAt                                                                                                                              <  7:168401/149‑1 (MQ=255)              aaTGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATc                                                                                                                             <  5:356372/149‑1 (MQ=255)                tGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACAACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATcag                                                                                                                           <  6:255305/149‑1 (MQ=255)                   tCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATcaggca                                                                                                                        <  5:3511/149‑1 (MQ=255)                     aGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATcaggcagg                                                                                                                      >  5:50342/1‑149 (MQ=255)                     aGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATcaggcag                                                                                                                       <  3:211568/148‑1 (MQ=255)                      gCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATcaggcag                                                                                                                       >  2:23018/1‑147 (MQ=255)                      gCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATcaggcag                                                                                                                       <  1:23018/147‑1 (MQ=255)                           gcagcaGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGAt                                                                                                                 >  6:2024/1‑148 (MQ=255)                           gcagcaGATATCGATGGGGCCGCATCACGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGAt                                                                                                                 <  2:324829/148‑1 (MQ=255)                              gcagATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGa                                                                                                             <  8:40181/149‑1 (MQ=255)                                     tCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGt                                                                                                       <  5:2024/148‑1 (MQ=255)                                               cGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCa                                                                                                                                  <  2:15678/111‑1 (MQ=255)                                               cGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCt                                                                                             >  7:200885/1‑148 (MQ=255)                                               cGCATCGCGGCCCTTTCGCGACGCAAGGACGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCa                                                                                                                                  >  1:15678/1‑111 (MQ=255)                                                gCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTa                                                                                           <  3:70496/149‑1 (MQ=255)                                                 cATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAg                                                                                          <  4:99610/149‑1 (MQ=255)                                                     gcgGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAGCCCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAg                                                                                          <  1:248474/146‑1 (MQ=255)                                                     gcgGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAg                                                                                          >  2:248474/1‑145 (MQ=255)                                                         cccTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGaaa                                                                                  >  6:175276/1‑149 (MQ=255)                                                           cTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATca                                                                                                                            >  4:55387/1‑105 (MQ=255)                                                           cTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATca                                                                                                                            <  3:55387/105‑1 (MQ=255)                                                           cTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATg                                                                                >  5:182102/1‑149 (MQ=255)                                                               cgcgACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGTGCTGCGCAGCAG‑CCCCAGGATca                                                                                                                            <  5:187700/101‑1 (MQ=255)                                                               cgcgACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGTGCTGCGCAGCAG‑CCCCAGGATca                                                                                                                            >  6:187700/1‑101 (MQ=255)                                                                gcgACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATgcgcg                                                                            >  1:147793/1‑148 (MQ=255)                                                                        aGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGc                                                                     >  6:295793/1‑147 (MQ=255)                                                                                cccACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGaa                                                           <  6:182102/149‑1 (MQ=255)                                                                                 ccACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCg                                                             <  2:74770/146‑1 (MQ=255)                                                                                 ccACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCg                                                             >  1:74770/1‑146 (MQ=255)                                                                                  cACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGTCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAgg                                                                                               <  5:216631/111‑1 (MQ=255)                                                                                  cACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGTCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAgg                                                                                               >  6:216631/1‑111 (MQ=255)                                                                                   aCGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATgcgcg                                                                            >  6:275952/1‑129 (MQ=255)                                                                                   aCGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATgcgcg                                                                            <  5:275952/129‑1 (MQ=255)                                                                                           gCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCAt                                                <  7:194187/149‑1 (MQ=255)                                                                                                cTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCa                                                                   >  4:149532/1‑125 (MQ=255)                                                                                                cTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCa                                                                   <  3:149532/125‑1 (MQ=255)                                                                                                   cAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCctg                                          <  5:175276/147‑1 (MQ=255)                                                                                                                  ccTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGaa                         <  7:319218/149‑1 (MQ=255)                                                                                                                  ccTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGaa                         >  8:18980/1‑149 (MQ=255)                                                                                                                   cTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCaa                                                    >  8:73608/1‑121 (MQ=255)                                                                                                                   cTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCaa                                                    <  7:73608/121‑1 (MQ=255)                                                                                                                     gtgGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAg                                                                                                              >  4:63537/1‑61 (MQ=255)                                                                                                                     gtgGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAg                                                                                                              <  3:63537/61‑1 (MQ=255)                                                                                                                             cTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAg              >  6:111856/1‑149 (MQ=255)                                                                                                                             cTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAg              >  2:314002/1‑149 (MQ=255)                                                                                                                                 tgtCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGcc           <  6:50342/148‑1 (MQ=255)                                                                                                                                      cgATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCa     <  7:18980/149‑1 (MQ=255)                                                                                                                                         tggggCTGCTGCGCAG‑CCCTAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTa  >  7:32222/3‑149 (MQ=255)                                                                                                                                         tggggCTGCTGCGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTa  >  5:66232/3‑149 (MQ=255)                                                                                                                                               tGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCctg                                          >  2:48259/1‑103 (MQ=255)                                                                                                                                                    |                                                                                                                                           GCACTGGACCCCGAATGATCAGCTTCGCAGCAGATATCGATGGGGCCGCATCGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACGAATACGCCACCTCCAGGAGAGAGGACGACCTGTGGGGGCCTTGTGTCGCGATAGGGCTGCGCAGCAG‑CCCCAGGATCAGGCAGGCAGATTAGATCTTGTCGAACAGGCTTAGCTGAGAAATGCGCGCAAACGCCAGTTGCGAAGCTTCAAGCATGGCCTGCTGCAGGGTCAGCTGAATGGTGGCCGAGGCCATGTCAGTA  >  NC_002947/4971809‑4972093

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