Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,794,488	C→A	S79Y (TCC→TAC)	gabP →	4‑aminobutyrate:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,794,488	0	C	A	100.0%	140.3 / NA	44	S79Y (TCC→TAC)	gabP	4‑aminobutyrate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (5/39);  total (5/39)

GTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTACTGGTGGTTTTGGGTACTGGTTATCCCGCTGGAAGCCAACATCGCCGCCATGATCCTGCACTCATGGGTTCCAGGCATT  >  NC_000913/2794352‑2794624                                                                                                                                           |                                                                                                                                          gTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTacac                                                                                                                                        <  1:1279467/140‑1 (MQ=255) gTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTacac                                                                                                                                        <  1:721859/140‑1 (MQ=255)      ttCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATg                                                                                                                                   <  1:1317223/140‑1 (MQ=255)      ttCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATg                                                                                                                                   <  2:1172410/140‑1 (MQ=255)           gCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGAt                                                                                                                              <  2:460867/140‑1 (MQ=255)            cGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATa                                                                                                                             <  1:550582/140‑1 (MQ=255)               cGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAg                                                                                                                          <  2:1342211/140‑1 (MQ=255)                gCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGGTG‑CCACCCCCCCTACCGGGTCGTTTTTCACCCATTCCCCTaaaac                                                                                                                         >  2:946716/1‑138 (MQ=255)                         gAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCg                                                                                                                <  1:692018/140‑1 (MQ=255)                           aGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGGTTC‑CCCCCCCCGCAACCGGTCGGTTTTTC‑CCCAACCCGCTAAAACCACTGGCCGCt                                                                                                              >  2:495990/1‑140 (MQ=255)                             cGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTgg                                                                                                            <  1:876921/140‑1 (MQ=255)                             cGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATG‑CGATAAAGCCATTGGTCGCTggg                                                                                                           <  1:188133/140‑1 (MQ=255)                              gggCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTggg                                                                                                           <  2:1022723/140‑1 (MQ=255)                              gggCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTggg                                                                                                           <  2:1345830/140‑1 (MQ=255)                              gggCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTggg                                                                                                           <  1:1204631/140‑1 (MQ=255)                              gggCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTggg                                                                                                           <  1:3362/140‑1 (MQ=255)                                gCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTgggcg                                                                                                         <  1:495990/140‑1 (MQ=255)                                gCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGGTT‑CCCCCCCCCCTAACGGGTTGTTTTTTACCCATGCCCATAAAGACAATTGTCGCTggggg                                                                                                         >  2:1050583/1‑138 (MQ=255)                                    ggcggTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCt                                                                                                     <  2:677036/140‑1 (MQ=255)                                          aTTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑GCCCGCCCCCGACCGGGTCTTTTTTTACCCATTGCGGTAAAAACCTTTGTCGCGGGGGGGGGtataacc                                                                                               >  2:205802/1‑137 (MQ=255)                                             aCTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCg                                                                                            <  1:1050583/140‑1 (MQ=255)                                             aCTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCg                                                                                            <  2:445421/140‑1 (MQ=255)                                             aCTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCg                                                                                            <  2:1102542/140‑1 (MQ=255)                                              cTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCgg                                                                                           <  2:1573/140‑1 (MQ=255)                                              cTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCgg                                                                                           <  1:1024136/140‑1 (MQ=255)                                              cTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCgg                                                                                           <  1:92370/140‑1 (MQ=255)                                                     aTCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCCCTACCGGTTTGTTTTTCACCTATTCCGATAAAACCATTTGTCGCTGGGGGGGCTATACCATCGGCT‑GGGTg                                                                                   >  2:222095/1‑140 (MQ=255)                                                      tCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCCCCCCCGGTACCGGCTCGTTTTATACCTATGCCGGAAAAAGCATTTGGCGCTGGGCGGGGTATACCATCGGCT‑gggttt                                                                                  >  2:1231005/1‑136 (MQ=255)                                                        tGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTAc                                                                                <  1:222095/140‑1 (MQ=255)                                                           tCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCCCGCCCCATAACGGGTCGTTTTTCACCCTTTCCCATAAAACCATTGGTCGCTGGGCGGGGCATACCATCCGCT‑GGgtgttttgg                                                                             >  2:857459/1‑135 (MQ=255)                                                           tCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTACtgg                                                                             <  1:967199/140‑1 (MQ=255)                                                           tCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTACtgg                                                                             <  1:205802/140‑1 (MQ=255)                                                           tCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGGTCGTTTTTCACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTACtgg                                                                             >  2:847826/1‑140 (MQ=255)                                                            cGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGGTCGTTTTTCACCTATGCCGATAAAACCATTTGTCGCTGGGCGGGCTATACCATCGGCT‑GGGTGTACtggg                                                                            >  2:139444/1‑139 (MQ=255)                                                             gCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTACtggtg                                                                           <  2:1167035/140‑1 (MQ=255)                                                                          gggtTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTACTGGTGGTTTTGGGTACTg                                                              <  1:208940/139‑1 (MQ=255)                                                                          gggtTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTACTGGTGGTTTTGGGTACTg                                                              <  1:1210561/139‑1 (MQ=255)                                                                               ttGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTACTGGTGGTTTTGGGTACTGGTTAt                                                         <  2:1129921/139‑1 (MQ=255)                                                                                tGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCCCGCCCGATACCGGTTCGTTTTTCACCTATGCCGATAAAACCATTGGGCGCTGGGCGGGGTATACCATCGGCT‑GGGTGTACTGGGGGGTTTTGGGGCTGGTTAtt                                                        >  2:680303/1‑139 (MQ=255)                                                                                      tGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGGTACCGGGTCGTTTTTCACCTATGCCGATAAAACCATTGGGCGCTGGGCGGGGTATACCATCGGCT‑GGGTGTACTGGGGGGTTTGGGGGCTGGTTTTTCCGccc                                                  >  2:1210561/1‑138 (MQ=255)                                                                                      tGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTACTGGTGGTTTTGGGTACTGGTTATCCCGCTg                                                  <  1:753986/140‑1 (MQ=255)                                                                                      tGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACCCCCGATAACGGTTTGTTTTACACCTAAGCCGGTAAAACCATTGGGCCCTGGGCGGGGCATACCATCGGCTGGGTTTTCTTGGGGGTTTTGGG‑GCTGGTTATCCCGccg                                                  >  2:1359539/1‑138 (MQ=255)                                                                                           aTGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCCATACCGGTTCGTTTTTCACCCATGCCGGTAAAACCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTACTGGGGGGTTTTGGTGCTGGGTATTCCCCTGGGAGc                                             >  2:69430/1‑140 (MQ=255)                                                                                                ggCGGAAATGGC‑GGTTGCCCCCCCCCGTAACCGGTT‑GTTTTTCACCCATGCCGGTAAAAGCATTGGGTGCTTGGGGGGCTTTTCCATCGGCT‑GGTTGTACTGGGGGGTTTTGGGGCTGGGTATCCCCCCGGGAGGCAACa                                        >  2:577260/1‑140 (MQ=255)                                                                                                         ggC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTACTGGTGGTTTTGGGTACTGGTTATCCCGCTGGAAGCCAACATCGCCGCCa                               <  1:903264/140‑1 (MQ=255)                                                                                                         ggC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTACTGGTGGTTTTGGGTACTGGTTATCCCGCTGGAAGCCAACATCGCCGCCa                               <  1:506146/140‑1 (MQ=255)                                                                                                         ggC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTACTGGTGGTTTTGGGTACTGGTTATCCCGCTGGAAGCCAACATCGCCGCCa                               <  1:139444/140‑1 (MQ=255)                                                                                                                 g‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTACTGGTGGTTTTGGGTACTGGTTATCCCGCTGGAAGCCAACATCGCCGCCATGATCCt                        <  1:577260/140‑1 (MQ=255)                                                                                                                 g‑CCACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTACTGGTGGTTTTGGGTACTGGTTATCCCGCTGGAAGCCAACATCGCCGCCATGATCCt                        <  2:59875/140‑1 (MQ=255)                                                                                                                 g‑CCACGCCCGATACCGGTTCGCTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTACTGGTGGTTTTGGGTACTGGTTATCCCGCTGGAAGCCAACATCGCCGCCATGATCCt                        <  1:1359539/140‑1 (MQ=255)                                                                                                                    cACGCCCGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTACTGGTGGTTTTGGGTACTGGTTATCCCGCTGGAAGCCAACATCGCCGCCATGATCCTGc                      <  1:1189799/140‑1 (MQ=255)                                                                                                                          cGATACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTACTGGTGGTTTTGGGTACTGGTTATCCCGCTGGAAGCCAACATCGCCGCCATGATCCTGCACTCAt                <  1:69430/140‑1 (MQ=255)                                                                                                                            aTACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTACTGGTGGTTTTGGGTACTGGTTATCCCGCTGGAAGCCAACATCGCCGCCa                               <  2:237925/123‑1 (MQ=255)                                                                                                                            aTACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTACTGGTGGTTTTGGGTACTGGTTATCCCGCTGGAAGCCAACATCGCCGCCa                               >  1:237925/1‑123 (MQ=255)                                                                                                                            aTACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTACTGGTGGTTTTGGGTACTGGTTATCCCGCTGGAAGCCAACATCGCCGCCATGATCCTGCACTCATgg              >  2:890135/1‑140 (MQ=255)                                                                                                                             tACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTACTGGTGGTTTTGGGGACTGGTTATCCCGCTGGAAGCCAACATCGCCGCCATGATCCTGCACTCATggg             >  1:696637/1‑140 (MQ=255)                                                                                                                             tACCGGTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTACTGGGGGTTTTGGGTACTGGTTATCCCGCTGGAAGCCAACATCGCCGCCATGATCCTGCACTCATggg             >  2:522996/1‑140 (MQ=255)                                                                                                                               ccGGTTCATTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTACTGGTGGTTTTGGGTACTGGTTATCCCGCTGGAAGCCAACATCGCCGCCATGATCCTGCACTCATGGGtt           <  1:415361/140‑1 (MQ=255)                                                                                                                                 ggTTCGTTTTACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTACTGGTGGTTTTGGGTACTGGTTATCCCGCTGGAAGCCAACATCGCCGCCATGATCCTGCACTCATGGGTTcc         >  2:525851/1‑140 (MQ=255)                                                                                                                                        tttACACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTACTGGTGGTTTTGGGTACTGGTTATCCCGCTGGAAGCCAACATCGCCGCCATGATCCTGCACTCATGGGTTCCAGGCAtt  <  1:890135/140‑1 (MQ=255)                                                                                                                                           |                                                                                                                                          GTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGC‑GGTTG‑CCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCT‑GGTTGTACTGGTGGTTTTGGGTACTGGTTATCCCGCTGGAAGCCAACATCGCCGCCATGATCCTGCACTCATGGGTTCCAGGCATT  >  NC_000913/2794352‑2794624

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