Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	31,025	A→C	7.6%	T199P (ACG→CCG)	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,025	0	A	C	7.6%	194.4 / 6.9	66	T199P (ACG→CCG)	PP_0026	cobalt/cadmium/zinc exporter
Reads supporting (aligned to +/- strand):  ref base A (32/29);  new base C (2/3);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

ACGTCAAAGGAGCCTACCTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCATGACTTGCACGTGTGGGCTGTCAGCAGTGG  >  NC_002947/30879‑31170                                                                                                                                                    |                                                                                                                                                   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                                                                                                                                        <  1:157349/147‑5 (MQ=255)             ccTACCTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCACGTCCGCGGCAgaacc                                                                                                                                        >  2:157349/1‑143 (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                                                                                                                   >  1:32592/1‑127 (MQ=255)                                                      gTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATcc                                                                                                                   <  2:32592/127‑1 (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)                                                                                                                                                 cgGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGCCGTGACGGACGTTCATGACTTGCACGTGTGGGCTGTCAGCAGTgg  <  1:232093/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   ACGTCAAAGGAGCCTACCTCGAAGTTTGGAGCGACATGCTCGGCTCGCTTGGCGTAATCATCGCGGCACTTATCATCCGCTTCACCGGCTGGACCTGGGTCG‑ACACGATTGTCGCTGTGGCGATCGGCCTTTGGGTTCTGCCGCGGACGTGGCAGTTGCTTCGCGAAAGCCTGGGGATCCTCATGGAGGGTGTCCCG‑AGGGGGCTCGACGTTACGGCAATCGAGGCCACCATCCTCGGAGTAGATGGCGTGACGGACGTTCATGACTTGCACGTGTGGGCTGTCAGCAGTGG  >  NC_002947/30879‑31170

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