Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	339,979	0	G	A	69.0%	24.0 / 19.6	29	V71M (GTG→ATG)	yahH	putative uncharacterized protein YahH
Reads supporting (aligned to +/- strand):  ref base G (4/5);  new base A (7/13);  total (11/18)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCC  >  NC_000913/339681‑340260                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                          cATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCAACATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGt                                                                                                                                                                                                                                                                                          >  1:188190/1‑300 (MQ=255)       ccTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTANGCCTGATAAGCGTAGCGCATCAGGCANTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCAGGATGCGGCGTGAACGCCTTATCAGGCCTAAGAGTGGcgcg                                                                                                                                                                                                                                                                                    >  2:551891/1‑300 (MQ=255)                   ttCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCTGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGc                                                                                                                                                                                                                                                                       <  1:425070/300‑1 (MQ=255)                           aGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAg                                                                                                                                                                                                                                                               >  1:229092/1‑300 (MQ=255)                              cATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGt                                                                                                                                                                                                                                                            <  1:110953/300‑1 (MQ=255)                                          gTGCGAACGTTACGGCGTCTGACCTACATGTTCATTCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATNAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                                                                                                                                                                                >  1:570516/1‑300 (MQ=255)                                           tGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCNCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                                                                                                                                                                               <  2:744075/300‑1 (MQ=255)                                               aaCGTTACTGCGTCTTACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCTGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAGCGCATCAGGCagtct                                                                                                                                                                                                                                           <  2:707887/300‑5 (MQ=11)                                                    tACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAGCGCATCAGGCagtctggcgt                                                                                                                                                                                                                                      <  1:367158/300‑10 (MQ=255)                                                                                                                                                 tatgtcgtcgtacctgcatgctcatGCCGGATGCGGCATGAACGCCTTATCCGGCCTGCGAGTGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCTGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCTAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGt                                                                                                                                          <  2:172893/276‑1 (MQ=11)                                                                                                                                                           gacctacatgttcatGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAc                                                                                                                                >  2:821251/15‑300 (MQ=11)                                                                                                                                                            acctacatgtccatGCCGGATGCGGCGTGAACGCCTTATCCGTCCTACAAATGGCGCAAGAATCTGTATGCCTGATGAGCGTAGCGCATCAGGCATTTTTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCGTTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAcc                                                                                                                               <  2:227838/287‑1 (MQ=11)                                                                                                                                                             cctacatgttcatGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAAccc                                                                                                                              <  1:354733/288‑1 (MQ=11)                                                                                                                                                                     ttcatGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCAccc                                                                                                                      <  1:154228/296‑1 (MQ=255)                                                                                                                                                                              gATGCGTCGTGAACGCCTTATCCTGCCTACGAGTGGCGCAAGAATCTGGAGGCCTTATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCTAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTTGCGTTGGTCATAACCCCATCAccccctctaccc                                                                                                             <  2:243820/300‑10 (MQ=11)                                                                                                                                                                               aTGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCgg                                                                                                            <  1:545668/300‑1 (MQ=255)                                                                                                                                                                                      gTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGNAGCGCATCAGGCATTTGTCACCATTGNCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATaac                                                                                                     <  2:561384/300‑1 (MQ=255)                                                                                                                                                                                      gTGAACGCCCTATCCTGCCTACGAGTGTTGTAAGAGTTTTGAGGCCTGATAAGCTTAGCGCATCAGGCATTGGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGTCGCAACAATCTGTAGGCCTGATAAACGTAGCGCATCAGGCATTTGTCACAATTGCCGGATGCGGCAGGAACACCTTATCCGGCCTACGATTTGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTTATAACCCCATCACCCTCTGTAGCGGACATaac                                                                                                     <  2:16999/300‑1 (MQ=14)                                                                                                                                                                                                tATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCAACATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGGATCTGTAGGCCTGATAAGCGTAACGCCTCAGGCCTTTGTCAACATTGCCGGATTCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGGATTGTAGGCCTGATAAGCGTAGTGCATCAGGCAATCCGGCGTTGGTCATAACCCAATCACAATCTGTAGCGGCAATAACAAACTGAAag                                                                                           >  2:454016/1‑300 (MQ=11)                                                                                                                                                                                                          aCGAACGGCGCAATAATCTTTAGGCCGGATAAGCGTAGCGAATCAGGCATTGGTCACCATTGCCGGAGGCGGAGGGAACGCCTCATCCGGCCTACGAATGGCACAAGAACCTTTATGCCTGATAAGCGTAGTGCATCAGGCATGTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGTGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGt                                                                                 <  2:225962/300‑1 (MQ=21)                                                                                                                                                                                                                 gcgcAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGcc                                                                          >  1:356987/1‑300 (MQ=255)                                                                                                                                                                                                                        aaTCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAAGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCACATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATtggtg                                                                   >  1:243872/1‑300 (MQ=255)                                                                                                                                                                                                                                ggCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCGTTGGTGGCAACAGc                                                           >  1:124386/1‑300 (MQ=255)                                                                                                                                                                                                                                            gTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGAt                                               <  1:17984/300‑1 (MQ=255)                                                                                                                                                                                                                                                 gcATCANGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTNCGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACgc                                          <  2:689798/300‑1 (MQ=255)                                                                                                                                                                                                                                                      aGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCTGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACgccagcc                                     <  1:363772/300‑1 (MQ=255)                                                                                                                                                                                                                                                                              gcggcgTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGAGAAGCGTAGCGCATCAGGCGTTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCTGCCTACGAGGGGCGCGAGAATTGTATGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGTTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCaggcggagg             <  2:546309/300‑1 (MQ=25)                                                                                                                                                                                                                                                                                   gtgaaCGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGGTAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTg        >  1:632280/5‑300 (MQ=255)                                                                                                                                                                                                                                                                                         gCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCATACGAGTGGCGGGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCAATCACCCTATGTAGCGGACATAACAAAATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGAGTGGGCATCAGGCGGAGGCGGTGAAGGcc  >  2:9268/1‑300 (MQ=25)                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                          CATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCGATTGAGCATCAGGCGGAGGCGGTGAAAGCC  >  NC_000913/339681‑340260

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