Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,740,795	T→G	10.3%	intergenic (+32/+721)	gltA → / ← yeiW	citrate synthase/putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,740,795	0	T	G	10.3%	139.2 / 6.5	56	intergenic (+32/+721)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (24/26);  new base G (3/3);  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 = 7.94e-01

GATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACAC‑CGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACT  >  NC_002947/4740652‑4740942                                                                                                                                                 |                                                                                                                                                    gATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACAC‑CGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCgg                                                                                                                                                  <  6:174363/147‑1 (MQ=255)       gCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACAC‑CGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCa                                                                                                                                               <  7:285584/144‑1 (MQ=255)       gCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACAC‑CGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtt                                                                                                                                          <  4:49036/149‑1 (MQ=255)        gaCTGGAAAGAAATGCTCTCCAGCCGGTACAAGATCGGCCGTCCGGGCGAGCTGTAGAC‑CGGCGAGCAGAAGCGCGTCATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAAAAGGCTGCCCGAGGGCAGCCttc                                                                                                                                         <  6:293594/148‑2 (MQ=255)         aCTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACAC‑CGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttt                                                                                                                                        <  7:194449/149‑1 (MQ=255)             gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACAC‑CGGCGAGCAGAAGCGCGCCATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                    >  3:225903/1‑147 (MQ=255)             gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACAC‑CGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttttta                                                                                                                                     <  4:46930/148‑2 (MQ=255)              aaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACAC‑CGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                    <  8:153130/148‑3 (MQ=255)                  aaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACAC‑CGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                    <  2:90845/144‑3 (MQ=255)                  aaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACAC‑CGGCAAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                    >  1:90845/1‑142 (MQ=255)                   aaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACAC‑CGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatttgtc                                                                                                                               <  3:125796/148‑8 (MQ=255)                       ctctCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACAC‑CGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                       >  7:293857/1‑136 (MQ=255)                       ctctCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACAC‑CGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                       <  8:293857/136‑1 (MQ=255)                          tcCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACAC‑CGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCg                                                                                                                      >  6:342455/1‑149 (MQ=255)                               cccGTACAAGATCGGCCGTCCGCGCCAGCTGTACAC‑CGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCT‑GTGCCGGCctct                                                                                                                >  6:144511/1‑149 (MQ=255)                                ccGTACAAGATCGGCCGTCCGCGCCAGCTGTACAC‑CGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctct                                                                                                                <  6:13883/149‑1 (MQ=255)                                    aCAAGATCGGCCGTCCGCGCCAGCTGTACAC‑CGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                            >  7:102713/1‑149 (MQ=255)                                     cAAGATCGGCCGTCCGCGCCAGCTGTACACGCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGCCCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGATGCCGGCGGGCCGCC‑TTTTTTATTTTTCCTGTGCCGGCCTCTTcgc                                                                                                            >  1:205637/1‑149 (MQ=255)                                        gATCGGCCGTCCGCGCCAGCTGTACAC‑CGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGcc                                                                                                        <  6:162239/149‑1 (MQ=255)                                              ccGTCCGCGCCAGCTGTACAC‑CGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGccc                                                                                                  >  7:216544/1‑149 (MQ=255)                                              ccGTCCGCGCCAGCTGTACAC‑CGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAG‑CGTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGccc                                                                                                  >  5:298903/1‑149 (MQ=255)                                               cGTCCGCGCCAGCTGTACAC‑CGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCcg                                                                                                 <  8:102713/149‑1 (MQ=255)                                                   cgcgCCAGCTGTACAC‑CGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgcg                                                                                                           <  5:105525/135‑1 (MQ=255)                                                   cgcgCCAGCTGTACAC‑CGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgcg                                                                                                           >  6:105525/1‑135 (MQ=255)                                                   cgcgCCAGCTGTACAC‑CGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                             <  8:206861/149‑1 (MQ=255)                                                   cgcgCCAGCTGTACAC‑CGGCGAGAAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTACTGTGCCGGCCACATCGCGGCCATGACCACGa                                                                                              >  7:216171/1‑148 (MQ=255)                                                     agcCAGCTGTACAC‑CGTCTAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTACGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                           <  8:216544/148‑1 (MQ=255)                                                        cAGCTGTACAC‑CGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCGCTTAGCGGGGATGGCCGCGaaccggg                                                                                        >  1:336142/1‑144 (MQ=255)                                                              tACAC‑GGGCTACCATAAGCGGCACATCGTTGCCCTTAAACACCAAAATTATGTTGCATTCCCCTAACCAAATACGGCGGCCCTCGGGCTGCTTTTTTTTTTTGTCCTGG‑GCCGGCCTCTTCGCGGCCTTGCCCGCTATGTGGCCGGca                                                                                   <  6:298903/148‑1 (MQ=255)                                                              tACAC‑CGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGca                                                                                   <  5:342455/148‑1 (MQ=255)                                                                 ac‑CGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAgg                                                                                         <  8:37483/139‑1 (MQ=255)                                                                 ac‑CGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCAggtcggaa                                                                               >  7:37483/1‑143 (MQ=255)                                                                     ggCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGacac                                                                            <  5:285397/149‑1 (MQ=255)                                                                     ggCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGacac                                                                            >  6:285397/1‑149 (MQ=255)                                                                        gAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGAcacca                                                                          >  3:294558/1‑148 (MQ=255)                                                                           cAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAAAAGGCTGCCCGAGGGCAGcc                                                                                                                                            >  6:175377/1‑80 (MQ=39)                                                                           cAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAAAAGGCTGCCCGAGGGCAGcc                                                                                                                                            <  5:175377/80‑1 (MQ=39)                                                                                    gACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGtt                                                             >  1:33037/1‑149 (MQ=255)                                                                                              cccTGAAAGACCGCAAGTACGCTGCAGTAGCCGACCGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTGCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgcc                                                   <  2:205637/149‑1 (MQ=255)                                                                                                       aCCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGc                                          >  8:186216/1‑149 (MQ=255)                                                                                                           cAAGTAAGCTGCAGTAGCCTAACGCAAGAAGGCTGCCGTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGTCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGcc                                      <  7:84080/149‑1 (MQ=255)                                                                                                              gTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACCGACACCACAAAACCACAGTTTTTCTCCGCCCGCTCACGCAccccgt                                   >  6:282645/1‑147 (MQ=255)                                                                                                                        gTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAtt                         >  2:226838/1‑149 (MQ=255)                                                                                                                           gCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGTCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                              <  1:162278/141‑1 (MQ=255)                                                                                                                           gCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCATGCCCCCGAAGAGGGGGGCACACACACCACAAAACCACACTTTTTATACGCCCCCTCACCCAGCCCATTAAgggagatcgg                      >  2:162278/1‑142 (MQ=255)                                                                                                                                    aagaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATc             >  1:210796/1‑149 (MQ=255)                                                                                                                                    aagaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGGTTTGAACGGCGCATc             >  5:81386/1‑149 (MQ=255)                                                                                                                                    aacacggCTGCCGTC‑GGCAAGCTTTTTTTATTTGTCCTGGCCCGCCCTCTCCGCGCCCTGGCCCGCTATGTGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATc             <  5:43071/144‑1 (MQ=255)                                                                                                                                     ccaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTCCCGCCCCCTTCGCGGCCTTCCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTTAACGGCGCATcc            <  4:225903/147‑1 (MQ=255)                                                                                                                                       aagGCTGCCCGAGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCa                                                                >  2:203670/1‑95 (MQ=25)                                                                                                                                       aagGCTGCCCGAGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCa                                                                <  1:203670/95‑1 (MQ=25)                                                                                                                                          gCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCcaccac       <  4:294558/149‑1 (MQ=255)                                                                                                                                            tGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCACACGCAGCCCTTTCAGGGTATTGAAAGGCGCATCAACAACg      >  4:56901/1‑148 (MQ=255)                                                                                                                                              cccTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGaa    <  3:56901/148‑1 (MQ=255)                                                                                                                                               ccTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACGTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCAACAACGAAca  >  3:187322/1‑148 (MQ=255)                                                                                                                                               ccTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACt  >  6:337407/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                    GATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACAC‑CGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACT  >  NC_002947/4740652‑4740942

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