Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,431,883	G→A	8.7%	intergenic (‑12/‑38)	PP_1253 ← / → xenA	ArsR family transcriptional regulator/protein XenA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,431,883	0	G	A	8.7%	157.5 / 2.8	57	intergenic (‑12/‑38)	PP_1253/xenA	ArsR family transcriptional regulator/protein XenA
Reads supporting (aligned to +/- strand):  ref base G (22/30);  new base A (3/2);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

CCAGGTGGCGGACGATTTCCAGGCGCACGGGGTCGCTGAGCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATGGCCGAAGACGGCATGATCAACGA  >  NC_002947/1431736‑1432027                                                                                                                                                    |                                                                                                                                                   ccAGGTGGCGGACGATTTCCAGGCGCACGGGGTCGCTGATCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGa                                                                                                                                                   <  8:175074/149‑1 (MQ=255)  cAGGTGGCGGACGATTTCCAGGCGCACGGGGTCGCTGAGCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGAt                                                                                                                                                  <  8:268875/149‑1 (MQ=255)  cAGGTGGCGGACGATTTCCAGGCGCACGGGGTCGCTGAGCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGAt                                                                                                                                                  <  5:302776/149‑1 (MQ=255)  cAGGTGGCGGACGATTTCCAGGCGCACGGGGTCGCTGAGCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGAt                                                                                                                                                  <  6:95382/149‑1 (MQ=255)      tGGCGGAAGATTTCGAGGCGCACGGGGTCGCTGAGCGCATAGAGCAGGCGTTCAAGAATGACGGCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGaa                                                                                                                                              <  2:345502/149‑1 (MQ=255)           gACGATTTCCAGGCGCATGGGGTCGCTGAGCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAacaaca                                                                                                                                         <  7:104975/149‑1 (MQ=255)           gACGATTTCCAGGCGCACGGGGTCGCTGAGCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATCATAATCGaactatg                                                                                                                                         <  2:189643/149‑5 (MQ=255)            aCGATTTCCAGGCGCACGGGGTCGCTGAGCTCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCTCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAacaaca                                                                                                                                         <  1:134023/148‑1 (MQ=255)            aCGATTTCCAGGCGCACGCGGTCGCAGAGCGCATAGAGCAGGCGTTCAAGAATCAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCAAGGTTGACAAGAAAACGCGGGTTTCATTAATTGCCATAGTTCGATAAAGATCGaacaaaa                                                                                                                                         >  2:134023/1‑146 (MQ=255)             cGATTTCCAGGCGCACGGGGTCGCTGAGCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATCATAATCGaactat                                                                                                                                          <  4:249214/146‑4 (MQ=255)             cGATTTCCAGGCGCACGGGGTCGCTGAGCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATCATAATCGaactat                                                                                                                                          >  3:249214/1‑143 (MQ=255)                 ttCCAGGCGCACGGGGTCGCTGAGCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTAtt                                                                                                                                    >  7:300413/1‑148 (MQ=255)                           aCGGGGTCGCTGAGCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAAccc                                                                                                                       <  5:75936/149‑1 (MQ=255)                                gTCGCTGAGCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAAACCGGAGt                                                                                                                  >  6:221530/1‑149 (MQ=255)                                     tGAGCGCATAGAGCAGGCGTTCAAGAATGGGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTcc                                                                                                              >  2:262634/1‑148 (MQ=255)                                     tGAGCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTcc                                                                                                              >  4:45051/1‑148 (MQ=255)                                        gcgcATAGAGCAGGCGTTCAATAATGATCTCGTCGGGGTGTGGATGTTTATTGGCTCGCATTGTGGACATGATAACGGGGGTTTCATTAATTTCCATTGTTCGATTAGGATCGTACAACATTATTCAA‑A‑TGAACCCGGAGTTTCCCATg                                                                                                          <  6:173179/149‑1 (MQ=255)                                          gcATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTccc                                                                                                             >  5:212720/1‑144 (MQ=255)                                          gcATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTccc                                                                                                             <  6:212720/144‑1 (MQ=255)                                          gcATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTc                                                                                                        <  2:131340/149‑1 (MQ=255)                                            aTAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATTGAACAACAGTATTCAA‑A‑TGAACCCg                                                                                                                      <  5:168920/133‑1 (MQ=255)                                            aTAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATTGAACAACAGTATTCAA‑A‑TGAACCCg                                                                                                                      >  6:168920/1‑133 (MQ=255)                                               gagCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCa                                                                                                    <  7:267786/148‑1 (MQ=255)                                                       gTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAACA‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGa                                                                                            <  1:51600/149‑1 (MQ=255)                                                        ttCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAAc                                                                                          <  5:302549/149‑1 (MQ=255)                                                        ttCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAAc                                                                                          >  4:68860/1‑149 (MQ=255)                                                             gAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTa                                                                                      >  1:266119/1‑148 (MQ=255)                                                               aTGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTacac                                                                                   <  1:89002/149‑1 (MQ=255)                                                                          ggggTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑ACTGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGAc                                                                         <  1:262634/149‑1 (MQ=255)                                                                          ggggTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACCCCCTCAAATACg                                                                        >  4:201824/1‑149 (MQ=255)                                                                                 ggATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCt                                                                 >  5:148918/1‑149 (MQ=255)                                                                                   aTGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCTAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGc                                                               <  3:68860/149‑1 (MQ=255)                                                                                          taGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAAcc                                                         <  3:45051/148‑1 (MQ=255)                                                                                           aGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCg                                                        >  5:86519/1‑148 (MQ=255)                                                                                             gCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGAc                                                                         <  1:144646/129‑1 (MQ=255)                                                                                             gCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGAc                                                                         >  2:144646/1‑129 (MQ=255)                                                                                                     ggTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGc                                                  >  8:196122/1‑144 (MQ=255)                                                                                                     ggTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGc                                                  <  7:196122/144‑1 (MQ=255)                                                                                                       ttGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTccg                                           >  1:21601/1‑149 (MQ=255)                                                                                                            aTGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGcc                                                 <  6:143175/138‑1 (MQ=255)                                                                                                            aTGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGcc                                                 >  5:143175/1‑138 (MQ=255)                                                                                                              gATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGAtgt                                    >  5:270037/1‑149 (MQ=255)                                                                                                               aTAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGAtgt                                    <  5:221530/148‑1 (MQ=255)                                                                                                                tAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGc                                  >  4:144742/1‑149 (MQ=255)                                                                                                                            cATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATGGc                       <  6:148918/148‑1 (MQ=255)                                                                                                                             aTTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATGGcc                      <  8:300413/148‑1 (MQ=255)                                                                                                                             aTTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATGGcc                      >  2:163698/1‑148 (MQ=255)                                                                                                                              ttAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATGGCCg                     >  2:282277/1‑148 (MQ=255)                                                                                                                                     ccATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATGGcc                      >  6:18329/1‑140 (MQ=255)                                                                                                                                     ccATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATGGcc                      <  5:18329/140‑1 (MQ=255)                                                                                                                                     ccATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATGGCCGAAGACgg              <  3:144742/148‑1 (MQ=255)                                                                                                                                     ccATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATGGCCGAAGACgg              <  6:270037/148‑1 (MQ=255)                                                                                                                                      ttttgTTCGATCATAATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACACTCAAATACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATGGCCGAAGACGGCa            >  7:220338/5‑149 (MQ=255)                                                                                                                                      tgtAGTTCGAACATAATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGATTCCCATGTCCGCTCTGTTCGAACCCTACACCCTGAAAGACGTCACCCTCCGTAACCGTATCGCCATTCCGCCCCTGTGCCAGTACATGGCCTAAGACTGCa            >  4:28354/3‑149 (MQ=255)                                                                                                                                             cGATTATGATCGCACAACAGTATTCAT‑A‑TGAACCCGGTGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGCCGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATg                         <  4:77222/129‑1 (MQ=255)                                                                                                                                             cGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGAGCCAGTACATg                         >  3:77222/1‑129 (MQ=255)                                                                                                                                              gATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATGGCCGAAGACGGCATGATCaa     <  6:86519/148‑1 (MQ=255)                                                                                                                                                    gATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATGGCCGAAGACGGCATGATCAACGa  >  7:143868/1‑145 (MQ=255)                                                                                                                                                    gATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATGGCCGAAGACGGCATGATCAACGa  <  8:143868/145‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CCAGGTGGCGGACGATTTCCAGGCGCACGGGGTCGCTGAGCGCATAGAGCAGGCGTTCAAGAATGAGGTCGTCGGGGTTGGGATGTTTATAGGCTCGCATGGTTGACATGATAACGGGGGTTTCATTAATTGCCATAGTTCGATTATGATCGAACAACAGTATTCAA‑A‑TGAACCCGGAGTTTCCCATGTCCGCACTGTTCGAACCCTACACCCTCAAAGACGTCACCCTGCGTAACCGTATCGCCATTCCGCCGATGTGCCAGTACATGGCCGAAGACGGCATGATCAACGA  >  NC_002947/1431736‑1432027

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