Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,680,490	G→A	C752C (TGC→TGT)	yphG ←	DUF5107 domain‑containing protein YphG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,680,490	0	G	A	97.0%	101.8 / NA	33	C752C (TGC→TGT)	yphG	DUF5107 domain‑containing protein YphG
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (13/19);  minor base C (1/0);  total (14/19)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GTCGTTTCTCCGGTGTGGCTCCACTTAACTTATCCAGCAAATCCCGTTCGAAAAGCAGACGTGCATCCTGCGGCGCAAGCTGATAAGCATTATCAAGATAACGCGCGGCCAGCTCATAATCGTGTTGCTTATTCCACGCATGGATCGCTAACCCGCGCCAGCCGTCGGCAAACTCCGGCGACATCTCTACGCAACGTTGCCAAAAGGCAATGGCTTTGTTGTAGCTACGTTTGTTGTAGT‑AGAAGCAAGCTAGTAAA‑TG‑GCGAGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAA‑CTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCC  >  NC_000913/2680246‑2680716                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                      gTCGTTTCTCCGGTGTGGCTCCACTTAACTTATCCAGCAAAGCCCGTTCGAAAAGCAGACGTGCATCCTGCGGCGCAAGCTGATAAGCATTATCAAGATAACGCGCGGCCAGCTCATAATCGTGTTGCTTATTCCACGCATGGATCGCTAACCCGCGCCAGCCGTCGGCAAACTCCGGCGACATCTCTACGCAACGTTGCCAAAAGGCAATGGCTTTGTTGTAGCTACGTTTGTTGTAGT‑AGAAACAAGCt                                                                                                                                                                                                                                 <  1:26134/251‑1 (MQ=255)              ggtgtgtccgcttaacttaTCCAGCAAACCCCTTTCGAAAAGAAGACGGGCATCCGGCGGCGCAAGCTGATAAGCATTATCAAGATAACGCGCGGCCAGCTCATAAGCTTGTTGCTTAGTCCACGCATGGATCGCTAACCCGCGCCAGCTGTCGGCAAACTCCGGCGACATCTCTACTCAACGTTGCCAAAAGGCAATGGCTTTGTTGTAGCTACGTTTGTTGTAGT‑AGAAACAAGCTAGTAAA‑TG‑GCGAg                                                                                                                                                                                                                  <  2:128982/241‑1 (MQ=255)                              ttaTCCAGCAAATCCCGTTCGAAAAGCAGCCGTGCATCCTGCGGCGCAAGCTGATAAGCATTATCAAGATAACGCGCGGCCAGCTCATAATCGTGTTGCTTATTCCACGCATGGATCGCTAACCCGCGCCAGCCGTCGGCAAACTCCGGCGACATCTCTACGCAACGTTGCCAAAAGGCAATGGCTTTGTTGTCGCTACGTTTGTTGTAGT‑AGAAACAAGCTAGTAAA‑TG‑GCGAGCAAACCAGCACTCTTc                                                                                                                                                                                                  >  1:318898/1‑251 (MQ=255)                                          ccccGTTCGAAAAGCAGACGTGCATCCTGCGGCGCAAGCTGATAAGCATTATCAAGATAACGCGCGGCCAGCTCATAATCGTGTTGCTTATTCCACGCATGGATCGCTAACCCGCGCCAGCCGTCGGCAAACTCCGGCGACATCTCTACGCAACGTTGCCAAAAGGCAATGGCTTTGTTGTAGCTACGTTTGTTGTAGT‑AGAAACAAGCTAGTAAA‑TG‑GCGAGCAAACCAGCACTCTTCAATACTCTCCAg                                                                                                                                                                                      <  1:305180/250‑1 (MQ=255)                                                  gAAAAGCAGACGTGCATCCTGCGGCGCAAGCTGATAAGCATTATCAAGATAACGCGCGGCCAGCTCATAATCGTGTTGCTTATTCCACGCATGGATCGCTAACCCGCGCCAGCCGTCGGCAAACTCCGGCGACATCTCTACGCAACGTTGCCAAAAGGCAATGGCTTTGTTGTAGCTACGTTTGTTGTAGT‑AGAAACAAGCTAGTAAA‑TG‑GCGAGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCa                                                                                                                                                                              <  2:37530/251‑1 (MQ=255)                                                   aaaaGCAGACGTGCATCATGCGGCGCAAGCTGATAAGCATTATCAAGATAACGCGCGGCCAGCTCATAATCGTGTTGCTTCTTCCACGCATGGATCGCTAACCCGCGCCAGCCGTCGGCAAACTCCGGCGACATCTCTACGCAACGTTGCCAAAAGGCAATGGCTTTGTTGTCGCTACGTTTGTTGTAGT‑AGAAACAAGCTAGAAAA‑TG‑GCGAGCAAACCCGCACTCTTCAATACTCTCCACCGCCGCCAc                                                                                                                                                                             >  2:274157/1‑251 (MQ=255)                                                            cGTGCATCCTGCGGCGCAAGCTGATAAGCATTATCAAGATAACGCGCGGCCAGCTCATAATCGTGTTGCTTATTCCACGCATGGATCGCTAACCCGCGCCAGCCGTCGGCAAACTCCGGCGACATCTCTACGCAACGTTGCCAAAAGGCAATGGCTTTGTTGTAGCTACGTTTGTTGTAGT‑AGAAACAAGCTAGTAAA‑TG‑GCGAGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTcc                                                                                                                                                                      >  1:108874/1‑249 (MQ=255)                                                                          cgcAAGCTGATAAGCATTATCAAGATATCGCGCGGCCAGCTCATAATCGGGTTGCTTATTCCTCGCATGGATCGCTAACCCGCGCCAGCCGTCGGCAAACTCCGGCGACATCCCTACGCAACGTTGCGAAAAGGCAATGGCTTTGTTGTAGCTACGTTTGTTGTAGT‑AGAAACAAGCTAGTAAA‑TG‑GCGAGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAAc                                                                                                                                                      <  1:246116/251‑1 (MQ=255)                                                                                      aGCATTATAAAGATAACGCGCGGCCAGCTCATAATCGTGTTGCTTATTCCACGCATGGATCGCTCACCCGCGCCAGCCGTCGGCAAACTCCGGCGACATATCTACGCAACGATGCCAAAAGGCAATGGCTTTGTTGTAGCTACGTTTGTTGTAGT‑AGAACCAAGCTAGTAAA‑TG‑GCGAGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAA‑CTGCg                                                                                                                                         >  2:375115/1‑251 (MQ=255)                                                                                             acAAGATAACGCGCGGCCAGCTCATAATCGTGTTGCTTATTCCACGCATGTATCGCTACCCCGCGCCAGCCGTCGGCAAACTCCGGCGACATCTCTACGCAACGTTGCCAAAAGGCAAGGGCTTTGTTGTAGCTACGTTTGTTGTAGT‑AGAAACAAGCTAGTAAT‑TG‑GCGAGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAA‑CTGCGGGAAGAc                                                                                                                                  <  1:38426/250‑1 (MQ=255)                                                                                               aaGATAACGCGCGGCCAGCTCATAAGCGTGTTGCTTATTCCACGCATGGATCGCTAACCCGCGCCAGCCGTCGGCAAACCCCGGCGACATCTCTACGCAACGTTGCCAAAAGGCAATGGCTTTGTTGTAGCTACGTTTGTTGTAGT‑AGAAACAAGCTAGTAAA‑TG‑GCGAGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAA‑CTGCGGGAAGACAt                                                                                                                                <  1:286397/251‑1 (MQ=255)                                                                                                              cAGCTCATAAGCGTTTTGCTTATTCCACGCTGGGACCGCTAACCCGCGCCAGCGGGCGGCAAACTCCGGCGACATCTCTACGCAACGTTGCAAAAAGGCAATGGCTTTGTTGTAGCTACGTTTGTTGTAGT‑AGAAACAAGCTAGTAAA‑TG‑GCGAGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAA‑CTGCGGGAAGACATCAATGGCTTTTGCGa                                                                                                                 <  2:108874/251‑1 (MQ=255)                                                                                                               aGCTCATAACCGCGTTGCTTATACCACGCTTGGAGCGCAATCCCGCGCCAGCCGTCGGCAAACTCCGGCTACATCTCTACGAAATGTTGCCAAAAGGCAATGGCTTTGTTGTATCTACGTTTGTTGTAGT‑AGAAACAAGCTAGTAAA‑GG‑GCGATCAAACCAGCACTCTCCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAA‑CTGCGGGAAGACATCAATGGCTTTTGCGAc                                                                                                                <  2:112294/251‑1 (MQ=255)                                                                                                                gCTCATAATCGTGTTGCTTATTCCACGCATGGATCGCTAACCCGCGCCAGCCGTCGGCAAACTCCGGCGACATCTCTACGCAACGTTGCCAAAAGGCAATGGCTTTGTTGTAGCTACGTTTGTTGTAGT‑AGAAACAAGCTAGTAAA‑TG‑GCGAGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAA‑CTGCGGGAAGCCCTCAATGGCTTTTGCGAcc                                                                                                               >  2:334082/1‑251 (MQ=255)                                                                                                                   cagaaTCGTGTGGCTTATCCCACGCTTGGATCGCTACCCCGCGCCAGCCGTCGGCAAACTCCGGCGACTTCTCTACGCAACGTTGCCAAAAGGCAATGGCTTTGTTGTAGCTACGTTTGTTGTAGT‑AGAAACAAGCTAGTAAA‑TG‑GCGAGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAA‑CTGCGGGAAGACATCAATGGCTTTTGCGACCAg                                                                                                             <  2:46238/247‑1 (MQ=255)                                                                                                                         cGTGTTGCTTATTCCACGCATGGATCGCTAACCCGCGCCAGCCGTCGGCAAACTCCGGCGACATCTCTACGCAACGTTGCCAACAGGCTATGGCTTTGTTGTAGCTACGTTTGTTGTAGT‑AGAAACAAGCTAGTAAA‑TG‑GCGAGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGCAACGGACAAA‑CTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGcc                                                                                                       >  2:286841/1‑250 (MQ=255)                                                                                                                                                aTCGCTCACCCGCGCCAGCCGTCGGCAAACTCCGGCGACAACTCTACGCAACGTGGCCAAAAGGCAATGGCTTTGTTGTAGCTACGTTTGTTGTAGT‑AGAAACAAGCTAGTAAA‑TG‑GCGAGCAAACCAGCCCTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAA‑CTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCt                                                                               <  2:79957/251‑1 (MQ=255)                                                                                                                                                            cgcCAGCCGTCGGCAAACTCCGGCGACATCTCTACGCAACGTTGCCAAAAGGCAATGGCTTTGTTGTAGCTACGTTTGTTGTAGT‑AGAAACAAGCTAGTAAA‑TG‑GCGAGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAA‑CTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTa                                                                   <  1:59335/251‑1 (MQ=255)                                                                                                                                                                    gTCGGCAAACGCCGGCGACATCTCTACGCACCGTTGCCAAAAGGCAATGGCTTTGTTTAAGCTACGTTTGTTGTAGT‑AGAAACAAGTTACTAAAGTGCGCGA‑CAAAC‑AGCACTCTTCAATACTCTCCAGCGCCGCCACTTCCTCCAGCGTATTCGGGAACCGGCCAAA‑CAGCGGGAAGACATCAATGGTTCTTGCGACCAGTTCTCCAAGTACGGCTTTCGGCAGCAAgatcatttgcaggtaacgcgccc                                                           >  2:53079/1‑229 (MQ=255)                                                                                                                                                                       ggCAAACTCCGGCGACATCTCTACGCAACGTTGCCAAAAGGCAATGGCTTTGTTGTAGCTACGTTTGTTGTAGT‑AGAAACAAGCTAGTAAA‑TG‑GCGAGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAA‑CTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCg                                                        <  1:119819/251‑1 (MQ=255)                                                                                                                                                                         cAAACTCCGGCGACCTCTCTACGCAACGTTGCCAAAAGGCAATGGCTTTGTTGTAGCTACGTTTGTTGTAGT‑AGAAACAAGCTAGTAAA‑TG‑GCGAGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAA‑CTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTg                                                      >  2:240944/1‑251 (MQ=255)                                                                                                                                                                                  gCGACATCTCTACGCAACGTTGCAAAAAGGCAATGGCTTTGTTGTAGCTACGTTTGTGGTAGT‑AGAAACAAGCTAGTAAA‑TG‑GCGAGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAA‑CTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAg                                             <  1:334082/251‑1 (MQ=255)                                                                                                                                                                                    gACATGTCTAAGCAACGTTGCCAAAAGGCAATGGCTTTGTTGTAGCTACGTTTGTTGTAGT‑AGAAACAAGCTAGTAAA‑TG‑GCGAGCAAACCAGAACTCTTCAATACTATCCAGCGCCGCAACTTCTTCCAGCGTATTCGGGAAACGGACAAA‑CTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGACGAAAGCTGGCTTGAAGGTAAAGGGGCAGCGGGGGCAGGCAGTc                                           >  2:168297/1‑251 (MQ=255)                                                                                                                                                                                             aCGCAACGTTGCCAAAAGGCAATGGCTTTGTTGTAGCTACGTTTGTTGTAGT‑AGAAACAAGCTAGTAAA‑TG‑GCGAGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAA‑CTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGtt                                  <  1:80479/251‑1 (MQ=255)                                                                                                                                                                                              gtcAACGCTGCCAAAAGGCAATGGCTTTGTTGTAGCTACGTTTGTTATAGT‑AGAAACAAGATAGTAAA‑TG‑GCGCGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGAATTCGGGAAACG‑ACAAAACTGAGGGAAGACATCAATGTCTTTTGCGACAAGTTCACCAGGTTCGGCTTTCGGAAGCAAGCTGGCTTGCATGTAAAGCGGCAGGGTGAGCTggagaacccgttcgttc                                 >  2:206626/3‑236 (MQ=255)                                                                                                                                                                                                 aacgtggaaaaaaGGAAATGGCTTTTTTGTAGTTACGTTTTTTGT‑GTGAGAAACAAGCTAGTAAT‑TG‑GCGAGCAACCCAGCACTCTCCAATACCCCCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAA‑CTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGc                              <  2:241432/242‑1 (MQ=255)                                                                                                                                                                                                      tGCCAAAAGGCAATGGCTTTGTTGTAGCTACGTTTGTTGTAGT‑AGAAACAAGCTAGTAAA‑TG‑GCGAGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAA‑CTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCAtctc                         >  1:292534/1‑251 (MQ=255)                                                                                                                                                                                                       gCCAAAAGGCAATGGCTTTGTTGTAGCTACGTTTGTTGTAGT‑AGAAACAAGCTAGTAAA‑TG‑GCGAGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAA‑CTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCAtctc                         <  1:73040/250‑1 (MQ=255)                                                                                                                                                                                                         cAAAAGGCAATGGCTTTGTTGTAGCTACGTTTGTTGTAGT‑AGAAACAAGCTAGTAAA‑TG‑GCGAGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAA‑CTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTgc                      >  1:65281/1‑251 (MQ=255)                                                                                                                                                                                                           aaaGGCAATGGCTTTGTTGTAGCTACGTTTGTTGTAGT‑AGAAACAAGCTAGTAAA‑TG‑GCGAGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAA‑CTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTg                                                 >  1:347804/1‑222 (MQ=255)                                                                                                                                                                                                             aGGCAATGGCTTTGTTGTAGCTACGTTTGTTGTAGT‑AGAAACAAGCTAGTAAA‑TG‑GCGAGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAA‑CTGCGGGCAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGACTTGAAGGTAAAGCGGCACAGGGCGCTGGCGGTcacgggcggtccgcatctttgacccc                  >  2:339060/1‑226 (MQ=255)                                                                                                                                                                                                                     gCTTTGTTGTAGCTACGTTTGTTGTAGT‑AGAAACAAGCTAGTAAA‑TG‑GCGAGCAAACCAGCACTCTTCAATACGCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAA‑CTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTggg          <  2:261659/251‑1 (MQ=255)                                                                                                                                                                                                                             gtatcTACGTTTGTTGTAGT‑AGAAACAAGCTAGTAAA‑TG‑GCGAGCAAACCAGCCCTCTTCAATACCCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAA‑CTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTcccc  <  1:240944/247‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                      GTCGTTTCTCCGGTGTGGCTCCACTTAACTTATCCAGCAAATCCCGTTCGAAAAGCAGACGTGCATCCTGCGGCGCAAGCTGATAAGCATTATCAAGATAACGCGCGGCCAGCTCATAATCGTGTTGCTTATTCCACGCATGGATCGCTAACCCGCGCCAGCCGTCGGCAAACTCCGGCGACATCTCTACGCAACGTTGCCAAAAGGCAATGGCTTTGTTGTAGCTACGTTTGTTGTAGT‑AGAAGCAAGCTAGTAAA‑TG‑GCGAGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAA‑CTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCC  >  NC_000913/2680246‑2680716

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