Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	31,021	G→A	7.9%	P197P (CCG→CCA)	PP_0026 →	cobalt/cadmium/zinc exporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	31,021	0	G	A	7.9%	176.4 / 4.3	63	P197P (CCG→CCA)	PP_0026	cobalt/cadmium/zinc exporter
Reads supporting (aligned to +/- strand):  ref base G (31/27);  new base A (2/3);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

CCTCAACGTCAAAGGAGCCTACCTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCATGACTTGCACGTGTGGGCTGT  >  NC_002947/30874‑31161                                                                                                                                                     |                                                                                                                                              ccTCAACGTCAAAGGAGCCTACCTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCcgc                                                                                                                                              >  3:213827/1‑149 (MQ=255)  cTCAACGTCAAAGGAGCCTACCTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCcgcg                                                                                                                                             <  6:22127/149‑1 (MQ=255)      aCGTCAAAGGAGCCTACCTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACg                                                                                                                                         >  5:96064/1‑149 (MQ=255)           aaaGGAGCCTACCTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCa                                                                                                                                    <  5:39641/149‑1 (MQ=255)              ggAGCCTACCTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCACGTCCGCGGCAgaa                                                                                                                                 <  3:56424/149‑3 (MQ=255)               gAGCCTACCTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCACGTCCGCGGCAgaac                                                                                                                                <  7:47106/149‑4 (MQ=255)                aGCCTACCTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGc                                                                                                                               <  8:150851/149‑1 (MQ=255)                  ccTACCTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCACGTCCGCGGCAgaacc                                                                                                                               >  2:157349/1‑143 (MQ=255)                  ccTACCTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCACGTCCGCGGCAgaacc                                                                                                                               <  1:157349/147‑5 (MQ=255)                   cTACCTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTc                                                                                                                            <  4:213827/149‑1 (MQ=255)                    tACCTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTcg                                                                                                                           >  1:160979/1‑149 (MQ=255)                       cTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTcgcg                                                                                                                         >  4:308703/1‑148 (MQ=255)                       cTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTcgcg                                                                                                                         >  3:155782/1‑148 (MQ=255)                        tCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGt                                                                                                                                  >  1:174425/1‑138 (MQ=255)                        tCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGt                                                                                                                                  <  2:174425/138‑1 (MQ=255)                             gTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCt                                                                                                                  >  1:213219/1‑149 (MQ=255)                              tttGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTg                                                                                                                 >  1:167660/1‑149 (MQ=255)                                tGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTgg                                                                                                                <  6:96064/148‑1 (MQ=255)                                          tGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATg                                                                                                     >  3:142931/1‑149 (MQ=255)                                              cggctcgCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAgg                                                                                                 >  5:267550/1‑149 (MQ=255)                                               ggctcgCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAggg                                                                                                >  4:136004/1‑149 (MQ=255)                                               ggctcgCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAg                                                                                                  <  8:83528/147‑1 (MQ=255)                                                     cTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTccc                                                                                          <  3:308703/149‑1 (MQ=255)                                                           gTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATcc                                                                                                          <  2:32592/127‑1 (MQ=255)                                                           gTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATcc                                                                                                          >  1:32592/1‑127 (MQ=255)                                                           gTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑Agggg                                                                                   >  2:123347/1‑149 (MQ=255)                                                               tcatcGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTc                                                                               >  8:58731/1‑149 (MQ=255)                                                               tcatcGCGGCACTTATCAACCGCTTCACCGGCTGGACCTGGGACG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTTGCATTTGCTTCGCGAAAGACTGGGGATCCTCATTGAGGGTGTCCCG‑AGGGGGCTc                                                                               >  5:290661/1‑149 (MQ=255)                                                                       gCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCTGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTAc                                                                       <  5:63920/149‑1 (MQ=255)                                                                        cACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACg                                                                      <  2:167660/149‑1 (MQ=255)                                                                         aCTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACgg                                                                     <  6:122972/149‑1 (MQ=255)                                                                           ttATCATCCGCTTCACCGGCTGGACCTGGGTCGTACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACgg                                                                     >  1:44373/1‑148 (MQ=255)                                                                             atcatcCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAAt                                                                 >  3:126393/1‑149 (MQ=255)                                                                             atcatcCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAAt                                                                 >  4:240577/1‑149 (MQ=255)                                                                              tcatcCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAAt                                                                 <  1:123347/148‑1 (MQ=255)                                                                                atcCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATc                                                                >  3:218995/1‑147 (MQ=255)                                                                                atcCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGa                                                              <  2:213219/149‑1 (MQ=255)                                                                                   cGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGc                                                           <  2:160979/149‑1 (MQ=255)                                                                                          ccGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑Aggg                                                                                    >  6:202580/1‑117 (MQ=255)                                                                                          ccGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑Aggg                                                                                    <  5:202580/117‑1 (MQ=255)                                                                                          ccGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATc                                                    >  2:320664/1‑149 (MQ=255)                                                                                             gCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTc                                                 >  1:277166/1‑149 (MQ=255)                                                                                              cTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCg                                                <  4:218995/149‑1 (MQ=255)                                                                                               tGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCAGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCg                                                >  6:73502/1‑148 (MQ=255)                                                                                                 gACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAg                                             <  4:126393/149‑1 (MQ=255)                                                                                                  aCCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGt                                            <  7:58731/149‑1 (MQ=255)                                                                                                   ccTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAg                                             >  2:232093/1‑147 (MQ=255)                                                                                                      gggTCG‑ACACGATTGTCGCTGTGGAGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGAt                                        >  1:324335/1‑149 (MQ=255)                                                                                                       ggTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCGTAGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGa                                         <  1:72304/148‑1 (MQ=255)                                                                                                             acacGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGt                                   >  1:146615/1‑148 (MQ=255)                                                                                                              cacGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGc                                     >  8:65592/1‑145 (MQ=255)                                                                                                              cacGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGc                                     <  7:65592/145‑1 (MQ=255)                                                                                                               acGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGc                                                           <  6:54627/122‑1 (MQ=255)                                                                                                               acGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGc                                                           >  5:54627/1‑122 (MQ=255)                                                                                                               acGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTgac                                <  3:136004/149‑1 (MQ=255)                                                                                                                 gATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTgacg                               >  1:129770/1‑148 (MQ=255)                                                                                                                         cTGTGGCGATCGGCCTTTGGGTTCTGCGGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCAt                      <  5:73502/149‑1 (MQ=255)                                                                                                                         cTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCAt                      >  7:298902/1‑149 (MQ=255)                                                                                                                            tgGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCATGAc                   <  6:267550/149‑1 (MQ=255)                                                                                                                              gCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCATGACtt                 <  3:240577/149‑1 (MQ=255)                                                                                                                              gCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCATGACtt                 >  2:242130/1‑149 (MQ=255)                                                                                                                                        ttttGGTTCTGCCGCGTAGGTGGCAGTTGCTTCGCGATAGCCTGGGGATCCTATTGGAGGGTGTCCCG‑TGGTGGCTCTACGTTAAGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCATGACTTGCACGTGTgg       <  6:290661/149‑1 (MQ=255)                                                                                                                                        tttGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCATGACTTGCACGTGTgg       >  2:113937/1‑149 (MQ=255)                                                                                                                                         gaagcaacTGCCACGTCCGCGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCATGACTTGCACGTGTggg      >  5:45373/8‑149 (MQ=255)                                                                                                                                          tGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCATGACTTGCACGTGTggg      <  2:277166/148‑1 (MQ=255)                                                                                                                                                    cgcgGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCATGACTTGCACGTGTGGGCTGt  >  6:87130/1‑142 (MQ=255)                                                                                                                                                    cgcgGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCATGACTTGCACGTGTGGGCTGt  <  5:87130/142‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                              CCTCAACGTCAAAGGAGCCTACCTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCATGACTTGCACGTGTGGGCTGT  >  NC_002947/30874‑31161

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