Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,731,970	A→T	8.7%	L196M (TTG→ATG) ‡	hutI ←	imidazolonepropionase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,731,970	0	A	T	8.7%	148.3 / 3.9	54	L196M (TTG→ATG) ‡	hutI	imidazolonepropionase
Reads supporting (aligned to +/- strand):  ref base A (23/25);  new base T (2/3);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCGCAGCATCTTG  >  NC_002947/5731831‑5732118                                                                                                                                            |                                                                                                                                                     gATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCaa                                                                                                                                             <  5:224997/149‑1 (MQ=255) gATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCaa                                                                                                                                             <  5:241008/149‑1 (MQ=255)  aTGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCTTGGCCGGcgaa                                                                                                                                            <  4:121636/149‑4 (MQ=255)  aTGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAAc                                                                                                                                            <  8:101525/149‑1 (MQ=255)    gACAGTTGTTCGTCATGCAGCTTGACCGGCAGGCCCAGTTCGCTCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCATCGCCGGCCTTGGCCGGcgaag                                                                                                                                           <  7:318152/148‑5 (MQ=255)    gACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCGGTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGGATCCACCGCCACCACCAGGCCTTCGCCGGCCACGGCCGGCAACAt                                                                                                                                          >  4:270243/1‑149 (MQ=255)       aGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACAtctc                                                                                                                                       <  5:145014/149‑1 (MQ=255)        gttgttCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCg                                                                                                                                      >  5:308195/1‑149 (MQ=255)        gttgttCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCg                                                                                                                                      <  3:270243/149‑1 (MQ=255)         ttgttCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGt                                                                                                                                     <  4:50327/149‑1 (MQ=255)            ttCGGGATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAAGAGGCGCTCGACCTGGGCCGGCGAGAAGGCGAGGTGTTCGCAGAAGTCATGCACCCGATCCACCACGCCTTCGCCGGCCATGGCCGTCAACATCTCGTca                                                                                                                                   <  2:139822/148‑1 (MQ=255)               ggCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTca                                                                                                                                   >  7:55002/1‑145 (MQ=255)                gCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGAt                                                                                                                              <  4:228079/149‑1 (MQ=255)                  aTGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGAtgt                                                                                                                            <  6:308195/149‑1 (MQ=255)                  aTGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGAtgt                                                                                                                            <  8:314606/149‑1 (MQ=255)                           gACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGAt                                                                                                                              <  8:73161/138‑1 (MQ=255)                           gACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGAt                                                                                                                              >  7:73161/1‑138 (MQ=255)                           gACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATg                                                                                                                    <  2:108459/148‑1 (MQ=255)                           gACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATg                                                                                                                    >  5:311983/1‑148 (MQ=255)                           gACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGCCAACATCTCGTCACAGATGTGCGCGATg                                                                                                                    <  1:209420/148‑1 (MQ=255)                               ggCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAgtc                                                                                                               >  3:26457/1‑149 (MQ=255)                                     cccAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGa                                                                                                                      >  1:145261/1‑136 (MQ=255)                                       cAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGccc                                                                                                       >  8:2236/1‑149 (MQ=255)                                                cgcCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCACGTCACACATGTGCGCGATGTAGTCGTCGGCCAGACCcgcg                                                                                               >  4:153288/1‑148 (MQ=255)                                                          aaCAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCgggg                                                                                     <  8:55002/148‑1 (MQ=255)                                                                cGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTTGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGc                                                                              >  8:239530/1‑149 (MQ=255)                                                                cGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGc                                                                              >  8:239529/1‑149 (MQ=255)                                                                 gCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCg                                                                             >  3:172870/1‑149 (MQ=255)                                                                           ggCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCa                                                                   <  7:239530/149‑1 (MQ=255)                                                                           ggCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCa                                                                   <  7:239529/149‑1 (MQ=255)                                                                           ggCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCa                                                                   <  7:2236/149‑1 (MQ=255)                                                                            gCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAg                                                                  <  4:172870/149‑1 (MQ=255)                                                                                  gagaACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGt                                                            >  1:240415/1‑149 (MQ=255)                                                                                      aCGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGgc                                                         <  4:26457/148‑1 (MQ=255)                                                                                        gcgAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTgc                                                                          <  7:201878/129‑1 (MQ=255)                                                                                        gcgAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTgc                                                                          >  8:201878/1‑129 (MQ=255)                                                                                          gAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGgcacgca                                                    >  5:365640/1‑149 (MQ=255)                                                                                              tgtTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGt                                                >  3:297256/1‑149 (MQ=255)                                                                                                     aGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTg                                                           <  2:309620/131‑1 (MQ=255)                                                                                                     aGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTg                                                           >  1:309620/1‑131 (MQ=255)                                                                                                         ggCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCt                                     >  5:314586/1‑149 (MQ=255)                                                                                                              ccaccgcatccaccAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCtcgtcg                                >  6:39399/1‑149 (MQ=255)                                                                                                                accgcatccaccAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGc                              >  6:22033/1‑149 (MQ=255)                                                                                                                    catccaccAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCtcgt                                  >  4:201498/1‑141 (MQ=255)                                                                                                                    catccaccAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCtcgt                                  <  3:201498/141‑1 (MQ=255)                                                                                                                         accaGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGgcggcg                     >  4:76048/1‑149 (MQ=255)                                                                                                                          ccaGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGcggcgg                    <  6:71538/149‑1 (MQ=255)                                                                                                                              gatgttgCCGGCCTTGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGGCGGCCAGGCGGCGGGCCa                >  7:293451/7‑149 (MQ=255)                                                                                                                                cTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACt              >  1:285731/1‑149 (MQ=255)                                                                                                                                  ttgCCGGCCTTGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCg            >  6:193576/3‑149 (MQ=255)                                                                                                                                    gCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCgca          <  3:76048/149‑1 (MQ=255)                                                                                                                                    gCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCgc           <  5:22033/148‑1 (MQ=255)                                                                                                                                        gCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCGCAGCAt      >  7:178659/1‑149 (MQ=255)                                                                                                                                         ccAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCGCAGCATc     <  5:250596/149‑1 (MQ=255)                                                                                                                                           aaGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTCCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCGCAGCATCtt   >  3:122813/1‑149 (MQ=255)                                                                                                                                            aGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCGCAGCATCTTg  <  8:178659/149‑1 (MQ=255)                                                                                                                                            aGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCGCAGCATCTTg  <  2:285731/149‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                     GATGACAGTTGTTCGGCATGCAGCTTGACCGGCAGGCCCAGTTCGCGCGCCTTGATGAACAGCCGCTCGACCTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGCACGTGGCACGCACCGTCAGCGGCAGCTCGTCGGCCAGGCGGCGGGCCACTCGCAGCATCTTG  >  NC_002947/5731831‑5732118

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