Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	4,635,606	T→G	T66P (ACC→CCC)	uvrY ←	BarA/UvrY two‑componenent system response regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,635,606	0	T	G	100.0%	171.6 / NA	54	T66P (ACC→CCC)	uvrY	BarA/UvrY two‑componenent system response regulator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (27/27);  total (27/27)

CCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAA  >  NC_002947/4635462‑4635752                                                                                                                                                 |                                                                                                                                                   ccTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCt                                                                                                                                               >  2:915450/1‑150 (MQ=255) ccTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCt                                                                                                                                               >  1:683062/1‑150 (MQ=255)  cTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTc                                                                                                                                              >  1:347198/1‑150 (MQ=255)  cTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTc                                                                                                                                              <  1:158316/150‑1 (MQ=255)      aCCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTTACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAgg                                                                                                                                          <  2:1336322/150‑1 (MQ=255)      aCCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAgg                                                                                                                                          <  1:556155/150‑1 (MQ=255)       ccATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGgc                                                                                                                                         <  2:502453/150‑1 (MQ=255)          tttCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAACGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGAGCAGCAGTTTGCGGGGGACCTCCAGgccgc                                                                                                                                      >  1:580082/1‑150 (MQ=255)                   ggCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCgg                                                                                                                             <  2:22805/150‑1 (MQ=255)                    gCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTAAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCggg                                                                                                                            >  1:81640/1‑150 (MQ=255)                     ccTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGc                                                                                                                           <  2:158744/150‑1 (MQ=255)                      cTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCa                                                                                                                          >  1:40593/1‑150 (MQ=255)                          aCCCTTGGTCAGGTACCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCTTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCAtttt                                                                                                                      <  1:424683/150‑1 (MQ=255)                            ccTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTCACTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGACCGCCGATGCCGGGCATTTTca                                                                                                                    >  1:875480/1‑150 (MQ=255)                            ccTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTca                                                                                                                    <  2:798563/150‑1 (MQ=255)                             cTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTcac                                                                                                                   >  2:622009/1‑150 (MQ=255)                             cTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTcac                                                                                                                   >  2:715591/1‑150 (MQ=255)                             cTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTcac                                                                                                                   >  1:586290/1‑150 (MQ=255)                                     ggTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCa                                                                                                           <  1:622009/150‑1 (MQ=255)                                       tAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGc                                                                                                         >  1:98288/1‑150 (MQ=255)                                        aGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGca                                                                                                        <  1:243181/150‑1 (MQ=255)                                            ggCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAAc                                                                                                    <  1:1235006/150‑1 (MQ=255)                                               aGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCTTCCAGGCCGCCGATGCCGACCATTTTCAAATCCATCAGCACAACGTc                                                                                                 <  1:1319043/150‑1 (MQ=255)                                                   ccGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGCCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGt                                                                                             >  1:791068/1‑150 (MQ=255)                                                        cTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAg                                                                                        >  1:860197/1‑150 (MQ=255)                                                               aaGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATTCATCGACACAACGTCCGGTTTCAGCTCTCgg                                                                                 >  2:563147/1‑150 (MQ=255)                                                                         ggAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCCGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCa                                                                       >  2:207917/1‑150 (MQ=255)                                                                          gAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCACCGATGCCGACCACAATATCAGCCATCAACACAACGTCCGGTTTCAGCTCCCGGGCCCGTTTCAg                                                                      >  2:249984/1‑150 (MQ=255)                                                                             cGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGc                                                                   <  1:922615/150‑1 (MQ=255)                                                                                      cTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGcc                                                          >  1:59272/1‑150 (MQ=255)                                                                                      cTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGcc                                                          <  2:683062/150‑1 (MQ=255)                                                                                       tCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGTCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCg                                                         <  1:797582/150‑1 (MQ=255)                                                                                       tCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCa                                                                                                           >  1:1036768/1‑100 (MQ=255)                                                                                       tCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCa                                                                                                           <  2:1036768/100‑1 (MQ=255)                                                                                       tCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCg                                                         <  2:1184630/150‑1 (MQ=255)                                                                                           acacCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAAt                                                     >  1:481623/1‑150 (MQ=255)                                                                                             acCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAAtct                                                   >  1:362357/1‑150 (MQ=255)                                                                                               cGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTcc                                                 >  1:828403/1‑150 (MQ=255)                                                                                               cGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTcc                                                 <  1:119637/150‑1 (MQ=255)                                                                                                            aCCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCaccacc                                    <  1:775587/150‑1 (MQ=255)                                                                                                             ccTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCt                                   >  1:222730/1‑150 (MQ=255)                                                                                                              cTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCAACACCTg                                  >  1:538760/1‑150 (MQ=255)                                                                                                              cTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTAAGCATCCGTTCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTg                                  >  2:1045037/1‑150 (MQ=255)                                                                                                                 gATGTCAGGGTGGCTGCGCAGCAGTTTGCTGGGGGCCTCCAGGCCGCCGGTGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCAAGTTTAACGGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAg                               <  1:1175635/150‑1 (MQ=255)                                                                                                                 gATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAg                               <  2:59272/150‑1 (MQ=255)                                                                                                                    gTCAGGGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGACGGAATTTCCCTCACCCACCACCTGCAGGcc                            <  1:1045037/150‑1 (MQ=255)                                                                                                                     tCAGTGTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCg                           <  2:811761/150‑1 (MQ=255)                                                                                                                         ggTGGCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTCCTGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGa                       >  1:428200/1‑150 (MQ=255)                                                                                                                          gTGGCTGCGCAGCAGTTTGCGGGGCACGATCATGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGAt                      >  1:1043988/1‑150 (MQ=255)                                                                                                                           acaatGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATg                     >  1:264289/5‑150 (MQ=255)                                                                                                                            ggCTGCGCAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGt                    <  1:915450/150‑1 (MQ=255)                                                                                                                                gcgcAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACGACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGc                <  2:40593/150‑1 (MQ=255)                                                                                                                                 cgcAGCAGTTTGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGcc               <  1:11634/150‑1 (MQ=255)                                                                                                                                          ttGCGGGGGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCTGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCgt      <  2:860197/150‑1 (MQ=255)                                                                                                                                              ggggggCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCATGTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTaa  <  1:164934/150‑1 (MQ=255)                                                                                                                                              ggggggCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTaa  >  1:490539/1‑150 (MQ=255)                                                                                                                                              ggggggCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTaa  >  1:599558/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                                   CCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAA  >  NC_002947/4635462‑4635752

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