Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,536,133	G→C	23.1%	R150R (CGC→CGG)	ompR ←	response regulator in two‑component regulatory system with EnvZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,536,133	0	G	C	23.1%	103.9 / 5.9	64	R150R (CGC→CGG)	ompR	response regulator in two‑component regulatory system with EnvZ
Reads supporting (aligned to +/- strand):  ref base G (21/25);  new base C (11/4);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 7.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCCCAGACGGTCTGAATGTAACGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GA‑ACA‑TTTC‑GCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAGCGCCAATCTCCAGGCCTACGATACGGTCCACTTCT‑TCCCCTTTCGCCGTCACCATAATGATCGGCATCGGGTTGCTCTGAC  >  NC_000913/3535901‑3536368                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  gcccccagacGGTCTGAATGTAACGCGGATGCGCTGGATCTTCGTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GA‑ACA‑TTTC‑GCGC‑G‑Ta                                                                                                                                                                                                                                           <  2:593985/237‑1 (MQ=255) gaccccagacGGTCTGAATGTAACGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCTCG‑GA‑CCA‑TTTC‑CCCC‑G‑t                                                                                                                                                                                                                                            >  2:758164/1‑238 (MQ=255)    cccagacGGTCTGAATGTAACGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GA‑ACA‑TTTC‑GCGC‑G‑TAcc                                                                                                                                                                                                                                         <  2:195338/238‑1 (MQ=255)       agacGGTCTGAATGTAACGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGCCCAGTGCCTTCAGTACCGCAAACTCACCGCTGGGGAGCGGCATCGGCC‑CGTC‑TTTCCG‑GG‑AAA‑TTTC‑TCCC‑G‑TACCCAgg                                                                                                                                                                                                                                     >  2:636031/1‑239 (MQ=255)       agacGGTCTGAATGTAACGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCCCG‑GG‑AAA‑TTTT‑GCGC‑G‑TAccccgg                                                                                                                                                                                                                                     >  2:944147/1‑235 (MQ=255)       agacGGTCTGAATGTAACGCGGATGCACTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGCA‑TTTCCG‑GG‑ACA‑ATTC‑CGCC‑G‑TACCGAgg                                                                                                                                                                                                                                     >  2:1011166/1‑239 (MQ=255)        ggcGGTCTGAAGTTAACGCGGATGCGCTGGATCTTCTCCCCCCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GA‑ACA‑TTTC‑GCGC‑G‑TACCGAGGt                                                                                                                                                                                                                                    <  1:758164/237‑1 (MQ=255)        gccgggcctgagtgaaccgcgATGCCGTGGATCTTCTGCCCCCATGCGGGGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATAGTCACGACCACGGGCAAGGGTCAGCAGCTTACCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTTAGTACCGCAAACTCGCCGCTGGTGAGCGGCATCGGCT‑CGTC‑TGCGGG‑GC‑ACA‑TATC‑GCGC‑G‑TAACTAGGt                                                                                                                                                                                                                                    <  1:485027/219‑1 (MQ=255)        gccGGGCTGAATGGAACGGGGAGTGGCTGGGTCTTTTGCCCCCAGGCGGGGCAAGCGGGAAATCTGCACGTCGATGGAGGGTTCCA‑TGGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GA‑ACA‑TTTC‑GCGC‑G‑TACCGAGGt                                                                                                                                                                                                                                    <  1:593501/237‑1 (MQ=255)        gacGGTCTGAATGTAACGCGGATGCCCTGGATCTTCTGCCACCCTGCGGCGCAAACGCGAAATCTGCCCGGCGATTGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GA‑ACA‑TTTC‑GCGC‑G‑TACCGAGGt                                                                                                                                                                                                                                    <  2:364124/239‑1 (MQ=255)        cgcggccgtgatttaacggtggtGCGCCTGATCTTGCTCCCTCAAGAGGCGCTGACGCGAAACCTGCGCGGCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTAGCGCGGGAGAGCGGCTCGCGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCTCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GA‑ACA‑TTTC‑GCGC‑G‑TACCGAGGt                                                                                                                                                                                                                                    <  1:566821/217‑1 (MQ=255)           gggcgtaaTGTTAGGCGGATGCGCTCAGCCTTCTTCCCCCATTCGGCGCAGACGAGAAAGCGGCACGGCGATGGAGCGTTCCA‑TTGCGTAATATTCACGACCCCGGGCAAGGTTCATCAGATTACCGCGGGAGAGCGGCTCACGC‑GGCTGGCTGACCAGTGCCTTCAGGACCGCAAACCCCCCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GA‑ACA‑TTTC‑GCGC‑G‑TACCGAGGt                                                                                                                                                                                                                                    <  1:1251754/230‑1 (MQ=255)            gTCTGAATGTAACGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GA‑ACA‑TTTT‑CCGC‑G‑TACCGAGGG‑TAA‑g                                                                                                                                                                                                                              >  2:933875/1‑239 (MQ=255)            gTCTGAATGTAACGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGGC‑TTCTCG‑GG‑ACATTTTC‑ACGC‑G‑ACCCGAGTT‑Tag                                                                                                                                                                                                                                >  2:778218/1‑238 (MQ=255)             tCTGAATGTAACGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCCCG‑GA‑ACA‑TTTC‑CCGC‑G‑TACCCAGGT‑TAC‑Gt                                                                                                                                                                                                                             >  2:236974/1‑239 (MQ=255)              cTGAATGTAACGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GA‑CCA‑TTTC‑GCGC‑G‑TACCCGGGT‑TAA‑Gtt                                                                                                                                                                                                                            >  2:831107/1‑239 (MQ=255)              cTGAATGTAACGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGGC‑CTCTCG‑GA‑ACA‑ATTC‑CCCC‑G‑GACCCACGG‑TAA‑Gtt                                                                                                                                                                                                                            >  2:1023332/1‑239 (MQ=255)              cTGAATGTAACGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGGGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCCCG‑GA‑ACA‑TTTT‑GCCC‑G‑GAccacgggtaatgt                                                                                                                                                                                                                              >  2:316672/1‑228 (MQ=255)               tGAATGTAACGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGGC‑TTCCCG‑GG‑AAA‑TTTCTGCG‑‑G‑TACCGAGGT‑TAA‑Ggtt                                                                                                                                                                                                                           >  2:939117/1‑239 (MQ=255)                       aCGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGGCCAGTGCCTTCAGTCCCGCAAACTCCCCGCTGGTGAGCGGCATCGGCT‑CGTG‑TTCGCGCGA‑CCA‑TTTC‑CCCG‑G‑GCCCGAGGT‑GAA‑ATTTGagatca                                                                                                                                                                                                                    >  1:258792/1‑234 (MQ=255)                       aCGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCCAACTCCCCGCTGGTGAGCGGCATCGGCTCCGTC‑‑TCGCG‑GG‑ACA‑CTTT‑GCGC‑G‑ACGCCAGGTATAA‑GTTT‑AACTTa                                                                                                                                                                                                                    >  1:859428/1‑238 (MQ=255)                       aCGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GA‑ACA‑TTTC‑GCGC‑G‑GACCCAGGG‑TAA‑GATTGTACTTa                                                                                                                                                                                                                    >  2:492500/1‑238 (MQ=255)                       aCGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGGGAGCGGCATCGGCT‑CGGC‑TTCGGG‑GG‑ACA‑TTCA‑CCGG‑G‑TACCGCGGtcttttataaacaat                                                                                                                                                                                                                      >  2:167483/1‑224 (MQ=255)                           ggATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGTCCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GA‑ACA‑TTTC‑GCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACCGaa                                                                                                                                                                                                               <  2:1274838/239‑1 (MQ=255)                           ggATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGGCTTTCCGG‑GA‑CCA‑TTTC‑GCGC‑G‑‑ACCCGGGT‑TAA‑GGTTTAACTTTcccca                                                                                                                                                                                                               >  2:581061/1‑236 (MQ=255)                             aTGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GG‑ACA‑TTTC‑GCCC‑G‑TACCCAGGT‑TAA‑GTTTGAACTTACCCAAAg                                                                                                                                                                                                             >  2:104073/1‑239 (MQ=255)                             aTGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCCCG‑GC‑ACA‑TTTT‑GCGC‑G‑TACCCAGGG‑TAA‑Gtttttaatatccgaaac                                                                                                                                                                                                             >  2:1194858/1‑225 (MQ=255)                             aTGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGGC‑TTCCCG‑GA‑ACA‑TTTC‑GCGC‑G‑GACCGAGGT‑TAA‑GGTTTAACTTACAGAAAg                                                                                                                                                                                                             >  2:1210867/1‑239 (MQ=255)                                          tcttcCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGGC‑CTCCCC‑GA‑AAC‑TTTC‑CCCC‑G‑TACCAAGGG‑TAA‑GTTTGAACTTTCCGAAAGCAATTACCGCct                                                                                                                                                                                                 >  2:534467/1‑238 (MQ=255)                                             tcCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGGC‑TTTGCC‑GG‑ACA‑TTTC‑CCAC‑T‑GACCGAGGG‑TTA‑GTTTGAACtttaccacacgaattaccccccctt                                                                                                                                                                                              >  2:532483/1‑215 (MQ=255)                                                    tGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGGC‑TTCGCG‑GG‑ACC‑TTTC‑GCCCGG‑AACAGAGG‑‑TAA‑GTTTGGACTTAACGCAAGCAATTACCCCCTCTTccctggac                                                                                                                                                                                      >  2:516605/1‑233 (MQ=255)                                                              aCGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GA‑ACA‑TTTC‑CCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACCGAAAGCAATTACCCCCTCTTCCTCTGACGGTGCGCCTg                                                                                                                                                                            >  2:908669/1‑239 (MQ=255)                                                              aCGCAAAATTGGCGCGGCGATGGAGCGTTCCA‑TTGCGGAATATCCCCCACCACGGGCAAGGTTCATCAGCTTACCGCGGGAGAGCGGCCCACGC‑GGATGGCTGACCAGTGCCTTCAGTCCCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GA‑ACA‑TTTC‑GCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGTGCGCCt                                                                                                                                                                             <  2:292098/238‑1 (MQ=255)                                                               cgcgAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGGC‑TTCGCG‑GA‑ACA‑TTTT‑GCCG‑G‑TCCCGAGGT‑TAA‑CGTTGAACTTACCGAAAACAATTACCGCCCCTTCACGTGGCGCTGGGccgg                                                                                                                                                                            >  2:1074465/1‑236 (MQ=255)                                                                   aaaTCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GG‑ACA‑TTTC‑GCGC‑G‑TACCCAGGG‑TAAGGTTTGAA‑TTTCCGCAAGCAATTACCCCCCCTTTCCTTGACGGGTCGCCCGGCAgg                                                                                                                                                                       >  2:27278/1‑238 (MQ=255)                                                                                ttggtgcgtttCCATTTGCGGAATTTTCGCGCCCCCGGGCAGGGTTAAACATCTTACGGCGGG‑GAGCGCCCGCCGCGGGATGGCTGACCAGGGCTTCCGGACCCGAAACCTCCCCGCTGGTGGGCGTCACGGCCT‑CGTC‑TCCGCG‑AA‑ACT‑TTTC‑GCGC‑G‑TACCGAGGT‑AAA‑GTTTGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTga                                                                                                                                                           <  1:27278/230‑1 (MQ=255)                                                                                    ggtgtgtCCA‑TTGGGGAAGTTTCCCGGCCAGGGCGAAGGGTCATCCGCTTACCGGGGGGGAGCGGCCCCGGC‑GGGTGGCTGACCTGTGCCTTCAGTACCGCAAACTCCCCGCTGGGGGGGGGCATCGGCT‑CGTC‑TTCGCG‑GA‑ACA‑TTCC‑GCGC‑G‑TACCGGGGT‑TAG‑GTTTGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGacgacg                                                                                                                                                      <  1:1023332/233‑1 (MQ=255)                                                                                         ttCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTTCCG‑GG‑ACA‑TTTT‑CCGC‑G‑TACCCGGGG‑TAA‑GTTTGAACTTACCGGAAGCAATTACCGCCCCTTCCTGTGACGGTGCGCCTGGCGGTTCGTTTGCCTGCCCAcgccgc                                                                                                                                                  >  2:17536/1‑238 (MQ=255)                                                                                         ttCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GG‑ACA‑TTTC‑GCGC‑G‑TACCGAGGG‑TAA‑GTTCGAAATTCCCGAACGCAATTTCCGCCCCTTTCTGTGACGGGGCGCCTGTCAGGTTGTTCGCCTcacgacccgaa                                                                                                                                                  >  1:411193/1‑233 (MQ=255)                                                                                                tGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CTGC‑TTTCCG‑GACACA‑CTTC‑ACCC‑G‑T‑CCGCAGG‑TCA‑GGTTGAACTTACCGGACGCAAATACCGCCCCTCCCCCTAGCGGGTCGCCGCGCAGTTGGTTccccgaccgcccccgcccccac                                                                                                                                            >  1:246103/1‑220 (MQ=255)                                                                                                             cgcccAGGGGCAAGGTTCATCAGCTTATCGCGGGAGCGGGGCTCACGG‑GGATGGCTGACCAGTGCCTTCAGTCCCGCAAACCCACCGCTGGTGAGCGGAATCGGCT‑CGTC‑TTCGCG‑GA‑CCA‑TTTC‑GCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCcag                                                                                                                              <  2:877783/236‑1 (MQ=255)                                                                                                                           gTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GA‑ACA‑TTTT‑GCGC‑G‑TAcc                                                                                                                                                                                                                                         >  2:1110951/1‑120 (MQ=255)                                                                                                                           gTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GA‑ACA‑TTTC‑GCGC‑G‑TAcc                                                                                                                                                                                                                                         <  1:1110951/120‑1 (MQ=255)                                                                                                                             tcaacaGTTTATGGGGGGAGGGCGGCCCACGC‑GGATGGCTGACCAGTGCCTTCAGTCCCGCAACCCCACCGCGGGTGAGCGGCAGCGGCT‑CGTC‑TTCGCG‑GA‑ACA‑TTTC‑GCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTa                                                                                                             <  2:1047253/239‑1 (MQ=255)                                                                                                                               aacgggttACCGCGGGAGAGCGGCTCACGC‑GGTGGGCTGACCGGTGCCTTCAGTCCCGCAAACTCCCCGCGGGTGGGCGGCATCGGCT‑CTTC‑TTCGCG‑GA‑ACA‑TTTC‑GCGC‑G‑TACCGAGGT‑TAA‑GTTGGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTaaa                                                                                                           <  2:54438/233‑1 (MQ=255)                                                                                                                                             gaaatgtgtataggggccggtGGCTGACCAGTGCCTTCAGTCCCGAAAACTCCCCGCTGGTGAGCGGCATCGTCT‑CGTC‑TTCGCG‑GA‑ACA‑TTTC‑GCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAAt                                                                                               <  2:1295057/218‑1 (MQ=255)                                                                                                                                                          cgggggAGGGCTGACAAGTGCCTTCAGGCCCGCAAACTCCCCGCGGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GA‑ACA‑TTTC‑GCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAgcgc                                                                                 <  2:1226985/235‑1 (MQ=255)                                                                                                                                                                  gggcgcccAGTCTCTTCAGGCCCGAAACCTCCCCGCGGGGGGGCGTACTCGGCG‑CTTC‑TGCGTG‑AA‑ACA‑TTCC‑CGGC‑T‑TACCGGGGT‑TAA‑GTTTGAACTGACCGAAAGCAATTAGCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAGCGCCAAtctc                                                                          <  1:167962/232‑1 (MQ=255)                                                                                                                                                                                          cAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GA‑ACA‑TTTC‑GCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGGGCGCCTTGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGGT‑CACGCGGGGTAAACGGGTTTGGAATGTAGTCGTCAGCGCCAATCTCCAGGCCTACGAAACGGTCCActtct                                                 >  2:1095995/1‑239 (MQ=255)                                                                                                                                                                                          cAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GA‑ACA‑ATTC‑GCGC‑G‑GACCGCGGT‑TAA‑GTTTGAACTTACCGCAAGCAATTACCGCCCCCTCCTGTGACGGTGCGCCCGGCAGGTCGTTCGCCTGCGTGCCCCACACCGCCCCGCTACGCGCCCCCAGTTCCGCCGGGGTAAAACGGGTTTGGTA‑GTAGTCGGCAGCCCCCCTCTCCACGCCTCCGCAACggccgccctc                                                  >  1:948613/1‑230 (MQ=255)                                                                                                                                                                                               tCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GA‑ACA‑TTTC‑GCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGGTTTGGAATGTAGTCGTCAGGGCCAATCTCCAGGCCTCCCATACGGTCCCcttct‑ttc                                             >  2:878889/1‑235 (MQ=255)                                                                                                                                                                                                 aCCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GA‑ACA‑TTTC‑GCGC‑G‑TACCGAGGC‑TAA‑GTTTGTACTTACCGAAAAACATTACTGCCCCTTCCTGGGGTGGGGCGCCTGGCCGTTCGTTCGCCATACGCCCCCGCCCCGCCCGGAGACGGGCCCGACGTT‑TACGCGGGGTCAACGGTTTTGGAGTGAGGGCGGCCAGGCGAATCCCCAGGGCTCCGAAACGGGCCACTTCTCTCCCCt                                          >  1:136064/1‑239 (MQ=255)                                                                                                                                                                                                      gggTGAGCGGATCCGTCT‑CGTC‑TTCGGG‑GA‑ACT‑TTGC‑CGGG‑T‑TACCGGGTT‑TAA‑GTTTGAACTTCCCGAAAGCAATTCCCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAGCGCCAATCTCCAGGCCTACGATACGGTCCACTTCT‑TCCCCTTTCGc                                     <  2:258792/237‑1 (MQ=255)                                                                                                                                                                                                         tGAGCGGCCTCGGCT‑GCTT‑CTCGGG‑GA‑ACA‑TTGC‑CGGC‑G‑TACCGAGGT‑TAA‑GTTTAAACTTACCGAAAGAAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAGCGCCAATCTCCAGGCCTACGATACGGTCCACTTCT‑TCCCCTTTCGCCGt                                  <  2:478793/238‑1 (MQ=255)                                                                                                                                                                                                         ggAGCGGCATGGGCT‑CGCT‑TCCGCG‑AA‑CCA‑TTTC‑GCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACAGAAAGCAATTCCCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAGCGCCAATCTCCAGGCCTACGATACGGTCCACTTCT‑TCCCCTTTCGCCGt                                  <  2:246103/237‑1 (MQ=255)                                                                                                                                                                                                                ccttagcgtctgtctggggg‑gA‑ACC‑TTTC‑GGGT‑C‑TACCGAGGG‑TAA‑GTTTGACCTGACCAAAAGCAATTTCCGCCCCTTCCTGGGAGGGGGCGCCGGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAGCGCCAATCTCCAGGCCTACGATACGGTCCACTTCT‑TCCCCTTTCGCCGTCACCa                             <  1:831107/219‑1 (MQ=255)                                                                                                                                                                                                                    ggCT‑GGTG‑TTCGGG‑GA‑ACA‑TTGG‑GCGC‑G‑ACCCGGGTT‑TAA‑TTTTGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAGCGCCAATCTCCAGGCCTACGATACGGTCCACTTCT‑TCCCCTTTCGCCGTCACCATAATGAt                      <  2:37404/239‑1 (MQ=255)                                                                                                                                                                                                                     ggctgttc‑tcCGGG‑GA‑ACA‑TTTC‑GCGA‑G‑TCCGGAGGT‑TAA‑GTTTGAACTTCCCGAAAGCAATTCCCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAGCGCCAATCTCCAGGCCTACGATACGGTCCACTTCT‑TCCCCTTTCGCCGTCACCATAATGAt                      <  2:7831/233‑1 (MQ=255)                                                                                                                                                                                                                           gttctCGGG‑GA‑ACA‑TTTC‑GCGG‑C‑GACAGGGGG‑TAA‑TTTTGAACTTACCGAAGGCAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAGCGCCAATCTCCAGGCCTACGATACGGTCCACTTCT‑TCCCCTTTCGCCGTCACCATAATGATCGGCa                 <  2:735274/235‑1 (MQ=255)                                                                                                                                                                                                                           gcctccgcg‑GA‑CCA‑TTCG‑CCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACCGAAGCCATTTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAGCGCCAATCTCCAGGCCTACGATACGGTCCACTTCT‑TCCCCTTTCGCCGTCACCATAATGATCGGCa                 <  2:594212/234‑1 (MQ=255)                                                                                                                                                                                                                            tgttCGCG‑GA‑AAC‑TTTT‑GCGC‑GTTTCCGGGGT‑TAA‑GTTTGAACTTCCCAAAATGAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAGCGCCAATCTCCAGGCCTACGATACGGTCCACTTCT‑TCCCCTTTCGCCGTCACCATAATGATCGGCa                 <  2:411193/237‑1 (MQ=255)                                                                                                                                                                                                                                        aCA‑TTCC‑GCGC‑G‑TACCGAGGT‑TAG‑TGTTGACCTTACCGAAAGCATTTCCCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAGCGCCAATCTCCAGGCCTACGATACGGTCCACTTCT‑TCCCCTTTCGCCGTCACCATAATGATCGGCATCGGGTTGct       <  2:66316/238‑1 (MQ=255)                                                                                                                                                                                                                                          a‑TTTC‑GCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAGCGCCAATCTCCAGGCCTACGATACGGTCCACTTCT‑TCCCCTTTCGCCGTCACCATAATGATCGGCATCGGGTTGCTCTg    >  1:403718/1‑239 (MQ=255)                                                                                                                                                                                                                                             ttC‑GCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTCAACGGTTGTGGACTGTAGTCGTCAACCCCACTCTCCAGGCCTACGATACGGTCCACTTCT‑TTCCCATTACCCGTACCCATAAAGATCGGCAACGGGTTGCTCTgca  >  1:1051400/1‑237 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  GACCCCAGACGGTCTGAATGTAACGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GA‑ACA‑TTTC‑GCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAGCGCCAATCTCCAGGCCTACGATACGGTCCACTTCT‑TCCCCTTTCGCCGTCACCATAATGATCGGCATCGGGTTGCTCTGAC  >  NC_000913/3535901‑3536368

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