Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,911,417	C→T	*745* (TAG→TAA)	relA ←	(p)ppGpp synthetase I/GTP pyrophosphokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,911,417	0	C	T	100.0%	248.4 / NA	75	*745* (TAG→TAA)	relA	(p)ppGpp synthetase I/GTP pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (48/27);  total (48/27)

GCCGGGATCCGCACCGCCGGTGAATTTCCCCAACGCTTTACGCTACTGTGGATCATAACCCTTTCCTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTA  >  NC_000913/2911284‑2911552                                                                                                                                      |                                                                                                                                        gCCGGGATCCGCACCGCCGGTGAATTTCCCCAACGCTTTACGCTACTGTGGATCATAACCCTTTCCTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTCAACt                                                                                                                                    >  2:528405/1‑139 (MQ=255)                 ccgGTGAATTTCCCCAACGCTTTACGCTACTGTGGATCATAACCCTTTCCTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGAcgc                                                                                                                    >  2:494772/1‑139 (MQ=255)                  cgGTGAATTTCCCCAACGCTTTACGCTACTGTGGATCATAGCCCTTTCCTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGAcgcg                                                                                                                   <  2:767659/139‑1 (MQ=255)                                             aCTGTGGATCATAACCCTTTCCTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGa                                                                                        >  2:624230/1‑139 (MQ=255)                                               tgtgGATCATAACCCTTTCCTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGt                                                                                      <  1:494772/139‑1 (MQ=255)                                                        tAACCCTTTCCTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGt                                                                                           >  2:653068/1‑125 (MQ=255)                                                        tAACCCTTTCCTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGt                                                                                           <  1:653068/125‑1 (MQ=255)                                                           cccTTTCCTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCAcg                                                                          <  2:503174/139‑1 (MQ=255)                                                             cTTTCCTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCAcgcg                                                                        >  1:785254/1‑139 (MQ=255)                                                                  cTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCa                                                                   <  2:249004/139‑1 (MQ=255)                                                                   tCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAg                                                                  >  2:608285/1‑139 (MQ=255)                                                                          gCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGc                                                           >  2:256461/1‑139 (MQ=255)                                                                                   tGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGACCCAGCACTTGCAGGTTGTaa                                                  >  2:264720/1‑139 (MQ=255)                                                                                        aTACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGt                                                                                           >  2:484212/1‑93 (MQ=255)                                                                                        aTACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGt                                                                                           <  1:484212/93‑1 (MQ=255)                                                                                        aTACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCa                                                                            <  2:259164/108‑1 (MQ=255)                                                                                        aTACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCa                                                                            >  1:259164/1‑108 (MQ=255)                                                                                        aTACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCAcg                                                                          <  1:291128/110‑1 (MQ=255)                                                                                        aTACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCAcg                                                                          >  2:291128/1‑110 (MQ=255)                                                                                        aTACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGcccc                                                                    <  1:289590/116‑1 (MQ=255)                                                                                        aTACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGcccc                                                                    >  2:289590/1‑116 (MQ=255)                                                                                        aTACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCa                                                          <  2:7530/126‑1 (MQ=255)                                                                                        aTACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCa                                                          >  1:7530/1‑126 (MQ=255)                                                                                        aTACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAAtctc                                             >  1:384994/1‑139 (MQ=255)                                                                                        aTACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAAtctc                                             >  1:117750/1‑139 (MQ=255)                                                                                        aTACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAAtctc                                             >  1:698605/1‑139 (MQ=255)                                                                                        aTACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAAtctc                                             >  2:16054/1‑139 (MQ=255)                                                                                        aTACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAAtctc                                             >  1:726224/1‑139 (MQ=255)                                                                                          aCAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGcccc                                                                    <  2:363501/114‑1 (MQ=255)                                                                                          aCAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGcccc                                                                    >  1:363501/1‑114 (MQ=255)                                                                                          aCAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAATCTCaa                                           >  1:399680/1‑139 (MQ=255)                                                                                              tatataTCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTg                                                     >  1:474535/1‑125 (MQ=255)                                                                                              tatataTCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTg                                                     <  2:474535/125‑1 (MQ=255)                                                                                              tatataTCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAATCTCAATGGt                                       >  2:230240/1‑139 (MQ=255)                                                                                                tatatCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAATCTCAATGGTCa                                     >  1:435339/1‑139 (MQ=255)                                                                                                 atatCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAATCTCAATGGTCAt                                    >  2:574264/1‑139 (MQ=255)                                                                                                   atCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGcccc                                                                    >  1:696986/1‑105 (MQ=255)                                                                                                   atCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGcccc                                                                    <  2:696986/105‑1 (MQ=255)                                                                                                   atCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAATCTCAATGGTCATGt                                  >  2:305425/1‑139 (MQ=255)                                                                                                   atCAATCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAATCTCAATGGTCATGt                                  >  1:375311/1‑139 (MQ=255)                                                                                                      aaTCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGcc                                                                      <  2:642058/100‑1 (MQ=255)                                                                                                      aaTCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGcc                                                                      >  1:642058/1‑100 (MQ=255)                                                                                                        tCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCaa                                                                                                                           <  1:354322/45‑1 (MQ=255)                                                                                                        tCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCaa                                                                                                                           >  2:354322/1‑45 (MQ=255)                                                                                                        tCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCa                                                                                                 <  2:558515/71‑1 (MQ=255)                                                                                                        tCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCa                                                                                                 >  1:558515/1‑71 (MQ=255)                                                                                                        tCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAg                                                                  >  2:38800/1‑102 (MQ=255)                                                                                                        tCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAg                                                                  <  1:269341/102‑1 (MQ=255)                                                                                                        tCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAg                                                                  <  1:38800/102‑1 (MQ=255)                                                                                                        tCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAg                                                                  <  2:753147/102‑1 (MQ=255)                                                                                                        tCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAg                                                                  >  1:753147/1‑102 (MQ=255)                                                                                                        tCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAg                                                                  >  2:269341/1‑102 (MQ=255)                                                                                                        tCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCa                                                          <  2:26066/110‑1 (MQ=255)                                                                                                        tCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCa                                                          >  1:26066/1‑110 (MQ=255)                                                                                                        tCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAATCTCAATGGTCATGTCGATg                             >  2:46996/1‑139 (MQ=255)                                                                                                        tCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAATCTCAATGGTCATGTCGATg                             >  1:702211/1‑139 (MQ=255)                                                                                                        tCTACATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAATCTCAATGGTCATGTCGATg                             >  2:732580/1‑139 (MQ=255)                                                                                                            cATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCAcg                                                                          <  2:35059/90‑1 (MQ=255)                                                                                                            cATTGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCAcg                                                                          >  1:35059/1‑90 (MQ=255)                                                                                                               tGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGc                                                                                                  >  1:793795/1‑63 (MQ=255)                                                                                                               tGTAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGc                                                                                                  <  2:793795/63‑1 (MQ=255)                                                                                                                 tAGATACGAGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGt                    >  2:585019/1‑139 (MQ=255)                                                                                                                      aCGAGCAAATTTCGGCTTAACTCCCGTGCAACCGAcg                                                                                                                     <  2:226148/37‑1 (MQ=255)                                                                                                                      aCGAGCAAATTTCGGCTTAACTCCCGTGCAACCGAcg                                                                                                                     >  1:226148/1‑37 (MQ=255)                                                                                                                         aGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTgg                                                                                            <  2:245050/59‑1 (MQ=255)                                                                                                                         aGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTgg                                                                                            >  1:245050/1‑59 (MQ=255)                                                                                                                         aGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGt                                                       <  2:546446/96‑1 (MQ=255)                                                                                                                         aGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGt                                                       >  1:546446/1‑96 (MQ=255)                                                                                                                         aGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGttt            >  2:547661/1‑139 (MQ=255)                                                                                                                         aGCAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGttt            >  1:732064/1‑139 (MQ=255)                                                                                                                           cAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTgg          >  2:618532/1‑139 (MQ=255)                                                                                                                           cAAATTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTgg          >  2:418754/1‑139 (MQ=255)                                                                                                                              aTTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCAcg                                                                          >  2:213642/1‑72 (MQ=255)                                                                                                                              aTTTCGGCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCAcg                                                                          <  1:213642/72‑1 (MQ=255)                                                                                                                                   ggCTTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTa  <  2:732064/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                        GCCGGGATCCGCACCGCCGGTGAATTTCCCCAACGCTTTACGCTACTGTGGATCATAACCCTTTCCTCAAACCGCTATCATATGTAGATACAGTATATATCAATCTACATTGTAGATACGAGCAAATTTCGGCCTAACTCCCGTGCAACCGACGCGCGTCGATAACATCCGGCACCTGGTTGAGTTTACCCAGCACGCGCCCCAGCACTTGCAGGTTGTAAATCTCAATGGTCATGTCGATGGTCGCCAGTTGCTGTTTGGTGTCGCTA  >  NC_000913/2911284‑2911552

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