Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,434,882	0	T	G	48.2%	60.7 / 15.3	56	Y252S (TAC→TCC)	truA	tRNA pseudouridine(38‑40) synthase
Reads supporting (aligned to +/- strand):  ref base T (10/19);  new base G (27/0);  total (37/19)
Fisher's exact test for biased strand distribution p-value = 5.24e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCA  >  NC_000913/2434736‑2435029                                                                                                                                                   |                                                                                                                                                    gTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCAGTCAGGGtga                                                                                                                                                   >  1:669107/1‑147 (MQ=255) gTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATCATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGGTCATACCGGTCAGGGGaa                                                                                                                                                   >  1:872047/1‑149 (MQ=255)   ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAACCGTCAGGGGAATcc                                                                                                                                                >  2:339310/1‑150 (MQ=255)   ccATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATcc                                                                                                                                                <  1:986239/150‑1 (MQ=255)      tatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGGCAGGGGAATCCAcc                                                                                                                                             >  2:1021965/1‑150 (MQ=255)       atTTGCCTATCCCGTGTACTCTTCATCTTGCAGACTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGGCAGAAATGGCGGGACGATTGGCGGTTTTGTAAGATCATACCTGTCTGGGTAATCCAc                                                                                                                                              <  1:857286/148‑1 (MQ=255)              tATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGGAATCCACCGCGAcc                                                                                                                                       >  2:17691/1‑148 (MQ=255)               aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAg                                                                                                                                     <  2:799214/149‑1 (MQ=255)               aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGAt                                                                                                                                   <  1:397494/151‑1 (MQ=255)                 gCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGCAGGTCATACCCGGCAGGGGAATCCAACGCGACCAGa                                                                                                                                    >  2:821800/1‑148 (MQ=255)                 gCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAACGGGCAGGGGAATCCACCCCGGCCAGATa                                                                                                                                  >  2:642969/1‑150 (MQ=255)                        aCTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCtt                                                                                                                          <  2:702089/151‑1 (MQ=255)                         cTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGTCAGGGGAAATCACCCCGACCAGATACAGCCCtt                                                                                                                          >  1:964882/1‑150 (MQ=255)                         cTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATACAGCCCtt                                                                                                                          >  1:189204/1‑150 (MQ=255)                                tCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAAATCATACCGGTGAGGGGAATCCAACCCGACCAGATACAGCCCTTCCCCtttt                                                                                                                  >  1:255312/1‑151 (MQ=255)                                 cTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATAACGGGCAGGGGAATCCACCGCGACCAGATACAACCCTTCCGCtttt                                                                                                                  >  2:824538/1‑150 (MQ=255)                                   ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTTcc                                                                                                               >  2:204628/1‑151 (MQ=255)                                    tCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGGCAGGGGAATCCACCGCGACCCGATACAGCCCCTCCGCTTTTggc                                                                                                               >  2:723871/1‑148 (MQ=255)                                      aGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGtt                                                                                                            <  2:430901/151‑1 (MQ=255)                                         cTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGTAATCCACCCCGGCCAGATACAGCCCCTCCGCTTTTGCCGtt                                                                                                            >  1:275275/1‑148 (MQ=255)                                         cTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGGCCGGGGAATCCACCGCGACCAGATACAACCCTTCCGCTTTTTCCggt                                                                                                            >  2:373233/1‑146 (MQ=255)                                                   gTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATAACGGTCAGGGGAATCCACCGCGACCAGATACAACCCCTCCGCTTTTGCCGTTGCTGCCCCCAGc                                                                                                >  2:958548/1‑150 (MQ=255)                                                   gTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGGAATCCACCGCGACCCGATACAACCCTTCCCCTTTTTCCGTTGCTgccccc                                                                                                   >  2:313857/1‑147 (MQ=255)                                                            aGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGTAGGTTATACCCGTCCAGGGAATCCACCCCGCCCCGATACAGCCCCTCCTCCTTTGCCCTTGCTGccgcccgccgtttgt                                                                                         >  1:74757/1‑141 (MQ=255)                                                            aGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAATCCACCGCGGCCAGATACAGCCCCTCCGCCTTTTCCGTTGCTGCCGCCAGCGGTCTTTcc                                                                                       >  2:649070/1‑150 (MQ=255)                                                               cTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCGGCAGGGGAATCCACCCCGACCAGATACAACCCTTTCCCCTTTTCCGGTGCTTCCCCCCCCGGTCTTTCCttt                                                                                    >  1:408812/1‑150 (MQ=255)                                                                tGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAACCGGCAGGGGAATCCACCGCGACCAGATACAACCCCTCCGCTTTTTGCGTTTCTGCCCCCCGCGTTCTGTCCttt                                                                                    >  2:761343/1‑149 (MQ=255)                                                                   aaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTgccgc                                                                               <  2:677771/151‑1 (MQ=255)                                                                      gTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCa                                                                             <  2:282620/150‑1 (MQ=255)                                                                              cATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAATCCACCCCGACCAGATACAGCCCCTCCCCTTTTTCCGTTGCTGCCGCCAGCGTTCTGTCCTTTTCCGCCCGCaa                                                                         >  2:600860/1‑146 (MQ=255)                                                                               attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGACAGCAACTCTGc                                                                   <  1:834862/151‑1 (MQ=255)                                                                               attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGGCAGGGGAATCCACCGCGACCAGATACAGCCCTTCCGCCTTTGCCGTTGCTGCCGCCCGCGTTCTTTCCTTTTCCGCCAGCAACTCTGc                                                                   >  2:299436/1‑151 (MQ=255)                                                                                  attCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTAt                                                                <  2:1004347/151‑1 (MQ=255)                                                                                        gTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAACGGTCTGGGGGATCCACCGCGACCAGAGATAGCCCCTCCGCCTTTGCCCGTGCTGCCCCCCGCGGTCTGCCCTTTTCCGCCCGCCACCCTGCTATCCAg                                                            >  2:49330/1‑149 (MQ=255)                                                                                          tAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCCGGGGAATCCACCGCGACCAGATACAGCCCCTCCCCCTTTGCCGTTGCTGCCGCCCGCGTTCTGTCCTTTTTCGCCAGCAACTCTGCTATCCAGctc                                                         >  2:268572/1‑150 (MQ=255)                                                                                          tAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATAATACCCGTCAGGGGAATCCACCGCGACCAGAAAACGCCCCTCCCCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCCGctc                                                         >  1:692650/1‑150 (MQ=255)                                                                                           aGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGTCAGGGGAATTCACCGCGACCAGATAAAACCCCTCCGCTTTTTCCGTTTCTGCCCCCAGCGTTCTGTCCTTTTCCGCCCGCAAATCTTCTCTCCAGctc                                                         >  2:838233/1‑149 (MQ=255)                                                                                           aGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAAATCATACCGGGCAGGGTAATCCACCCCGACCAGATAAAACCCCTCCGCTTTTTCCGGTTCTGCCGCCCGCGTTCTGTCCTTTGCCCCCCGCCACTCCGCTATACCGctccc                                                       >  1:1021365/1‑149 (MQ=255)                                                                                            gTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTcc                                                      <  1:157891/151‑1 (MQ=255)                                                                                                gCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATGATAACCGTCAGGGGAAACCACCGCGACCAGATACAGCCCTTCCCCTTTTTCCGTTGCTGCCGCCCGCGTTCTGTTCTTTGCCCCCAGCCACTCTGCTATCCCGCTcccccg                                                    >  2:306147/1‑147 (MQ=255)                                                                                                  cAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGGTCATACCGGGCAGGGTAATCCCCCCCCACCCGATACAGCCCTTCCGCTTTTTCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCCCCCGCAACTCTGCTATTCCGCTCTCCgg                                                    >  1:948791/1‑147 (MQ=255)                                                                                                              gCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCCCCg                                    <  1:91868/151‑1 (MQ=255)                                                                                                                 cATTGGCGGTTTTTGAAGATCAAAACGGTCAGGGGAAAACACCACGGCCAAAGACAGACCTTCCGCTTTTGTCGTTGCTGCCCCCCGCGTTCTTGCCCTTGCCGCCCGCAACTCTGCTATCCAGCTCTCCGGCCGGTTGTGGGCgacgac                                  >  1:137608/1‑150 (MQ=255)                                                                                                                 cATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACt                                 <  2:683105/151‑1 (MQ=255)                                                                                                                      gCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCAt                            >  2:1015621/1‑151 (MQ=255)                                                                                                                          ttttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAgg                        >  2:543164/1‑151 (MQ=255)                                                                                                                            ttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGc                       <  2:685281/150‑1 (MQ=255)                                                                                                                              ggAAGATAATACCGGTCTGGCTAATCCACCGCGACCGAATACAGCCGTTCCGCTTTTGCCGTTGCTTCCCGCAGCGTTCTGTCCTTTGCCGCCAGCAACGCAGCTAGCCAGCTTGCCGGCTGGTTGTGGGCTTCGACTTccacca                          <  1:862545/145‑1 (MQ=255)                                                                                                                               gAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTACCGTTGCTCCCGCCAGCGTTCTCTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGTCTGGTTGTAGGCACCGACTTCCATCat                         <  1:480434/145‑2 (MQ=255)                                                                                                                                aaGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCg                  >  1:730352/1‑151 (MQ=255)                                                                                                                                aaGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCg                  >  1:391711/1‑151 (MQ=255)                                                                                                                                     cATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAat             >  2:341525/1‑151 (MQ=255)                                                                                                                                      aTACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAat             >  1:425930/1‑150 (MQ=255)                                                                                                                                       tcccggtcgtcgTAATCCACCGCGGCCAGATTCAGCCCTGCCGCTTTTGCCGTTGCTGCGGCCCGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCGCCGGCTGGTTGTGGGGACCGACTGCCAGCGGGCTGCAGACaaga            <  1:876444/139‑3 (MQ=255)                                                                                                                                         ccGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTc         <  1:766249/151‑1 (MQ=255)                                                                                                                                               aGGGTAATCCACCGCGACCAGATCCAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGCTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAcc   >  1:43406/1‑151 (MQ=255)                                                                                                                                                gggTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCa  <  2:970521/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTGCTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCA  >  NC_000913/2434736‑2435029

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