Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,904,604	T→A	7.1%	C287S (TGC→AGC)	nirB →	nitrite reductase large subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,904,604	0	T	A	7.1%	163.2 / 2.6	56	C287S (TGC→AGC)	nirB	nitrite reductase large subunit
Reads supporting (aligned to +/- strand):  ref base T (25/27);  new base A (2/2);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

AAACCGACCTGATCGTGTTCTCCGCCGGTATCCGCCCACAGGATGCACTGGGCCGGGGTTGCGGCCTGGACATCGCCGCGCGCGGTGGCGTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAGGCGGCGGTGGCGTTCACCGGCGCCGACATGTCGACCAAGCTCAAGCTGCTGGGTGTGGATGT  >  NC_002947/1904458‑1904752                                                                                                                                                   |                                                                                                                                                       aatccGACCTGATCGTGTTCTCCGCCGGTATCCGCCCACAGGATTCACTGGTCCGGGGTTGCTGCCTGGACATCGCCGCGCGCGGTGGCGTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTgc                                                                                                                                                      <  2:284446/146‑1 (MQ=255)         cTGATCGTGTTCTCCGCCGGTATCCGCCCACAGGATGCACTGGGCCGGGGTTGCGGCCTGGACATCGCCGCGCGCGGTGGCGTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑g                                                                                                                                             <  1:213711/149‑1 (MQ=255)         cTGATCGTGTTCTCCGCCGGTATCCGCCCACAGGATGCACTGGGCCGGGGTTGCGGCCTGGACATCGCCGCGCGCGGTGGCGTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑g                                                                                                                                             <  7:131239/149‑1 (MQ=255)                tgtTCTCCGCCGGTATCCGCCCACAGGATGCACTGGGCCGGGGTTGCGGCCTGGACATCGCCGCGCGCGGTGGCGTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACgg                                                                                                                                       >  5:11788/1‑148 (MQ=255)                tgtTCTCCGCCGGTATCCGCCCACAGGATGCACTGGGCCGGCGTTGCGGCCTGGACATCGCCGCGCGCGGTGGCGTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACgg                                                                                                                                       >  7:144850/1‑148 (MQ=255)                 gtTCTCCGCCGGTATCCGCCCACAGGATGCACTGGGCCGGGGTTGCGGCCTGGACATCGCCGCGCGCGGTGGCGTGGTCATCGATAACCACTGCCGCCGCAGCGACCCCCGCCTCTTCGCCATCGGGGAGTGCGCCTCGG‑GGAAAgg                                                                                                                                       >  8:275696/1‑147 (MQ=255)                  ttCTCCGCCGGTATCCGCCCACAGGATGCACTGGGCCGGGGTTGCGGCCTGGACATCGCCGCGCGCGGTGGCGTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAg                                                                                                                                    >  7:36402/1‑149 (MQ=255)                     tcCGCCGGTATCCGCCCACAGGATGCACTGGGCCGGGGTTGCGGCCTGGACATCGCCGCGCGCGGTGGCGTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAAc                                                                                                                                         >  5:224600/1‑141 (MQ=255)                     tcCGCCGGTATCCGCCCACAGGATGCACTGGGCCGGGGTTGCGGCCTGGACATCGCCGCGCGCGGTGGCGTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAAc                                                                                                                                         <  6:224600/141‑1 (MQ=255)                         ccgGTATCCGCCCACAGGATGCACTGGGCCGGGGTTGCGGCCTGGACATCGCCGCGCGCGGTGGCGTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTc                                                                                                                             <  3:231732/149‑1 (MQ=255)                         ccgGTATCCGCCCACAGGATGCACTGGGCCGGGGTTGCGGCCTGGACATCGCCGCGCGCGGTGGCGTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTc                                                                                                                             <  3:44879/149‑1 (MQ=255)                             tATCCGCCCACAGGATGCACTGGGCCGGGGTTGCGGCCTGGACATCGCCGCGCGCGGTGGCGTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGcc                                                                                                                         <  1:354287/149‑1 (MQ=255)                               tCCGCCCACAGGATGCACTGGGCCGGGGTTGCGGCCTGGACATCGCCGCGCGCGGTGGCGTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTg                                                                                                                       <  2:103341/149‑1 (MQ=255)                                   cccACAGGATGCACTGGGCCGGGGTTGCGGCCTGGACATCGCCGCGCGCGGTGGCGTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCg                                                                                                                   <  2:59274/149‑1 (MQ=255)                                    ccACAGGATGCACTGGGCCGGGGTTGCGGCCTGGACATCGCCGCGCGCGGTGGCGTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGTCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGc                                                                                                                  >  7:255096/1‑149 (MQ=255)                                    ccACAGGATGCACTGGGCCGGGGTTGCGGCCTGGACATCGCCGCGCGCGGTGGCGTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGc                                                                                                                  <  8:35707/149‑1 (MQ=255)                                          gATGCACTGGGCCGGGGTTGCGGCCTGGACATCGCCGCGCGCGGTGGCGTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCggg                                                                                                            >  3:249100/1‑149 (MQ=255)                                           aTGCACTGGGCCGGGGTTGCGGCCTGGACATCGCCGCGCGCGGTGGCGTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGAGCGCCTCGTCGGAACGGCAGCGTGTTCGGCCTGGTCGCCCCg                                                                                                              <  1:40049/147‑1 (MQ=255)                                           aTGCACTGGGCCGGGGTTGCGGCCTGGACATCGCCGCGCGCGGTGGCGTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGAGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCg                                                                                                              >  2:40049/1‑146 (MQ=255)                                           aTGCACTGGGCCGGGGTTGCGGCCTGGACATCGCCGCGCGCGGTGGCGTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGAGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCg                                                                                                              >  2:45805/1‑146 (MQ=255)                                           aTGCACTGGGCCGGGGTTGCGGCCTGGACATCGCCGCGCGCGGGGGCGTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGAGCGCCTCGTCGGAACGGCAGCGTGTTCGGCCTGGTCGCCCCg                                                                                                              <  1:45805/147‑1 (MQ=255)                                                        gggTTGCGGCCTGGACATCGCCGCGCGCGGTGGCGTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGACCTGGTCGCCCCGGGCTACAGCATGgcgc                                                                                              >  2:113106/1‑149 (MQ=255)                                                          gTTGCGGCCTGGACATCGCCGCTCGCGGTGGCGTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGTCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGgcgc‑gc                                                                                           <  8:255096/149‑1 (MQ=255)                                                                         tCGCCGCGCGCGGTGGCGTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGgctg                                                                            <  8:144850/149‑1 (MQ=255)                                                                            ccgcGCGCGGTGGCGTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGgctgct                                                                          >  2:75430/1‑148 (MQ=255)                                                                                  gcgGTGGCGTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGc                                                                   <  7:275696/149‑1 (MQ=255)                                                                                          gTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCg                                                                                                                            <  5:274746/85‑1 (MQ=255)                                                                                          gTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCg                                                                                                                            >  6:274746/1‑85 (MQ=255)                                                                                            ggTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAGGcggcgg                                                         >  1:220146/1‑149 (MQ=255)                                                                                               cATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAGGCGGCGGTgg                                                      >  1:241199/1‑149 (MQ=255)                                                                                                aTCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAGGCGGCGGTGGc                                                     >  2:219630/1‑149 (MQ=255)                                                                                                  cGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAGGCGGCGGTGGCGt                                                   >  8:354945/1‑149 (MQ=255)                                                                                                   gATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATg                                                                                                  <  8:159463/102‑1 (MQ=255)                                                                                                   gATAACCACTGCCGCAGCAGCGACCCACGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATg                                                                                                  >  7:159463/1‑102 (MQ=255)                                                                                                        ccACTGCCGCAGCAGCGACCCCCGCATCCTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAGGCGGCGGTGGCGTTCACCg                                             >  6:246416/1‑149 (MQ=255)                                                                                                         cACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAggcgg                                                            >  8:327124/1‑133 (MQ=255)                                                                                                         cACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAggcgg                                                            <  7:327124/133‑1 (MQ=255)                                                                                                            tGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAGGCGGCGGTGGCGTTCACCGgcgc                                         >  8:46290/1‑149 (MQ=255)                                                                                                              ccGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAGGCGGCGGTGGCGTTCACCGGCGcc                                        >  1:185909/1‑148 (MQ=255)                                                                                                                gcagcaGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAGGCGGCGGTGGCGTTCACCGGCGCCGAc                                     <  8:170933/149‑1 (MQ=255)                                                                                                                       cGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAGGCGGCGGTGGCGTTCACCGGCGCCGACATGTCGa                              <  1:75430/149‑1 (MQ=255)                                                                                                                       cGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAGGCGGCGGTGGCGTTCACCGGCGCCGACATGTCGa                              <  6:11788/149‑1 (MQ=255)                                                                                                                         aCCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAGGCGGCGGTGGCGTTCACCGGCGCCGACATGTCGAcc                            <  2:241199/149‑1 (MQ=255)                                                                                                                         aCCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAGGCGGCGGTGGCGTTCACCGGCGCCGACATGTCGAcc                            <  7:46290/149‑1 (MQ=255)                                                                                                                         aCCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGACCTGGTCGCCCCGGGCTACAGCATGGCGCAGCAACCTGGCCGGGCTGCTGATGGGCGAGGCGGCGGTGGCGTTCACCGGCGCCGACATGTCGAc                             <  1:113106/149‑1 (MQ=255)                                                                                                                             ccGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAGGCGGCGGTGGCGTTCACCGGCGCCGACATGTCGACcaagc                        >  3:142620/1‑149 (MQ=255)                                                                                                                             ccGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAGGCGGCGGTGGCGTTCACCGGCGCCGACATGTCGACcaag                         >  8:3140/1‑148 (MQ=255)                                                                                                                               gCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAGGCGGCGGTGGc                                                     >  8:101851/1‑118 (MQ=255)                                                                                                                               gCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAGGCGGCGGTGGc                                                     <  7:101851/118‑1 (MQ=255)                                                                                                                                   cTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAGGCGGCGGTGGCGTTCAcc                                              >  8:350491/1‑121 (MQ=255)                                                                                                                                   cTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAGGCGGCGGTGGCGTTCAcc                                              <  7:350491/121‑1 (MQ=255)                                                                                                                                    ttCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAGGCGGCGGTGGCGTTCACCGGCGCCGACATGTCGACcaagctcaagct                 >  3:68671/1‑149 (MQ=255)                                                                                                                                      cGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAGGCGGCGGTGGCGTTCACCGGCGCCGACATGTCGACcaagctcaag                   <  7:354945/145‑1 (MQ=255)                                                                                                                                         cATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAGGCGGCGGTGGCGTTCACCGGCGCCGACATGTCGACcaagctcaag                   >  6:212063/1‑142 (MQ=255)                                                                                                                                         cATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAGGCGGCGGTGGCGTTCACCGGCGCCGACATGTCGACcaagctcaag                   <  5:212063/142‑1 (MQ=255)                                                                                                                                                 aGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAGGCGGCGGTGGCGTTCACCGGCGCCGACATGTCGACCAAGCTCAAGCTGCTGGGTGTGGAt    <  5:246416/149‑1 (MQ=255)                                                                                                                                                   tGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAGGCGGCGGTGGCGTTCACCGGCGCCGACATGTCGACCAAGCTCAAGCTGCTGGGTGTGGATGt  >  7:353114/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                       AAACCGACCTGATCGTGTTCTCCGCCGGTATCCGCCCACAGGATGCACTGGGCCGGGGTTGCGGCCTGGACATCGCCGCGCGCGGTGGCGTGGTCATCGATAACCACTGCCGCAGCAGCGACCCCCGCATCTTCGCCATCGGCGAGTGCGCCTCGT‑GGAACGGCAGCGTGTTCGGCCTGGTCGCCCCGGGCTACAGCATGGCGC‑GCAACCTGGCCGGGCTGCTGATGGGCGAGGCGGCGGTGGCGTTCACCGGCGCCGACATGTCGACCAAGCTCAAGCTGCTGGGTGTGGATGT  >  NC_002947/1904458‑1904752

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: ‑