Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	71,555	A→C	M69L (ATG→CTG)	yabI →	DedA family protein YabI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	71,555	0	A	C	97.1%	112.8 / ‑4.3	34	coding (205/765 nt)	yabI	DedA family protein YabI
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (5/28);  minor base . (1/0);  total (6/28)
Fisher's exact test for biased strand distribution p-value = 1.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTATTGGCTG‑CTTGATGGGCGACTGGATTTCTTTCTGGCTGGGTTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCG  >  NC_000913/71316‑71779                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                  cAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTgg                                                                                                                                                                                                                        <  1:125986/251‑1 (MQ=255)           gCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATctccttct                                                                                                                                                                                                              <  1:56336/251‑1 (MQ=255)               cGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTctggct                                                                                                                                                                                                          <  1:227799/251‑1 (MQ=255)                       tGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGgctgggctgg                                                                                                                                                                                                   <  1:178322/250‑1 (MQ=255)                                               gCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGACTGGCG‑TTTTAAAAAGCCGTTGCATCGCt                                                                                                                                                                         <  2:65731/251‑1 (MQ=255)                                                  ggAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGt                                                                                                                                                                      <  2:336469/251‑1 (MQ=255)                                                   gAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTc                                                                                                                                                                     <  2:247733/251‑1 (MQ=255)                                                     acacTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCAt                                                                                                                                                                   <  1:278090/251‑1 (MQ=255)                                                       acTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGAATGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCAttt                                                                                                                                                                 <  2:251761/251‑1 (MQ=255)                                                          ttattaCCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGGTTCTGGGACGCCTGGGTGGGCGGGGTTGTCGGGTGTCTTG‑TGGGGGACTGGGTCTCCCTCTTGCTGGGCTGGCGGTTTTAAAAA‑CCGGTGCATCGCTGGGCAtttttt                                                                                                                                                              >  1:95307/1‑248 (MQ=255)                                                                        cccGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGc                                                                                                                                                <  1:3238/250‑1 (MQ=255)                                                                        aCCGTGTATGCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAg                                                                                                                                                 <  2:52009/250‑1 (MQ=255)                                                                                     tGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGGTGGGAGCGCTGATTGGCAGCGGCGAGGTAAGGTTTTGGCACGCCTGGGTGGGAGGGGTTGTCGGGTG‑TCTGCTGGGGGACTGGGTCTCCTTTTTGCTGGGCGGGCGGTTTTAAAAA‑CCGTTGCATCGCTGGTCATTTCTGAAGAaaaaaaaaacattaatcgataa                                                                                                                                    >  1:50887/1‑236 (MQ=255)                                                                                       aTGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTTCTCGATAAAAc                                                                                                                                 <  2:29931/251‑1 (MQ=255)                                                                                               ggggtTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGGTTGAATCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCGGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCAGGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGTTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTTCTCGATAAAACTGAACATg                                                                                                                         <  2:189758/248‑1 (MQ=255)                                                                                                    gggtggCCGTTCTGGGGTCGCGGGTGCTGGTCGGTTTGATTCGACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCGCGTCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGTGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGTTCAGTGCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACAGGCGTTg                                                                                                                    <  2:42702/250‑1 (MQ=255)                                                                                                        ggCCTTTCTGGAGTCGCTGGCGCTGGTCGGTGTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATc                                                                                                                <  1:4794/251‑1 (MQ=255)                                                                                                          ccTTTCTGGAGTCGCTGGCGCGGGTCGGTTGGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCATGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTGCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCaa                                                                                                              <  2:239647/251‑1 (MQ=255)                                                                                                                      cgctggcgctggTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGTATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACTTGCGTTGCATCAACACAGCATGtt                                                                                                   <  2:183737/250‑1 (MQ=255)                                                                                                                           gcgctggTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCAt                                                                                             <  1:13744/251‑1 (MQ=255)                                                                                                                            cgctggTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCAtt                                                                                            <  1:322979/251‑1 (MQ=255)                                                                                                                                      tttGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTgggcggtc                                                                                  >  1:40386/1‑251 (MQ=255)                                                                                                                                               ccccGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTttgttg                                                                         <  2:50887/250‑1 (MQ=255)                                                                                                                                               aCCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGATGGCG‑TTTTAAAAAGCCGGTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTttg                                                                            <  1:238961/248‑1 (MQ=255)                                                                                                                                                   ggTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAAc                                                                                                                                                     <  1:291859/171‑1 (MQ=255)                                                                                                                                                              gATGGCGGGGCTGGGGGCGCTGATTGGCAGCGGCGATTTAAGTTTCTTGCACGCCTGGCTGGCATGGATTGTCGGTTT‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTTGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCTATAAAACTGACCATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGGTGGCCCGActcatct                                                            <  2:231897/249‑7 (MQ=255)                                                                                                                                                                  gCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGt                                                      <  1:70493/251‑1 (MQ=255)                                                                                                                                                                     gggCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGcc                                                   <  1:186286/251‑1 (MQ=255)                                                                                                                                                                                             ggcgAGTTAAGTTTCTGGCACGCCTGGCTGGCCGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCt                           >  2:167354/1‑251 (MQ=255)                                                                                                                                                                                                    tAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTg                    >  1:91042/1‑251 (MQ=255)                                                                                                                                                                                                             ggCACGCCTGGCTGGCAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAAt             <  1:11003/249‑1 (MQ=255)                                                                                                                                                                                                               cACGCCTGGCTGGCAGGGATGGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAAttcat         <  2:330176/251‑4 (MQ=255)                                                                                                                                                                                                                     tggctggcAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACACGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGcc   >  1:324192/1‑251 (MQ=255)                                                                                                                                                                                                                      ggctggcAGGGATTGTCGGTTG‑TCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCg  <  1:45257/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                  CAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTATTGGCTG‑CTTGATGGGCGACTGGATTTCTTTCTGGCTGGGTTGGCG‑TTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCG  >  NC_000913/71316‑71779

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