Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	684,841	Δ1 bp	6.9%	coding (1406/2400 nt)	cadA‑II ←	cadmium translocating P‑type ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	684,841	0	T	.	6.9%	181.6 / 7.0	58	coding (1406/2400 nt)	cadA‑II	cadmium translocating P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base T (25/29);  new base . (2/2);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

CAGGGCGCCGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGCCAGCCAC  >  NC_002947/684693‑684987                                                                                                                                                     |                                                                                                                                                    cAGGGCGCCGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCt                                                                                                                                                     >  1:93591/1‑149 (MQ=255) cAGGGCGCCGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCt                                                                                                                                                     >  3:164745/1‑149 (MQ=255)   gggCGCCGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCTCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGTTGAGGGTGCCGGTTTTGTCGAACACTACGGGGTTGACGGCATGGGCACGCTCCAGGGCTTCGGCATCCTTg                                                                                                                                                   <  2:9656/149‑1 (MQ=255)    ggCGCCGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGa                                                                                                                                                  <  2:93591/149‑1 (MQ=255)       gcCGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑a                                                                                                                                              <  6:277914/149‑1 (MQ=255)         cGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑Agg                                                                                                                                            <  4:118830/149‑1 (MQ=255)            ccAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATg                                                                                                                                         <  3:220186/149‑1 (MQ=255)             cAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGc                                                                                                                                        >  5:177905/1‑149 (MQ=255)                  ggTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCC‑TGATCAAGGATGCCgaag                                                                                                                                   <  1:131782/149‑4 (MQ=255)                       aGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTg                                                                                                                                    <  3:185284/143‑1 (MQ=255)                       aGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTg                                                                                                                                    >  4:185284/1‑143 (MQ=255)                         gTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGgcgggcgg                                                                                                                            <  2:152861/149‑1 (MQ=255)                                   ttGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCgg                                                                                                                   >  7:173110/1‑148 (MQ=255)                                   ttGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCgg                                                                                                                   >  4:166195/1‑148 (MQ=255)                                        gTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCgg                                                                                                             >  6:264395/1‑149 (MQ=255)                                         tCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGc                                                                                                            >  5:64079/1‑149 (MQ=255)                                         tCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGc                                                                                                            >  7:150154/1‑149 (MQ=255)                                           cAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGgcgggcgg                                                                                                                            >  6:3192/1‑131 (MQ=255)                                           cAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGgcgggcgg                                                                                                                            <  5:3192/131‑1 (MQ=255)                                                ccTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATc                                                                                                     >  3:63301/1‑149 (MQ=255)                                                            ccaccTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTg                                                                                         >  5:197995/1‑149 (MQ=255)                                                             caccTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGc                                                                                                            <  6:184594/129‑1 (MQ=255)                                                             caccTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGc                                                                                                            >  5:184594/1‑129 (MQ=255)                                                                   cGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCTAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGc                                                                                  <  6:177905/149‑1 (MQ=255)                                                                               ggTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAgg                                                                      <  4:164745/149‑1 (MQ=255)                                                                                gTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAggg                                                                     <  6:197995/149‑1 (MQ=255)                                                                                gTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAggg                                                                     >  5:39496/1‑149 (MQ=255)                                                                                 tGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAgggc                                                                    >  1:262556/1‑148 (MQ=255)                                                                                  gAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAgggg                                                                    >  8:130525/1‑148 (MQ=255)                                                                                     ggTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAgg                                                                <  7:130525/149‑1 (MQ=255)                                                                                      gTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGc                                                               <  4:63301/149‑1 (MQ=255)                                                                                       tGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCg                                                              <  4:23671/149‑1 (MQ=255)                                                                                        gCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGa                                                             >  6:136270/1‑149 (MQ=255)                                                                                         ccGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGAt                                                            <  5:264395/149‑1 (MQ=255)                                                                                          cGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATg                                                           <  3:166195/149‑1 (MQ=255)                                                                                                tGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAg                                                      <  6:39496/148‑1 (MQ=255)                                                                                                    gAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGGCGAGGGCGCAGGGGCAGGCGATGACCAGCacgg                                                 <  8:173110/149‑1 (MQ=255)                                                                                                    gAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATCATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacgg                                                 >  5:224633/1‑149 (MQ=255)                                                                                                     aaCACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacgg                                                 <  1:120302/148‑1 (MQ=255)                                                                                                          tACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggcgacgg                                           >  4:239119/1‑149 (MQ=255)                                                                                                          tACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggcgacg                                            >  7:71144/1‑148 (MQ=255)                                                                                                            cgcgGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggcgacg                                            <  8:71144/146‑1 (MQ=255)                                                                                                             gcgGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGt                                        >  4:73599/1‑149 (MQ=255)                                                                                                               ggTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGAcc                                                        <  2:98045/131‑1 (MQ=255)                                                                                                               ggTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGAcc                                                        >  1:98045/1‑131 (MQ=255)                                                                                                                     ccGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAgc                                >  5:38240/1‑149 (MQ=255)                                                                                                                           gggCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGa                                                                               <  2:25809/96‑1 (MQ=255)                                                                                                                           gggCACGCTCCAGGGCTTCGGCATCCTTGAAC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGa                                                                               >  1:25809/1‑96 (MQ=255)                                                                                                                            ggCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGttt                         >  2:25349/1‑149 (MQ=255)                                                                                                                             gCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGTCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGTGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTc                        <  3:239119/149‑1 (MQ=255)                                                                                                                                 gCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGc                    <  6:64079/149‑1 (MQ=255)                                                                                                                                 gCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGc                    <  8:150154/149‑1 (MQ=255)                                                                                                                                     cAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTa                >  4:281766/1‑149 (MQ=255)                                                                                                                                      aGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGCTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTAc               >  5:162498/1‑149 (MQ=255)                                                                                                                                          cTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCcca           <  5:136270/149‑1 (MQ=255)                                                                                                                                           cacGGCATCC‑TGATCAAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCccag          >  7:257593/3‑149 (MQ=255)                                                                                                                                            acGGCATCC‑TGATCAAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGCGGCAGGCGATGACCAGCACGGCGACGGCGATGATCAGCGCGGTTTCCAGCGCTACAccaga         >  8:196393/2‑148 (MQ=255)                                                                                                                                             cGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCccagc         <  7:200776/148‑1 (MQ=255)                                                                                                                                             cGGCATCCTTGATC‑AGGATGCAGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGAAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCccagc         <  7:197225/148‑1 (MQ=255)                                                                                                                                               gCATCC‑TGATCAAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCccagccag      <  8:257593/149‑1 (MQ=255)                                                                                                                                                     ttGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGccagccac  >  7:63634/1‑147 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CAGGGCGCCGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGCCAGCCAC  >  NC_002947/684693‑684987

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