Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	1,964,559	G→A	N226N (AAC→AAT)	LY180_09345 ←	transposase IS609

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	1,964,559	0	G	A	88.5%	38.9 / ‑6.9	29	N226N (AAC→AAT)	LY180_09345	transposase IS609
Reads supporting (aligned to +/- strand):  ref base G (1/1);  major base A (9/14);  minor base C (2/1);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 5.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCTGCTTGTACTCAAGCTGGCGGCGCATTTCATACCAGCCCTGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTG‑CGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTGTTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGACCATTGATGCTGAAGGGTGAACCGGAGTGGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACAGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA  >  CP006584/1964290‑1964842                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                            agttgctTGTACGCAATGTGGCGGCGCATTTCATACCTGCCCTGATCCAGTATCGACCGGTTTAACCCGTATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCGGACTTTGACCTGTTGCTGACCTTCAAATCCGCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTGTGAAGGTAGTCCCTGCGGATATTTG‑CGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTGTTGCTGAATTTGACCTTGCGGCTTAACt                                                                                                                                                                                                                                                                 <  2:235119/295‑1 (MQ=1)   ctgctTGTACTCAAGCTGGCGGCGCATTTCATACCAGCCCTGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACCTGTTGCTGACCTTCCAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTG‑CGTTACGGGAATGCAGTCGCTGTATTTTTGGTTTCTGCTTCTGCCAGTTTTTGCTGCATTTGACCTTggcgttttacttt                                                                                                                                                                                                                                                               >  2:54072/1‑286 (MQ=11)          tACTCAAGCTGGCGGCGCATTTCATACCAGCCCTGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTG‑CGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGaa                                                                                                                                                                                                                                                        <  1:97248/298‑1 (MQ=14)                           cATTTCATACCAGCCCTGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTG‑CGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGtcttc                                                                                                                                                                                                                                       <  1:54072/298‑1 (MQ=14)                                  tACCAGCCCTGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTG‑CGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGtttt                                                                                                                                                                                                                                >  1:111891/1‑298 (MQ=14)                                    ccAGCCCTGATCCAGTATCGAACGGTTTAAACCTGATATTGCCCGGACATTGCGCCACGGCTGGCTGACAGTAACCGCTGATGCCTTTGACATGTTGCTGCCCTTCAAATCCTCAATGAATATCCTTGCGTGGTTTTTGCTGACGGTCGTTGTCACTATATGCAGGTAGTCCCTGCGGATATTTG‑CGATACGGTAATGAATGCGCTGTGCTTTTCGTTTCTGCTTCTGCCAGTTCATGCTAGATTTGACTTTGCGGGTAAACTGTCCCTGGAGGCTACCCACAGTCTTGGCGTTTTTc                                                                                                                                                                                                                              >  2:268462/1‑298 (MQ=11)                                      aGCCCTGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTG‑CGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTg                                                                                                                                                                                                                            >  1:144345/1‑298 (MQ=14)                                        cccTGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTG‑CGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGaa                                                                                                                                                                                                                          <  1:300394/298‑1 (MQ=14)                                               ccAGTATTGAACGGTTTAAACCGTAGTTTGCCCGGACATTGGGCCCCGGCGGGCTGACCGTACCCGCTGGTGACTTTGGTATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTG‑CGATACGGGAATGGAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGtt                                                                                                                                                                                                                   <  2:268722/298‑1 (MQ=11)                                                 aGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGCCTATCATTGCGTGGTTTTTGCTGGCGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTG‑CGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCTGTTCTTGCTGACTTTGTCCTTGCGGCTTACCTTTTGCTATCGTCTCCCCCCCTCTTTTTGGTCTTTTTTAAAACTGtttt                                                                                                                                                                                                                 >  2:222510/1‑297 (MQ=11)                                                                        tgttGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTG‑CGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCt                                                                                                                                                                                          <  1:266751/296‑1 (MQ=17)                                                                                        cgcCCGGGCGGGCTGACCGTTCCCGGTTCTGACTTTGGCCTGTTGCTGACCCTCAAATCCTCAAGGACTATCATTGCGTTGTTGTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTG‑CGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGtt                                                                                                                                                                          <  2:65259/298‑1 (MQ=11)                                                                                             cGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTG‑CGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGcc                                                                                                                                                                     >  1:274487/1‑298 (MQ=17)                                                                                                  ggctgACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGCCTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTGTGAAGGTAGTCCCTGCGGATATTTG‑CGATACGGGAATGCAGTCGCTGTATCTTGCGTTTCTGCTTCTGCCAGTTGTTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCCCCAGCGTCTTCTGGCTTTTCTGAAAACCTTTTCACGGCCCCAAAACCGTGCCATCTGCCATCGTTGCGAGTTTAGCCAccac                                                                                                                                                                >  2:89030/1‑295 (MQ=11)                                                                                                   gctgACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTG‑CGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCa                                                                                                                                                               >  1:244057/1‑298 (MQ=17)                                                                                                         cGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTG‑CGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATcc                                                                                                                                                         >  1:96088/1‑298 (MQ=17)                                                                                                                                          cTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTG‑CGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGa                                                                                                                        <  1:303655/298‑1 (MQ=17)                                                                                                                                               aaTCCTCAATGACAATCATTGCGTGGTTTTTGCTGACGGACGTTGTGACTATATGAAGGTAGTCACTGCGGATATTTG‑CGATACGGGAATGCAGTCGCTGCATCTTGCGTTTCTGCTTATGCCAGTTATTGCGCAATTGGACCTAGCTGCTTCACTGTCGATGCAGTATCGCCAGAGCCTTATAGTTGTGCTGAATACAGGTTCAAGGTCCAAAGAATGTGCCCTATAACAGCGCGGCGAGTTTAGTCACGCACGCATGCACCCAGATCACTACGGCTGAAGGTGGAACCGGAggcga                                                                                                                   >  2:128234/1‑294 (MQ=2)                                                                                                                                                   cTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTG‑CGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACt                                                                                                               >  1:198917/1‑298 (MQ=17)                                                                                                                                                             ttgcttgggggggttttGTTGACGTGCTGTGTAACTGTATGGAGGTAGTCCATTCCGGTATTTGCCGGTA‑GGGGGTGCTGTCGGTTTATTTTGCGTTGCTGCTTTCGGCAGTTCTTGCTGAATTTGTCCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAGCGGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTtctg                                                                                                     <  2:285433/285‑1 (MQ=11)                                                                                                                                                              tCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTG‑CGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTtctgt                                                                                                    <  1:78928/298‑1 (MQ=17)                                                                                                                                                               cATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTG‑CGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTtctgtc                                                                                                   <  1:163828/298‑1 (MQ=17)                                                                                                                                                                      tGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTG‑CGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTTAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGATTTTAGCCACCCCAGCACCCCGCCCCGCCATTCATCCCGAAGGGTGCACCGGGGTGGATACCTCACTTTGTGTCTTAAttc                                                                                            >  2:72888/1‑296 (MQ=11)                                                                                                                                                                              gCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTG‑CGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTAc                                                                                    <  1:72888/298‑1 (MQ=17)                                                                                                                                                                                            tgACTTTATGAAGGTAGTCCCTGCGGATATTTG‑CGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGa                                                                      <  1:222510/298‑1 (MQ=17)                                                                                                                                                                                                tttATGAAGGTAGTCCCTGCGGATATTTG‑CGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTgg                                                                  >  1:309521/1‑298 (MQ=17)                                                                                                                                                                                                                              cGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGt                                     <  1:176736/298‑1 (MQ=17)                                                                                                                                                                                                                                                 cTGTATTTTGCGTTTCTGCTTCTGCCAGTTATGGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCgg                  <  1:212069/298‑1 (MQ=9)                                                                                                                                                                                                                                                                 tgcttctgcCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCa  <  1:243872/298‑1 (MQ=14)                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                            AGCTGCTTGTACTCAAGCTGGCGGCGCATTTCATACCAGCCCTGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTG‑CGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTGTTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGACCATTGATGCTGAAGGGTGAACCGGAGTGGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACAGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCA  >  CP006584/1964290‑1964842

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑