Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,083,620	G→A	R60C (CGT→TGT)	speB ←	agmatinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,083,620	0	G	A	100.0%	245.6 / NA	75	R60C (CGT→TGT)	speB	agmatinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (38/37);  total (38/37)

TCGCTCATCTCACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACGGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTAAAAAACCAAAGG  >  NC_000913/3083490‑3083751                                                                                                                                   |                                                                                                                                    tCGCTCATCTCACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTg                                                                                                                             >  1:1095757/1‑139 (MQ=255)                gCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGac                                                                                                              <  1:1757003/139‑1 (MQ=255)                   tCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGaccac                                                                                                           >  1:1148491/1‑139 (MQ=255)                   tCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGaccac                                                                                                           >  2:916660/1‑139 (MQ=255)                   tCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGaccac                                                                                                           >  2:1528831/1‑139 (MQ=255)                          aGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTg                                                                                                                   >  2:839433/1‑124 (MQ=255)                          aGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTg                                                                                                                   <  1:839433/124‑1 (MQ=255)                               tataCCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGc                                                                                                       >  1:278692/1‑131 (MQ=255)                               tataCCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGc                                                                                                       <  2:278692/131‑1 (MQ=255)                                      gATCGCCGCAGTCCATGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGAt                                                                                                                                  <  1:618480/97‑1 (MQ=255)                                      gATCGCCGCAGTCCATGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGAt                                                                                                                                  >  2:618480/1‑97 (MQ=255)                                      gATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAg                                                                                            <  1:1004350/135‑1 (MQ=255)                                      gATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAg                                                                                            >  2:1004350/1‑135 (MQ=255)                                      gATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGc                                                                                        >  2:1475896/1‑139 (MQ=255)                                       aTCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGc                                                                                                       >  2:171412/1‑123 (MQ=255)                                       aTCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGc                                                                                                       <  1:171412/123‑1 (MQ=255)                                           ccgcAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGAt                                                                                                                                  <  2:1951513/92‑1 (MQ=255)                                           ccgcAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGAt                                                                                                                                  >  1:1951513/1‑92 (MQ=255)                                           ccgcAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGacc                                                                                                             <  2:1938492/113‑1 (MQ=255)                                           ccgcAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGacc                                                                                                             >  1:1938492/1‑113 (MQ=255)                                           ccgcAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCatat                                                                                   >  1:1376288/1‑139 (MQ=255)                                           ccgcAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCatat                                                                                   >  2:724896/1‑139 (MQ=255)                                            cgcAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGAc                                                                                                  >  2:1911723/1‑123 (MQ=255)                                            cgcAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGAc                                                                                                  <  1:1911723/123‑1 (MQ=255)                                            cgcAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATc                                                                                  >  1:427344/1‑139 (MQ=255)                                               aGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGaa                                                                               <  1:724896/139‑1 (MQ=255)                                                 tCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGAt                                                                                                                                  >  2:949784/1‑86 (MQ=255)                                                 tCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGAt                                                                                                                                  <  1:949784/86‑1 (MQ=255)                                                       acgTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGc                                                                       <  2:427344/139‑1 (MQ=255)                                                           tCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGcc                                                                      >  1:1227868/1‑136 (MQ=255)                                                           tCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGcc                                                                      <  2:1227868/136‑1 (MQ=255)                                                           tCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGt                                                                   >  1:1552756/1‑139 (MQ=255)                                                            cAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTa                                                                  <  2:593371/139‑1 (MQ=255)                                                                  cTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGt                                                                   <  1:214356/132‑1 (MQ=255)                                                                  cTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGt                                                                   >  2:214356/1‑132 (MQ=255)                                                                   tCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGacca                                                                                                            >  1:2018688/1‑90 (MQ=255)                                                                   tCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGacca                                                                                                            <  2:2018688/90‑1 (MQ=255)                                                                   tCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCAccc                                                           <  2:1880872/139‑1 (MQ=255)                                                                   tCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCAccc                                                           <  2:1376288/139‑1 (MQ=255)                                                                    cACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCa                                                          >  1:1734695/1‑139 (MQ=255)                                                                           gTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCa                                                   >  2:1169385/1‑139 (MQ=255)                                                                               aaaTTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGc                                               <  1:1169385/139‑1 (MQ=255)                                                                                     cACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGAGAGATCGCTGCCGGACCGTGGCGACCACGCGCACGACCAGa                                                                                              >  1:1187976/1‑86 (MQ=255)                                                                                     cACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGa                                                                                              <  2:1187976/86‑1 (MQ=255)                                                                                      aCGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATa                                        >  1:1715484/1‑139 (MQ=255)                                                                                       cGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATAc                                       <  2:1552756/139‑1 (MQ=255)                                                                                            aGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGa                                                                                                                        <  2:759714/53‑1 (MQ=255)                                                                                            aGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGa                                                                                                                        <  2:797089/53‑1 (MQ=255)                                                                                            aGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGa                                                                                                                        >  1:797089/1‑53 (MQ=255)                                                                                            aGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGa                                                                                                                        >  1:759714/1‑53 (MQ=255)                                                                                            aGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCg                                                                                                                <  1:1922285/61‑1 (MQ=255)                                                                                            aGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCg                                                                                                                >  2:1922285/1‑61 (MQ=255)                                                                                              cGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCg                                                                                 <  2:5077/90‑1 (MQ=255)                                                                                              cGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCg                                                                                 >  1:5077/1‑90 (MQ=255)                                                                                                gTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGt                                                        <  2:1735421/113‑1 (MQ=255)                                                                                                gTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGt                                                        >  1:1735421/1‑113 (MQ=255)                                                                                                    tgtTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGAc                                                                                                  <  1:549823/67‑1 (MQ=255)                                                                                                    tgtTCCCAGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGAc                                                                                                  >  2:549823/1‑67 (MQ=255)                                                                                                           aGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGaa                                                                               <  1:1695654/79‑1 (MQ=255)                                                                                                           aGGCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGaa                                                                               >  2:1695654/1‑79 (MQ=255)                                                                                                            ggCCAGATTCGTCGAAACCTGACATATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAgg                  <  2:1148491/139‑1 (MQ=255)                                                                                                             gCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGc                                                                                                       >  1:1368810/1‑53 (MQ=255)                                                                                                             gCCAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGc                                                                                                       <  2:1368810/53‑1 (MQ=255)                                                                                                              ccAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGa                                                                                              <  2:1862275/61‑1 (MQ=255)                                                                                                              ccAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGa                                                                                              >  1:1862275/1‑61 (MQ=255)                                                                                                              ccAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCAcc                                                            >  1:1636753/1‑95 (MQ=255)                                                                                                              ccAGATTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCAcc                                                            <  2:1636753/95‑1 (MQ=255)                                                                                                                  aTTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACTCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGc                     >  2:516370/1‑130 (MQ=255)                                                                                                                  aTTCGTCGAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACTCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGc                     <  1:516370/130‑1 (MQ=255)                                                                                                                       tcgAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGAc                                                                                                  >  1:193850/1‑48 (MQ=255)                                                                                                                       tcgAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGAc                                                                                                  <  2:193850/48‑1 (MQ=255)                                                                                                                       tcgAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTAAAAAAcc       >  2:918152/1‑139 (MQ=255)                                                                                                                         gAAACCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTAAAAAACCaa     <  1:918152/139‑1 (MQ=255)                                                                                                                            aCCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTAAAAAACCAAAgg  <  1:1475896/139‑1 (MQ=255)                                                                                                                            aCCTGACAGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTAAAAAACCAAAgg  <  1:1528831/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                                    TCGCTCATCTCACGGGCATCGCCAAAGGCATATACCAGATCGCCGCAGTCCACGACGTTCAGACGCTCACGCATGTCGAAATTCCACGGGAAGCGGTTGTGTTCCCAGGCCAGATTCGTCGAAACCTGACGGATCGCTGCCGGACCGTGGCGACCACCCGCACGACCAGAAGTGGCCATATCGAACGGCACGCCAGTAATCACCCAGTCTGCATCGCTGTCATACGGCTGGAAGTTCATCGGCAGGCGTAAAAAACCAAAGG  >  NC_000913/3083490‑3083751

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