Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	76,358	C→A	G314G (GGG→GGT)	sgrR ←	DNA‑binding transcriptional dual regulator SgrR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	76,358	0	C	A	97.1%	108.6 / NA	35	G314G (GGG→GGT)	sgrR	DNA‑binding transcriptional dual regulator SgrR
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base A (22/11);  minor base G (0/1);  total (23/12)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAACAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCGCCCTGTGGACCTTTTAGCATCAGCCCTCCGGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGAATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAGTAACGGTAACCAGCGGTCCGGTTG  >  NC_000913/76114‑76600                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   cGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAGTAATAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCCCCCTGTGGACCTTTTAGCATCAGACCTcca                                                                                                                                                                                                                                              <  1:29876/251‑2 (MQ=255) cGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTcca                                                                                                                                                                                                                                              <  1:285608/251‑2 (MQ=255)  gCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTccag                                                                                                                                                                                                                                             <  1:284359/251‑3 (MQ=255)              gCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACTTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCtcgtcg                                                                                                                                                                                                                                  <  2:190445/250‑1 (MQ=255)                            aaCGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGc                                                                                                                                                                                                                      >  2:168278/1‑248 (MQ=255)                                    aGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCTCCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGAc                                                                                                                                                                                                           <  2:161583/251‑1 (MQ=255)                                        gTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGAAGtt                                                                                                                                                                                                        >  1:161066/1‑250 (MQ=255)                                             tAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGCTTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCGCACCCTGTGGACCTTTTAGCATCAGACCTCCCGCTGGCTCGTCGGCAATTTGCGGCAGAACCCAGACGTTAACTtc                                                                                                                                                                                                  >  2:91931/1‑251 (MQ=255)                                                  tgggAACCACAGTTGCTGGTACTGTTCCTCAGCGAACTAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTGGGCAATTTCCGGCAGAACCCAGCCGTTAACTtcgtcg                                                                                                                                                                                              >  2:93273/2‑250 (MQ=255)                                                                   ggTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCg                                                                                                                                                                            >  2:266105/1‑251 (MQ=255)                                                                                   gAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGTTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCCAATg                                                                                                                                                            >  2:86124/1‑251 (MQ=255)                                                                                       tAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTg                                                                                                                                                        >  1:48965/1‑251 (MQ=255)                                                                                                      ggAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGAtt                                                                                                                                         >  1:253515/1‑251 (MQ=255)                                                                                                               acaTAGCTTACCCAGACCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGtt                                                                                                                                >  1:268919/1‑251 (MQ=255)                                                                                                                                             ttCGCCCCGCGATGGGTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCgat                                                                                                  >  2:293626/1‑251 (MQ=255)                                                                                                                                               cGCCCCGCGATGGGTGAGGCTGTCGAATCGTAAATAGTAGCAACCTTCCTCCAGGCGGATTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCgatgg                                                                                                >  2:39901/1‑251 (MQ=255)                                                                                                                                                          gggTGCGGCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGAtggctggc                                                                                      <  1:205588/250‑1 (MQ=255)                                                                                                                                                                ggCTGTCGAATAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGtt                                                                                >  1:10392/1‑250 (MQ=255)                                                                                                                                                                          taGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGttt                                                                     >  2:22555/2‑251 (MQ=255)                                                                                                                                                                           aGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATg                                                                                                                                                            <  1:304695/163‑1 (MQ=255)                                                                                                                                                                           aGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATg                                                                                                                                                            >  2:304695/1‑163 (MQ=255)                                                                                                                                                                              aaaTAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTGATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGGTTCCCa                                                                 <  2:252920/251‑1 (MQ=255)                                                                                                                                                                                aTAGGAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCAtt                                                               >  2:90978/1‑251 (MQ=255)                                                                                                                                                                                  agtagCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTc                                                              <  2:221560/250‑1 (MQ=255)                                                                                                                                                                                                  cAGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCg                                             >  1:245332/1‑251 (MQ=255)                                                                                                                                                                                                   aGGCGGCTTTCAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGc                                            >  2:316577/1‑251 (MQ=255)                                                                                                                                                                                                             cAATCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGc                                  >  1:43401/1‑251 (MQ=255)                                                                                                                                                                                                               aTCTCTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGCCCTCCAGCTGGTTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGccc                                >  2:39463/1‑251 (MQ=255)                                                                                                                                                                                                                tctcTTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGcccc                               >  1:27453/1‑251 (MQ=255)                                                                                                                                                                                                                   ctTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCcagc                            >  1:147930/1‑251 (MQ=255)                                                                                                                                                                                                                   ctTTTTCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCcagc                            >  1:204672/1‑251 (MQ=255)                                                                                                                                                                                                                     ttttCCTCACCCTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCcagca                           <  1:352205/250‑1 (MQ=255)                                                                                                                                                                                                                         ccTCACCCTGTGGTACTCTTAGCATCAGGCCTCCAGCTGGCTCGTCGGCAATTTCCGGCCGAACCCAGATGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACAGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTCCTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCACCAGCAATAAc                      <  2:196393/251‑1 (MQ=255)                                                                                                                                                                                                                          ctcacccTGTGGACCTTTTAGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACg                     >  1:5650/5‑251 (MQ=255)                                                                                                                                                                                                                                             aGCATCAGACCTCCAGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGGATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAATAACGGTAACCAGCGGTCCGGTTg  <  2:147930/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   CGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTGACTTGCTGATTCGCCCCGCGATGGGTGCGGCTGTCGAACAGTAAATAGTAGCAACCTTCCTCCAGGCGGCTTTCAATCTCTTTTTCCTCGCCCTGTGGACCTTTTAGCATCAGCCCTCCGGCTGGCTCGTCGGCAATTTCCGGCAGAACCCAGACGTTAACTTCGTCGATTAATGCCCGGTAACCGAAGAAGTCATCGAATGCCTGAATTTTCAGTTGATTGGTGCTGTTGCGAATCACCGCATACGGACCGGTGCCGATGGGATGGCTGGCAAAGTTACTGAGGGTTTCCCATTCGCGCGGCAGGATCATCGCCGGAACTTGCCCCAGCAGTAACGGTAACCAGCGGTCCGGTTG  >  NC_000913/76114‑76600

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