Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,658,893	G→A	*176* (TAG→TAA)	gadE →	gad regulon transcriptional activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,658,893	0	G	A	100.0%	249.2 / NA	77	*176* (TAG→TAA)	gadE	gad regulon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (37/40);  total (37/40)

CGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAGTCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTATCCACCAGTA  >  NC_000913/3658781‑3659023                                                                                                                 |                                                                                                                                   cGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGAt                                                                                                          >  2:740410/1‑139 (MQ=255)    tCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGAtggc                                                                                                       <  1:138369/139‑1 (MQ=255)     cAATATTTCGCAACAAACATTAAAGATGCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCCTATTTTTAATCAGGACATAAGCAACTGAAATTGAtggct                                                                                                      <  1:136666/139‑1 (MQ=255)        tatTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGAtggctggc                                                                                                   <  1:912966/139‑1 (MQ=255)                  caaacaTTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGAGATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTg                                                                                                                  >  1:833612/1‑114 (MQ=255)                  caaacaTTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGAGATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTg                                                                                                                  <  2:833612/114‑1 (MQ=255)                     acaTTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATg                                                                                      >  2:388919/1‑139 (MQ=255)                             gATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCt                                                                              <  2:452093/139‑1 (MQ=255)                               ttCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGa                                                                            >  2:958252/1‑139 (MQ=255)                                tCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                             >  2:843045/1‑121 (MQ=255)                                tCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                             <  1:843045/121‑1 (MQ=255)                                tCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGAt                                                                           >  1:1180908/1‑139 (MQ=255)                                         cAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGAtgg                                                                                                        >  2:447853/1‑101 (MQ=255)                                         cAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGAtgg                                                                                                        <  1:447853/101‑1 (MQ=255)                                                aTGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                             >  1:13441/1‑105 (MQ=255)                                                aTGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                             <  2:13441/105‑1 (MQ=255)                                                 tGTACAATCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCa                                                                                    >  2:938120/1‑113 (MQ=255)                                                 tGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCa                                                                                    <  1:938120/113‑1 (MQ=255)                                                  gTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGAt                                                         <  1:740410/139‑1 (MQ=255)                                                  gTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGAt                                                         <  1:882765/139‑1 (MQ=255)                                                       aGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTa                                                    >  2:265794/1‑139 (MQ=255)                                                       aGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTa                                                    <  1:388919/139‑1 (MQ=255)                                                        gCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                             >  2:100826/1‑97 (MQ=255)                                                        gCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                             <  1:100826/97‑1 (MQ=255)                                                        gCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTg                                                                             <  1:576566/113‑1 (MQ=255)                                                        gCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTg                                                                             >  2:576566/1‑113 (MQ=255)                                                        gCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTaa                                                   >  1:148725/1‑139 (MQ=255)                                                        gCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTaa                                                   >  1:1145365/1‑139 (MQ=255)                                                           gAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAgg                                                >  2:1275565/1‑139 (MQ=255)                                                             aaCTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATg                                                                                                <  2:1287596/89‑1 (MQ=255)                                                             aaCTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATg                                                                                                >  1:1287596/1‑89 (MQ=255)                                                             aaCTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGAtt                                                                          <  2:7181/111‑1 (MQ=255)                                                             aaCTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGAtt                                                                          >  1:7181/1‑111 (MQ=255)                                                             aaCTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGgag                                              >  2:496405/1‑139 (MQ=255)                                                                  aGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGAt                                                         <  2:218193/123‑1 (MQ=255)                                                                  aGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGAt                                                         >  1:218193/1‑123 (MQ=255)                                                                  aGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGga                                               >  2:835874/1‑133 (MQ=255)                                                                  aGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGga                                               <  1:835874/133‑1 (MQ=255)                                                                        aTGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACg                                   >  1:624773/1‑139 (MQ=255)                                                                         tGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGa                                  <  2:1180908/139‑1 (MQ=255)                                                                          gAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGAtggctggc                                                                                                   >  2:319185/1‑73 (MQ=255)                                                                          gAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGAtggctggc                                                                                                   <  1:319185/73‑1 (MQ=255)                                                                          gAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGAc                                 >  1:459360/1‑139 (MQ=255)                                                                           aGCGACATCGTCACCCTGGGAATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACt                                >  1:60034/1‑139 (MQ=255)                                                                             cGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTAccc                                                                 >  1:636314/1‑104 (MQ=255)                                                                             cGACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTAccc                                                                 <  2:636314/104‑1 (MQ=255)                                                                              gACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTa                             <  1:738943/139‑1 (MQ=255)                                                                              gACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTa                             <  2:624773/139‑1 (MQ=255)                                                                              gACATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTa                             <  1:265794/139‑1 (MQ=255)                                                                               aCATCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTAt                            >  2:524789/1‑139 (MQ=255)                                                                                 aTCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGAt                                                                                 >  2:1129151/1‑84 (MQ=255)                                                                                 aTCGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGAt                                                                                 <  1:1129151/84‑1 (MQ=255)                                                                                   cGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                             <  2:1040561/70‑1 (MQ=255)                                                                                   cGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                             >  1:1040561/1‑70 (MQ=255)                                                                                   cGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTccc                        <  2:459360/139‑1 (MQ=255)                                                                                   cGTCACCCTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTccc                        <  1:958252/139‑1 (MQ=255)                                                                                    gTCACCCTGGGTATCACATCTTATTTGTGATCAGGACATAAGCAACTGAAATTGATGGCGTGCATGACGAGGGAGGCAGATGCTGATTGCATTACCCCCGGTGATTACTAAAGGAGAGTCTAAAACGACTTTATTcccc                       <  2:60034/139‑1 (MQ=255)                                                                                          cTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTAt                 <  1:1399307/139‑1 (MQ=255)                                                                                          cTGGGTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTAt                 <  1:1275565/139‑1 (MQ=255)                                                                                            gggTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAAtt                                                                                                             >  1:1210500/1‑45 (MQ=255)                                                                                            gggTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAAtt                                                                                                             <  2:1210500/45‑1 (MQ=255)                                                                                             ggTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGga                                               <  1:1255105/106‑1 (MQ=255)                                                                                             ggTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGga                                               >  2:1255105/1‑106 (MQ=255)                                                                                              gTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                             <  2:820587/59‑1 (MQ=255)                                                                                              gTATCACATCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACg                                                                                             >  1:820587/1‑59 (MQ=255)                                                                                                   acaTCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTATcca        >  2:673631/1‑139 (MQ=255)                                                                                                    caTCTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTATccac       <  1:524789/139‑1 (MQ=255)                                                                                                       cTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTAcccc                                                                <  1:1009138/79‑1 (MQ=255)                                                                                                       cTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTAcccc                                                                >  2:1009138/1‑79 (MQ=255)                                                                                                       cTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTcccc                       >  2:965851/1‑120 (MQ=255)                                                                                                       cTTATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGAGTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTcccc                       <  1:965851/120‑1 (MQ=255)                                                                                                         tATTTTTAATCAGGACATAAGCAACTGAAATTGAtgg                                                                                                        >  2:1129582/1‑37 (MQ=255)                                                                                                         tATTTTTAATCAGGACATAAGCAACTGAAATTGAtgg                                                                                                        <  1:1129582/37‑1 (MQ=255)                                                                                                         tATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTa                             <  2:1258795/112‑1 (MQ=255)                                                                                                         tATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTa                             >  1:1258795/1‑112 (MQ=255)                                                                                                         tATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTATCCACCAGTa  <  2:1145365/139‑1 (MQ=255)                                                                                                         tATTTTTAATCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTATCCACCAGTa  >  1:126628/1‑139 (MQ=255)                                                                                                                 |                                                                                                                                   CGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCACATCTTATTTTTAGTCAGGACATAAGCAACTGAAATTGATGGCTGGCATGACGAGGGATGCAGATGCTGATTTCATTACCCCCGGTGATTACTAAAGGAGAGGCTAAAACGACTTTATTCCCCTGGTATGTGTATCCACCAGTA  >  NC_000913/3658781‑3659023

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