Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,175,843	C→T	intergenic (+14/‑101)	thrT → / → tufB	tRNA‑Thr/translation elongation factor EF‑Tu 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,175,843	0	C	T	100.0%	218.4 / NA	68	intergenic (+14/‑101)	thrT/tufB	tRNA‑Thr/translation elongation factor EF‑Tu 2
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/38);  total (30/38)

TGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAA  >  NC_000913/4175713‑4175971                                                                                                                                   |                                                                                                                                 tgatgaTGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGtttt                                                                                                                          <  1:915435/139‑1 (MQ=255) tgatgaTGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGtttt                                                                                                                          >  1:1048042/1‑139 (MQ=255)       tgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCtt                                                                                                                    <  1:1096142/139‑1 (MQ=255)       tgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCtt                                                                                                                    <  1:265631/139‑1 (MQ=255)           ggTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTc                                                                                                                       <  2:52969/132‑1 (MQ=255)           ggTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTc                                                                                                                       >  1:52969/1‑132 (MQ=255)                                tCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCgg                                                                                           <  2:174251/139‑1 (MQ=255)                                tCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCgg                                                                                           >  1:1534606/1‑139 (MQ=255)                                        gtgCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTg                                                                                   <  2:1039910/139‑1 (MQ=255)                                         tgCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGc                                                                                  >  2:241741/1‑139 (MQ=255)                                          gCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTCTTGCATTCAACAAGTCGGGCATGTTGcc                                                                                 >  2:967854/1‑139 (MQ=255)                                          gCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTAt                                                                                                            >  2:474116/1‑112 (MQ=255)                                          gCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTAt                                                                                                            <  1:474116/112‑1 (MQ=255)                                          gCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGcc                                                                                 >  1:932861/1‑139 (MQ=255)                                          gCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGcc                                                                                 >  2:1324380/1‑139 (MQ=255)                                            tGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTg                                                                               <  1:697604/139‑1 (MQ=255)                                             gATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTgg                                                                              <  1:241741/139‑1 (MQ=255)                                             gATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTgg                                                                              <  2:1577193/139‑1 (MQ=255)                                                ataGCTCAGTTGGTAGAGTGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTg                                                                           <  2:456718/139‑1 (MQ=255)                                                  aGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTgct                                                                         <  1:967854/139‑3 (MQ=255)                                                  aGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGAt                                                                         <  1:1324380/139‑1 (MQ=255)                                                   gCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTg                                                                           <  2:740796/136‑1 (MQ=255)                                                   gCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTg                                                                           >  1:740796/1‑136 (MQ=255)                                                    cTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATg                                                                                      <  1:981582/124‑1 (MQ=255)                                                    cTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATg                                                                                      >  2:981582/1‑124 (MQ=255)                                                         ttGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGt                                                                                                                >  2:1151952/1‑93 (MQ=255)                                                         ttGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGt                                                                                                                <  1:1151952/93‑1 (MQ=255)                                                         ttGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGt                                                                                     >  2:126037/1‑120 (MQ=255)                                                         ttGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGt                                                                                     <  1:126037/120‑1 (MQ=255)                                                          tGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTGTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCAGTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGat                                                                 <  1:762576/139‑1 (MQ=255)                                                            gTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTCTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCTGGCATGTTGCCTGGTTGATGTGGTGatat                                                               <  2:1077965/139‑1 (MQ=255)                                                             tAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATc                                                              <  2:831173/139‑1 (MQ=255)                                                                 gcgcACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCg                                                          >  2:5939/1‑139 (MQ=255)                                                                    cACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGAtt                                                       <  1:5939/139‑1 (MQ=255)                                                                     aCCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGAttt                                                      >  2:1233705/1‑139 (MQ=255)                                                                         ttGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATcc                                                  <  2:932861/139‑1 (MQ=255)                                                                                 ggTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCt                                                                                                                     <  2:324751/64‑1 (MQ=255)                                                                                 ggTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCt                                                                                                                     >  1:324751/1‑64 (MQ=255)                                                                                 ggTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATc                                                   <  1:897830/130‑1 (MQ=255)                                                                                 ggTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATc                                                   >  2:897830/1‑130 (MQ=255)                                                                                     aGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCa                                                                                                        <  1:763869/73‑1 (MQ=255)                                                                                     aGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCa                                                                                                        >  2:763869/1‑73 (MQ=255)                                                                                     aGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGt                                                                                              >  2:428973/1‑83 (MQ=255)                                                                                     aGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGt                                                                                              <  1:428973/83‑1 (MQ=255)                                                                                           gCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCAt                                                                                       <  2:814078/84‑1 (MQ=255)                                                                                           gCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCAt                                                                                       >  1:814078/1‑84 (MQ=255)                                                                                           gCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGAt                                                                         >  1:260611/1‑98 (MQ=255)                                                                                           gCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGAt                                                                         <  2:260611/98‑1 (MQ=255)                                                                                           gCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATgtgg                                                                     <  2:863578/102‑1 (MQ=255)                                                                                           gCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATgtgg                                                                     >  1:863578/1‑102 (MQ=255)                                                                                            cAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATc                               >  2:922677/1‑139 (MQ=255)                                                                                               ttCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAgg                                      >  1:137976/1‑129 (MQ=255)                                                                                               ttCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAgg                                      <  2:137976/129‑1 (MQ=255)                                                                                                 cGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGt                          >  1:158592/1‑139 (MQ=255)                                                                                                 cGAATCTGCCTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGt                          >  1:1517968/1‑139 (MQ=255)                                                                                                         ccTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAGCAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATc                                                              <  2:1085593/95‑1 (MQ=255)                                                                                                         ccTATCAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATc                                                              >  1:1085593/1‑95 (MQ=255)                                                                                                              cAGCACCACTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGtt             <  2:1517968/139‑1 (MQ=255)                                                                                                                 caccacTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGtt                                                                            >  2:613879/1‑73 (MQ=255)                                                                                                                 caccacTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGtt                                                                            <  1:613879/73‑1 (MQ=255)                                                                                                                 caccacTTCTTCTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGAGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGa          <  2:1048042/139‑1 (MQ=255)                                                                                                                   ccacTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGAt                                                                         >  2:838719/1‑74 (MQ=255)                                                                                                                   ccacTTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGAt                                                                         <  1:838719/74‑1 (MQ=255)                                                                                                                      cttcttTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTa     >  2:1453008/1‑139 (MQ=255)                                                                                                                       ttcttTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTg                                                                                   >  2:590353/1‑60 (MQ=255)                                                                                                                       ttcttTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTg                                                                                   <  1:590353/60‑1 (MQ=255)                                                                                                                       ttcttTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTAc    <  2:1534606/139‑1 (MQ=255)                                                                                                                         cttTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACaa  <  2:158592/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                                 TGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAA  >  NC_000913/4175713‑4175971

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