Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	14,232	G→A	R22K (AGA→AAA)	dnaJ →	chaperone protein DnaJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	14,232	0	G	A	97.1%	104.7 / NA	34	R22K (AGA→AAA)	dnaJ	chaperone protein DnaJ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (21/12);  minor base C (1/0);  total (22/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGG  >  NC_000913/13985‑14475                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    agcagcagACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGAAGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCTTTTCCAAAACAGCGGAAGAGCAGGAAATCAaaaa                                                                                                                                                                                                                                                  <  2:47066/251‑1 (MQ=255)               ggttcTGATGCTTCTGCAAACAACTCGAAATATGACGATGTTGACGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATGATGACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGcc                                                                                                                                                                                                                                    <  2:178066/247‑1 (MQ=255)                gTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCt                                                                                                                                                                                                                                   <  2:339277/251‑1 (MQ=255)                                         gAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCCAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGt                                                                                                                                                                                                          >  1:172091/1‑251 (MQ=255)                                                             aCGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGc                                                                                                                                                                                      >  2:255191/1‑251 (MQ=255)                                                                         aataagTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCTGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATGACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAAtg                                                                                                                                                                          <  1:313170/248‑2 (MQ=255)                                                                         aagaagTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAAtt                                                                                                                                                                          >  1:197633/1‑251 (MQ=255)                                                                                                        taAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTg                                                                                                                                           >  1:15607/1‑251 (MQ=255)                                                                                                             gggTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCAGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCg                                                                                                                                       >  1:464669/1‑250 (MQ=255)                                                                                                                           acacGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTg                                                                                                                        <  1:255191/251‑1 (MQ=255)                                                                                                                           aaacGGGCGAAGGGGAATTGCCTCCCCGCCCGTTCATGCAGCAAGGGGCAATTTAAAAAAGATGGCTAAGCATGATTATTTCGAGATTCTACGCGTTACCAAAACAGCCGAAGAGCACGAAACCAACAAGGCCTACAAACCCCTGACGATGAAATGCCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGATATCAAGGAAGCCTATGGAGTACTGACCGACTCCCACAAACGTg                                                                                                                        <  2:457724/249‑1 (MQ=255)                                                                                                                                  cGAAGGGGAATTTCCTCTCCGCCCGCGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCAGAGGCGAACTTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCggaat                                                                                                                  >  2:502216/1‑247 (MQ=255)                                                                                                                                          aaTTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGt                                                                                                         <  2:15607/251‑1 (MQ=255)                                                                                                                                           aTTTCCTCTCCGCCCGTGCATTCATCTAGGGGCATTTTAACAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGATGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCACGGAAGCTTATGAAGTTCTGACCGACTAGCAAAAACGTGCGGCATACGATCAGt                                                                                                         >  2:350193/1‑250 (MQ=255)                                                                                                                                                        ccGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCg                                                                                           >  1:62699/1‑251 (MQ=255)                                                                                                                                                               ttCATCTAGGGGCCATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGc                                                                                    >  2:452497/1‑251 (MQ=255)                                                                                                                                                                             aTTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATgggg                                                                      >  2:242948/1‑250 (MQ=255)                                                                                                                                                                                tAAAAAAGATGGATAAGCAAGATTATTACGAGATTTTAGGCGATTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGggcggc                                                                   >  1:114824/1‑251 (MQ=255)                                                                                                                                                                                       gctGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAACACAGCGGAAGAGCATGAGATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTCACCAGGGTGACAAAGAGGCCGCGGCGAATTTTAACGAGATCAAGTAAGCTTATGAAGTTCTGACCGACTCGCATAAACGTGCGGCATACGATCAGTCTGGTAATGCTGCATTTGAGCAAGGTgcgattggagagtgcggtt                                                            >  2:128996/3‑233 (MQ=255)                                                                                                                                                                                           gCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTgg                                                        <  1:35942/251‑1 (MQ=255)                                                                                                                                                                                              aaGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACCCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGGATGGGCGGCGGCGGGTTTggcgg                                                     >  2:498596/1‑251 (MQ=255)                                                                                                                                                                                                 cAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTggcggcgg                                                  >  1:252898/1‑251 (MQ=255)                                                                                                                                                                                                     attattACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGCGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGAAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGgcgc                                               >  2:180289/1‑250 (MQ=255)                                                                                                                                                                                                       tagtaCGAGATTTTAGGCGTTTACAAAACAGCGGAAGAGCAAGAAATAAAAAAGGCCTACAAACGACTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGGTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTAGGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAg                                             <  2:332964/247‑1 (MQ=255)                                                                                                                                                                                                        attaCGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTAACAAAGAGGCCGAGGCGAACTTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGCCCGACTCGCACAAACGTGCGGCATCCGATCAGTAGGGTCTTACTGCGTTTGAGAGATGTGGTATGGTCGGCGGCGGTTTTgggggcggggccgac                                           >  2:497889/1‑244 (MQ=255)                                                                                                                                                                                                            cGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGACTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGCAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACTATCAGTATGGTAATGCTGAGTTTGAGCAAGGTGGTCGCGGCGGCGGCGGATTTGGCGGCGgcgccgtattct                                       >  2:441951/1‑243 (MQ=255)                                                                                                                                                                                                            cGAGATTGTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCa                                       <  1:452497/251‑1 (MQ=255)                                                                                                                                                                                                             gagaTTTTAGGCGTTTCCAAAACAGCGGAAGAGCATTAAAACAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGATGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAg                                      <  2:473907/251‑1 (MQ=255)                                                                                                                                                                                                                    tAGGCGGTTCCAAACCAGCGGAAGCGCATGAACTCACAAAGGCCTACAAACGCCTGGCCATGAAATCCCACCCGGAACGTACCCAGGGTGACAAAGAGGCCGTGGCGAAATTTAAGGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCACAAACGTGCGGCATAAGATCAGTATGGTCATGCTGCGATTGAGCACGGAGGAATGGGCGGCGGCGGTGTTggaggcgggccacaattacgcgtgatt                               >  2:450241/1‑232 (MQ=255)                                                                                                                                                                                                                                   cAGCGGAAGAGCATGAAAGCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCTCCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTTATGCTGTGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTc                <  1:51813/251‑1 (MQ=255)                                                                                                                                                                                                                                    aGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCg               >  1:460792/1‑251 (MQ=255)                                                                                                                                                                                                                                              catGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATAtttt     >  1:346059/3‑251 (MQ=255)                                                                                                                                                                                                                                               atGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATAttttt    <  1:242948/250‑1 (MQ=255)                                                                                                                                                                                                                                                tGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGAGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTg   >  1:73375/1‑251 (MQ=255)                                                                                                                                                                                                                                                 gAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAACGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGCGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATCTTTTTGGTGACGTTTTCGGCGTTATTTTTgg  >  1:230009/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    AGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGG  >  NC_000913/13985‑14475

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