Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	306,996	G→A	N779N (AAC→AAT)	ecpC ←	ECP production outer membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	306,996	0	G	A	100.0%	225.0 / NA	58	N779N (AAC→AAT)	ecpC	ECP production outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/36);  total (22/36)

GATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATAGAGAGTC  >  NC_000913/306850‑307143                                                                                                                                                   |                                                                                                                                                    gATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCAtttt                                                                                                                                                 <  2:740387/151‑1 (MQ=255)  aTGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTc                                                                                                                                                >  1:354947/1‑151 (MQ=255)    gagaGGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCAt                                                                                                                                              <  1:745405/151‑1 (MQ=255)             cTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCCGGATTTATTGCCACTGTAGCGAAAATCGATAGAGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCggttcgg                                                                                                                                      >  2:87912/1‑150 (MQ=255)               gCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCggttcgggtt                                                                                                                                   <  2:345687/151‑1 (MQ=255)                  gACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCggttcgggttcg                                                                                                                                 <  1:476984/150‑1 (MQ=255)                        aTCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGata                                                                                                                          <  2:409846/151‑1 (MQ=255)                          cACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGatat                                                                                                                         <  2:707604/150‑1 (MQ=255)                          cACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATg                                                                                                                        <  2:7415/151‑1 (MQ=255)                           aCCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATg                                                                                                                        <  2:509250/150‑1 (MQ=255)                                cccAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGgttgtt                                                                                                                  <  1:701685/151‑1 (MQ=255)                                cccAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGgttgtt                                                                                                                  >  2:355320/1‑151 (MQ=255)                                cccAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGgttgt                                                                                                                   <  1:372307/150‑1 (MQ=255)                                     aCGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATcc                                                                                                             <  2:178435/151‑1 (MQ=255)                                           ccGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGt                                                                                                       <  1:199256/151‑1 (MQ=255)                                              cgcgcCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAg                                                                                                    <  2:455055/151‑1 (MQ=255)                                               gcgcCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGc                                                                                                   >  2:272792/1‑151 (MQ=255)                                                 gcCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGcc                                                                                                  >  2:232836/1‑150 (MQ=255)                                                 gcCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCca                                                                                                 <  1:380991/151‑1 (MQ=255)                                                   cTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCcagc                                                                                               >  1:4209/1‑151 (MQ=255)                                                     ggTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCcagcag                                                                                             <  2:456939/151‑1 (MQ=255)                                                     ggTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCcagca                                                                                              <  2:749600/150‑1 (MQ=255)                                                     ggTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCcagc                                                                                               >  1:618441/1‑149 (MQ=255)                                                     ggTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCcagc                                                                                               >  2:317120/1‑149 (MQ=255)                                                          aGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCAGGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGc                                                                                        <  1:81318/151‑1 (MQ=255)                                                            ttCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCg                                                                                      >  2:348171/1‑151 (MQ=255)                                                              ccAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTc                                                                                    >  2:202104/1‑151 (MQ=255)                                                                     cACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGc                                                                              >  1:532285/1‑150 (MQ=255)                                                                           gCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGa                                                                         <  1:326781/149‑1 (MQ=255)                                                                           gCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGAt                                                                        >  1:117356/1‑150 (MQ=255)                                                                           gCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATa                                                                       <  1:684249/151‑1 (MQ=255)                                                                           gCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATa                                                                       <  2:745206/151‑1 (MQ=255)                                                                                 tttATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGccc                                                                  <  1:88237/150‑1 (MQ=255)                                                                                    aTTGCCGCTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCTGGTTCGGCCGATATGGTTGTTAAGCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGgag                                                              <  1:672467/151‑1 (MQ=255)                                                                                       gCCACTGTAGCGAAAATCGATAGTGGGGAATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACg                                                           >  1:734901/1‑151 (MQ=255)                                                                                        ccACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACgg                                                          >  1:338130/1‑151 (MQ=255)                                                                                            tGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTaa                                                       >  1:599344/1‑150 (MQ=255)                                                                                               aGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCAt                                                   >  1:419352/1‑151 (MQ=255)                                                                                                 cGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCt                                                 <  2:21550/151‑1 (MQ=255)                                                                                                  gAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTg                                                >  1:130186/1‑151 (MQ=255)                                                                                                          aTAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACct                                        >  2:177690/1‑151 (MQ=255)                                                                                                            aGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACctct                                      <  1:504026/151‑1 (MQ=255)                                                                                                               ggggTATTTCTTATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTgg                                    <  2:181979/150‑1 (MQ=255)                                                                                                                          tATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGAc                           <  2:690389/148‑1 (MQ=255)                                                                                                                          tATCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACa                          <  1:664614/149‑1 (MQ=255)                                                                                                                           aTCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGc                        <  2:245403/150‑1 (MQ=255)                                                                                                                            tCCACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCg                       >  2:184982/1‑150 (MQ=255)                                                                                                                             ccACGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGAc                     <  1:274889/151‑1 (MQ=255)                                                                                                                                cGTCCATGACAAACTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACAtt                  <  2:361933/151‑1 (MQ=255)                                                                                                                                          aaaCTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGAt          >  2:569959/1‑149 (MQ=255)                                                                                                                                          aaaCTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATag        <  2:485258/151‑1 (MQ=255)                                                                                                                                          aaaCTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATag        <  2:386020/151‑1 (MQ=255)                                                                                                                                             cTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATag        >  2:176025/1‑148 (MQ=255)                                                                                                                                             cTCGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATag        >  1:317385/1‑148 (MQ=255)                                                                                                                                               cGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGTGACATTGCCTGGATagagag    <  2:130219/150‑1 (MQ=255)                                                                                                                                               cGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATagagag    <  2:381558/150‑1 (MQ=255)                                                                                                                                               cGCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATAGAGAGt   <  2:125500/151‑1 (MQ=255)                                                                                                                                                gCCATTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATAGAGAGTc  <  1:501977/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GATGAGAGGCCGCTGCAGACCACATCACCGACCCAGACGGCACCGCGCGCCTGGTTGAGTTCCAGAGCCACTTCGCAGGTTTTATTGCCACTGTAGCGAAAATCGATAGTGGGGTATTTCTTATCCACGTCCATGACAAACTCGCCGTTTTCATCGGTTCGGGTTCGGCCGATATGGTTGTTAATCCGTGCGTTAGCCAGCAGTGTGCCGTCTTCCGCACGGATACGACCGGAGACGGTAACCATCTGCTTCACCTCTGGCTCAATGACAGCGACATTGCCTGGATAGAGAGTC  >  NC_000913/306850‑307143

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