Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,648,455	A→T	6.7%	I247F (ATC→TTC)	PP_3214 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,648,455	0	A	T	6.7%	181.1 / 2.7	60	I247F (ATC→TTC)	PP_3214	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/29);  new base T (2/2);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

GTGGGGACGAAGGGTACCTCGAACTATGTAGCCAAGAGTTTCCACTTCAACCCCCAGATCAAGTACGACGTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTCTAATGTGATGTTTCAGACGTACCACTGAGTATCTGGTACAGGCTCCGTTGATCATCACTCATCAAGCGTTGTGTAT  >  NC_002947/3648308‑3648596                                                                                                                                                     |                                                                                                                                              gtgGGGACGAAGGGTACCTCGAACTATGTAGCCAAGAGTTTCCACTTCAACCCCCAGATCAAGTACGACGTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCa                                                                                                                                               <  2:319991/148‑1 (MQ=255) gtgGGGACGAAGGGTACCTCGAACTATGTAGCCAAGAGTTTCCACTTCAACCCCCAGATCAAGTACGACGTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCAt                                                                                                                                              >  6:211923/1‑149 (MQ=255)         gAAGGGTACCTCGAACTATGTAGCCAAGAGTTTCCACTTCAACCCCCAGATCAAGTACGACGTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACa                                                                                                                                      <  3:282185/149‑1 (MQ=255)         gAAGGGTACCTCGAACTATGTAGCCAAGAGTTTCCACTTCAACCCCCAGATCAAGTACGACGTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCAGACAAATATGGCATCGAAGACa                                                                                                                                      >  8:286692/1‑149 (MQ=255)                 ccTCGAACTATGTAGCCAAGAGTTTCCACTTCAACCCCCAGATCAAGTACGACGTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGc                                                                                                                              <  4:287467/149‑1 (MQ=255)                  cTCGAACTATGTAGCCAAGAGTTTCCACTTCAACCCCCAGATCAAGTACGACGTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCt                                                                                                                             >  2:289172/1‑149 (MQ=255)                  cTCGAACTATGTAGCCAAGAGTTTCCACTTCAACCCCCAGATCAAGTACGACGTCGGTAAGGCGAT‑GGAGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCt                                                                                                                             >  1:300402/1‑149 (MQ=255)                              aGCCAAGAGTTTCCACTTCAACCCCCAGATCAAGTACGACGTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTccgac                                                                                                                      >  1:223032/1‑144 (MQ=255)                              aGCCAAGAGTTTCCACTTCAACCCCCAGATCAAGTACGACGTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTccgac                                                                                                                      <  2:223032/144‑1 (MQ=255)                                  aaGAGTTTCCACTTCAACCCCCAGATCAAGTACGACGTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGaa                                                                                                             >  1:264624/1‑149 (MQ=255)                                    gagTTTCCACTTCAACCCCCAGATCAAGTACGACGTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAACACAGCCACGGCTTTCCGACCGATCAGAAcg                                                                                                           <  1:29124/149‑1 (MQ=255)                                         tCCACTTCAACCCCCAGATCAAGTACGACGTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGtt                                                                                                       <  2:69066/148‑1 (MQ=255)                                         tCCACTTCAACCCCCAGATCAAGTACGACGTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTAGTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGtt                                                                                                       <  4:308942/148‑1 (MQ=255)                                          ccACTTCAACCCCCAGATCAAGTACGACGTCGGTAAGGCGATGGGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTccgaccg                                                                                                                    >  1:245387/1‑135 (MQ=255)                                          ccACTTCAACCCCCAGATCAAGTACGACGTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTccgaccg                                                                                                                    <  2:245387/134‑1 (MQ=255)                                           cACTTCAACCCCCAGATCAAGTACGACGTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGa                                                                                                     <  1:289172/148‑1 (MQ=255)                                              ttCAACCCCCAGATCAAGTACGACGTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGtt                                                                                                       >  5:101755/1‑143 (MQ=255)                                              ttCAACCCCCAGATCAAGTACGACGTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGtt                                                                                                       <  6:101755/143‑1 (MQ=255)                                                  aCCCCCAGATCAAGTACGACGTCGGTAAGGCGAT‑GGGGCATGCGCCGGGGCGTTTCTATGCGGGGATCGAGTGGGACTATTGGTCCGACAATTATGGCATGGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTg                                                                                             <  2:151560/149‑1 (MQ=255)                                                   ccccaAGATCAAGTACGACGTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTg                                                                                             >  3:67957/1‑148 (MQ=255)                                                    ccccAGATCAAGTACGACGTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGAt                                                                                           >  5:252020/1‑149 (MQ=255)                                                     cccAGATCAAGTACGACGTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGAt                                                                                           >  7:206387/1‑148 (MQ=255)                                                       cAGATCAAGTACGACGTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCaa                                                                                        <  7:286692/149‑1 (MQ=255)                                                        aGATCAAGTACGACGTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGt                                                                                                        >  8:135542/1‑132 (MQ=255)                                                        aGATCAAGTACGACGTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGt                                                                                                        <  7:135542/132‑1 (MQ=255)                                                             aaGTACGACGTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTcc                                                                                                                         <  6:176878/110‑1 (MQ=255)                                                             aaGTACGACGTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTcc                                                                                                                         >  5:176878/1‑110 (MQ=255)                                                               gTACGACGTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGt                                                                                     <  4:67957/144‑1 (MQ=255)                                                                 acgacgTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCAt                                                                                 >  6:42481/1‑146 (MQ=255)                                                                 acgacgTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCAt                                                                                 <  5:42481/146‑1 (MQ=255)                                                                  cgacgTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTc                                                                              >  8:286364/1‑148 (MQ=255)                                                                     cgTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTCTAAt                                                                          <  2:264624/149‑1 (MQ=255)                                                                      gTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTccgac                                                                                                                      <  7:214640/104‑1 (MQ=255)                                                                      gTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTccgac                                                                                                                      >  8:214640/1‑104 (MQ=255)                                                                          gTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTCTAATGTGAt                                                                     >  6:65183/1‑149 (MQ=255)                                                                           tAAGGCGAT‑GGCGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTCTAATGTGATg                                                                    <  2:300402/149‑1 (MQ=255)                                                                                        gCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGTGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTCTAATGTGATGTTTCAGACGTAc                                                        <  5:65183/149‑1 (MQ=255)                                                                                        gCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCa                                                                                         >  8:261443/1‑116 (MQ=255)                                                                                        gCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTCTAATGTGATGTTTCAGACGTAc                                                        <  8:206387/149‑1 (MQ=255)                                                                                        gCAAGCGCCGGCGCGTTTCTATGCCGGGATAGAGTGGGACTATTGGTCCGACAAATATGGCATCGAACACAGCCACCGCTTTCCGACCGATCAGAACGCGTTCACTGCAATCATCAAGGTGCATTTCTAATGTGATGTTTCAGACGTAc                                                        >  2:303712/1‑149 (MQ=255)                                                                                        gCAAGCGCCGGCGCGTTGCTATGCCGGGATCGAGTGGGACTATTGGTCCGAGAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCa                                                                                         <  7:261443/116‑1 (MQ=255)                                                                                                   cgcgTTTCTATGCCGGGATCGAGTGGGACTATTGGTACGACAAATTTGTCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTCTAATGTGATGTTTCAGACGTACCACTGAGTATc                                             <  1:303712/149‑1 (MQ=255)                                                                                                             tGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTCTAATGTGATGTTTCAGACGTACCACTGAGTATCTGGTACAg                                     <  6:252020/147‑1 (MQ=255)                                                                                                              gCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTCTAATGTGATGTTTCAGACGTACCACTGAGTATCTGGTACAGGCt                                  >  6:339589/1‑149 (MQ=255)                                                                                                                  ggATCGAGTGGGACTATTGGTCCGACAAATATGGCTTCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTCTAATGTGATGTTTCAGACGTACCACTGAGTATCTGGTa                                        >  5:236957/1‑139 (MQ=255)                                                                                                                  ggATCGAGTGGGACTATTGGTCCGACAAATATGGCTTCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTCTAATGTGATGTTTCAGACGTACCACTGAGTATCTGGTa                                        >  5:231813/1‑139 (MQ=255)                                                                                                                  ggATCGAGTGGGACTATTGGTCCGACAAATATGGCTTCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTCTAATGTGATGTTTCAGACGTACCACTGAGTATCTGGTa                                        <  6:231813/139‑1 (MQ=255)                                                                                                                  ggATCGAGTGGGACTATTGGTCCGACAAATATGGCTTCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTCTAATGTGATGTTTCAGACGTACCACTGAGTATCTGGTa                                        <  6:236957/139‑1 (MQ=255)                                                                                                                     tCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTCTAATGTGATGTTTCAGACGTACCACTGAGTATCTGGTACAGGCTCCGTTGa                           <  7:286364/149‑1 (MQ=255)                                                                                                                        aGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTCTAATGTGATGttt                                                                 <  3:12944/108‑1 (MQ=255)                                                                                                                        aGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTCTAATGTGATGttt                                                                 >  4:12944/1‑108 (MQ=255)                                                                                                                              aCTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTCTAATGTGATGTTTCAGACGTACCACTGAGTATCTGGTACAGGCTCCGTTGATCATCACt                   <  5:211923/148‑1 (MQ=255)                                                                                                                                    ggTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTCTAATGTGATGTTTCAGACGTACCACTGAGTATCTGGTACAGGCTCCGTTGATCATCACTCATCAAg            >  7:10761/1‑149 (MQ=255)                                                                                                                                        cGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTCTAATGTGATGTTTCAGACGTACCACTGAGTATCTGGTACAGGCTCCGTTGATCATCACTCATCAAGCGtt        <  1:4601/149‑1 (MQ=255)                                                                                                                                         gACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTCTAATGTGATGTTTCAGACGTACCACTGAGTATCTGGTACAGGCTCCGTTGATCATCACCCATCAAGCGTtg       >  6:166496/1‑149 (MQ=255)                                                                                                                                         gACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTCTAATGTGATGTTTAAGACGTACCACTGAGTATCTGGTACAGGCTCCGTTGATCATCACTCATCAAGCGTtg       >  7:96883/1‑149 (MQ=255)                                                                                                                                              atatGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTCTAATGTGATGTTTCAGACGTAc                                                        >  8:219880/1‑95 (MQ=255)                                                                                                                                              atatGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTCTAATGTGATGTTTCAGACGTAc                                                        <  7:219880/95‑1 (MQ=255)                                                                                                                                              atatGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTCTAATGTGATGTTTCAGACGTACCACTGAGTATCTGGTACAGGCTCCGTTGATCATCACTCATCAAGCGTTGTGTAt  >  3:239483/1‑149 (MQ=255)                                                                                                                                                 tGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTCTAATGTGATGTTTCAGACGTACCACTGAGTATCTGGTACAGGCTCCGtt                             >  5:308479/1‑119 (MQ=255)                                                                                                                                                 tGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTCTAATGTGATGTTTCAGACGTACCACTGAGTATCTGGTACAGGCTCCGtt                             <  6:308479/119‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                              GTGGGGACGAAGGGTACCTCGAACTATGTAGCCAAGAGTTTCCACTTCAACCCCCAGATCAAGTACGACGTCGGTAAGGCGAT‑GGGGCAAGCGCCGGCGCGTTTCTATGCCGGGATCGAGTGGGACTATTGGTCCGACAAATATGGCATCGAAGACAGCCACGGCTTTCCGACCGATCAGAACGCGTTGAGTGCATTGATCAAGGTGCATTTCTAATGTGATGTTTCAGACGTACCACTGAGTATCTGGTACAGGCTCCGTTGATCATCACTCATCAAGCGTTGTGTAT  >  NC_002947/3648308‑3648596

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