Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	409,497	T→C	R571R (CGA→CGG)	sbcC ←	exonuclease subunit SbcC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	409,497	0	T	C	96.8%	208.2 / ‑1.0	62	R571R (CGA→CGG)	sbcC	exonuclease subunit SbcC
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (28/32);  total (29/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.78e-01

GCTGAGTAACCGCAGCTGGCGTTCGTGCTCATCTTGTGCATCCAGCCACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGCACCTTCTTCACCGAGCTTTTTAACTTCGTTTTCCAGCGCCAGTAATCGAGACTGAT  >  NZ_CP009273/409362‑409627                                                                                                                                        |                                                                                                                                   gCTGAGTAACCGCAGCTGGCGTTCGTGCTCATCTTGTGCATCCAGCCACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCgg                                                                                                                                 <  1:205954/139‑1 (MQ=255)               gCTGGCGTTCGTGCTCATCTTGTGCATCCAGCCACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTc                                                                                                                   >  2:883110/1‑139 (MQ=255)               gCTGGCGTTCGTGCTCATCTTGTGCATCCAGCCACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTc                                                                                                                   >  1:173981/1‑139 (MQ=255)                        cGTGCTCATCTTGTGCATCCAGCCACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTtcat                                                                                                           <  2:868089/138‑1 (MQ=255)                        cGTGCTCATCTTGTGCATCCAGCCACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTtcat                                                                                                           >  1:309789/1‑138 (MQ=255)                        cGTGCTCATCTTGTGAATCCAGCCACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTtcat                                                                                                           <  2:455555/138‑1 (MQ=255)                                     tgCATCCAGCCACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAgctgc                                                                                             >  1:475680/1‑139 (MQ=255)                                     tgCATCCAGCCACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAgctg                                                                                              <  2:309789/138‑1 (MQ=255)                                       cATCCAGCCACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAgctgct                                                                                            <  2:225228/138‑1 (MQ=255)                                       cATCCAGCCACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAgctgct                                                                                            <  2:527217/138‑1 (MQ=255)                                         tccatccaCTGTTTAATATCGTCCATTGGCTGCATGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAATTTCTTGCTCATCTTGCCGGTTTCTTTGCGCTTCGTTTTCTTCACGCTGCAGCTGCTTTg                                                                                         >  2:133720/1‑139 (MQ=255)                                         tCCAGCCACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCgg                                                                                                                                 >  2:734224/1‑99 (MQ=255)                                         tCCAGCCACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCgg                                                                                                                                 <  1:734224/99‑1 (MQ=255)                                         tCCAGCCACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTg                                                                                         >  1:503862/1‑139 (MQ=255)                                          ccAGCCACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGt                                                                                        >  2:1282851/1‑139 (MQ=255)                                            aGCCACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTa                                                                                      <  2:155974/139‑1 (MQ=255)                                              ccACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAgctg                                                                                              <  2:122278/129‑1 (MQ=255)                                              ccACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAgctg                                                                                              >  1:122278/1‑129 (MQ=255)                                              ccACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATg                                                                                    <  2:1031977/139‑1 (MQ=255)                                                   gTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTc                                                                               <  1:725655/139‑1 (MQ=255)                                                      gAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCa                                                                             <  2:890725/138‑1 (MQ=255)                                                      gAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCa                                                                             <  2:1085995/138‑1 (MQ=255)                                                       aaTATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGt                                                                           >  1:437483/1‑139 (MQ=255)                                                           tCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGc                                                                        <  2:161602/138‑1 (MQ=255)                                                              tCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTa                                                                                      >  1:652712/1‑121 (MQ=255)                                                              tCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTa                                                                                      <  2:652712/121‑1 (MQ=255)                                                              tCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCa                                                                             >  2:1182325/1‑130 (MQ=255)                                                              tCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCa                                                                             <  1:1182325/130‑1 (MQ=255)                                                              tCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTg                                                                         <  1:1174483/134‑1 (MQ=255)                                                              tCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTg                                                                         >  2:1174483/1‑134 (MQ=255)                                                              tCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCAc                                                                    >  2:1138727/1‑139 (MQ=255)                                                              tCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCAc                                                                    >  2:1290369/1‑139 (MQ=255)                                                                   tGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGc                                                               >  1:88870/1‑139 (MQ=255)                                                                   tGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCAAGTaac                                                               >  2:685088/1‑137 (MQ=255)                                                                      cTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTc                                                            <  2:56532/139‑1 (MQ=255)                                                                        gCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTa                                                                 <  1:281476/132‑1 (MQ=255)                                                                        gCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTa                                                                 >  2:281476/1‑132 (MQ=255)                                                                         cAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGCa                                                         <  2:475680/139‑1 (MQ=255)                                                                            ggTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGCACct                                                      <  2:88870/139‑1 (MQ=255)                                                                                        gCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGc                                                                        <  1:930400/109‑1 (MQ=255)                                                                                        gCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGc                                                                        >  2:930400/1‑109 (MQ=255)                                                                                        gCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGCACCTTCTTCACCGAGc                                          >  2:378210/1‑139 (MQ=255)                                                                                           ggCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCg                                                           >  1:796361/1‑119 (MQ=255)                                                                                           ggCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCg                                                           <  2:796361/119‑1 (MQ=255)                                                                                            gCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGc                                                          <  1:795934/119‑1 (MQ=255)                                                                                            gCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGc                                                          >  2:795934/1‑119 (MQ=255)                                                                                                      cTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTc                                                                               <  1:1148740/88‑1 (MQ=255)                                                                                                      cTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTc                                                                               >  2:1148740/1‑88 (MQ=255)                                                                                                      cTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGCACCTTCTTCACCGAGCTTTTTAACTTCGtt                            >  1:990120/1‑139 (MQ=255)                                                                                                             ttgttgAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGCACCTTCTTCACCGAGCTTTTTAACTTCGTTTTCCAgc                     <  2:495594/139‑1 (MQ=255)                                                                                                                ttgAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGCACCTTCTTCACCGAGCTTTTTAACTTCGTTTTCCAgcgc                   <  1:883110/138‑1 (MQ=255)                                                                                                                ttgAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGCACCTTCTTCACCGAGCTTTTTAACTTCGTTTTCCAgcgc                   <  1:951680/138‑1 (MQ=255)                                                                                                                  gAGTAAGTGCTTGCTCATCTTGCTGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGCACCTTCTTCACCGAGCTTTTTAACTTCGTTTTCCAGCGCCa                 <  2:512041/138‑1 (MQ=255)                                                                                                                  gAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGCACCTTCTTCACCGAGCTTTTTAACTTCGTTTTCCAGCGCCa                 <  2:510671/138‑1 (MQ=255)                                                                                                                  gAGTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGCACCTTCTTCACCGAGCTTTTTAACTTCGTTTTCCAGCGCCAg                >  1:181709/1‑139 (MQ=255)                                                                                                                    gTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGCAc                                                        >  1:95484/1‑97 (MQ=255)                                                                                                                    gTAAGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGCAc                                                        <  2:95484/97‑1 (MQ=255)                                                                                                                       aGTGCTTGCTCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGCACCTTCTTCACCGAGCTTTTTAACTTCGTTTTCCAGCGCCAGTAAt            <  1:1276132/138‑1 (MQ=255)                                                                                                                          gcttgctCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGCACCTTCTTCACCGAGCTTTTTAACTTCGTTTTCCAGCGCCAGTAATc           >  1:697208/1‑136 (MQ=255)                                                                                                                          gcttgctCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGCACCTTCTTCACCGAGC‑TTTTAACTTCGTTTTCCAGCGCCAGTAANNnnga       >  2:641081/1‑139 (MQ=255)                                                                                                                             tgctCATCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGCACCTTCTTCACCGAGCTTTTTAACTTCGTTTTCCAGCGCCAGTAATCGAGACt     >  1:40649/1‑139 (MQ=255)                                                                                                                                  aTCTTGCCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGCACCTTCTTCACCGAGCTTTTTAACTTCGTTTTCCAGCGCCAGTAATCGAGACTGAt  <  1:378210/137‑1 (MQ=255)                                                                                                                                        |                                                                                                                                   GCTGAGTAACCGCAGCTGGCGTTCGTGCTCATCTTGTGCATCCAGCCACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGCACCTTCTTCACCGAGCTTTTTAACTTCGTTTTCCAGCGCCAGTAATCGAGACTGAT  >  NZ_CP009273/409362‑409627

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