Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,200,596	T→G	7.0%	Q152P (CAG→CCG)	yafJ ←	amidotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,200,596	0	T	G	7.0%	156.6 / 4.9	57	Q152P (CAG→CCG)	yafJ	amidotransferase
Reads supporting (aligned to +/- strand):  ref base T (30/23);  new base G (2/2);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

CCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGTTGTGCGCGA  >  NC_002947/5200448‑5200737                                                                                                                                                     |                                                                                                                                              ccGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTg                                                                                                                                              <  1:844024/150‑1 (MQ=255)                          cggGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGacg                                                                                                                     >  1:125983/1‑148 (MQ=255)                          cggGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTGCCACCGGCACCGGTTCGGGGAAg                                                                                                                     <  2:125983/150‑1 (MQ=255)                                   tgCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCgg                                                                                                            >  1:1087850/1‑150 (MQ=255)                                       ccAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTATCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGc                                                                                                        <  2:170170/150‑1 (MQ=255)                                       ccAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGc                                                                                                        >  2:429966/1‑150 (MQ=255)                                          gCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGt                                                                                                     <  1:461794/150‑1 (MQ=255)                                           cTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACAAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGtt                                                                                                    <  2:740229/150‑1 (MQ=255)                                            tttGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTg                                                                                                   <  1:429966/150‑1 (MQ=255)                                             ttGTCGAGCAGAAGCAGAACAGCGAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTGGa                                                                                                  >  1:249945/1‑150 (MQ=255)                                             ttGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCTGATGTGGTTGa                                                                                                  <  2:403769/150‑1 (MQ=255)                                               gTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGc                                                                                                >  2:1144827/1‑150 (MQ=255)                                               gTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGc                                                                                                <  1:101052/150‑1 (MQ=255)                                                         agcagaaCAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGa                                                                                      >  1:170790/1‑150 (MQ=255)                                                          gcagaaCAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGa                                                                                                  <  2:529137/137‑1 (MQ=255)                                                          gcagaaCAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGa                                                                                                  >  1:529137/1‑137 (MQ=255)                                                          gcagaaCAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGaa                                                                                     >  2:719552/1‑150 (MQ=255)                                                          gcagaaCAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGAAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGaa                                                                                     >  2:1142591/1‑150 (MQ=255)                                                           cagaaCAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAg                                                                                    <  2:1087850/150‑1 (MQ=255)                                                             gaaCAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGc                                                                                  >  2:949063/1‑150 (MQ=255)                                                               aCAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCg                                                                                >  2:853855/1‑150 (MQ=255)                                                                cAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGc                                                                               >  1:126216/1‑150 (MQ=255)                                                                 aGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCt                                                                              <  1:1144827/150‑1 (MQ=255)                                                                    cAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCg                                                                           >  2:453824/1‑150 (MQ=255)                                                                         gCCATCGCTGAGCATGCAGTTGAGCACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTc                                                                      <  1:607972/150‑1 (MQ=255)                                                                          ccATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCg                                                                     <  1:949063/150‑1 (MQ=255)                                                                                         cAGTCGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCCGCTGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGc                                                      <  2:508881/150‑1 (MQ=255)                                                                                           gTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCg                                                    >  1:577866/1‑150 (MQ=255)                                                                                           gTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCGGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCgg                                                                                                      >  1:337095/1‑100 (MQ=255)                                                                                           gTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCGGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCgg                                                                                                      >  1:336877/1‑100 (MQ=255)                                                                                           gTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCGGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCgg                                                                                                      <  2:337095/100‑1 (MQ=255)                                                                                           gTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCGGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCgg                                                                                                      <  2:336877/100‑1 (MQ=255)                                                                                                 cacCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCGACTGGCCGGTAGaa                                              >  2:696763/1‑150 (MQ=255)                                                                                                   cccccAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATTCGGTTGAGCAAGTCGCAGAAGGCGGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAAc                                            <  1:247471/150‑1 (MQ=255)                                                                                                     cccAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTg                                          <  2:126216/150‑1 (MQ=255)                                                                                                     cccAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTg                                          >  1:46956/1‑150 (MQ=255)                                                                                                      ccAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTgg                                         >  2:220227/1‑150 (MQ=255)                                                                                                      ccAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTgg                                         >  2:856883/1‑150 (MQ=255)                                                                                                            cccGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGc                                   <  1:754715/150‑1 (MQ=255)                                                                                                                   gCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGaa                            <  1:453824/150‑1 (MQ=255)                                                                                                                   gCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGaa                            <  2:46956/150‑1 (MQ=255)                                                                                                                   gCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGaa                            <  2:624265/150‑1 (MQ=255)                                                                                                                         acaGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCAcc                      >  2:747264/1‑150 (MQ=255)                                                                                                                           aGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCa                    >  1:851451/1‑150 (MQ=255)                                                                                                                             gCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCGGAAGGCGGCTTCGCTGTCGGTGTCGCCCACTTGGCGGTAGAAACTGGCCAGGCCCTTGAAGTCGCCCAGc                  >  2:796091/1‑150 (MQ=255)                                                                                                                             gCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCAGc                  >  2:931165/1‑150 (MQ=255)                                                                                                                                 cGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTTCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGc              <  1:719552/150‑1 (MQ=255)                                                                                                                                 cGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGc              >  1:26671/1‑150 (MQ=255)                                                                                                                                 cGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGc              >  2:668312/1‑150 (MQ=255)                                                                                                                                    ccAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGt           >  1:325792/1‑150 (MQ=255)                                                                                                                                     cAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGc                                        <  2:643660/120‑1 (MQ=255)                                                                                                                                     cAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGc                                        >  1:643660/1‑120 (MQ=255)                                                                                                                                     cAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGtt          >  2:715857/1‑150 (MQ=255)                                                                                                                                             gcagcaGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTTACCCAGCTGGCCGTTGTGCGCGa  >  1:1262904/1‑150 (MQ=255)                                                                                                                                             gcagcaGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGTTGTGCGCGa  >  2:864518/1‑150 (MQ=255)                                                                                                                                                gcagcTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCAc                       <  2:57970/126‑1 (MQ=255)                                                                                                                                                gcagcTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCAc                       >  1:57970/1‑126 (MQ=255)                                                                                                                                                     |                                                                                                                                              CCGTGCGGCACCGAATGGCGCACGGCGGGTAATGTGCACCAGCTTTGTCGAGCAGAAGCAGAACAGCCAGTCGCCATCGCTGAGCATGCAGTTGAACACCCCCAGATCCCGGTAGCCGGCACAGGCTTCGACCAGCACTGGCAGCAGCTGTTCCACCGGCACCGGTTCGGGGAAGGCGCTGCGGATGCGGTTGAGCAAGTCGCAGAAGGCCGCTTCGCTGTCGGTGTCGCCCACTGGCCGGTAGAAACTGGCCAGGCCCTTGAAGTCACCCAGCTGGCCGTTGTGCGCGA  >  NC_002947/5200448‑5200737

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