Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	685,191	G→A	7.7%	S352S (AGC→AGT) ‡	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	685,191	0	G	A	7.7%	190.4 / 4.3	65	S352S (AGC→AGT) ‡	cadA‑II	cadmium translocating P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base G (30/30);  new base A (2/3);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.89e-01

TGACACGGTCGACCAGCTTCTGGATCGGCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGG  >  NC_002947/685043‑685336                                                                                                                                                     |                                                                                                                                                    tGACACGGTCGACCAGCTTCTGGATCGGCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACg                                                                                                                                                     >  6:131879/1‑149 (MQ=255)     acGGTCGACCAGCTTCTGGATCGGCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGt                                                                                                                                                 >  6:169750/1‑149 (MQ=255)       ggTCGACCAGCTTCTGGATCGGCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCa                                                                                                                                               >  2:71060/1‑149 (MQ=255)          cGACCAGCTTCTGGATCGGCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑g                                                                                                                                           <  7:243599/149‑1 (MQ=255)           gACCAGCTTCTGGATCGGCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACAGCGTCACC‑gg                                                                                                                                          <  6:81742/149‑1 (MQ=255)           gACCAGCTTCTGGATCGGCGCGTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACAGCGTCACC‑gg                                                                                                                                          <  5:168179/149‑1 (MQ=255)             ccAGCTTCTGGATCGGCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACAGCGTCACC‑ggtg                                                                                                                                        >  1:189031/1‑147 (MQ=255)             ccAGCTTCTGGATCGGCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACAGCGTCACC‑ggtg                                                                                                                                        <  2:189031/149‑3 (MQ=255)                 cTTCTGGATCGGCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCt                                                                                                                                    >  8:98492/1‑149 (MQ=255)                 cTTCTGGATCGGCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCt                                                                                                                                    >  8:1651/1‑149 (MQ=255)                    cTGGATCGGCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTg                                                                                                                                  <  5:169750/148‑1 (MQ=255)                    cTGGATCGGCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTg                                                                                                                                  <  8:141415/148‑1 (MQ=255)                          cGGCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTg                                                                                                                           >  5:95701/1‑149 (MQ=255)                           ggCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTgg                                                                                                                          >  3:189010/1‑149 (MQ=255)                            gcgcTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTgg                                                                                                                          >  6:2370/1‑148 (MQ=255)                            gcgcTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTgg                                                                                                                          >  1:270896/1‑148 (MQ=255)                             cgcTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCa                                                                                                                        <  1:236365/149‑1 (MQ=255)                                 ttGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGct                                                                                                                    >  6:225993/1‑149 (MQ=255)                                               gTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCTCCGGTTACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCa                                                                                                      <  7:167131/149‑1 (MQ=255)                                               gTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCt                                                                                                                                                   <  5:81256/105‑1 (MQ=255)                                               gTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCt                                                                                                                                                   >  6:81256/1‑105 (MQ=255)                                                  ttCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑ggg                                                                                                  >  6:75175/1‑149 (MQ=255)                                                      accaGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTg                                                                                                                           >  2:127741/1‑121 (MQ=255)                                                      accaGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTg                                                                                                                           <  1:127741/121‑1 (MQ=255)                                                              gatgatGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCAt                                                                                      <  8:142736/149‑1 (MQ=255)                                                                  atgCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGgct                                                                                  <  5:131879/149‑1 (MQ=255)                                                                    gcgcgcCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGc                                                                                                 >  3:55735/1‑132 (MQ=255)                                                                    gcgcgcCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGc                                                                                                 <  4:55735/132‑1 (MQ=255)                                                                    gcgcgcCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGgctgc                                                                                >  8:248972/1‑149 (MQ=255)                                                                        gcCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGcc                                                                            <  3:107158/149‑1 (MQ=255)                                                                        gcCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGcc                                                                            <  6:259126/149‑1 (MQ=255)                                                                            gCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCt                                                                        >  7:78369/1‑149 (MQ=255)                                                                             cACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCtt                                                                       <  4:189010/149‑1 (MQ=255)                                                                               cGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCg                                                                     <  3:171664/149‑1 (MQ=255)                                                                                 gTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGAc                                                                   <  5:2370/149‑1 (MQ=255)                                                                                  tctcGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGAcc                                                                  >  8:245750/1‑149 (MQ=255)                                                                                  tctcGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGAc                                                                   >  3:131342/1‑148 (MQ=255)                                                                                           ccAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTTTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCATAGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTc                                                          <  1:188398/149‑1 (MQ=255)                                                                                               cgcCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCg                                                     <  6:95701/149‑1 (MQ=255)                                                                                                gcCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑gg                                                                                                   >  3:33973/1‑102 (MQ=255)                                                                                                gcCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑gg                                                                                                   <  4:33973/102‑1 (MQ=255)                                                                                                       tACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGt                                             >  3:73259/1‑149 (MQ=255)                                                                                                        aCGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGtt                                            >  7:44073/1‑149 (MQ=255)                                                                                                         cGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGtt                                            <  8:78369/148‑1 (MQ=255)                                                                                                          gCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCg                                          >  6:67744/1‑149 (MQ=255)                                                                                                               agcagcCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCggg                                     <  4:73259/149‑1 (MQ=255)                                                                                                               agcagcCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCggg                                     <  5:67744/149‑1 (MQ=255)                                                                                                                    ccGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACCGGGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGATTg                                 >  1:51950/1‑149 (MQ=255)                                                                                                                    ccGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGa                                >  4:141322/1‑149 (MQ=255)                                                                                                                       gCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCtt                                  <  7:163275/144‑1 (MQ=255)                                                                                                                       gCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCtt                                  >  8:163275/1‑144 (MQ=255)                                                                                                                             gCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTAc                         >  6:239676/1‑147 (MQ=255)                                                                                                                               cATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGc                                                                                        >  4:2913/1‑82 (MQ=255)                                                                                                                               cATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGc                                                                                        <  3:2913/82‑1 (MQ=255)                                                                                                                               cATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGc                                                                                        >  4:6767/1‑82 (MQ=255)                                                                                                                               cATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGc                                                                                        <  3:6767/82‑1 (MQ=255)                                                                                                                                aTTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCGTGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGt                    >  2:206712/1‑149 (MQ=255)                                                                                                                                aTTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCg                                                                     <  1:231849/100‑1 (MQ=255)                                                                                                                                aTTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCg                                                                     >  2:231849/1‑100 (MQ=255)                                                                                                                                    aTGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGcc                <  2:51950/149‑1 (MQ=255)                                                                                                                                       gcgcCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCa              <  7:1651/148‑1 (MQ=255)                                                                                                                                          ccACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCAcgcg          <  7:245750/149‑1 (MQ=255)                                                                                                                                          aacccGGTGACAGCGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCAcgcg          >  6:189571/4‑149 (MQ=255)                                                                                                                                                gTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCCCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGCGAAGCGTTCGCCCGCCTTGACCAGTACCAGGTCGCCCACGCCCAAGtc    >  8:125467/1‑148 (MQ=255)                                                                                                                                                 tGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTgg   <  4:131342/149‑1 (MQ=255)                                                                                                                                                   aCGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGGCGACTGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTggg  <  1:71060/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TGACACGGTCGACCAGCTTCTGGATCGGCGCTTTGGCGGCCTGGGCGTCTTCCACCAGGCGGATGATGCGCGCCAGCACGGTCTCGGTGCCCAGCGCCTGGGTACGCACCAGCAGCCGGCCTTCGCCATTGATGGCGCCACCGGTGACGCTGTCACC‑GGGTTGCTTGGGCACTGGCAGGCTCTCGCCGCTGATCA‑GGGCCTCATCGGCATGGCTGCTGCCGTCTTCGACCACGCCGTCGACCGGGAAGCGTTCGCCGGGCTTGACCAGTACCAGGTCGCCCACGCGCAAGTGGG  >  NC_002947/685043‑685336

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