Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	372,649	T→A	7.0%	intergenic (+69/‑110)	dgcA → / → dgcB	dimethylglycine dehydrogenase subunit/dimethylglycine dehydrogenase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	372,649	0	T	A	7.0%	175.1 / 3.2	57	intergenic (+69/‑110)	dgcA/dgcB	dimethylglycine dehydrogenase subunit/dimethylglycine dehydrogenase subunit
Reads supporting (aligned to +/- strand):  ref base T (29/24);  new base A (2/2);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

CCGCATCGGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCC  >  NC_002947/372502‑372797                                                                                                                                                      |                                                                                                                                                     ccGCATCGGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGA‑TCAACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACt                                                                                                                                                      >  8:76201/1‑148 (MQ=255)     aTCGGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGacaa                                                                                                                                                 <  4:74662/149‑1 (MQ=255)         gCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACt                                                                                                                                             >  6:72097/1‑149 (MQ=255)           gACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATTCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTg                                                                                                                                            >  1:53231/1‑149 (MQ=255)            aCTGCGTGGCCCAGCGCAACGTGCATGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGtt                                                                                                                                          <  7:314314/149‑1 (MQ=255)                    gCCCAGCGCAACGTGCATGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGAcc                                                                                                                                  >  4:20998/1‑149 (MQ=255)                    gCCCAGCGCAACGTGCATGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGAc                                                                                                                                   <  7:76201/148‑1 (MQ=255)                    gCCCAGCGCAACGTGCATGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGAc                                                                                                                                   <  3:123836/148‑1 (MQ=255)                      ccAGCGCAACGTGCATGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCgg                                                                                                                                <  5:72097/149‑1 (MQ=255)                             aaCGTGCATGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCAtt                                                                                                                         >  3:299486/1‑149 (MQ=25)                                gTGCATGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCATTcgc                                                                                                                      <  4:299486/149‑1 (MQ=21)                                   cATGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCgg                                                                                                                    <  1:326093/148‑1 (MQ=255)                                    aTGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTcgc                                                                                                                      <  4:339315/145‑1 (MQ=255)                                    aTGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTcgc                                                                                                                      >  3:339315/1‑145 (MQ=255)                                    aTGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGt                                                                                                                  <  2:53231/149‑1 (MQ=255)                                        ggcgA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGt                                                                                                                  >  4:18382/1‑145 (MQ=255)                                        ggcgA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGt                                                                                                                  <  3:18382/145‑1 (MQ=255)                                            a‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGc                                                                                                          >  1:293771/1‑149 (MQ=255)                                                  cgacgCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTa                                                                                                     <  4:84422/149‑1 (MQ=255)                                                     cgCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACAt                                                                                                  <  3:20998/149‑1 (MQ=255)                                                        cTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCATTc                                                                                                                        >  7:264508/1‑124 (MQ=18)                                                        cTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCATTc                                                                                                                        <  8:264508/124‑1 (MQ=18)                                                               cTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGtc                                                                                        >  6:79031/1‑149 (MQ=255)                                                                  tgCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTAACCCCCAATTCCGAAACCAATCAAAACAATGTTGCATGACCGGACCAATACCCGGTTGAACCCGCTGTTACATGGGCCCataaac                                                                                      >  1:269908/1‑144 (MQ=255)                                                                    cAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGAc                                                                                                                              <  1:153639/106‑1 (MQ=255)                                                                    cAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGAc                                                                                                                              >  2:153639/1‑106 (MQ=255)                                                                    cAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGaa                                                                                                               >  3:168827/1‑121 (MQ=255)                                                                    cAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGaa                                                                                                               <  4:168827/121‑1 (MQ=255)                                                                           ttCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGt                                                                            >  5:339444/1‑149 (MQ=255)                                                                              tGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCgg                                                                         >  8:162435/1‑149 (MQ=255)                                                                              tGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCgg                                                                         >  6:165445/1‑149 (MQ=255)                                                                                    ttAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATggg                                                                                               >  2:305663/1‑121 (MQ=255)                                                                                      aaCAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGaa                                                                 >  2:206526/1‑149 (MQ=255)                                                                                           tCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAg                                                            >  5:314226/1‑149 (MQ=255)                                                                                                 ccATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGttt                                                      <  1:206526/149‑1 (MQ=255)                                                                                                    tCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATct                                                              >  2:253000/1‑137 (MQ=255)                                                                                                    tCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATct                                                              <  1:253000/138‑2 (MQ=255)                                                                                                    tCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTtgtg                                                   <  2:293771/149‑1 (MQ=255)                                                                                                       ggCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCaa                                                                    <  2:7863/129‑1 (MQ=255)                                                                                                       ggCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCaa                                                                    >  1:7863/1‑129 (MQ=255)                                                                                                              tCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATg                                                                                                 <  1:265497/93‑1 (MQ=255)                                                                                                              tCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATg                                                                                                 >  2:265497/1‑93 (MQ=255)                                                                                                                 cTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATg                                      <  7:162435/149‑1 (MQ=255)                                                                                                                    tAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTg                                   <  2:269908/149‑1 (MQ=255)                                                                                                                       gcgCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGACCCCGCACTTACATGGGGCCAGTCTCCCCCACTGGTCGGCGCAAGAAGCCAGCTGGTAGTCGGAGCCTCCCATGTAGAac                                >  3:346063/3‑149 (MQ=255)                                                                                                                       gAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGc                                                                                                          >  7:94749/1‑75 (MQ=255)                                                                                                                       gAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGc                                                                                                          <  8:94749/75‑1 (MQ=255)                                                                                                                       gAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAac                                >  3:53333/1‑149 (MQ=255)                                                                                                                         gCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAacac                              >  1:206488/1‑149 (MQ=255)                                                                                                                         gCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAac                                >  7:208354/1‑147 (MQ=255)                                                                                                                         gCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAac                                <  8:208354/147‑1 (MQ=255)                                                                                                                             gTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCtt                          >  8:254805/1‑149 (MQ=255)                                                                                                                                    ccGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTtgt                                                    <  6:150935/117‑1 (MQ=255)                                                                                                                                    ccGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTtgt                                                    >  5:150935/1‑117 (MQ=255)                                                                                                                                    ccGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTAcc                     >  8:10445/1‑148 (MQ=255)                                                                                                                                              gTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGt          >  8:289488/1‑149 (MQ=255)                                                                                                                                               tCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGtt         >  7:102913/1‑149 (MQ=255)                                                                                                                                                 aGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTc        <  6:314226/148‑1 (MQ=255)                                                                                                                                                      tGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGcc  >  8:252601/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CCGCATCGGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCC  >  NC_002947/372502‑372797

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