Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	760,032	A→C	17.5%	T443P (ACC→CCC)	sucA →	2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	760,032	0	A	C	17.5%	178.5 / 2.1	74	T443P (ACC→CCC)	sucA	2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base A (46/12);  new base C (2/14);  total (48/26)
Fisher's exact test for biased strand distribution p-value = 1.68e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GGTTCAGGAAACCCTGAACATGTCGAAAGCGCGTGGTTATGAAGTTGGCGGTACGGTACGTATCGTTATCAACAACCAGGTTGGTTTCACCACCTCTAATCCGCTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCA  >  NC_000913/759800‑760268                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               ggTTCAGGAAACCCTGAACATGGCGAAAGCGCGTGGTTATGAAGTTGGCGGTACGGTACGTATCGTTATCAACAACCAGGTTGGTTTCACCACCTCTAATCCGCTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTa                                                                                                                                                                                                                                          <  1:107948/239‑1 (MQ=255)      aGGAAACCCTGAACATGTCGAAAGCGCGTGGTTATGAAGTTGGCGGTACGGTACGTATCGTTATCAACAACCAGGTTGGTTTCACCACCTCTAATCCGCTGGATGCCCGTTCTACGCCGTACTGTACTGCTATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TGCCCTTTGT‑CCCCGCCGTG‑CGCCCGATTTCCCTAAAACCTTTAAACgg                                                                                                                                                                                                                                     >  2:611772/1‑238 (MQ=255)          accccTGATCATTTCGAACGCGTGTTGTTTGGAATTTGGCGTTCTGGTCCGGATCGTTATCAACAACCAGGTTGGTTTCACCACCTCTAATCCGCTGGGTGCCCGTTCTACGCCGTACTGTACGGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCCCGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATg                                                                                                                                                                                                                                 <  2:644942/237‑1 (MQ=255)            cccTGAACATGTCGAAAGCGCGTGGTTATGAAGTTGGCGGTACGGTACGTATCGTTATCAACAACCAGGTTGGTTTCACCACCTCTAATCCGCTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATg                                                                                                                                                                                                                                 >  1:964514/1‑237 (MQ=255)             ccTGAACATGTCGAAAGCGCGTGGTTATGAAGTTGGCGGTACGGTACGTATCGTTATCACCAACCAGGTTGGTTTCACCACCTCTAATCCGCTGGATGCCAGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCAGGTTAACGC‑GGACGATCCGGAAGCCGT‑GG‑CCTTTGTGACCCGCATGG‑CGCTCGATTTCCGTAACCCCTTTAAACGTGATGGCt                                                                                                                                                                                                                              >  1:1227438/1‑239 (MQ=255)             ccTGAACATGTCGAAAGCGCGTGGTTATGAAGTTGGCGGTACGGTACGTATCGTTATCAACAACCAGGTTGGTTTCACCACCTCTAATCCGCTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGCCGATCCGGAGGCCGT‑TG‑CCTTTGTGACCCGTCTGC‑CGCTCGTCTTCCGTAACAACTTTCAACGTGCTGt                                                                                                                                                                                                                                >  1:462681/1‑237 (MQ=255)             ccTGAACATGTCGAAAGCGCGTGGTTATGAAGTTGGCGGTACGGTACGTATCGTTATCAACAACCAGGTTGGATTCACCACCTCTAATCCGCTGGATGCCCGTTCTACGCCGAACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTCTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGATCCGCCTGG‑CGCCTGATTTCCGTAAAACCTTTAACCGTGATGGCt                                                                                                                                                                                                                              >  1:48351/1‑239 (MQ=255)               tGAACATGTCGAAAGCGCGTGGTTATGAAGTTGGCGGTACGGTACGTATCGTTATCAACAACCAGGTTGGTTTCACCACCTCTAATCCGCTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTc                                                                                                                                                                                                                            >  1:185588/1‑239 (MQ=255)                gAACATGTCGAAAGCGGGTGGTTAAGAAGTTGGCGGTACGGTACGTATCGTTATCAACAACCAGGTTGGTTTCACCACCTCTAATCCGCTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCa                                                                                                                                                                                                                           <  1:1254955/239‑1 (MQ=255)                      gTCGAAGGGGCGTTGTTATGTGTTGGGGGGTGCGGTCCGATCCGGTCTCAACAACCATGGTGGTTGCACCACCTCTAATCCGCTGGAGCCCCGTTTTACGCCGGACTGTACTGATAGCGGTAAGATGGTTCAGGCCCCGATTTTCCCCGTTAACGC‑GGCCGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGAcc                                                                                                                                                                                                                     <  2:1172442/239‑1 (MQ=255)                      gTCGAAAGCGCGTGGTTATGAAGTTGGCGGTACGGTACGTATCGTTATCAACAACCAGGTTGGTTTCACCACCTCTAATCCGCTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGa                                                                                                                                                                                                                       >  1:1050367/1‑237 (MQ=255)                      gTCGAAAGCGCGTGGTTATGAAGTTGGCGGTACGGTACGTATCGTTATCAACAACCAGGTTGGTTTCACCACCTCTAATCCGCTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGAGGTCTTCATCGa                                                                                                                                                                                                                       >  1:39445/1‑237 (MQ=255)                      gTCGAAAGCGCGTGGTTATGAAGTTGGCGGTACGGTACGTATCGTTATCAACAACCAGGTTGGTTTCACCACCTCTAATCCGCTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCAGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCGGGAACACCTTTAAAAGTGATGTGTTCACCgcc                                                                                                                                                                                                                      >  2:1107273/1‑236 (MQ=255)                                cgTGGTTATGAAGTTGGCGGTACGGTACGTATCGTTATCAACAACCAGGTTGGTTTCACCACCTCTAATACGCTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCCGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGC‑TG‑CCTTTGTGACCCGTCTGG‑CGCCCGATTTCCGAAAAACCTTTAAACGGGATTGCTTCCTCTACCT‑GGTGT‑GCt                                                                                                                                                                                                          >  2:706743/1‑238 (MQ=255)                                               ggCGGTACGGTACGTATCGTTATCAACAACCAGGTTGGTTTCACCACCTCTAATCCGCTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGCAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTACCACCTTTAAACGGAATGACTTCATAGACCG‑GGTGG‑GCGAccaccgcacaggcc                                                                                                                                                                                           >  2:860857/1‑230 (MQ=255)                                               ggCGGTACGGTACGTATCGTTATCAACAACCAGGTTGGTTTCACCACCTCTAATCCGCTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCa                                                                                                                                                                                                                           >  2:278422/1‑208 (MQ=255)                                               ggCGGTACGGTACGTATCGTTATCAACAACCAGGTTGGTTTCACCACCTCTAATCCGCTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCa                                                                                                                                                                                                                           <  1:278422/208‑1 (MQ=255)                                               ggCGGTACGGTACGTATCGTTATCAACAACCAGGTTGGTTTCACCACCTCTAATCCGCTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCca                                                                                                                                                                                          >  1:588231/1‑239 (MQ=255)                                               ggCGGTACGGTACGTATCGTTATCAACAACCAGGTTGGTTTCACCACCTCTAATCCGCTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCca                                                                                                                                                                                          >  1:620793/1‑239 (MQ=255)                                               ggCGGTACGGTACGTATCGTTATCAACAACCAGGTTGGTTTCACCACCTCTAATCCGCTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCCCGATTTCCGTAACACCTTTAAACGCGATGGCTTCATCCCCCTCGGTG‑‑GCCACCGCCGTCACCGccc                                                                                                                                                                                          >  2:580436/1‑238 (MQ=255)                                               ggCGGTACGGTACGTATCGTTATCAACAACCAGGTTGGTTTCACCACCTCTAATCCGCTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CCCTCGATTTCCGTAACACCTTTAAACGTGATTTCTTCATCACCCT‑TGTGG‑GCTACCGCCGTCCCGcccc                                                                                                                                                                                          >  1:1179044/1‑238 (MQ=255)                                                                cGTTATCAACAACCAGGTTGGTTTCACCACCTCTAATCCGCTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCCTCTGG‑CGCTCGATTTCCGTTACACCTTTAAACTTGATGCCTTCATCGACCT‑GGTGG‑GCTCCCCCCGTCCCGGCCACAACGAAACCACcgag                                                                                                                                                                          >  2:71729/1‑238 (MQ=255)                                                                            ccAGGTTGTTTTCACCCCCTCTAACCCGCTGGATGCCCGTTCTCGGCCGTACTGTACTGATACCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAAccc                                                                                                                                                             <  1:1048107/239‑1 (MQ=255)                                                                                  tggttTCACCTCCTTTAATCCGCTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGGTGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGCCGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGc                                                                                                                                                       <  1:1276766/239‑1 (MQ=255)                                                                                                        cTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTCAACCCTTGAAACGTGCTGTCCTCATCGACCT‑GTTGT‑GCTACCGCCGGCACGTCCACAACGAAGCCGCCCAGCCGAGCGCAACCCAGCCGCCAGTGTATCAGAAAAGaaaaaa                                                                                                                                 >  2:901165/1‑239 (MQ=255)                                                                                                        cTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCCCAACGAAGCCGCCGAGCCGAGCGCAACCCACCCGCTGATGCATCAGAAAATCaaaaa                                                                                                                                 >  1:710621/1‑239 (MQ=255)                                                                                                        cTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCaaaaa                                                                                                                                 >  1:243504/1‑239 (MQ=255)                                                                                                        cTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGG‑GCTACCGCCGTCACGGGCCCAACGAAGCCGACGAGCCGGACGCAACCCACCCGCTGATGTATCCGAAAATCaaaaa                                                                                                                                 >  1:925828/1‑239 (MQ=255)                                                                                                                     cTACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACAcc                                                                                                                     >  1:1034062/1‑238 (MQ=255)                                                                                                                             tactgtactgATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTa                                                                                                                                                                                                                                          >  2:401417/1‑115 (MQ=255)                                                                                                                             tactgtactgATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTa                                                                                                                                                                                                                                          <  1:401417/115‑1 (MQ=255)                                                                                                                             tactgtactgATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACCCCGCGCaaaa                                                                                                             >  1:510286/1‑238 (MQ=255)                                                                                                                             tactgtactgATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCaaaa                                                                                                             >  1:856750/1‑238 (MQ=255)                                                                                                                                  tactgATCGCGGTAAGGTGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAACCCTT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACg                                                                                                       <  1:860857/239‑1 (MQ=255)                                                                                                                                  tactgATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCAAGCGCAACCAAGCCCCTGATGGATCAGAAAATCAACAAACATCCGACACCACGCAAAATCTag                                                                                                        >  1:970151/1‑237 (MQ=255)                                                                                                                                  tactgATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACAAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCAAc                                                                                                        >  1:542938/1‑238 (MQ=255)                                                                                                                                  tactgATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGCGT‑GCTACCGCCGTCCCGGCCACAACGAAGCCGCCGAGCCGAGCCGAACCCCACCGCGGTTGCATCCCAAAACCAAAAAACACCCCACCCCCCCCaaaacccac                                                                                                        >  2:978120/1‑233 (MQ=255)                                                                                                                                     tgATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGCCGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTg                                                                                                    >  1:761556/1‑239 (MQ=255)                                                                                                                                       atatCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGGTGGCTCAGAAAATCAAAAAACATCCGACACCACCCAAAAACTAAGCTAAc                                                                                                  >  2:993845/1‑239 (MQ=255)                                                                                                                                       atatCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGAc                                                                                                  >  1:113710/1‑239 (MQ=255)                                                                                                                                       atatCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCCGATGTATCCGAAAATCAAAAAACAACCGAAACCGCGCAAAAACtcagcggac                                                                                                  >  2:244200/1‑231 (MQ=255)                                                                                                                                       atatCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAACCGCGGATGTATCGGAAAATCAAAAAAAATCCAAAACCGCGAAAAAACCACGcaaaa                                                                                                  >  2:247654/1‑235 (MQ=255)                                                                                                                                       atatCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGA‑CGCTCGATTTCCGTAACACCTTTCAACGTGATGTCTTCATCGACCT‑GGTGG‑GCTACCGCCTTCACGGCCCAAACGAAACGCACGACCCACACGCCAACCAGCCGCGGATGGATCAGAAAATAAAAAAACAACCCACAACGCCAaaaaaccacatgcaa                                                                                                  >  2:755581/1‑228 (MQ=255)                                                                                                                                            ggTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGCGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACaa                                                                                                >  1:222320/1‑236 (MQ=255)                                                                                                                                            ggTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACaa                                                                                                >  1:370368/1‑236 (MQ=255)                                                                                                                                            ggTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCAACGCTgcaaa                                                                                                >  1:916395/1‑232 (MQ=255)                                                                                                                                            ggTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTAGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACCAGCCGAGCGCAACCCAGCCGCTGATGAATCAGAACATCAAACAACATCCGACCCCGCCCAAAATCTACGATAACaaccg                                                                                             >  1:1174143/1‑236 (MQ=255)                                                                                                                                            cggtaGATGTTTCTCGGCCCGATTTTCCCCGTTTACGG‑GGGCAGACCGGGAGCCGTGTC‑CCTTTGGACCCCTTCTGGTCGC‑GGTTTCCCGTAACCCCTGAAAAGGTGTTGCCTTCACCGACCG‑GGTGG‑CCTCCCGCCGTCAGGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGc                                                                                              <  2:856750/235‑1 (MQ=255)                                                                                                                                                   tGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGGGCAg                                                                                      >  1:204760/1‑239 (MQ=255)                                                                                                                                                              cccttttcCCCCGTTAGCGC‑GGCGGTCCCGGACCCCGT‑GG‑CCTTTGGGCCCCTTCTGC‑GCCGCGTTTTCCGAAACACTTTTAAACGTGATGTCTTCAGCGCCCT‑GGTGT‑GCTCCCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGc                                                                           <  2:964514/235‑1 (MQ=255)                                                                                                                                                              ccGTTTTTCC‑CCTTTACGCGGGACCATCCGGAAGCCGT‑TG‑CCTTTGGGACCCGTTTGG‑CGCCCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGc                                                                           <  1:483281/239‑1 (MQ=255)                                                                                                                                                                gttttcCCCCGTTACCGC‑GGCCGATCGGGAAGCCGT‑TG‑CCTTTGGGACCCGTCTGG‑CGCTCGATTTCCGTAACCCCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGa                                                                         <  1:244200/237‑1 (MQ=255)                                                                                                                                                                gATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCt                                                                                             >  1:617003/1‑219 (MQ=255)                                                                                                                                                                   tttCCCCGTTACCGG‑GGACCTTCCGGAAGCCGT‑TG‑CCTTTGTTACCCGTCTGG‑CGCCCGATTTCCGTAACCCCTTTAAACGTGATGTCTTCACCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGc                                                                      <  1:223106/239‑1 (MQ=255)                                                                                                                                                                               cgc‑gGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGGGGCGACGCTGGAAGATGca                                                           >  1:1045421/1‑237 (MQ=255)                                                                                                                                                                                cggtgACGAGCCGGAAGCCGT‑TC‑CCTTTGTGCCCCGTCTGG‑CGCTCGATTTCCGTAACCCCTTTAAACGTGAGGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCAAAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCa                                                          <  1:926974/235‑1 (MQ=255)                                                                                                                                                                                  gggACGATCCGGAAGCCTT‑TG‑CCTTTGTGGCCCGTCTGG‑CGCTCGATTGCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCa                                                          <  1:611772/236‑1 (MQ=255)                                                                                                                                                                                     ccGATCGGGAACCCGT‑TG‑CCTTTGGGACCCGTCTGG‑CGCTCGATTTCCGTAACCCCTTTAAACGTGATGTCTTCATCGACCT‑GGTGG‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCga                                                      <  1:580436/237‑1 (MQ=255)                                                                                                                                                                                     aggttCCGGAAGCCGT‑TG‑CCTTTGGGACCCGTTTGG‑CGCTCGATTTCCTTAACCCCTTTAAACGTGATGCCTTCATCGACCT‑GGTGG‑GGTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCga                                                      <  1:1198038/234‑1 (MQ=255)                                                                                                                                                                                     aCGATCCGGAAGCCGT‑TT‑CCTTTGTGCCCCGTCTGG‑CGCTCGATTTCCGTAACCCCTTTAAACGGGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCga                                                      <  1:993845/238‑1 (MQ=255)                                                                                                                                                                                      cGATCCGGAAGCCGT‑TG‑CCTTTGGCACCCGTCTGG‑CGCTCGATTTCCGTAACCCCTTTAACCGTTATGTCTTCATCGCCCT‑GGGGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCgag                                                     <  1:1107273/238‑1 (MQ=255)                                                                                                                                                                                         tCCGGAAGCCTG‑TG‑CCTTTTGGCCCCGTTCGG‑CGCTCGATTTCCGTAACCCCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATgg                                                 <  1:247654/239‑1 (MQ=255)                                                                                                                                                                                         tCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GTGGT‑GCTACCGCCGTCCCGGCCACAACGAAGCCACCGAGCCGAGCGCAACCCAGCCGCTCATGTATCAGAAAATCCAACACAATCCGCCACCCCCCAAAACCTACGCTGCCAACCTGGAGACGGACAAAGTGGCGAGCCGGGAAGAAGCAACCGAAATgg                                                 >  2:1002315/1‑239 (MQ=255)                                                                                                                                                                                         tCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCACAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATgg                                                 >  1:1204204/1‑239 (MQ=255)                                                                                                                                                                                            gcgcgccgtgtG‑CCTTTGTGCCCCGTGGGG‑ACCTATATTTCCTCAACCCCTTTAAACGGGATGTCTCGATCTACCT‑GCTGT‑GCTGCCGCCCGCACGGCCACAACGAAGCCGACGAGCCGTGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAAACTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGtt                                               <  2:1045421/229‑1 (MQ=255)                                                                                                                                                                                               aGCCGT‑TC‑CCTTGGTGCCCCGTCTGG‑CGCCCGATTGCCGAAACCCCTTTAAACGTGATGTCTTCACCGCCCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAAc                                            <  1:706743/238‑1 (MQ=255)                                                                                                                                                                                                    t‑tG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCAACCCTGGAAGATGCCACCGAGATGGTTAACCTGTa                                       >  1:479185/1‑238 (MQ=255)                                                                                                                                                                                                                  cccGTCGGG‑CCCTCGATTTCCGTAACCCCTTTAAACGTGGTGTTTTCGTCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCt                           <  1:937364/238‑1 (MQ=255)                                                                                                                                                                                                                      tCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGCCAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGAt                       >  1:661299/1‑238 (MQ=255)                                                                                                                                                                                                                      tCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGTTGGAt                       >  1:888357/1‑238 (MQ=255)                                                                                                                                                                                                                        tGG‑CGCTCGATTTCCGTACCCCTTTTAAACGTGATGTCTTCATCGACCG‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCt                    <  1:805472/239‑1 (MQ=255)                                                                                                                                                                                                                              cTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCCAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGAAGGCCTGGATGCTgggca               >  2:110705/1‑236 (MQ=255)                                                                                                                                                                                                                               ttgttttCCGTAACCCTTTTAACCTTGGTG‑CTTCAACGGCCC‑GGTGTGGCTCCCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGAt              <  2:113710/235‑1 (MQ=255)                                                                                                                                                                                                                                 gATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACg                                                                                                       >  1:64065/1‑148 (MQ=255)                                                                                                                                                                                                                                           acacCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGAGCCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGAGGCTGCCGATTGCGTAGTGGCa  >  1:687914/1‑239 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               GGTTCAGGAAACCCTGAACATGTCGAAAGCGCGTGGTTATGAAGTTGGCGGTACGGTACGTATCGTTATCAACAACCAGGTTGGTTTCACCACCTCTAATCCGCTGGATGCCCGTTCTACGCCGTACTGTACTGATATCGGTAAGATGGTTCAGGCCCCGATTTTCCACGTTAACGC‑GGACGATCCGGAAGCCGT‑TG‑CCTTTGTGACCCGTCTGG‑CGCTCGATTTCCGTAACACCTTTAAACGTGATGTCTTCATCGACCT‑GGTGT‑GCTACCGCCGTCACGGCCACAACGAAGCCGACGAGCCGAGCGCAACCCAGCCGCTGATGTATCAGAAAATCAAAAAACATCCGACACCGCGCAAAATCTACGCTGACAAGCTGGAGCAGGAAAAAGTGGCGACGCTGGAAGATGCCACCGAGATGGTTAACCTGTACCGCGATGCGCTGGATGCTGGCGATTGCGTAGTGGCA  >  NC_000913/759800‑760268

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