Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_012971	767,863	A→C	A51A (GCA→GCC)	ybcW →	protein YbcW

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_012971	767,863	0	A	C	100.0%	206.6 / NA	56	A51A (GCA→GCC)	ybcW	protein YbcW
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (19/37);  total (19/37)

TGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCAGGCTT  >  NC_012971/767731‑768149                                                                                                                                     |                                                                                                                                                                                                                                                                                               tgatgccccctcccaggagcaggttcgtgtaaaaactatgcttaatagCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACTAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtat                                                                                                                         <  2:404036‑M1/277‑92 (MQ=255)   atgccccctcccaggagcaggttcgtgtaaaaactatgcttaatAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgcattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtattt                                                                                                                       <  1:501038‑M1/279‑94 (MQ=255)   atgccccctcccaggagcaggttcgtgtaaaaacgatgcttaatAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaaggcctcgtattt                                                                                                                       <  1:499836‑M1/279‑94 (MQ=255)           tcccatgagcaggttcgtgtaaaaactaTGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGTTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttagg                                                                                                                    <  2:148296‑M1/282‑97 (MQ=255)           tcccatgagcaggttcgtgtaaaaactaTGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttagg                                                                                                                    <  1:21447‑M1/282‑97 (MQ=255)           tcccaggagcaggttcgtgtaaaaactaTGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttagg                                                                                                                    <  1:18761‑M1/282‑97 (MQ=255)            cccaggagcaggttcgtgtaaaaactATGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaatt                                                                                                              >  2:460623‑M1/13‑198 (MQ=255)               aggagcaggttcgtgtaaaaACTATGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACTAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttagg                                                                                                                    <  1:205602‑M1/282‑97 (MQ=255)               aggagcaagttcgtgtaaaaACTATGCTTAATAGCACCACTTCGATGTGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGATGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttagg                                                                                                                    <  2:148005‑M1/282‑97 (MQ=255)                 gagcaggttcgtgtaaAAACTATGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcc                                                                                                         <  1:596828‑M1/293‑108 (MQ=255)                      ggttcgTTTAAAAACTATTCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCNAGTGTCTCTGGAAGCATTCAGTGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACTAAGGATTATTCCCTGGTGGTTGACTGATCACCAGAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggggaattgatcctacc                                                                                                     <  1:568971‑M1/297‑112 (MQ=255)                      ggttcgTTTAAAAACTATGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTTAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctaccc                                                                                                    <  2:437921‑M1/298‑113 (MQ=255)                                aaaacTATGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGNAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatg                                                                                          <  1:174434‑M1/300‑123 (MQ=255)                                     tatGCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacttaatatggacac                                                                                     <  2:209808‑M1/300‑128 (MQ=255)                                     tatGCTTAATAGCACCATCTCGATGGGTTACCCGGGTGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagttacagagtcaatatcgtattctgtcactggcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacac                                                                                     <  2:207919‑M1/300‑128 (MQ=255)                                     tatGCTGAGTAGCACCATTTCGATGGGTTACCCGGAGTTTGTAATTGCATGTATAAAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacac                                                                                     <  2:208052‑M1/300‑128 (MQ=255)                                        gCTTAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacagg                                                                                  <  1:279725‑M1/300‑131 (MQ=255)                                          ttAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcc                                                                                <  2:544866‑M1/300‑133 (MQ=255)                                          ttAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTAAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggc                                                                                 <  1:479984‑M1/299‑132 (MQ=255)                                           tAATAGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTATGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCATGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacataggc                                                                                 >  2:421832‑M1/1‑167 (MQ=255)                                               agcgccatttCGATGGGTTACCCGGATGGTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGAGAAGAATACGAAGGATTATTCCCTGGTGGGTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagc                                                                           <  2:645857‑M1/296‑138 (MQ=255)                                               aGCACCATTTTGATGTGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCGGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAACACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagc                                                                           <  2:148873‑M1/300‑138 (MQ=255)                                               aGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACTAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagc                                                                           <  2:383109‑M1/300‑138 (MQ=255)                                               aGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACTAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagc                                                                           <  1:624132‑M1/300‑138 (MQ=255)                                               aGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagc                                                                           <  2:298042‑M1/300‑138 (MQ=255)                                               aGCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAAAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtnaccatttcacaagtgacagagtcaatatcgtantctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagc                                                                           <  2:592078‑M1/300‑138 (MQ=255)                                               aGCACCATTTCGATGGGTTACCCGGAGTTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagc                                                                           <  2:756355‑M1/300‑138 (MQ=255)                                                gCACCATTTCGATGGGTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAANTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcg                                                                          <  1:171360‑M1/300‑139 (MQ=255)                                                             gggttACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgtttt                                                             >  1:494028‑M1/3‑149 (MQ=255)                                                               gTTACCCGGATGTTGTAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACNGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcataagcgaggttcttgttttca                                                           >  2:572251‑M1/1‑147 (MQ=255)                                                                         tgttgtAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcataagcgaggttcttgttttcaaattgttcc                                                  >  2:456343‑M1/1‑137 (MQ=255)                                                                         tgttgtAACTGCATGTATAGAACACCAAGTGTCTATGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaacatttcacaagtgacagagtcaatatcgtattctgtcactgccaggataatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttc                                                   >  2:9357‑M1/1‑137 (MQ=255)                                                                          gttgtAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttcc                                                  >  1:600610‑M1/1‑136 (MQ=255)                                                                           ttgtAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctagtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccgga                                               >  1:144966‑M1/1‑135 (MQ=255)                                                                           ttgtAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCngaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccgga                                               >  2:568542‑M1/1‑135 (MQ=255)                                                                           ttgtAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCATGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactaaagcaatgcctagtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccgga                                               >  2:145225‑M1/1‑135 (MQ=255)                                                                            tgtAATTGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTNtgaccacttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggac                                              <  2:259080‑M1/300‑167 (MQ=255)                                                                                 ttGCATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATANGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgagg                                         <  1:149703‑M1/300‑172 (MQ=255)                                                                                    cATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaattcaattcaactacagcaatgcctagtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggc                                        >  1:317393‑M1/1‑126 (MQ=255)                                                                                    cATGTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaattcaattcaactaaagcaatgcctcgtatttaggtgaattgatcctaaccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttcaggac                                              >  1:317398‑M1/1‑126 (MQ=255)                                                                                       gTATAGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctagtatttaggtgaattgatcctaaccacgtaatatggacacaggcctaagcgaggtttttgttttcaaattgttccggactgaggacgcca                                   >  2:177427‑M1/1‑123 (MQ=255)                                                                                       gTATAGAACATCAAGTGTCTATGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAAATGCTAAACATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatctta                                                                                                       >  2:34168‑M1/1‑123 (MQ=255)                                                                                        tataGAACATCAAGTGTCTTTGGAAGCATTCAGGGCAATTGAGGCCGCATTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccac                                  <  2:88922‑M1/300‑179 (MQ=255)                                                                                        tataGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccac                                  <  2:88292‑M1/300‑179 (MQ=255)                                                                                        tataGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccac                                  <  1:541273‑M1/300‑179 (MQ=255)                                                                                          taGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACTAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacac                                <  1:381305‑M1/300‑181 (MQ=255)                                                                                          taGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaattgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacac                                <  1:382937‑M1/300‑181 (MQ=255)                                                                                          taGAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacac                                <  1:381221‑M1/300‑181 (MQ=255)                                                                                            gAACATCAAGTGTCTCTGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacacca                              >  2:689352‑M1/1‑118 (MQ=255)                                                                                              acatcaagtgtctctgtaAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgncagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaac                            <  2:260792‑M1/293‑185 (MQ=255)                                                                                              acatcaagtgtctctggnAGCATTCAGGGCAATTGANGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactnaggccgccacaccaac                            <  1:678106‑M1/293‑185 (MQ=255)                                                                                                atcaggtGTCTCTGGAAGCAGTTAGGGCAATTTAGGCCGCGTTAGTAAAGCACGATAAGAATACTAAGGATTATTCCTTGGTGGTTGACTGATCACTATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactg                          <  2:521313‑M1/297‑187 (MQ=255)                                                                                                         ctctGGAAGCATTCAGGGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcaaattgttccggactgaggccgccacaccaactgtgccgccgc                 >  1:24949‑M1/1‑105 (MQ=255)                                                                                                                 gCATTCAGAGCAATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagcaatgcctagtatttaggtgaattgatcctacccacgtaatatggacacaggactaagcgagggtcttgttttcaagttgttcaggaatgaggccgccacaccaaatgtgccgccgcaaccgat          >  2:221390‑M1/1‑97 (MQ=255)                                                                                                                        gggcaATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaatttaactacagcaatgcatcgtatttaggtggattgatcctacccacgtaatatggacacaggcctacgcgaggttcttgttttcaaattgttccggactgaggccgcccaaccaactgtgccgacgccacagattgtaacc   >  2:75070‑M1/3‑90 (MQ=255)                                                                                                                        gggcaATTGAGGCCGCGTTAGTAAAGCACGATAAGAATACGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAACCATTCtgaccatttcacaagtgacagagtcaatatcgtattctgtcactgtcaggctaatacagcaatgcaattcaactacagccatgcctcgtatttaggtgaattgatcctacccacgtaatatggacacaggcctaagcgaggttcttgttttcacattgttccggactgaggccgacgcacccactgtgccgacgccaccgattgtaatca  >  2:150274‑M1/3‑90 (MQ=255)                                                                                                                                     |                                                                                                                                                                                                                                                                                               TGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGCAATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCAGGCTT  >  NC_012971/767731‑768149

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑