Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	646,777	0	G	T	40.4%	68.5 / 27.8	52	Q245K (CAA→AAA)	citG	2‑(5''‑triphosphoribosyl)‑3'‑dephosphocoenzyme‑A synthase
Reads supporting (aligned to +/- strand):  ref base G (9/22);  new base T (21/0);  total (30/22)
Fisher's exact test for biased strand distribution p-value = 7.80e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑TTCGAA‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCA  >  NC_000913/646627‑646926                                                                                                                                                         |                                                                                                                                                      ataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGACCTGCCGGAGATAATCGAGATCGGCGGGGG‑TTCGAA‑TGCCCCCTTTTTGCAATAATGTTTg                                                                                                                                                       <  1:804410/151‑1 (MQ=255)    atTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑TTCGAA‑TGCCCCCTTTTTGCAATAATGTTTgcgc                                                                                                                                                    <  1:467184/151‑1 (MQ=255)       aaaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGGCTGAGATTTCGTTCGATACACTCCCTTTCGAACTGCCGGAGATAATCGAGATCGGCGTGGG‑TTCGAA‑TGCACCCTTTTTGCAATAATGTTTGCGCCt                                                                                                                                                  <  1:653627/150‑1 (MQ=255)       aaaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑TTCGAA‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTc                                                                                                                                                 <  1:608234/151‑1 (MQ=255)       aaaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑GTCGAA‑TGCCCCCCTTTTTCAATAATGTTTTCGcccc                                                                                                                                                 >  1:893906/1‑149 (MQ=255)               gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑TTCGAA‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTa                                                                                                                                         <  2:298213/151‑1 (MQ=255)                 tAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑GTCGAA‑TGCCCCCCTTTTTCAATAATGTTTTCGCCCCccgc                                                                                                                                             >  1:499719/1‑145 (MQ=255)                 tAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑GTCGAA‑TGCCCCCCTTTTTCAATAATGTTTTCGCCCCCCGCCGTAg                                                                                                                                        >  1:224902/1‑150 (MQ=255)                 tAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑GTCGAA‑TGCCCCCCTTTTTCAATAAAGTTTTCGCCTCGCGCTGTAGc                                                                                                                                       >  1:587552/1‑151 (MQ=255)                  aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑TTCGAA‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAg                                                                                                                                        <  1:425358/149‑1 (MQ=255)                       ccAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑GTTGGA‑AGCCCCCCTTTTTCCATAATTTTTTCGcccccccccgttgcc                                                                                                                                      >  1:354248/1‑133 (MQ=255)                        cAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑GTTGGA‑AGCCCCCCCTTTCCCATAAACTTTTGGCCgcgcgcgccag                                                                                                                                        >  2:679754/1‑138 (MQ=255)                             aaGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGCGATCGGCGGGGG‑TTGGGA‑GGGGCCCTCTTTTTAATAATATGTTTTCCCCGGCCCGTAGCCAGCGCAGGcc                                                                                                                            >  2:629092/1‑150 (MQ=255)                             aaGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑GTCGGA‑TGCCCCCCTTTTTCAATAAAGTTTTCTCCCCCCGCCGGAGccccacccggcac                                                                                                                           >  2:722498/1‑140 (MQ=255)                             aaGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑GTCGAAGTGCCCCCTTTTT‑TAATAAAGTTTTCCCCTCCCGCCGTAGGCAGCCCAGGccc                                                                                                                           >  2:249506/1‑151 (MQ=255)                                gATCAGTAGGTCATCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGGCGGAGATAATCGAGATCGGCGGGGG‑TTCGAA‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGcccacc                                                                                                                        <  2:337022/151‑1 (MQ=255)                                        ggTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑TTCGAA‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCAcc                                                                                                                <  1:817981/151‑1 (MQ=255)                                        ggTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑TTCGAA‑TGACCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCAcc                                                                                                                <  1:495567/151‑1 (MQ=255)                                           cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGATTGGAG‑GGCCCCCTTTTT‑CAATAATTTTTTCGCCTCCCCCTGTTGGCCCCGCCGGcccccccccccccccc                                                                                                              >  1:740660/1‑140 (MQ=255)                                           cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑TTCGAA‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcgc                                                                                                              <  2:847618/150‑1 (MQ=255)                                              cACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑TTCGAA‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGa                                                                                                          <  1:877821/151‑1 (MQ=255)                                                   ccgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑GTCGAA‑TGCCCCCTTTTTGTAATAATATTTTTGCCCCCCGCCGTTGACAGCGCAGGCCCCCCCCTCCCCCCCCCGATGc                                                                                                       >  2:69061/1‑149 (MQ=255)                                                     gccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑GTCGAA‑TGCCCCCCTTTTTCAATAATGTTTTCGCCCCCCGCTGTAGCCCGCGCAGGCCCCCCCCCCCACCCCCCGATGCAAc                                                                                                    >  1:342902/1‑150 (MQ=255)                                                            cTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑GTCGAA‑TGCCCCCCTTTTTCAATAAAATTTTCGCCTCGCGCTGGAGCCCGCGCAGGCCCCCCTCGCCACCCCCCGCTGCAACGCTggga                                                                                             >  1:237258/1‑148 (MQ=255)                                                               agaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑GTCGAA‑TGCCCCCCTTTTTCCATAAATTTTTCGCCCCGCGCCGTAGCCCGCGCAGGGCCCCCCCCCCCCCCCGCGATGCAACGGTgg                                                                                               >  1:884702/1‑145 (MQ=255)                                                                 aTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑GTCGAA‑TGCCCCCCTTTTTCAATAATGTTTTCGCCCCGCGCTGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCGCGATGCAACGTTGGGATCGcc                                                                                        >  1:945304/1‑150 (MQ=255)                                                                  tttCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGAGG‑GTCGCA‑TTCCCCCTCTTTTCCATAATGTTTTCGCCTCGCGCCGTAAACCGCGCCGGACCCCCCCGCCCCCGCACGCAGCAACGCTgggttcccc                                                                                        >  1:224265/1‑142 (MQ=255)                                                                   ttcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑GTTGAA‑TGCCCCCCTTTTTCAATAATTTTTTCGCCCCCCGCTGTAGCCAGCGCAGGCCCCCCCCGCCACCCCCCGATGCAACGTTGGGATCGcc                                                                                        >  1:610625/1‑148 (MQ=255)                                                                   ttcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑GTCTAA‑TTCCCCCCTTTTTCAATAAATTTTTTGCCCCCCGCCGTAGCCCGCCCCCGCGCCCCCCCCCCCCCCGCCATGCAACGTTGGGGTCGCCggt                                                                                     >  2:8420/1‑149 (MQ=255)                                                                           aTACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑GTCCCA‑TGACCCCCTTTTTCAATAATGTTTTTGCCCCCCGCCGTAGCCAGCGCCGGGCCCCCCCCCCACCCCGCCGGGCAACGTTGGGATATCCGCTGATAGccc                                                                             >  2:916846/1‑150 (MQ=255)                                                                             acacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑TTCGAA‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCAt                                                                            <  1:540766/150‑1 (MQ=255)                                                                                 tCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑GTCGGA‑ATCCCCCCTTTTCTAATCAAATTTTCGCCTCGCGCCGTTGCCAGCGCCGGCCCCCCCCCCCACCCCGCGGGGCCACGCTGGGATCGccccggctcgcccccc                                                                          >  2:429336/1‑134 (MQ=255)                                                                                  cccTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑GTCGAA‑TGCCCCCCTTTTTCAATAAAGTTTTCGCCCCCCGCTGTAGCCCGCGCAGGCCCCCCCCCCCACCCCGCGGGGCCACGGTGGTATCGCCGTTGATCGcc                                                                              >  1:255148/1‑143 (MQ=255)                                                                                                  gagaTAATCGAGATCGGCGGGGG‑TTCGAA‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGc                                                      <  2:420263/151‑1 (MQ=255)                                                                                                   agaTAATCGAGATCGGCGGGGG‑TTCGAA‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGc                                                      <  2:128868/150‑1 (MQ=255)                                                                                                       aaTCGAGATCGGCGGGGG‑TTCGAA‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATg                                                  <  1:325817/150‑1 (MQ=255)                                                                                                            agaTCGGCGGGGG‑TTCGAA‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGtt                                            <  2:71945/151‑1 (MQ=255)                                                                                                             gaTCGGCGGGGG‑TTCGAA‑TGCCCCCTTTTTTCCATAATGTTTGCGCCCCGCGCTGTAACCAGCGCCGGCCCCCCTCGCCCCCGCGCGATGCAACGTTGGGATCGCCGTTGATCGCCATCCGTAGGGGCAAGGGATCGAGGAATggccggt                                            >  2:584099/1‑144 (MQ=255)                                                                                                                cggcggGGG‑TTCGAA‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAgg                                        <  1:99818/151‑1 (MQ=255)                                                                                                                  gcggGGG‑TTCGAA‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGAt                                      >  1:67939/1‑151 (MQ=255)                                                                                                                   cggGGG‑TTCGAA‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATc                                     >  1:160711/1‑151 (MQ=255)                                                                                                                    ggggg‑TTCGAA‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt                                    >  1:340605/1‑151 (MQ=255)                                                                                                                        g‑TTCGAA‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAcc                                <  1:426686/151‑1 (MQ=255)                                                                                                                            cGAA‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCt                             <  1:799014/151‑1 (MQ=255)                                                                                                                             gAA‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTg                            >  2:179529/1‑151 (MQ=255)                                                                                                                             gAA‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTg                            <  1:485308/151‑1 (MQ=255)                                                                                                                              aa‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGa                           <  1:344571/151‑1 (MQ=255)                                                                                                                                    ccccTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATcc                        >  1:716453/1‑149 (MQ=255)                                                                                                                                      ccTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCcagca                    <  1:96841/151‑1 (MQ=255)                                                                                                                                         ttttGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCcagcag                   <  1:304811/149‑1 (MQ=255)                                                                                                                                             gCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCCCCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCatag                 <  2:471617/147‑1 (MQ=255)                                                                                                                                                ataatGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCACTGATCCAGCAGAGTGa              >  2:336556/1‑147 (MQ=255)                                                                                                                                                 taatGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTa          >  2:526515/1‑150 (MQ=255)                                                                                                                                                  aatGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAAt        >  2:439190/1‑151 (MQ=255)                                                                                                                                                         gcgcCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCa  <  1:419812/150‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                      ATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGG‑TTCGAA‑TGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCA  >  NC_000913/646627‑646926

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