Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,636,925	G→C	R77P (CGC→CCC)	creC →	sensory histidine kinase CreC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,636,925	0	G	C	97.0%	219.3 / ‑6.5	66	R77P (CGC→CCC)	creC	sensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (30/34);  minor base A (2/0);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGC  >  NC_000913/4636679‑4637170                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       aGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGTTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACTCATGGGCAACTGGCGCAGGCGTTGAATCAGCTAC‑AACATCCCCCg                                                                                                                                                                                                                                                    <  1:408231/251‑1 (MQ=255)  gCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGGTTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCCGGCGTTTAATCAGCTAC‑AACATCCCCCGt                                                                                                                                                                                                                                                   >  2:45736/1‑251 (MQ=255)      gagtggCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCTGGCGTGCGAAGAGCAACGGATGGGACGTTGAGCGACACCGCAACGTTGCTGGCGGATCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTc                                                                                                                                                                                                                                                <  1:263059/246‑1 (MQ=255)        ggggCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCCACTGGCGCCGGCGTTTTATCCGCTAC‑AACATCCCCCGTTTtgcg                                                                                                                                                                                                                                             >  2:400282/1‑251 (MQ=255)               gTAATGCGTGTCGGCATGCGGTTGTTGCTGGGCTAGTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTGTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCTCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAAGTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAatat                                                                                                                                                                                                                                      <  1:226871/251‑1 (MQ=255)               gTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAat                                                                                                                                                                                                                                        <  1:68734/249‑1 (MQ=255)                 aaTGCGTATCGGCATGCGGCTGATGGAAGGCTATTTTTTACTGGTTGCGGTTGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTCAACCGGGCGAAAGAAGAGCAACGGAGGGGACGTTGTTCGACACCGCAACGTTGCTGGCGGAGCGTGCGCGTCCAGATTTGCTCTCTGGGGAACCAACACATTGGCAACTGGGGCAGGCGTTTGATAAGCTACCAAC‑TCACCCATTTCGCGCCCatattg                                                                                                                                                                                                                                    >  2:187599/1‑249 (MQ=255)                     cGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTgg                                                                                                                                                                                                                                <  1:165947/251‑1 (MQ=255)                     cGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGGTTCGCGCCCATATCggggg                                                                                                                                                                                                                                >  2:165947/1‑248 (MQ=255)                          cGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGTCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCAtt                                                                                                                                                                                                                            <  2:407915/250‑1 (MQ=255)                                gCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTTGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGGTTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTTATCCGATAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACaa                                                                                                                                                                                                                       >  2:121927/1‑249 (MQ=255)                                 cGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGt                                                                                                                                                                                                                    <  2:181354/251‑1 (MQ=255)                                     tgttgCTTGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACAGCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTgcg                                                                                                                                                                                                                 <  2:405548/250‑1 (MQ=255)                                     tgttgATGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGATTAGTACTGGCCATTTTTGTTAAAGACGTTAAACCGGGCGCGCGAAGAGCCACGGAGGGGACGCTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGAGTTTAATCAGCTAC‑AACATCCCCAGTTTCGCGCCAAAACCGGTGGCCTTAACAAAGTgcgg                                                                                                                                                                                                                >  2:28812/1‑250 (MQ=255)                                         gCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTT‑GTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATg                                                                                                                                                                                                            <  2:7003/250‑1 (MQ=255)                                                   tttttACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGt                                                                                                                                                                                                  >  1:212057/1‑251 (MQ=255)                                                    ttttACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGt                                                                                                                                                                                                  <  1:220163/250‑1 (MQ=255)                                                      ttACTGGTGGCGGTGGCCGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCCGGCGTTTAATCAGCTAC‑AACCTCCCCCGTTTCGCGCCAATATCGGGGGCATTAACAAAGTGCGCAATGAAAATCATGTCtt                                                                                                                                                                                               >  2:263338/1‑250 (MQ=255)                                                              ggcggtggcAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCg                                                                                                                                                                                       <  1:45414/251‑1 (MQ=255)                                                                             gTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAg                                                                                                                                                                        <  2:17157/251‑1 (MQ=255)                                                                                 gTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTAtt                                                                                                                                                                    <  1:80763/251‑1 (MQ=255)                                                                                   acgggCCATTCTTGTCAAAGAAGTTAAACCGGGCGTGCGCAGAGCAACGGAGGGGACGTGGAGCGGCACCGCAACGGTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCGCTGGGGACCCAACGCATGCGCAAGTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGACAATATCGGGGGCATTAACAAAGAGCGCAATGAATATCACGTCTATAAGCCCGATGCGCAGGGCAAAGTATTGt                                                                                                                                                                  <  2:87272/248‑1 (MQ=255)                                                                                      ggCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTc                                                                                                                                                                <  1:22904/250‑1 (MQ=255)                                                                                      ggCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTc                                                                                                                                                                <  2:135757/250‑1 (MQ=255)                                                                                      ggCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGGCCCGATTTGCTCTCTGGTGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTAACGCGCCAATATCTGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGtt                                                                                                                                                                 <  2:224160/249‑1 (MQ=255)                                                                                            ttttGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCg                                                                                                                                                          <  2:212057/250‑1 (MQ=255)                                                                                             tttGTCAAAGAAGTTAAAACGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑CACATCCCCCGTTTCGCGCCACTATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGAATTGTTCGATTCGGc                                                                                                                                                        >  1:357832/1‑251 (MQ=255)                                                                                                  cAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATa                                                                                                                                                   <  1:158708/251‑1 (MQ=255)                                                                                                   aaaGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACAGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCCACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAAAAAAGCGCGCAATGCATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGTTTCGGCAAATa                                                                                                                                                   >  2:390562/1‑250 (MQ=255)                                                                                                              aCCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAACTAAAGCCGTTGGAc                                                                                                                                       >  2:382029/1‑251 (MQ=255)                                                                                                              aCCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAAGCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGAGTCGGCAAATAAAGCCGTTGGAc                                                                                                                                       <  1:72241/251‑1 (MQ=255)                                                                                                              aCCGGGCGAGCGAAGAGCAACGAAGGGGACGTTGATCGGCACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCACCGCATGGGCAACTGGCGCAGGCGTTTAATAAGCTAC‑AACATCCACCGTTTCGCGCCAATATAGGGGGCATTAACAAAGGGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGTCAAATAAAGCCGTTGGAc                                                                                                                                       >  2:114600/1‑251 (MQ=255)                                                                                                                cGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑g                                                                                                                                    <  1:24008/251‑1 (MQ=255)                                                                                                                  ggCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑gg                                                                                                                                   >  1:253153/1‑250 (MQ=255)                                                                                                                  ggCGTGCGAAGAGCAACGGAGGGGACGGTGATCGACACCGCAACGTTGCCGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCTGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GGa                                                                                                                                  <  1:424197/251‑1 (MQ=255)                                                                                                                        ggAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GGATTATTc                                                                                                                            <  2:246654/250‑1 (MQ=255)                                                                                                                         gAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GGATTATTcg                                                                                                                           >  2:221641/1‑251 (MQ=255)                                                                                                                               cAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGCGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGGGCGAAATGCATATCATGTCTATATGACCGATGCGCAGGGCAACGGATTGTTCGATTCGGCAAATaaa                                                                                                                                                 >  2:56545/1‑224 (MQ=255)                                                                                                                               cAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCGCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATaaa                                                                                                                                                 <  1:56545/224‑1 (MQ=255)                                                                                                                                         gACGTTGATCGACACCGAAACGTTACTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑ACCATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GGCTTATTCGCGCTGGAATGACGtat                                                                                                           >  2:222023/1‑249 (MQ=255)                                                                                                                                             ttGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GGATTATTCGCGCTGGAATGACGTCTGGCt                                                                                                       >  1:134350/1‑251 (MQ=255)                                                                                                                                                       ccGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTT‑GGTTCGGGAAATAAAACCGGTGGACAAGGATT‑TTCGCGCTGGAATGAAGTCCGGCTAACGTTGGGTg                                                                                            >  2:267351/1‑251 (MQ=255)                                                                                                                                                              gTTGCTGGCGGAGCTGCAGCGTCCCGATTTGCTCTCTGGGGCCCCAACGCATGGGCCACTGGCCCAGGCGTTTAATAAGCTAC‑AAAATGACCCGTTTCGCGCCAACATAGGTGGCATTAACAAAGTGAGCAATGAATATCATGTCTATCTGCCCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAgg                                                                                      >  2:441258/1‑250 (MQ=255)                                                                                                                                                                  cTGGCGGAGCCGCCGCGGCCCGATTTCCTCTTCGGGGAGCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGCCATTAACAAAGTGCGCCATGAATATACTGTCTATATGACCGATGAGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GGATTATTTGCGCTGGAATGACGTCTGGCTAACGCTGCGTGGTCCGTATgg                                                                                  >  2:213116/1‑251 (MQ=255)                                                                                                                                                                    ggCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAAc                                                                                                                                                                                                                         <  1:85533/115‑1 (MQ=255)                                                                                                                                                                    ggCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAAc                                                                                                                                                                                                                         >  2:85533/1‑115 (MQ=255)                                                                                                                                                                        gAGCGGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTGTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GGATTATTCGCGCTGGAATGACGGCTGGCTAACGTTGCGTGGTCAGTATGGTgcgcg                                                                            <  2:290189/251‑1 (MQ=255)                                                                                                                                                                            tGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCgcagc                                                                        >  1:186485/1‑251 (MQ=255)                                                                                                                                                                              gcgcGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCAc                                                                      <  1:324948/251‑1 (MQ=255)                                                                                                                                                                                       gATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCaaaa                                                             <  1:100226/251‑1 (MQ=255)                                                                                                                                                                                            gCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAGAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGGGCGCGCAGCACGTTGCAAAATCCTg                                                        >  2:325587/1‑251 (MQ=255)                                                                                                                                                                                                   ggggACCCAACGCATGGGCAACTGGCGCCGGCGTTTAATCAGCAAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GGCTTATTCGCGCTGGAATGACGTCTGGCTAAAGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATcc                                                 >  1:315271/1‑251 (MQ=255)                                                                                                                                                                                                    gggACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATAGCGGGGGCATTAACAAAGTGCGTAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGTTGCGCGCAGCACGTTGCAAAATCCTGCCGATccc                                                <  2:315271/251‑1 (MQ=255)                                                                                                                                                                                                      gACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCg                                               >  1:429879/1‑250 (MQ=255)                                                                                                                                                                                                       aCCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGaa                                             <  1:400282/251‑1 (MQ=255)                                                                                                                                                                                                              gCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCtg                                      >  1:191388/1‑251 (MQ=255)                                                                                                                                                                                                               cATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGAACGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCtgt                                     >  1:404364/1‑251 (MQ=255)                                                                                                                                                                                                                 tGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GGATTATTCGCGCTGGAATGACGTCt                                                                                                           >  2:94416/1‑179 (MQ=255)                                                                                                                                                                                                                 tGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GGATTATTCGCGCTGGAATGACGTCt                                                                                                           <  1:94416/179‑1 (MQ=255)                                                                                                                                                                                                                       aCTGGCGCAGGCGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGatgtatg                             <  1:126287/251‑1 (MQ=255)                                                                                                                                                                                                                                  cGTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GTATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGTCTCACCg                  <  2:404364/251‑1 (MQ=255)                                                                                                                                                                                                                                   gTTTAATCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCg                  >  1:160167/1‑250 (MQ=255)                                                                                                                                                                                                                                      tAATCAGCTAC‑AACATCCCCCGTTTTGCGCCAATATCGGTGGCATTAACCAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGAtt               >  2:262993/1‑250 (MQ=255)                                                                                                                                                                                                                                       aaTCAGCTAC‑AACATCCCCCGTTTCACGCCAAAATCGGTGGCATTAACAAAGAGCGCAATGAATATCATGTCTATATGAACGATGCGCAGGGCAAAGTATTGTTCGATTCGGCCAATAAAGCCGTTGGACA‑GGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCCCCGATTAt             >  2:259854/1‑251 (MQ=255)                                                                                                                                                                                                                                         tCAGCTAC‑AACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATgg           >  1:354014/1‑251 (MQ=255)                                                                                                                                                                                                                                                   acaTCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAAGTAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTcgc  <  2:414984/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       AGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTAC‑AACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACA‑GGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGC  >  NC_000913/4636679‑4637170

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