Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,699,027	G→A	R90H (CGC→CAC)	yfhH →	putative DNA‑binding transcriptional regulator YfhH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,699,027	0	G	A	100.0%	211.5 / NA	60	R90H (CGC→CAC)	yfhH	putative DNA‑binding transcriptional regulator YfhH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/29);  total (31/29)

CGGAGTCAGTCAGTCCAGCGTCGTGAAGTTCGCGCAAAAACTCGGCTATAAAGGTTTTCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCGCGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAAGAGAAACTGCATGAATGCGTAACAATGTTGCGCTCTGCGCGGCGGAT  >  NC_000913/2698890‑2699162                                                                                                                                          |                                                                                                                                        cGGAGTCAGTCAGTCCAGCGTCGTGAAGTTCGCGCAAAAACTCGGCTATAAAGGTTTTCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCcacg                                                                                                                                       <  1:390352/140‑1 (MQ=255)  ggAGTCAGTCAGTCCAGCGTCGTGAAGTTCGCGCAAAAACTCGGCTATAAAGGTTTTCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACgg                                                                                                                                      >  2:832866/1‑140 (MQ=255)     gtcagtcagtcCAGCGTCGTGAAGTTCGCGCAAAAACTCGGCTATAAAGGTTTTCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGtga                                                                                                                                   >  2:873014/1‑140 (MQ=255)       cagtcagtcCAGCGTCGTGAAGTTCGCGCAAAAACTCGGCTATAAAGGTTTTCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGtgatg                                                                                                                                 <  2:139013/140‑1 (MQ=255)        agtcagtcCAGCGTCGTGAAGTTCGCGCAAAAACTCGGCTATAAAGGTTTTCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGtgatga                                                                                                                                <  1:491500/140‑1 (MQ=255)        agtcagtcCAGCGTCGTGAAGTTCGCGCAAAAACTCGGCTATAAAGGTTTTCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGtgatga                                                                                                                                <  1:340373/140‑1 (MQ=255)              tcCAGCGTCGTGAAGTTCGCGCAAAAACTCGGCTATAAAGGTTTTCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGtgatg                                                                                                                                 <  2:526327/133‑1 (MQ=255)              tcCAGCGTCGTGAAGTTCGCGCAAAAACTCGGCTATAAAGGTTTTCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGtgatg                                                                                                                                 >  1:526327/1‑133 (MQ=255)                        tGAAGTTCGCGCAAAAACTCGGCTATAAAGGTTTTCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTc                                                                                                                >  2:542411/1‑140 (MQ=255)                        tGAAGTTCGCGCAAAAACTCGGCTATAAAGGTTTTCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTc                                                                                                                >  2:1231799/1‑140 (MQ=255)                        tGAAGTTCGCGCAAAAACTCGGCTATAAAGGTTTTCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTc                                                                                                                >  2:1061988/1‑140 (MQ=255)                         gAAGTTCGCGCAAAAACTCGGCTATAAAGGTTTTCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCg                                                                                                               <  1:657642/140‑1 (MQ=255)                             ttCGCGCAAAAACTCGGCTATAAAGGTTTTCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGa                                                                                                           >  2:655126/1‑140 (MQ=255)                                gcgcAAAAACTCGGCTATAAAGGTTTTCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGaaaa                                                                                                        >  2:984115/1‑140 (MQ=255)                                 cgcAAAAACTCGGCTATAAAGGTTTTCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGaaaaa                                                                                                       >  1:906648/1‑140 (MQ=255)                                 cgcAAAAACTCGGCTATAAAGGTTTTCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGaaaaa                                                                                                       <  1:1231799/140‑1 (MQ=255)                                       aaCTCGGCTATAAAGGTTTTCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGAtt                                                                                                 <  1:873014/140‑1 (MQ=255)                                              cTATAAAGGTTTTCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGaaa                                                                                          >  1:594641/1‑140 (MQ=255)                                                 taAAGGTTTTCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATa                                                                                       <  1:1092934/140‑1 (MQ=255)                                                    aGGTTTTCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTg                                                                                    <  1:202548/140‑1 (MQ=255)                                                    aGGTTTTCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTg                                                                                    <  1:542411/140‑1 (MQ=255)                                                      gTTTTCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGcc                                                                                  <  2:594641/140‑1 (MQ=255)                                                       ttttCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCg                                                                                 <  1:744836/140‑1 (MQ=255)                                                        tttCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGc                                                                                >  1:572538/1‑140 (MQ=255)                                                        tttCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGc                                                                                >  1:671667/1‑140 (MQ=255)                                                        tttCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGc                                                                                >  2:265320/1‑140 (MQ=255)                                                         ttCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCt                                                                               >  2:1205243/1‑140 (MQ=255)                                                          tCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTa                                                                              >  2:860367/1‑140 (MQ=255)                                                           ccGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTAt                                                                             >  1:1089638/1‑140 (MQ=255)                                                           ccGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTAt                                                                             <  1:1213614/140‑1 (MQ=255)                                                            cGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATg                                                                            <  2:572538/140‑1 (MQ=255)                                                                 cTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGc                                                                       <  1:655126/140‑1 (MQ=255)                                                                        tGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTa                                                                <  1:984115/140‑1 (MQ=255)                                                                                aGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTaa                                                        >  1:239979/1‑140 (MQ=255)                                                                                     aGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTg                                                   >  2:615037/1‑140 (MQ=255)                                                                                          tGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAagag                                              >  2:207764/1‑140 (MQ=255)                                                                                          tGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAagag                                              >  2:1020408/1‑140 (MQ=255)                                                                                          tGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAagag                                              >  1:1035373/1‑140 (MQ=255)                                                                                           ggCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAAgaga                                             <  1:1205243/140‑1 (MQ=255)                                                                                           ggCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAAgaga                                             <  1:832866/140‑1 (MQ=255)                                                                                            gCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAAGAGaa                                            >  1:1183780/1‑140 (MQ=255)                                                                                               agccagccGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAAgaga                                             >  2:267719/1‑136 (MQ=255)                                                                                               agccagccGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAAgaga                                             <  1:267719/136‑1 (MQ=255)                                                                                               agccagccGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAAGAGAAACt                                         <  1:523418/140‑1 (MQ=255)                                                                                                     ccGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAAGAGAAACTGCATGa                                   <  2:239979/140‑1 (MQ=255)                                                                                                     ccGGAATCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAAGAGAAACTGCATGa                                   <  1:811325/140‑1 (MQ=255)                                                                                                          aTCACCCTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAAGAGAAACTGCATGAATGCg                              <  1:1061988/140‑1 (MQ=255)                                                                                                               ccTCCGTGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAAGAGAAACTGCATGAATGCGTAACa                         <  2:157878/140‑1 (MQ=255)                                                                                                                     tGCCCATTCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAAGAGAAACTGCATGAATGCGTAACAATGTTg                   >  2:291861/1‑140 (MQ=255)                                                                                                                           ttCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAAGAGAAACTGCATGAATGCGTAACAATGTTGCGctct             >  1:659689/1‑140 (MQ=255)                                                                                                                            tCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAAGAGAAACTGCATGAATGCGTAACAATGTTGCGCTCTg            >  1:1172731/1‑140 (MQ=255)                                                                                                                            tCATAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGCTAATAGTGAAGAGAAACTGCATGAATGCGTAACAATGTTGCGCTTTg            >  1:180483/1‑140 (MQ=255)                                                                                                                              aTAACCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAAGAGAAACTGCATGAATGCGTAACAATGTTGCGCTCTgcg          <  1:953942/140‑1 (MQ=255)                                                                                                                                aaCCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAAGAGAAACTGCATGAATGCGTAACAATGTTGCGCTCTGcgcg        <  2:1172731/140‑1 (MQ=255)                                                                                                                                aaCCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAAGAGAAACTGCATGAATGCGTAACAATGTTGCGCTCTGcgcg        <  2:1089638/140‑1 (MQ=255)                                                                                                                                aaCCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAAGAGAAACTGCATGAATGCGTAACAATGTTGCGCTCTGcgcg        <  2:180483/140‑1 (MQ=255)                                                                                                                                aaCCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAAGAGAAACTGCATGAATGCGTAACAATGTTGCGCTCTGcgcg        <  1:291861/140‑1 (MQ=255)                                                                                                                                 aCCAAATCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAAGAGAAACTGCATGAATGCGTAACAATGTTGCGCTCTGCgcgg       >  1:1124465/1‑140 (MQ=255)                                                                                                                                    aaaTCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAAGAGAAACTGCATGAATGCGTAACAATGTTGCGCTCTGCGcggcgg    >  2:482223/1‑140 (MQ=255)                                                                                                                                     aaTCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAAGAGAAACTGCATGAATGCGTAACAATGTTGCGCTCTGCGCGGCGGa   <  2:671667/140‑1 (MQ=255)                                                                                                                                      aTCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAAGAGAAACTGCATGAATGCGTAACAATGTTGCGCTCTGCGCGGCGGAt  >  2:826313/1‑140 (MQ=255)                                                                                                                                      aTCCACGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAAGAGAAACTGCATGAATGCGTAACAATGTTGCGCTCTGCGCGGCGGAt  >  2:885065/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                        CGGAGTCAGTCAGTCCAGCGTCGTGAAGTTCGCGCAAAAACTCGGCTATAAAGGTTTTCCGGCGCTTAAGTTGGCGTTGAGTGAAGCGCTGGCAAGCCAGCCGGAATCACCCTCCGTGCCCATTCATAACCAAATCCGCGGTGATGATCCGTTACGGCTGGTCGGCGAAAAACTGATTAAAGAAAATACTGCCGCTATGTACGCAACGCTAAACGTTAATAGTGAAGAGAAACTGCATGAATGCGTAACAATGTTGCGCTCTGCGCGGCGGAT  >  NC_000913/2698890‑2699162

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