Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,344,363	A→G	8.4%	intergenic (+50/‑106)	PP_3817 → / → PP_3818	polyamine ABC transporter ATP‑binding protein/phosphate‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,344,363	0	A	G	8.4%	156.4 / 2.5	57	intergenic (+50/‑106)	PP_3817/PP_3818	polyamine ABC transporter ATP‑binding protein/phosphate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (27/24);  new base G (3/2);  total (31/26)
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

GACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCTGCTACCACTC  >  NC_002947/4344221‑4344510                                                                                                                                               |                                                                                                                                                     gACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATc                                                                                                                                               <  2:33155/150‑1 (MQ=255)       gTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAg                                                                                                                                         <  2:685433/150‑1 (MQ=37)       gTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAg                                                                                                                                         >  1:146178/1‑150 (MQ=255)       gTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTGATAACTAATGGGCCGAGAGCGACCTGAGCAAGGCGTCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCCCACGATCTGGCAg                                                                                                                                         >  1:687947/1‑150 (MQ=255)           aGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGgcgc                                                                                                                                     >  2:821914/1‑150 (MQ=255)            ggTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCt                                                                                                                                    <  2:146178/150‑1 (MQ=255)             gTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTATGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTg                                                                                                                                   >  2:224711/1‑150 (MQ=255)                  ttGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAg                                                                                                                                         <  1:783595/139‑1 (MQ=25)                  ttGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAg                                                                                                                                         >  2:783595/1‑139 (MQ=25)                  ttGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTa                                                                                                                              <  2:13585/150‑1 (MQ=255)                     ggCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTg                                                                                                                           <  2:61188/150‑1 (MQ=255)                       cAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGa                                                                                                                         <  1:772331/150‑1 (MQ=255)                        aGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAg                                                                                                                        <  1:405462/150‑1 (MQ=255)                            gCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTgg                                                                                                                          <  1:301551/144‑1 (MQ=255)                            gCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTgg                                                                                                                          >  2:301551/1‑144 (MQ=255)                            gCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGt                                                                                                                    >  2:327664/1‑150 (MQ=255)                              aGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTATCGCTGTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTAc                                                                                                                  <  2:341884/150‑1 (MQ=255)                                        aTTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTcac                                                                                                        >  2:1185754/1‑150 (MQ=255)                                              aCCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTatc                                                                                                  >  1:680693/1‑150 (MQ=255)                                              aCCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGAGCATTAACGCTCTTGTCCGTTGCAGAGCCAGTCAGAGCGACCATCACCCTCGCGGCACACGATCTTGCAGGTGTTGCACTATTGGAGGCGTACGCCCCGTCACATTatc                                                                                                  <  2:1042113/150‑1 (MQ=255)                                                 gAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACCCGATCTTGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTatcatc                                                                                               >  2:938522/1‑150 (MQ=255)                                                       gCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTTGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCa                                                                                         <  1:755622/150‑1 (MQ=255)                                                       gCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCa                                                                                         >  1:333615/1‑150 (MQ=255)                                                         cGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCGCACGAGCTGGCGGGCGCGGCACTATTGGAGGCGTACGCCCCGTCACAGTATCATCACGCCACt                                                                                       >  1:664454/1‑150 (MQ=255)                                                         cGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACt                                                                                       >  2:897456/1‑150 (MQ=255)                                                              gCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGt                                                                                  >  2:663767/1‑150 (MQ=255)                                                                 aCCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑CAt                                                                              >  1:1140880/1‑150 (MQ=255)                                                                 aCCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AAt                                                                              <  2:246864/150‑1 (MQ=255)                                                                 aCCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AAt                                                                              <  1:897456/150‑1 (MQ=255)                                                                     gAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCGGGCTCTGCACTATTGGAGGCGCACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATCttc                                                                          >  1:203656/1‑148 (MQ=255)                                                                     gAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATCTGc                                                                          >  1:1040648/1‑150 (MQ=255)                                                                     gAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATCTGc                                                                          >  1:794297/1‑150 (MQ=255)                                                                       gCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATCTGCCg                                                                        >  2:678355/1‑150 (MQ=255)                                                                          aGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATCTGCCGGcc                                                                     >  1:910051/1‑150 (MQ=255)                                                                                       gcATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATCTGCCGGCCAACTCTATGCcac                                                        >  2:782406/1‑150 (MQ=255)                                                                                         aTTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATCTGCCGGCCAACTCTATGCcaca                                                       >  2:702457/1‑149 (MQ=255)                                                                                         aTTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATCTGCCGGCCAACTCTATGCcaca                                                       <  1:702457/149‑1 (MQ=255)                                                                                         aTTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATCTGCCGGCCAACTCTATGCCACAg                                                      <  1:663767/150‑1 (MQ=255)                                                                                         aTTAACGCTCTTGTACGTTGCACTGCATTTCAGAGCTGTCACCTCCCGCGCGGCACCCGATCTGGCTCGATCCTCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATCTGCCGGCCAACTCTATGCCACAg                                                      <  1:327664/150‑1 (MQ=255)                                                                                           tAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATCTGCCGGCCAACTCTATGCCACAGGt                                                    <  1:938522/150‑1 (MQ=255)                                                                                                         gTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATCTGCCGGCCAACTCTATGCCACAGGTCACCCgg                                             <  2:609009/143‑1 (MQ=255)                                                                                                         gTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATCTGCCGGCCAACTCTATGCCACAGGTCACCCgg                                             >  1:609009/1‑143 (MQ=255)                                                                                                         gTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGAGGCACCCGATCGGGGTGACCCTTCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATCTGCCGGCCAACTCTACGCCACAGGTCACCCGGCATGCTc                                      >  1:151174/1‑150 (MQ=255)                                                                                                             cAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACGATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACT‑TGGTGAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGcc                                  <  1:821914/150‑1 (MQ=255)                                                                                                                  ccTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCtt                             <  1:224711/150‑1 (MQ=255)                                                                                                                     gCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATTTAGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATCTGCCCGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGGTTGTc                          >  1:585386/1‑150 (MQ=37)                                                                                                                     gCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATCTGCCGGCCAACTCTATGCAACATCGCACCCGGCATGCTCCGCCTGCTTGTc                          <  1:337215/150‑1 (MQ=255)                                                                                                                      ccAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑CATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCt                         >  2:337215/1‑150 (MQ=255)                                                                                                                      ccAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCt                         >  1:706619/1‑150 (MQ=255)                                                                                                                       cAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCtt                        >  2:1142231/1‑150 (MQ=255)                                                                                                                        aGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtg                       <  2:794297/150‑1 (MQ=255)                                                                                                                               gCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATTTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCa                >  2:193339/1‑150 (MQ=255)                                                                                                                                     cccgcGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCTGc          <  2:687947/150‑1 (MQ=255)                                                                                                                                       cgcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCTGCTa        <  2:333615/150‑1 (MQ=255)                                                                                                                                        gcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCTGCTAc       >  1:375524/1‑150 (MQ=255)                                                                                                                                        gcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCTGCTAc       >  1:1074873/1‑150 (MQ=255)                                                                                                                                             gCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCTGCTACCACTc  <  2:1074873/150‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                     GACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGT‑AATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCTGCTACCACTC  >  NC_002947/4344221‑4344510

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