Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	355,222	G→T	16.1%	intergenic (‑51/+34)	cbcW ← / ← cbcX	choline/betaine/carnitine ABC transporter permease/choline/betaine/carnitine ABC transporter substrate‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	355,222	0	G	T	16.1%	127.8 / 18.8	56	intergenic (‑51/+34)	cbcW/cbcX	choline/betaine/carnitine ABC transporter permease/choline/betaine/carnitine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (26/21);  new base T (4/5);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.62e-01

TCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATCCATTTCTTGGC  >  NC_002947/355086‑355369                                                                                                                                          |                                                                                                                                                     tCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGACCGTGCCTGCCCCggg                                                                                                                                          <  5:329703/149‑3 (MQ=255)      gTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGCGGATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAg                                                                                                                                      >  1:249659/1‑149 (MQ=255)                   ttCTGGGTCAGCCATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGt                                                                                                                        >  7:277827/1‑149 (MQ=255)                        ggTCACCCATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCg                                                                                                                     <  3:305245/147‑1 (MQ=255)                         gTCAGCCATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTa                                                                                                                  >  8:154381/1‑149 (MQ=255)                         gTCAGCCATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTa                                                                                                                  >  1:321807/1‑149 (MQ=255)                          tCAGCCATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGACCGTGCCTGCCCCg                                                                                                                                            <  1:147633/122‑1 (MQ=255)                          tCAGCCATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGACCGTGCCTGCCCCg                                                                                                                                            <  1:47830/122‑1 (MQ=255)                          tCAGCCATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGACCGTGCCTGCCCCg                                                                                                                                            >  2:147633/1‑122 (MQ=255)                          tCAGCCATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGACCGTGCCTGCCCCg                                                                                                                                            >  2:47830/1‑122 (MQ=255)                            aGCCATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTg                                                                                                               >  8:218763/1‑149 (MQ=255)                             gCCATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTTGTTCGTTATTGc                                                                                                              <  7:150264/149‑1 (MQ=255)                              ccATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTAtt                                                                                                                <  5:132499/146‑1 (MQ=255)                                aTTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGt                                                                                                            >  7:208844/1‑148 (MQ=255)                                aTTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGt                                                                                                            >  2:107121/1‑148 (MQ=255)                                    gacgaACGCCGC‑GATGTACAGGACCAGGGGAAGTTTCTGAACGATAAGAATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCAGGGCAGGCACGGTCCTGCACCTAGGTAGCGCTACTCCGTAGATAGTTAATGCGTATGCt                                                                                                       >  8:155488/1‑149 (MQ=255)                                         acgacgC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAACGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGcc                                                                                                  >  6:296227/1‑149 (MQ=255)                                           gacgC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCTGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCtt                                                                                                <  8:177344/149‑1 (MQ=255)                                                        tGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTccag                                                                                    >  3:104080/1‑149 (MQ=255)                                                             cAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTccaggc‑cag                                                                              >  5:45442/1‑149 (MQ=255)                                                               ggggTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCcaggc‑caggc                                                                            <  3:226195/149‑1 (MQ=255)                                                                  gTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTg                                                                         >  6:9765/1‑149 (MQ=255)                                                                   tATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGc                                                                        <  7:154381/149‑1 (MQ=255)                                                                    aTTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTccaggc‑cag                                                                              >  2:252568/1‑142 (MQ=255)                                                                    aTTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCCTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTccaggc‑cag                                                                              <  1:252568/142‑1 (MQ=255)                                                                     ttttCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGcc                                                                       >  8:161844/1‑148 (MQ=255)                                                                      tttCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTcc                                                                                      >  5:178225/1‑133 (MQ=255)                                                                      tttCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTcc                                                                                      <  6:178225/133‑1 (MQ=255)                                                                          tGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAAc                                                                 >  2:287162/1‑149 (MQ=255)                                                                            aTCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCaa                                                                  <  2:249659/146‑1 (MQ=255)                                                                               gATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACggtg                                                             <  2:321807/148‑1 (MQ=255)                                                                                aTAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACggtggt                                                           <  5:9765/149‑1 (MQ=255)                                                                                aTAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACggtgg                                                            >  6:294227/1‑148 (MQ=255)                                                                                       tGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCg                                                    <  5:294227/149‑1 (MQ=255)                                                                                             cGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCTCGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGgccagc                                              >  5:204294/1‑149 (MQ=255)                                                                                              gAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCCTAGTTCCTTATTGCGTAAGCTTGGCCTAGACCGCCTCCAGGC‑CACGCTTGCCATCAACGGTGGTAACGCCGgccagc                                              >  2:4268/1‑148 (MQ=255)                                                                                               aGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGgccagcca                                            <  6:45442/149‑1 (MQ=255)                                                                                               aGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGgccagcca                                            <  4:104080/149‑1 (MQ=255)                                                                                                           cATAGATTGAAATCAACCCGGGGCCGGCACGGTCCTGCCCCGAGGTAGGGCTTCTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACggtg                                                             <  8:306166/120‑1 (MQ=255)                                                                                                                aTTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGc                                  >  8:16722/1‑142 (MQ=255)                                                                                                                aTTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGc                                  <  7:16722/142‑1 (MQ=255)                                                                                                                   gAAATCAACCCGGGGCAGGCTCGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTggg                         <  6:204294/148‑1 (MQ=255)                                                                                                                    aaaTCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGtt                       >  1:108063/1‑149 (MQ=255)                                                                                                                         caCCCGGTGCAGGCAGGGTCGTTCCCCGAGGTAGCGCTACTGGGTAGTTCGTTTTTGCGTAAGCTTGGCCGTGACCGCCTCCCGGC‑CAGGCTTGGGATGAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGtgc                   <  7:155488/147‑3 (MQ=255)                                                                                                                         aaCCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCt                  <  7:161844/149‑1 (MQ=255)                                                                                                                           ctcGGGGCAGGACCGTGCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTgg                          <  3:236551/137‑1 (MQ=255)                                                                                                                           ctcGGGGCAGGACCGTGCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTgg                          >  4:236551/3‑139 (MQ=255)                                                                                                                           cccGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCtt                 >  5:30601/1‑148 (MQ=255)                                                                                                                             cGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGCTCAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGa               <  1:107121/149‑1 (MQ=255)                                                                                                                              ggggCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATc             <  8:208844/149‑1 (MQ=255)                                                                                                                               gggCAGGACCGTGCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTcc                   >  6:56683/1‑142 (MQ=255)                                                                                                                               gggCAGGACCGTGCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTcc                   <  5:56683/142‑1 (MQ=255)                                                                                                                                   aggcaCGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTTCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCAACAAGCTTGCCTGGGTTTTCCTTGATCCaat         >  7:50115/1‑146 (MQ=255)                                                                                                                                   aggcaCGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATCCAttt        >  1:43217/1‑149 (MQ=255)                                                                                                                                     gcaCGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATCCATTTCt      >  4:248819/1‑149 (MQ=255)                                                                                                                                      caCGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATCCATTTCtt     >  1:145193/1‑149 (MQ=255)                                                                                                                                         ggTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATCCATTTCTTggc  >  8:62106/1‑149 (MQ=255)                                                                                                                                          |                                                                                                                                                     TCGAAGTAATTCGCGCCGTTCTGGGTCAGCCATTCGACGAACGACGC‑GATGTACTGGCCCAGGGGTATTTTCTGATCGATAAGCATGATAGCGAGCTTCCACCTGCATAGATTGAAATCAACCCGGGGCAGGCACGGTCCTGCCCCGAGGTAGCGCTACTGCGTAGTTCGTTATTGCGTAAGCTTGGCCTTGACCGCCTCCAGGC‑CAGGCTTGCCATCAACGGTGGTAACGCCGGCCAGCCAGGCATCCAGCTTGCCTGGGTTGTCCTTGATCCATTTCTTGGC  >  NC_002947/355086‑355369

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