Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,110,596	G→T	6.6%	K806N (AAG→AAT)	PP_3616 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,110,596	0	G	T	6.6%	181.8 / 2.9	61	K806N (AAG→AAT)	PP_3616	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (33/24);  new base T (2/2);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

GCTATTGGTTCCGCTGTCGTCCTTGAAGCAGGACTTGCAATTGCAGGCGTGACATCGGTTGTGCCCTTTGCCGGATGGGCCGCCGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑CCGACATCAACTTGGACTTTGACAAAAAACTGCCTGCTTTCACCTCACTGCAAACGGAAAT  >  NC_002947/4110453‑4110742                                                                                                                                                |                                                                                                                                                    gCTATTTGATCCGCGGTCGGCCTTGCAGCATGACTTGCAATTTCAGGCGTGACATCGGTTGTGCCATTTGCCGGATTGGCCGCCGATGCCCTTGTCCTCGTAGTTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGcac                                                                                                                                                 <  4:214598/148‑1 (MQ=255)         ttCCGCTGTCGTCCTTGAAGCGGGACTTGCAATTGCAGGCGTGACATCGGTTGTGCCCTTTGCCGGATGGGCCGCCGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGcc                                                                                                                                        >  3:61369/1‑149 (MQ=255)                      cTTGAAGCAGGACTTGCAATTGCAGGCGTGACATCGGTTGTGCCCTTTGCCGGATGGGCCGCCGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCaa                                                                                                                            >  3:338529/1‑148 (MQ=255)                      cTTGAAGCAGGACTTGCAATTGCAGGCGTGACATCGCTTGTGCCCTTTGCCGGATGGGCCGCCGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCaa                                                                                                                            >  4:288733/1‑148 (MQ=255)                         gAAGCAGGACTTGCAATTGCAGGCGTGACATCGGTTGTGCCCTTTGCCGGATGGGCCGCCGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATa                                                                                                                          <  4:338529/147‑1 (MQ=255)                         gAAGCAGGACTTGCAATTGCAGGCGTGACATCGGTTGTGCCCTTTGCCGGATGGGCCGCCGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATa                                                                                                                          <  1:54361/147‑1 (MQ=255)                            gCAGGACTTGCAATTGCAGGCGTGACATCGGTTGTGCCCTTTGCCGGATGGGCCGCCGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACt                                                                                                                     >  6:300152/1‑149 (MQ=255)                              aGGACTTGCAATTGCAGGCGTGACATCGGTTGTGCCCTTTGCCGGATGGGCCGCCGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATa                                                                                                                          <  8:109676/142‑1 (MQ=255)                              aGGACTTGCAATTGCAGGCGTGACATCGGTTGTGCCCTTTGCCGGATGGGCCGCCGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGAACATCAACGCCTTCACAGCCCAATa                                                                                                                          >  7:109676/1‑142 (MQ=255)                                     gCAATTGCAGGCGTGACATCGGTTGTGCCCTTTGCCGGATGGGCCGCCGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGgca                                                                                                              <  3:288733/147‑1 (MQ=255)                                      cAATTGCAGGCGTGACATCGGTTGTGCCCTTTGCCGGATGGGCCGCCGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGgcagca                                                                                                           <  6:265602/149‑1 (MQ=255)                                      cAATTGCAGGCGTGACATCGGTTGTGCCCTTTGCCGGATGGGCCGCCGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGgcagca                                                                                                           >  8:118705/1‑149 (MQ=255)                                         ttGCAGGCGTGACATCGGTTGTGCCCTTTGCCGGATGGGCCGCCGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACg                                                                                                         >  2:40554/1‑148 (MQ=255)                                           gCAGGCGTGACATCGGTTGTGCCCTTTGCCGGATGGGCCGCCGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTa                                                                                                       <  5:300152/148‑1 (MQ=255)                                            cAGGCGTGACATCGGTTGTGCCCTTTGCCGGATGGGCCGCCGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAg                                                                                                      >  2:265663/1‑148 (MQ=255)                                                    aCATCGGTTGTGCCCTTTGCCGGATGGGCCGCCGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCgg                                                                                             <  7:118705/149‑1 (MQ=255)                                                      aTCGGTTGTGCCCTTTGCCGGATGGGCCGCCGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGc                                                                                            >  5:138240/1‑148 (MQ=255)                                                              tgCCCTTTGCCGGATGGGCCGCCGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATc                                                                                   <  6:138240/149‑1 (MQ=255)                                                               gCCCTTTGCCGGATGGGCCGCCGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACt                                                                                                                     >  3:342217/1‑114 (MQ=255)                                                               gCCCTTTGCCGGATGGGCCGCCGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACt                                                                                                                     <  4:342217/114‑1 (MQ=255)                                                                  cTTTGCCGGATGGGCCGCCGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTa                                                                                                       >  3:219878/1‑125 (MQ=255)                                                                  cTTTGCCGGATGGGCCGCCGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTa                                                                                                       <  4:219878/125‑1 (MQ=255)                                                                   tttGCCGGATGGGCCGCCGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAatga                                                                              >  6:7776/1‑149 (MQ=255)                                                                       ccGGATGGGCCGCCGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGga                                                                          <  5:7776/149‑1 (MQ=255)                                                                        cGGATGGGCCGCCGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGTCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGga                                                                          <  1:265663/148‑1 (MQ=255)                                                                               gccgccGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCata                                                                                                  >  2:59056/1‑117 (MQ=255)                                                                               gccgccGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCata                                                                                                  <  1:59056/117‑1 (MQ=255)                                                                                  gccgcTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑c                                                              >  5:169787/1‑149 (MQ=255)                                                                                      cTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑CCGAc                                                          <  1:40554/149‑1 (MQ=255)                                                                                       tGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑CCGACa                                                         >  3:103761/1‑149 (MQ=255)                                                                                         cccTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑CCGACa                                                         <  4:61369/147‑1 (MQ=255)                                                                                           cTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑CCGACTTCaa                                                     >  6:281471/1‑149 (MQ=255)                                                                                                ccTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑CCGACATCAACTTg                                                 >  6:61722/1‑148 (MQ=255)                                                                                                  tCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCgg                                                                                             <  3:72124/103‑1 (MQ=255)                                                                                                  tCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCgg                                                                                             >  4:72124/1‑103 (MQ=255)                                                                                                  tCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑CCGACATCAACTTGGAc                                              >  3:167829/1‑149 (MQ=255)                                                                                                    gTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑CCGACATCAACTTGGACtt                                            >  7:74959/1‑149 (MQ=255)                                                                                                        gtgtAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑CCGACATCAACTTGGACTTTGa                                         <  8:24563/148‑1 (MQ=255)                                                                                                          gtAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGagag                                                                       >  4:54661/1‑117 (MQ=255)                                                                                                          gtAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGagag                                                                       <  3:54661/117‑1 (MQ=255)                                                                                                          gtAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑CCGACATCAACTTg                                                 >  7:63976/1‑138 (MQ=255)                                                                                                          gtAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑CCGACATCAACTTg                                                 <  8:63976/138‑1 (MQ=255)                                                                                                          gtAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑CCGACATCAACTTGGACTTTGACa                                       >  4:142372/1‑148 (MQ=255)                                                                                                           tAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑CCGACATCAACTTGGACTTTGACaa                                      >  2:264997/1‑148 (MQ=255)                                                                                                             gCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑CCGACATCAACTTGGACTTTGACaaaa                                    >  7:313910/1‑148 (MQ=255)                                                                                                               catcatCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑CCGACATCAACTTGGACTTTGACAAAAAAc                                 >  6:92481/1‑149 (MQ=255)                                                                                                                    tcGCAGGAGGCATCTACCTACACGCCAAGGCACATGAACGCCTTAACAGCCCAATAGAACTGTTGGCAGCACGTAGAATATTCGGCAATCGCATCAATGATGGAGAGGCGCACT‑CAGACATCAACTTGTAAGTAGAAAAAAAAatccct                            >  4:364740/1‑149 (MQ=255)                                                                                                                     cGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑CCGACATCAACTTGGACTTTGACAAAAAActgcctg                           >  8:162069/1‑149 (MQ=255)                                                                                                                      gCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑CCGACATCAACTTGGACTTTGACAAAAAActgcctgc                          <  8:74959/149‑1 (MQ=255)                                                                                                                       cAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACGCCCGACATCAACTTGGACTTTGACAAAAAActgcctgc                          <  1:264997/149‑1 (MQ=255)                                                                                                                            ggCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑CCGACATCAACTTGGACTTTGACAAAAAACTGCCTGCTTTCAc                    <  5:61722/149‑1 (MQ=255)                                                                                                                               ctctACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCACTCGCATCAATGATGGAGAGGCGCACG‑CCGACATCAACTTGGACTTTGACAAAAAACTGCCTGCTTTCACCt                  >  8:327824/1‑148 (MQ=255)                                                                                                                                  tacctacACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑CCGACATCAACTTGGACTTTGACAAAAAACTGCCTGCTTTCACCTCa                <  5:92481/147‑1 (MQ=255)                                                                                                                                      tacACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑CCGACATCAACTTGGACTTTGACAAAAAACTGCCTGCTTTCACCTCACTGCa           >  7:16494/1‑148 (MQ=255)                                                                                                                                       acacGCCAATGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑CCGACATCAACTTGGACTTTGACaaaaa                                   <  4:153966/123‑1 (MQ=255)                                                                                                                                       acacGCCAATGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑CCGACATCAACTTGGACTTTGACaaaaa                                   <  4:148275/123‑1 (MQ=255)                                                                                                                                       acacGCCAATGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑CCGACATCAACTTGGACTTTGACaaaaa                                   >  3:153966/1‑123 (MQ=255)                                                                                                                                       acacGCCAATGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑CCGACATCAACTTGGACTTTGACaaaaa                                   >  3:148275/1‑123 (MQ=255)                                                                                                                                           gCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCAACAATGATGGAGAGGCGCACG‑CCGACATCAACTTGGACTTTGACAAAAAACTGCCTGCTTTCACCTCACTGCAAACgg      >  2:10368/1‑148 (MQ=255)                                                                                                                                             cAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑CCGACATCAACTTGGACTTTGACAAAAAACTGCCTGCTTTCACCTCACTGCAAACGGaaa   >  7:231989/1‑149 (MQ=255)                                                                                                                                               aGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑CCGACATCAACTTGGACTTTGACAAAAAACTGCCTGCTTTCACCTCACTGCAAACGGAAAt  >  1:13800/1‑148 (MQ=255)                                                                                                                                                |                                                                                                                                                    GCTATTGGTTCCGCTGTCGTCCTTGAAGCAGGACTTGCAATTGCAGGCGTGACATCGGTTGTGCCCTTTGCCGGATGGGCCGCCGCTGCCCTTGTCCTCGTAGGTGTAGCCATCATCGCAGGAGGCCTCTACCTACACGCCAAGGCACATGAACGCCTTCACAGCCCAATAGAACTGTGGGCAGCACGTAGCATATTCGGCAATCGCATCAATGATGGAGAGGCGCACG‑CCGACATCAACTTGGACTTTGACAAAAAACTGCCTGCTTTCACCTCACTGCAAACGGAAAT  >  NC_002947/4110453‑4110742

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