Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,536,135	G→C	22.2%	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,135	0	G	C	22.2%	96.3 / 3.6	56	R150R (CGC→CGG)	ompR	response regulator in two‑component regulatory system with EnvZ
Reads supporting (aligned to +/- strand):  ref base G (20/21);  new base C (9/4);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCCCAGACGGTCTGAATGTAACGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GAACA‑TTTC‑GCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAGCGCCAATCTCCAGGCCTACGATACGGTCCACTTCT‑TCCCCTTTCGCCGTCACCATAATGATCGGCATCGGGTTGCTCTGAC  >  NC_000913/3535903‑3536370                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  gcccccagacGGTCTGAATGTAACGCGGATGCGCTGGATCTTCGTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GAACA‑TTTC‑GCGC‑G‑Ta                                                                                                                                                                                                                                           <  1:593985/237‑1 (MQ=255) gaccccagacGGTCTGAATGTAACGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCTCG‑GACCA‑TTTC‑CCCC‑G‑t                                                                                                                                                                                                                                            >  1:758164/1‑238 (MQ=255)    cccagacGGTCTGAATGTAACGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GAACA‑TTTC‑GCGC‑G‑TAcc                                                                                                                                                                                                                                         <  1:195338/238‑1 (MQ=255)       agacGGTCTGAATGTAACGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGCCCAGTGCCTTCAGTACCGCAAACTCACCGCTGGGGAGCGGCATCGGCC‑CGTC‑TTTCCG‑GGAAA‑TTTC‑TCCC‑G‑TACCCAgg                                                                                                                                                                                                                                     >  1:636031/1‑239 (MQ=255)       agacGGTCTGAATGTAACGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCCCG‑GGAAA‑TTTT‑GCGC‑G‑TAccccgg                                                                                                                                                                                                                                     >  1:944147/1‑235 (MQ=255)       agacGGTCTGAATGTAACGCGGATGCACTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGCA‑TTTCCG‑GGACA‑ATTC‑CGCC‑G‑TACCGAgg                                                                                                                                                                                                                                     >  1:1011166/1‑239 (MQ=255)        ggcGGTCTGAAGTTAACGCGGATGCGCTGGATCTTCTCCCCCCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GAACA‑TTTC‑GCGC‑G‑TACCGAGGt                                                                                                                                                                                                                                    <  2:758164/237‑1 (MQ=255)        gccGGGCTGAATGGAACGGGGAGTGGCTGGGTCTTTTGCCCCCAGGCGGGGCAAGCGGGAAATCTGCACGTCGATGGAGGGTTCCA‑TGGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GAACA‑TTTC‑GCGC‑G‑TACCGAGGt                                                                                                                                                                                                                                    <  2:593501/237‑1 (MQ=255)        gacGGTCTGAATGTAACGCGGATGCCCTGGATCTTCTGCCACCCTGCGGCGCAAACGCGAAATCTGCCCGGCGATTGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GAACA‑TTTC‑GCGC‑G‑TACCGAGGt                                                                                                                                                                                                                                    <  1:364124/239‑1 (MQ=255)           gggcgtaaTGTTAGGCGGATGCGCTCAGCCTTCTTCCCCCATTCGGCGCAGACGAGAAAGCGGCACGGCGATGGAGCGTTCCA‑TTGCGTAATATTCACGACCCCGGGCAAGGTTCATCAGATTACCGCGGGAGAGCGGCTCACGC‑GGCTGGCTGACCAGTGCCTTCAGGACCGCAAACCCCCCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GAACA‑TTTC‑GCGC‑G‑TACCGAGGt                                                                                                                                                                                                                                    <  2:1251754/230‑1 (MQ=255)            gTCTGAATGTAACGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GAACA‑TTTT‑CCGC‑G‑TACCGAGGG‑TAA‑g                                                                                                                                                                                                                              >  1:933875/1‑239 (MQ=255)            gTCTGAATGTAACGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGGC‑TTCTCG‑GGACATTTTC‑ACGC‑G‑ACCCGAGTT‑Tag                                                                                                                                                                                                                                >  1:778218/1‑238 (MQ=255)             tCTGAATGTAACGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCCCG‑GAACA‑TTTC‑CCGC‑G‑TACCCAGGT‑TAC‑Gt                                                                                                                                                                                                                             >  1:236974/1‑239 (MQ=255)              cTGAATGTAACGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GACCA‑TTTC‑GCGC‑G‑TACCCGGGT‑TAA‑Gtt                                                                                                                                                                                                                            >  1:831107/1‑239 (MQ=255)              cTGAATGTAACGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGGC‑CTCTCG‑GAACA‑ATTC‑CCCC‑G‑GACCCACGG‑TAA‑Gtt                                                                                                                                                                                                                            >  1:1023332/1‑239 (MQ=255)              cTGAATGTAACGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGGGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCCCG‑GAACA‑TTTT‑GCCC‑G‑GAccacgggtaatgt                                                                                                                                                                                                                              >  1:316672/1‑228 (MQ=255)               tGAATGTAACGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGGC‑TTCCCG‑GGAAA‑TTTCTGCG‑‑G‑TACCGAGGT‑TAA‑Ggtt                                                                                                                                                                                                                           >  1:939117/1‑239 (MQ=255)                       aCGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGGCCAGTGCCTTCAGTCCCGCAAACTCCCCGCTGGTGAGCGGCATCGGCT‑CGTG‑TTCGCGCGACCA‑TTTC‑CCCG‑G‑GCCCGAGGT‑GAA‑ATTTGagatca                                                                                                                                                                                                                    >  2:258792/1‑234 (MQ=255)                       aCGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCCAACTCCCCGCTGGTGAGCGGCATCGGCTCCGTC‑‑TCGCG‑GGACA‑CTTT‑GCGC‑G‑ACGCCAGGTATAA‑GTTT‑AACTTa                                                                                                                                                                                                                    >  2:859428/1‑238 (MQ=255)                       aCGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GAACA‑TTTC‑GCGC‑G‑GACCCAGGG‑TAA‑GATTGTACTTa                                                                                                                                                                                                                    >  1:492500/1‑238 (MQ=255)                           ggATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGTCCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GAACA‑TTTC‑GCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACCGaa                                                                                                                                                                                                               <  1:1274838/239‑1 (MQ=255)                           ggATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGGCTTTCCGG‑GACCA‑TTTC‑GCGC‑G‑‑ACCCGGGT‑TAA‑GGTTTAACTTTcccca                                                                                                                                                                                                               >  1:581061/1‑236 (MQ=255)                             aTGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GGACA‑TTTC‑GCCC‑G‑TACCCAGGT‑TAA‑GTTTGAACTTACCCAAAg                                                                                                                                                                                                             >  1:104073/1‑239 (MQ=255)                             aTGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGGC‑TTCCCG‑GAACA‑TTTC‑GCGC‑G‑GACCGAGGT‑TAA‑GGTTTAACTTACAGAAAg                                                                                                                                                                                                             >  1:1210867/1‑239 (MQ=255)                                          tcttcCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGGC‑CTCCCC‑GAAAC‑TTTC‑CCCC‑G‑TACCAAGGG‑TAA‑GTTTGAACTTTCCGAAAGCAATTACCGCct                                                                                                                                                                                                 >  1:534467/1‑238 (MQ=255)                                                    tGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGGC‑TTCGCG‑GGACC‑TTTC‑GCCCGG‑AACAGAGG‑‑TAA‑GTTTGGACTTAACGCAAGCAATTACCCCCTCTTccctggac                                                                                                                                                                                      >  1:516605/1‑233 (MQ=255)                                                              aCGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GAACA‑TTTC‑CCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACCGAAAGCAATTACCCCCTCTTCCTCTGACGGTGCGCCTg                                                                                                                                                                            >  1:908669/1‑239 (MQ=255)                                                              aCGCAAAATTGGCGCGGCGATGGAGCGTTCCA‑TTGCGGAATATCCCCCACCACGGGCAAGGTTCATCAGCTTACCGCGGGAGAGCGGCCCACGC‑GGATGGCTGACCAGTGCCTTCAGTCCCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GAACA‑TTTC‑GCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGTGCGCCt                                                                                                                                                                             <  1:292098/238‑1 (MQ=255)                                                               cgcgAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGGC‑TTCGCG‑GAACA‑TTTT‑GCCG‑G‑TCCCGAGGT‑TAA‑CGTTGAACTTACCGAAAACAATTACCGCCCCTTCACGTGGCGCTGGGccgg                                                                                                                                                                            >  1:1074465/1‑236 (MQ=255)                                                                   aaaTCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GGACA‑TTTC‑GCGC‑G‑TACCCAGGG‑TAAGGTTTGAA‑TTTCCGCAAGCAATTACCCCCCCTTTCCTTGACGGGTCGCCCGGCAgg                                                                                                                                                                       >  1:27278/1‑238 (MQ=255)                                                                                ttggtgcgtttCCATTTGCGGAATTTTCGCGCCCCCGGGCAGGGTTAAACATCTTACGGCGGG‑GAGCGCCCGCCGCGGGATGGCTGACCAGGGCTTCCGGACCCGAAACCTCCCCGCTGGTGGGCGTCACGGCCT‑CGTC‑TCCGCG‑AAACT‑TTTC‑GCGC‑G‑TACCGAGGT‑AAA‑GTTTGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTga                                                                                                                                                           <  2:27278/230‑1 (MQ=255)                                                                                    ggtgtgtCCA‑TTGGGGAAGTTTCCCGGCCAGGGCGAAGGGTCATCCGCTTACCGGGGGGGAGCGGCCCCGGC‑GGGTGGCTGACCTGTGCCTTCAGTACCGCAAACTCCCCGCTGGGGGGGGGCATCGGCT‑CGTC‑TTCGCG‑GAACA‑TTCC‑GCGC‑G‑TACCGGGGT‑TAG‑GTTTGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGacgacg                                                                                                                                                      <  2:1023332/233‑1 (MQ=255)                                                                                         ttCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTTCCG‑GGACA‑TTTT‑CCGC‑G‑TACCCGGGG‑TAA‑GTTTGAACTTACCGGAAGCAATTACCGCCCCTTCCTGTGACGGTGCGCCTGGCGGTTCGTTTGCCTGCCCAcgccgc                                                                                                                                                  >  1:17536/1‑238 (MQ=255)                                                                                         ttCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GGACA‑TTTC‑GCGC‑G‑TACCGAGGG‑TAA‑GTTCGAAATTCCCGAACGCAATTTCCGCCCCTTTCTGTGACGGGGCGCCTGTCAGGTTGTTCGCCTcacgacccgaa                                                                                                                                                  >  2:411193/1‑233 (MQ=255)                                                                                                             cgcccAGGGGCAAGGTTCATCAGCTTATCGCGGGAGCGGGGCTCACGG‑GGATGGCTGACCAGTGCCTTCAGTCCCGCAAACCCACCGCTGGTGAGCGGAATCGGCT‑CGTC‑TTCGCG‑GACCA‑TTTC‑GCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCcag                                                                                                                              <  1:877783/236‑1 (MQ=255)                                                                                                                           gTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GAACA‑TTTT‑GCGC‑G‑TAcc                                                                                                                                                                                                                                         >  1:1110951/1‑120 (MQ=255)                                                                                                                           gTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GAACA‑TTTC‑GCGC‑G‑TAcc                                                                                                                                                                                                                                         <  2:1110951/120‑1 (MQ=255)                                                                                                                             tcaacaGTTTATGGGGGGAGGGCGGCCCACGC‑GGATGGCTGACCAGTGCCTTCAGTCCCGCAACCCCACCGCGGGTGAGCGGCAGCGGCT‑CGTC‑TTCGCG‑GAACA‑TTTC‑GCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTa                                                                                                             <  1:1047253/239‑1 (MQ=255)                                                                                                                               aacgggttACCGCGGGAGAGCGGCTCACGC‑GGTGGGCTGACCGGTGCCTTCAGTCCCGCAAACTCCCCGCGGGTGGGCGGCATCGGCT‑CTTC‑TTCGCG‑GAACA‑TTTC‑GCGC‑G‑TACCGAGGT‑TAA‑GTTGGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTaaa                                                                                                           <  1:54438/233‑1 (MQ=255)                                                                                                                                                          cgggggAGGGCTGACAAGTGCCTTCAGGCCCGCAAACTCCCCGCGGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GAACA‑TTTC‑GCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAgcgc                                                                                 <  1:1226985/235‑1 (MQ=255)                                                                                                                                                                  gggcgcccAGTCTCTTCAGGCCCGAAACCTCCCCGCGGGGGGGCGTACTCGGCG‑CTTC‑TGCGTG‑AAACA‑TTCC‑CGGC‑T‑TACCGGGGT‑TAA‑GTTTGAACTGACCGAAAGCAATTAGCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAGCGCCAAtctc                                                                          <  2:167962/232‑1 (MQ=255)                                                                                                                                                                                          cAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GAACA‑TTTC‑GCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGGGCGCCTTGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGGT‑CACGCGGGGTAAACGGGTTTGGAATGTAGTCGTCAGCGCCAATCTCCAGGCCTACGAAACGGTCCActtct                                                 >  1:1095995/1‑239 (MQ=255)                                                                                                                                                                                          cAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GAACA‑ATTC‑GCGC‑G‑GACCGCGGT‑TAA‑GTTTGAACTTACCGCAAGCAATTACCGCCCCCTCCTGTGACGGTGCGCCCGGCAGGTCGTTCGCCTGCGTGCCCCACACCGCCCCGCTACGCGCCCCCAGTTCCGCCGGGGTAAAACGGGTTTGGTA‑GTAGTCGGCAGCCCCCCTCTCCACGCCTCCGCAACggccgccctc                                                  >  2:948613/1‑230 (MQ=255)                                                                                                                                                                                               tCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GAACA‑TTTC‑GCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGGTTTGGAATGTAGTCGTCAGGGCCAATCTCCAGGCCTCCCATACGGTCCCcttct‑ttc                                             >  1:878889/1‑235 (MQ=255)                                                                                                                                                                                                 aCCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GAACA‑TTTC‑GCGC‑G‑TACCGAGGC‑TAA‑GTTTGTACTTACCGAAAAACATTACTGCCCCTTCCTGGGGTGGGGCGCCTGGCCGTTCGTTCGCCATACGCCCCCGCCCCGCCCGGAGACGGGCCCGACGTT‑TACGCGGGGTCAACGGTTTTGGAGTGAGGGCGGCCAGGCGAATCCCCAGGGCTCCGAAACGGGCCACTTCTCTCCCCt                                          >  2:136064/1‑239 (MQ=255)                                                                                                                                                                                                      gggTGAGCGGATCCGTCT‑CGTC‑TTCGGG‑GAACT‑TTGC‑CGGG‑T‑TACCGGGTT‑TAA‑GTTTGAACTTCCCGAAAGCAATTCCCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAGCGCCAATCTCCAGGCCTACGATACGGTCCACTTCT‑TCCCCTTTCGc                                     <  1:258792/237‑1 (MQ=255)                                                                                                                                                                                                         tGAGCGGCCTCGGCT‑GCTT‑CTCGGG‑GAACA‑TTGC‑CGGC‑G‑TACCGAGGT‑TAA‑GTTTAAACTTACCGAAAGAAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAGCGCCAATCTCCAGGCCTACGATACGGTCCACTTCT‑TCCCCTTTCGCCGt                                  <  1:478793/238‑1 (MQ=255)                                                                                                                                                                                                         ggAGCGGCATGGGCT‑CGCT‑TCCGCG‑AACCA‑TTTC‑GCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACAGAAAGCAATTCCCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAGCGCCAATCTCCAGGCCTACGATACGGTCCACTTCT‑TCCCCTTTCGCCGt                                  <  1:246103/237‑1 (MQ=255)                                                                                                                                                                                                                    ggCT‑GGTG‑TTCGGG‑GAACA‑TTGG‑GCGC‑G‑ACCCGGGTT‑TAA‑TTTTGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAGCGCCAATCTCCAGGCCTACGATACGGTCCACTTCT‑TCCCCTTTCGCCGTCACCATAATGAt                      <  1:37404/239‑1 (MQ=255)                                                                                                                                                                                                                     ggctgttc‑tcCGGG‑GAACA‑TTTC‑GCGA‑G‑TCCGGAGGT‑TAA‑GTTTGAACTTCCCGAAAGCAATTCCCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAGCGCCAATCTCCAGGCCTACGATACGGTCCACTTCT‑TCCCCTTTCGCCGTCACCATAATGAt                      <  1:7831/233‑1 (MQ=255)                                                                                                                                                                                                                           gttctCGGG‑GAACA‑TTTC‑GCGG‑C‑GACAGGGGG‑TAA‑TTTTGAACTTACCGAAGGCAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAGCGCCAATCTCCAGGCCTACGATACGGTCCACTTCT‑TCCCCTTTCGCCGTCACCATAATGATCGGCa                 <  1:735274/235‑1 (MQ=255)                                                                                                                                                                                                                           gcctccgcg‑GACCA‑TTCG‑CCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACCGAAGCCATTTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAGCGCCAATCTCCAGGCCTACGATACGGTCCACTTCT‑TCCCCTTTCGCCGTCACCATAATGATCGGCa                 <  1:594212/234‑1 (MQ=255)                                                                                                                                                                                                                            tgttCGCG‑GAAAC‑TTTT‑GCGC‑GTTTCCGGGGT‑TAA‑GTTTGAACTTCCCAAAATGAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAGCGCCAATCTCCAGGCCTACGATACGGTCCACTTCT‑TCCCCTTTCGCCGTCACCATAATGATCGGCa                 <  1:411193/237‑1 (MQ=255)                                                                                                                                                                                                                                       aCA‑TTCC‑GCGC‑G‑TACCGAGGT‑TAG‑TGTTGACCTTACCGAAAGCATTTCCCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAGCGCCAATCTCCAGGCCTACGATACGGTCCACTTCT‑TCCCCTTTCGCCGTCACCATAATGATCGGCATCGGGTTGct       <  1:66316/238‑1 (MQ=255)                                                                                                                                                                                                                                         a‑TTTC‑GCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAGCGCCAATCTCCAGGCCTACGATACGGTCCACTTCT‑TCCCCTTTCGCCGTCACCATAATGATCGGCATCGGGTTGCTCTg    >  2:403718/1‑239 (MQ=255)                                                                                                                                                                                                                                            ttC‑GCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTCAACGGTTGTGGACTGTAGTCGTCAACCCCACTCTCCAGGCCTACGATACGGTCCACTTCT‑TTCCCATTACCCGTACCCATAAAGATCGGCAACGGGTTGCTCTgca  >  2:1051400/1‑237 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  GACCCCAGACGGTCTGAATGTAACGCGGATGCGCTGGATCTTCTTCCACCATGCGGCGCAGACGCGAAATCTGCACGTCGATGGAGCGTTCCA‑TTGCGGAATATTCACGACCACGGGCAAGGTTCATCAGCTTATCGCGGGAGAGCGGCTCACGC‑GGATGGCTGACCAGTGCCTTCAGTACCGCAAACTCACCGCTGGTGAGCGGCATCGGCT‑CGTC‑TTCGCG‑GAACA‑TTTC‑GCGC‑G‑TACCGAGGT‑TAA‑GTTTGAACTTACCGAAAGCAATTACCGCCTCTTCCTGTGACGGTGCGCCTGGCAGTTCGTTCGCCTGACGACGCAGCACCGCACGGATACGGGCCAGCAGTT‑CACGCGGGTTAAACGGTTTTGGAATGTAGTCGTCAGCGCCAATCTCCAGGCCTACGATACGGTCCACTTCT‑TCCCCTTTCGCCGTCACCATAATGATCGGCATCGGGTTGCTCTGAC  >  NC_000913/3535903‑3536370

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