Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	2,134,732	T→G	24.1%	L471V (TTG→GTG)	comEC →	ComE operon protein 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	2,134,732	0	T	G	24.1%	54.7 / 4.8	52	L471V (TTG→GTG)	comEC	ComE operon protein 3
Reads supporting (aligned to +/- strand):  ref base T (10/28);  new base G (12/2);  total (22/30)
Fisher's exact test for biased strand distribution p-value = 2.50e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGTTGGGAG‑GGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTACGCGGCATCACAGCCTGCT  >  KT2440/2134599‑2134862                                                                                                                                      |                                                                                                                                    gTTCCGGTGGCTGGCGCTGGGGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGtt                                                                                                                                    <  2:220179/135‑1 (MQ=255)  ttCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGTTgg                                                                                                                                  <  4:2060/136‑1 (MQ=255)  ttCCGTGGGCGGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGtt                                                                                                                                    <  3:189312/134‑1 (MQ=255)  ttCCGTGGGCGGGCGCTGGGGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGtt                                                                                                                                    <  3:294444/134‑1 (MQ=255)  gtCCGTGGGCTGGCGATGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGtt                                                                                                                                    <  4:258090/133‑1 (MQ=255)          gCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGGGGTGTGCCGCGGCGGGGGGTggggg‑ggg                                                                                                                           >  5:198088/1‑134 (MQ=255)              gcgcTGGGGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGTTGGGAG‑GGGTAAtg                                                                                                                      <  2:319409/135‑1 (MQ=255)               cgcTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTGGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGTTGGGAG‑GGGTTAtgc                                                                                                                     >  6:193431/1‑135 (MQ=255)               cgcTGGTGGCGCAGCAGCACCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGtt                                                                                                                                    <  4:32894/121‑1 (MQ=255)               cgcTGGTGGCGCAGCAGCACCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGtt                                                                                                                                    >  3:32894/1‑121 (MQ=255)                    gtggCGCAGCAGGGCCCGGCATGGGTCCCGCCGGCCCTGGCGTTGTGGGCCTGGCTGCTGGGATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGGGCCGCTTCGCGGGTGTGGAG‑GGGTAATGctgctg                                                                                                                <  5:273588/135‑1 (MQ=255)                       gcgcAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGGTGCGCGGGGTGGGAG‑GGGTAATGCTGCTGgcg                                                                                                             >  8:360530/1‑135 (MQ=255)                                 gcCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGTTGGGAG‑GGGTAATGCTGCTGGCGCTGTGGGTGcc                                                                                                  <  4:1277/136‑1 (MQ=255)                                 gcCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGGTGGGAG‑GGGTTATGCTGCTGGCGCTGTGGGTGcc                                                                                                  >  4:356549/1‑136 (MQ=255)                                     ggCATGGGTCCCGCCTGACCTGCCGTTGTGTGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCACCGTGGTGTGGCCCTTCGAGGGGTGGGAGTGGGTAAT‑CTTCTGGGGCGGTGGGttaccag                                                                                               >  6:156689/1‑129 (MQ=255)                                         tGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGTTGGGAG‑GGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAg                                                                                           <  8:164280/135‑1 (MQ=255)                                          gggTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTTTGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCAGCCCGGGGTGGGAG‑GGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGCGcc                                                                                         >  4:83076/1‑136 (MQ=255)                                          gggTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGTTGGGAG‑GGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGcc                                                                                         <  4:174296/136‑1 (MQ=255)                                           gggcccGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGGGGGCCGCTGCGCGGGTTGGGAG‑GGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGcc                                                                                         <  6:300961/132‑1 (MQ=255)                                                cgccgGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCGGCGGGGGGTGGGGG‑GGGGTATGCTGCTGGCGCGGTGGGTGCCCCGGGCGCCGGTGc                                                                                    >  8:233173/1‑135 (MQ=255)                                                    ggCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCGGGGCGGGGTGGGAG‑GGGGAATGCGGCTGGCGCGGGGGGTGCCCCGGGGGGCgggggcgg                                                                                 >  3:66788/1‑128 (MQ=255)                                                        cTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGGGCGGGTTGGGGG‑GGGTGATGCTGCGGGCGCTGTGGGTGCCCCGGGGGCCGGTGCCGtggggc                                                                            >  7:354814/1‑130 (MQ=255)                                                         tGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGTTGGGAG‑GGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGcc                                                                           <  7:350413/135‑1 (MQ=255)                                                              ttGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGGGCGGGTTGGGAG‑GGGGAATGCTGCTGGCGCGGTGGGTGCCCCGGGAGCCGGTGGCGTTCGGCCCgggc                                                                      >  1:42707/1‑133 (MQ=255)                                                                  gggCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTTCTGCCCCGTGGTGTGGCGCATCGCGGGGTGGGAG‑GGGTAAATCTTCTTGCGCTGTGGGTGCCCCCGGCTCCGGTGCCCTTCGGCCAGGTCgggg                                                                  >  6:330559/1‑135 (MQ=255)                                                                         gctgctGGTATGTCTGGGGGGGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGTTGGGAG‑GGGTAACGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGc                                                            <  4:147415/134‑1 (MQ=255)                                                                                     tCTGGGTGTGCTGTTGGTTTTGCTGCCCCGGGGTGTGCCGCTGCGCGGGTTGGGAG‑GGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTgg                                               <  4:229505/135‑1 (MQ=255)                                                                                      cTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGTTTGGAG‑GGGTAATGCTGCTGGAGCTGTTGGTGCCCCGGGAGCCGGTTCAGTTCGGACATGTCGTGGTCTTGAACCAGGATGTAGGc                                              >  7:315577/1‑135 (MQ=255)                                                                                        gggTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCCCGGGGGGGGGGGGGTG‑GGGTAAAGCTGGAGGCGGTGTGGGGGCCCCGGGAGCCGGTGCCGGTAGGGCAGGTCCAGGTCTGGGCGCTGGGTGTTTGCCAg                                           >  4:37687/1‑136 (MQ=255)                                                                                         ggTGTGCTGTTGGTTTTGCTGCCCCGGGGTGTGCCGCTGCGCGGGTTGGGAG‑GGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAg                                           <  7:53887/135‑1 (MQ=255)                                                                                            gtgCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCGGCGCGGGGTGGGTG‑GGGGTATGCTGCTGGGGCGGTGGGTGCCCCGGGCGCCCGTGCCGGTAGGGCAGGTCGCGGGGGGGGAGCGGGATTTTGGCCCCGGGc                                       >  3:8306/1‑136 (MQ=255)                                                                                             tgCTGTTGGGTTGGCTGCCCCGGGGTGTGCCGCTGCGCGGGTGGGGAG‑GGGTAATGCTGCTGGCGCGGTGGGTGCCCCGGGAGCCGGTGCCGGTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGc                                       <  5:140921/135‑1 (MQ=255)                                                                                              gCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGTTGGGAG‑GGGAAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTg                                     <  4:107265/136‑1 (MQ=255)                                                                                                tGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGGTGGGTC‑GGGTAATGCTTCGGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGGTAGGCCAGGGCCAGGGCGGGCAGCGGGATGTTGTCCAgggggggg                                    >  2:295962/1‑131 (MQ=255)                                                                                                tGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGGTGGGAG‑GGGTAATGCAGCAGGCGCAGTTGGTGCCCCGGGAGCCCGTGCCCTTCGGCCACGTCGAGGTCTGGCAGCAGGAAGTTGGGCAGGGGCTgg                                    >  2:248636/1‑135 (MQ=255)                                                                                                  ttGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGTTGGGAG‑GGGTAATGCTGGTGGCGCTGTGGGTGCCCCGGGCGCCGGTGCCGGTCGGCCAGGTCGAGGTCTGGGAGCTGGATGTTTGCCAGGGGCTGGcc                                  >  7:150838/1‑134 (MQ=255)                                                                                                  ttGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGTTGGGAG‑GGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGCGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCg                                  >  2:95733/1‑135 (MQ=255)                                                                                                  ttGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGTTGGGAG‑GGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGGCGAGGTCTTGCAGCTGGATGTTGGCCAGGGGCTGGCg                                  >  5:343322/1‑135 (MQ=255)                                                                                                    ggTTTTGCTGCCCCGTGGTGGGCCGCGGCGCGGGGTGGGAG‑GGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGt                                <  7:205378/135‑1 (MQ=255)                                                                                                    ggTGTTGCTGCCCCGTGGTGTGGCGCGGCGCGGGTTGTGAG‑GGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGGCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGt                                <  7:327239/135‑1 (MQ=255)                                                                                                           ctgcCCCGTGGTGTGCCGCTGCGCGGGTTGGGAG‑GGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGc                        <  4:68913/136‑1 (MQ=255)                                                                                                               cccGTGGTGTGCCGCTGCGGGGGTTGGGAG‑GGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTAc                    <  3:83076/136‑1 (MQ=255)                                                                                                               cccGTGGGGTGCCGCTGCGCGGGTGGGGAG‑GGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGGTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTaa                    <  3:234548/136‑2 (MQ=255)                                                                                                                ccGTGGTGTGCCGCTGCGCGGGTTGGGAG‑GGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTAc                    <  8:47280/135‑1 (MQ=255)                                                                                                                ccGTGGTGTGCCGCTGCGCGGGTTGGGAG‑GGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTAc                    <  2:174869/135‑1 (MQ=255)                                                                                                                ccGTGGTGGGCCGCTGCGCGGGTTGGGAG‑GGGTAATGCTGCTGGCGCTGTGGGGGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTAc                    <  1:61159/135‑1 (MQ=255)                                                                                                                        tgCCGCTGCGCGGGTTGGGAG‑GGGTAATGCGGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTACGCGGCATc            <  6:198088/135‑1 (MQ=255)                                                                                                                          ccGCTGCGCGGGTGGGGAG‑TGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGACGTGCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTTTTGCGTACGCGGCATcac          <  7:82193/135‑1 (MQ=255)                                                                                                                          ccGCGGCGCGGGTTGGGAG‑GGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTACGCGGCATcac          <  8:318840/135‑1 (MQ=255)                                                                                                                           cGCTGCGCGGGTTGGGAG‑GGGTAATGCTGCTGGCGCTGTGGGTTCCCCGGGACCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTTCGTACGAGGCATcaccg        >  3:117317/1‑134 (MQ=255)                                                                                                                           cGCTGCGCGGGTGGGGAG‑GGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTACGCGGCATcaca         <  5:306962/135‑1 (MQ=255)                                                                                                                            gcggcgcGGGTTGGGAG‑GGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTACGCGGCATCACAGc       <  3:37687/133‑1 (MQ=255)                                                                                                                            gCTGCGCGGGTTGGGAG‑GGGTAATGCTGCTGGCGCTGTGGGGGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTACGCGGCATCACAGc       <  3:165508/136‑1 (MQ=255)                                                                                                                             cTGCGCGGGTTTGGAG‑GGGTAATTCTGCTGGCGCTGTGGGTGCCCCCGGAGCCCGTTCCGATCGGCCACGTCGAGGTCTTGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTTCGTACGCGGCATCACAGc       >  3:217101/1‑135 (MQ=255)                                                                                                                                  cgGGTTGGGAG‑GGGTACTGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTACGCGGCATCACAGCCTGct  <  8:354814/135‑1 (MQ=255)                                                                                                                                  cgGGGTGGGCG‑GGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTACGCGGCATCACAGCCTGct  <  7:233173/135‑1 (MQ=255)                                                                                                                                     gTTGGGAG‑GGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTAc                    <  4:114247/114‑1 (MQ=255)                                                                                                                                     gTTGGGAG‑GGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTAc                    >  3:114247/1‑114 (MQ=255)                                                                                                                                      |                                                                                                                                    GTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGTTGGGAG‑GGGTAATGCTGCTGGCGCTGTGGGTGCCCCGGGAGCCGGTGCCGTTCGGCCAGGTCGAGGTCTGGCAGCTGGATGTTGGCCAGGGGCTGGCGGTGCTCTTGCGTACGCGGCATCACAGCCTGCT  >  KT2440/2134599‑2134862

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