Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	372,665	T→A	6.4%	intergenic (+85/‑94)	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,665	0	T	A	6.4%	192.6 / 3.1	62	intergenic (+85/‑94)	dgcA/dgcB	dimethylglycine dehydrogenase subunit/dimethylglycine dehydrogenase subunit
Reads supporting (aligned to +/- strand):  ref base T (33/25);  new base A (2/2);  total (35/27)
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.46e-01

GCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCTTGGCCTGGCAGT  >  NC_002947/372521‑372811                                                                                                                                                  |                                                                                                                                                   gCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGAcc                                                                                                                                                >  4:20998/1‑149 (MQ=255) gCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGAc                                                                                                                                                 <  3:123836/148‑1 (MQ=255) gCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGAc                                                                                                                                                 <  7:76201/148‑1 (MQ=255)   ccAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCgg                                                                                                                                              <  5:72097/149‑1 (MQ=255)          aaCGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCAtt                                                                                                                                       >  3:299486/1‑149 (MQ=25)             gTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCATTcgc                                                                                                                                    <  4:299486/149‑1 (MQ=21)                cATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCgg                                                                                                                                  <  1:326093/148‑1 (MQ=255)                 aTGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTcgc                                                                                                                                    >  3:339315/1‑145 (MQ=255)                 aTGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTcgc                                                                                                                                    <  4:339315/145‑1 (MQ=255)                 aTGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGt                                                                                                                                <  2:53231/149‑1 (MQ=255)                     ggcgATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGt                                                                                                                                <  3:18382/145‑1 (MQ=255)                     ggcgATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGt                                                                                                                                >  4:18382/1‑145 (MQ=255)                         aTCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑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                                                                                                                                      <  8:264508/124‑1 (MQ=18)                                    cTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCATTc                                                                                                                                      >  7:264508/1‑124 (MQ=18)                                           cTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGtc                                                                                                      >  6:79031/1‑149 (MQ=255)                                              tgCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTAACCCCCAATTCCGAAACCAATCAAAACAATGTTGCATGACCGGACCAATACCCGGTTGAACCCGCTGTTACATGGGCCCataaac                                                                                                    >  1:269908/1‑144 (MQ=255)                                                cAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGAc                                                                                                                                            >  2:153639/1‑106 (MQ=255)                                                cAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGAc                                                                                                                                            <  1:153639/106‑1 (MQ=255)                                                cAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGaa                                                                                                                             <  4:168827/121‑1 (MQ=255)                                                cAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGaa                                                                                                                             >  3:168827/1‑121 (MQ=255)                                                       ttCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGt                                                                                          >  5:339444/1‑149 (MQ=255)                                                          tGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCgg                                                                                       >  6:165445/1‑149 (MQ=255)                                                          tGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCgg                                                                                       >  8:162435/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                                                                                  >  1:7863/1‑129 (MQ=255)                                                                                   ggCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCaa                                                                                  <  2:7863/129‑1 (MQ=255)                                                                                          tCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATg                                                                                                               >  2:265497/1‑93 (MQ=255)                                                                                          tCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATg                                                                                                               <  1:265497/93‑1 (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                                                                                                                        <  8:94749/75‑1 (MQ=255)                                                                                                   gAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGc                                                                                                                        >  7:94749/1‑75 (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                                                                  >  5:150935/1‑117 (MQ=255)                                                                                                                ccGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTtgt                                                                  <  6:150935/117‑1 (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)                                                                                                                                    acaacaCTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATc                               >  6:215181/1‑132 (MQ=255)                                                                                                                                    acaacaCTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATc                               <  5:215181/132‑1 (MQ=255)                                                                                                                                     caacaCTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCtt             <  2:206488/149‑1 (MQ=255)                                                                                                                                         acTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCTTGGcc         >  4:104350/1‑149 (MQ=255)                                                                                                                                             ttGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGTCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCTTGGCCTgg      <  4:53333/148‑1 (MQ=255)                                                                                                                                                ccTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAGCACCCTTCTACCCATCCTGCTGTTCGCTGCCCTTGGCCTGGCAGt  >  6:121672/1‑149 (MQ=255)                                                                                                                                                ccTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCTTGGCCTGGCAGt  >  2:318143/1‑149 (MQ=255)                                                                                                                                                ccTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCTTGGCCTGGCAGt  >  6:220876/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCTTGGCCTGGCAGT  >  NC_002947/372521‑372811

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: ‑