Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,363,064	A→C	8.6%	D234E (GAT→GAG)	glmM ←	phosphoglucosamine mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,363,064	0	A	C	8.6%	173.5 / 7.6	58	D234E (GAT→GAG)	glmM	phosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base A (23/30);  new base C (2/3);  total (25/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

ACCGACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATACGCTGGGCGCAACCTTG  >  NC_002947/5362923‑5363208                                                                                                                                              |                                                                                                                                                  aCCGACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGATCTCGGCAt                                                                                                                                            >  1:259085/1‑149 (MQ=255)   cGACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGATCTCGGCAttg                                                                                                                                          <  2:259085/149‑3 (MQ=255)      ccaccCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGATCTCGGCAttgcct                                                                                                                                       <  2:270441/149‑6 (MQ=255)         ccccGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGAcc                                                                                                                                    <  3:258204/149‑1 (MQ=255)         ccccGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGAc                                                                                                                                     >  2:30870/1‑148 (MQ=255)           ccGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAg                                                                                                                                  <  5:46169/149‑1 (MQ=255)             gCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTa                                                                                                                                <  1:30870/149‑1 (MQ=255)               cTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGAc                                                                                                                                     <  3:95851/142‑1 (MQ=255)               cTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGAc                                                                                                                                     >  4:95851/1‑142 (MQ=255)                   cAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCg                                                                                                                           <  7:4894/148‑1 (MQ=255)                   cAGATTGCCATGGTCCTTCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGTATCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACATCg                                                                                                                           <  4:204907/148‑1 (MQ=255)                           cATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGc                                                                                                                 >  7:72431/1‑149 (MQ=255)                                  cTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGa                                                                                                          <  1:48439/149‑1 (MQ=255)                                  cTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGa                                                                                                          >  1:213784/1‑149 (MQ=255)                                         tCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGAt                                                                                                         <  3:263616/143‑1 (MQ=255)                                         tCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGAt                                                                                                         >  4:263616/1‑143 (MQ=255)                                           gcgAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTc                                                                                                  >  8:248180/1‑148 (MQ=255)                                            cgAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGa                                                                                                <  1:66887/149‑1 (MQ=255)                                                  atgatgAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGAtgtgt                                                                                                     >  6:191702/1‑138 (MQ=255)                                                  atgatgAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGAtgtgt                                                                                                     <  5:191702/138‑1 (MQ=255)                                                     atgaACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAAcc                                                                                       >  2:186842/1‑149 (MQ=255)                                                        aaCAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCtt                                                                                     <  8:21704/148‑1 (MQ=255)                                                           agcagcTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGAtt                                                                                                             >  2:207963/1‑121 (MQ=255)                                                           agcagcTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGAtt                                                                                                             <  1:207963/121‑1 (MQ=255)                                                             cagcTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGa                                                                               <  7:248180/149‑1 (MQ=255)                                                                      aCCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAgg                                                                       <  5:76246/148‑1 (MQ=255)                                                                      aCCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGc                                                                      >  7:167488/1‑149 (MQ=255)                                                                        cGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCg                                                                    >  2:6169/1‑149 (MQ=255)                                                                             acgaTGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCgg                                                               <  6:128194/149‑1 (MQ=255)                                                                              cgaTGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGAc                                                                                                                                     >  5:226205/1‑79 (MQ=255)                                                                              cgaTGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGAc                                                                                                                                     <  6:226205/79‑1 (MQ=255)                                                                               gaTGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGtt                                                             >  6:3234/1‑149 (MQ=255)                                                                                 tGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGTCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTgc                                                           <  7:200723/149‑1 (MQ=255)                                                                                 tGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGAGGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTgc                                                           <  5:294338/149‑1 (MQ=255)                                                                                      cccGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATg                                                      <  8:72431/149‑1 (MQ=255)                                                                                       ccGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATg                                                      >  1:141477/1‑148 (MQ=255)                                                                                             ggTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGtt                                                             <  1:310675/135‑1 (MQ=255)                                                                                             ggTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGtt                                                             >  2:310675/1‑135 (MQ=255)                                                                                                  aCCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCAc                                          <  1:62830/149‑1 (MQ=255)                                                                                                     aTCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTc                                       <  1:186842/149‑1 (MQ=255)                                                                                                         gCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCa                                   >  5:78627/1‑149 (MQ=255)                                                                                                                   gCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCGGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCAc                         >  1:60365/1‑149 (MQ=255)                                                                                                                   gCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCa                          <  5:3234/148‑1 (MQ=255)                                                                                                                    ccgtcaccgtcaAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACg                        <  6:78627/149‑1 (MQ=255)                                                                                                                     cgtcaccgtcaAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACg                        >  5:301667/1‑148 (MQ=255)                                                                                                                        caccgtcaAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAAt                    <  2:141477/149‑1 (MQ=255)                                                                                                                         accgtcaAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATa                   >  2:289566/1‑149 (MQ=255)                                                                                                                           cgtcaAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGATCAGGCCGTCCGGTAGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATa                   >  1:78823/1‑148 (MQ=255)                                                                                                                             tcaAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATACGCt               <  7:152498/149‑1 (MQ=255)                                                                                                                               aaaGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATACGCTg              <  8:167488/148‑1 (MQ=255)                                                                                                                                aaGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATACGCTggg            >  6:20769/1‑149 (MQ=255)                                                                                                                                 aGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATACGCTGGgc           <  2:213784/149‑1 (MQ=255)                                                                                                                                   gacaTGCCGATCTCGGCATGTCTGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACg                                               <  5:29033/108‑1 (MQ=255)                                                                                                                                   gacaTGCCGATCTCGGCATGTCTGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACg                                               >  6:29033/4‑111 (MQ=255)                                                                                                                                       tGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGt                                             <  4:74132/109‑1 (MQ=255)                                                                                                                                       tGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGt                                             >  3:74132/1‑109 (MQ=255)                                                                                                                                       gtgcGAGATCGGGATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGTAATACGCTGGGCGCTAcc     <  4:71594/146‑1 (MQ=255)                                                                                                                                          cGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATACGCTGGGCGCAACCTTg  >  1:9712/1‑149 (MQ=255)                                                                                                                                           gagaTCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATACGCTGGGCGCAACCTTg  >  5:41747/1‑148 (MQ=255)                                                                                                                                              |                                                                                                                                                  ACCGACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGC‑AGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATACGCTGGGCGCAACCTTG  >  NC_002947/5362923‑5363208

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