Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,642,313	G→C	13.2%	intergenic (‑99/+53)	PP_1442 ← / ← lon‑I	hypothetical protein/DNA‑binding ATP‑dependent protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,642,313	0	G	C	13.2%	154.5 / 12.5	53	intergenic (‑99/+53)	PP_1442/lon‑I	hypothetical protein/DNA‑binding ATP‑dependent protease
Reads supporting (aligned to +/- strand):  ref base G (23/23);  new base C (2/5);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCG  >  NC_002947/1642166‑1642459                                                                                                                                                    |                                                                                                                                                    cAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGATGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAggg                                                                                                                                                   <  1:191935/150‑1 (MQ=255)              gAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGt                                                                                                                                     >  2:387064/1‑150 (MQ=255)                    cAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGt                                                                                                                               <  2:351369/150‑1 (MQ=255)                       gggAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCg                                                                                                                            <  1:358418/150‑1 (MQ=255)                        ggAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCACGGGTACGC‑TACCCCTGCAGGCGCgg                                                                                                                           <  2:351045/150‑1 (MQ=38)                         gAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCggg                                                                                                                          <  2:322409/150‑1 (MQ=255)                             gAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTg                                                                                                                      >  2:52911/1‑150 (MQ=255)                              aGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGc                                                                                                                     <  2:603998/150‑1 (MQ=255)                                     gaggagAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAggg                                                                                                                                                   >  2:353624/1‑114 (MQ=255)                                     gaggagAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAggg                                                                                                                                                   <  1:353624/114‑1 (MQ=255)                                          gagGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAggg                                                                                                         <  1:241247/150‑1 (MQ=255)                                           agGTAATGGTGGCGCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGGGCAAGCCCGCACCTGCACGGGTA‑CGCTACCCCTGCAGGCGCGGGCTTACTCGCGaaaggacc                                                                                                        <  1:93337/150‑7 (MQ=17)                                              tCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGt                                                                                                     <  1:613479/150‑1 (MQ=255)                                                   gtggCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAgg                                                                                                <  2:373801/150‑1 (MQ=255)                                                       cTCCGAAGGCTGGCGATGGCCTGAGGATATCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAg                                                                                            >  1:314456/1‑150 (MQ=255)                                                       cTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGAGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAg                                                                                            >  1:292760/1‑150 (MQ=255)                                                          cGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAAc                                                                                         >  1:467473/1‑150 (MQ=255)                                                           gAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACa                                                                                        >  2:615288/1‑150 (MQ=255)                                                                 tGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACtt                                                                                  >  1:22112/1‑150 (MQ=255)                                                                 tGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACtt                                                                                  >  1:163945/1‑150 (MQ=255)                                                                     gATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGc                                                                              <  1:52911/150‑1 (MQ=255)                                                                                 aTACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCaaa                                                                  <  1:562949/150‑1 (MQ=255)                                                                                              gAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGaaa                                                     >  1:500876/1‑150 (MQ=255)                                                                                                 gAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCg                                                                                                       >  2:611965/1‑97 (MQ=255)                                                                                                 gAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCg                                                                                                       <  1:611965/97‑1 (MQ=255)                                                                                                 gAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGc                                                  >  2:375941/1‑150 (MQ=255)                                                                                                  aaaTTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCa                                                 <  1:387064/150‑1 (MQ=255)                                                                                                  aaaTTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCa                                                 <  2:22112/150‑1 (MQ=255)                                                                                                   aaTTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCAc                                                <  2:500876/150‑1 (MQ=255)                                                                                                                   gCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCgcgc                                >  1:253912/1‑150 (MQ=255)                                                                                                                    ccGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCa                               <  1:185899/150‑1 (MQ=255)                                                                                                                     cGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAaga                                                                                                            >  2:81387/1‑72 (MQ=255)                                                                                                                     cGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAaga                                                                                                            <  1:81387/72‑1 (MQ=255)                                                                                                                       gTCCTTTCGCGAGTCAGACCTCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGt                            >  2:458560/1‑150 (MQ=255)                                                                                                                          cTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTc                                                                                >  2:210491/1‑95 (MQ=255)                                                                                                                          cTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTc                                                                                <  1:210491/95‑1 (MQ=255)                                                                                                                          cTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGt                         >  2:542222/1‑150 (MQ=255)                                                                                                                          cTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGt                         >  1:247324/1‑150 (MQ=255)                                                                                                                           tttCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTc                        <  1:615288/150‑1 (MQ=255)                                                                                                                           tttCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTc                        >  1:527732/1‑150 (MQ=255)                                                                                                                           tttCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTc                        >  2:164800/1‑150 (MQ=255)                                                                                                                            ttCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTcc                       <  2:467473/150‑1 (MQ=255)                                                                                                                                cgAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCa                   <  2:543664/150‑1 (MQ=255)                                                                                                                                   gcaaGCCCGCACCTGCACGGGTACGC‑TACCCCTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCa                   <  1:34711/145‑1 (MQ=38)                                                                                                                                   gcaaGCCCGCACCTGCACGGGTA‑CGCTACCCCTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCa                   >  2:34711/3‑147 (MQ=38)                                                                                                                                       gCCCGCGCCTGCAGGGGTAGCG‑TACCCGAGCAGGTGCGGGCTTGCCCGCGAAGAAGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGttcttc            >  1:22197/1‑150 (MQ=255)                                                                                                                                        cccGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCg           <  1:358348/150‑1 (MQ=255)                                                                                                                                         ccGCGCCTGTAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTcgc                                  >  1:74019/1‑126 (MQ=255)                                                                                                                                         ccGCGCCTGTAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTcgc                                  <  2:74019/126‑1 (MQ=255)                                                                                                                                         ccGCACCTGCACGGGTACGC‑TACCCCTGCAGGCGCGGGCTTACTCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGttcttc            >  2:215455/1‑148 (MQ=25)                                                                                                                                         ccGCACCTGCACGGGTACGC‑TACCCCTGCAGGCGCGGGCTTACTCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGttcttc            <  1:215455/148‑1 (MQ=25)                                                                                                                                           gcgcCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAg        >  1:214064/1‑150 (MQ=255)                                                                                                                                               cTGCACGGGTACGC‑TACCCCTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGttcttc            <  2:429894/142‑1 (MQ=39)                                                                                                                                                tGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTcc                       >  2:112503/1‑130 (MQ=255)                                                                                                                                                tGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTcc                       <  1:112503/130‑1 (MQ=255)                                                                                                                                                 gCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTcgccg  <  2:163945/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCG‑TACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCG  >  NC_002947/1642166‑1642459

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