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	96.4%	188.2 / NA	56	R77P (CGC→CCC)	creC	sensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (27/27);  minor base A (2/0);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTG‑GTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTC  >  NC_000913/4636695‑4637168                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    aaTGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTG‑GTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATc                                                                                                                                                                                                                                  <  1:151275/250‑1 (MQ=255) aaTGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTG‑GTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGGTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGTACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCg                                                                                                                                                                                                                                 <  1:233096/251‑1 (MQ=255) aaTGCGTATCGGCATGCGGTTGTTGCTGGGCTATTGTTTACTG‑GTGGCGGTGTCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGTGCGTGCGAAGAGCAACGGAGGGGACGTTGAGCGACACCGCAACGTTGCTGGCGTAGCTGGCGCGTTCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCTCGCCAATATCg                                                                                                                                                                                                                                 <  1:326040/251‑1 (MQ=255)  aTGCGTATCGGCATGCGGTTATTGCTGGGCTATTTTTTACTG‑GTGGCGGTGGCAGCCTGGGTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGGGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGGGACGGCGTTTAATCCGCTACAACATCCCCCGTTTCGCGCCAATATCgg                                                                                                                                                                                                                                >  2:365111/1‑251 (MQ=255)  aTGCGTATCGGCATGAGGCTGTTGATGGGCTATTTTTTACTG‑GTGGCGGTGCCAGGCTGGGTCGGACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGGGCAACGGAGGGGACGTTGATAGACACCGCAACGTTGCTGGCGGAGCGGGCGAGTCCCGATTTGCTCTCTGGGGACCCACCGCACGGGGACCTGGCGCCGGCGTTTAATCAGGCACAACATCCCCCGTTTCGCGCCCatatagg                                                                                                                                                                                                                                >  2:47779/1‑248 (MQ=255)    gcgtttcggcattagtttgttgCTGGGCTATTTTTTA‑TGTGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGTGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGTCGCAGGCGTGTAATCAGCTACAACATCCCCCTTTTCGAGCCAATATCg                                                                                                                                                                                                                                 <  2:425752/232‑1 (MQ=255)    gCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTGTCCTG‑GTGGCGGTGGCAGCCTGGTTCGTACTGGCCATGTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTGGCTGGCGGAGCTGGCGCGTCCTGATTTGTTCGCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTGTAATCAGCTACATCATCCCCCGTTGCGCGCTAATAACGGt                                                                                                                                                                                                                               <  2:489984/250‑1 (MQ=255)     ggTATCGGCATGCGGTTGGTGCTGGGCTATTTTTTACTG‑GTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTgg                                                                                                                                                                                                                             <  2:374269/250‑1 (MQ=255)       tATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTG‑GTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGAAGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGGGGc                                                                                                                                                                                                                            >  1:110283/1‑250 (MQ=255)           gggaTGCGGTTGTTGCTGGGCTATTATTTACTG‑GTGGCGGTGGCAGCCTGGTTCGTACTGGCCAGTTTTGTCAAAGAAGTGAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAACTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTaa                                                                                                                                                                                                                       <  2:313884/248‑1 (MQ=255)             caggCGGTTGTGGCTGGGCTATTTTTTACTG‑GGGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGTACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCTGTGGCAGTAACa                                                                                                                                                                                                                     <  2:234113/248‑1 (MQ=255)              aTGCGGTTGTTGCTGGGCTATTTTTTACTG‑GTGGCGGTGGCAGCCTGGTTCGTAATGGCCATTTTTGTCAAAGAAGTTAAACCGGGAGGGGGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGGTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCCTGGGCAACTGGCGCCGGCGTTTAATCAGCTCAAACATCACCCGTTTCGCACCCATATAGGTAGCATTagca                                                                                                                                                                                                                     >  1:218145/1‑247 (MQ=255)                gCGGTTGTTGCTGCGCTATGTTTTCCGG‑GTGGCGGTGGCAGCATAGTTCGTACTGGCCATGTTTGTCAAAGACGTGAAACCGGTCGTGCGAAGAGCAACGGAGGGGACGTTGATCGTCACCGCAACGTTGCTGGCGGATCTGGCGCGTCCCGATTTGCGCTCGTTGGTCCCAACGCATTTGCAGCTGGCGTATGTTTTTAATCAGCTACAACATCCCCCCTTTCGCGGCAATGTCGGTGTCATTAACaaa                                                                                                                                                                                                                   <  2:336535/250‑1 (MQ=255)                   gttgttGCTGGGCTATTTTTTACTG‑GTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCACCTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGt                                                                                                                                                                                                                 >  2:10453/1‑249 (MQ=255)                   gttgttGCTGGGCTATTTTTTACTG‑GTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGt                                                                                                                                                                                                                 <  1:10453/249‑1 (MQ=255)                         gCTGGGCTATTTTTTTCTG‑GTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTTCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATTTCGGTGGCATTAACAAAGTGCGCAAt                                                                                                                                                                                                          <  2:265748/250‑1 (MQ=255)                         gCTGGGCTATTTTTTACTG‑GTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTGAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGATCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTTGCATTAACAAAGGGCGCAATg                                                                                                                                                                                                         <  2:230268/251‑1 (MQ=255)                            gggCTATTTTTTACTG‑GTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAat                                                                                                                                                                                                      <  1:365393/251‑1 (MQ=255)                                    ttttACTG‑GTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTc                                                                                                                                                                                              >  1:459405/1‑251 (MQ=255)                                      ttACTG‑GTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCt                                                                                                                                                                                             <  1:422317/250‑1 (MQ=255)                                                  ggtggcAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATg                                                                                                                                                                                 >  1:78214/1‑251 (MQ=255)                                                    gggcAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATgcg                                                                                                                                                                               <  1:345309/250‑1 (MQ=255)                                                     ggcAGCCTGGCTCGTACTGGCAATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGCCGTTGATCGACACCGCAACGTTGCTGGCGGAGATGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCCCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATgcgc                                                                                                                                                                              >  2:202544/1‑251 (MQ=255)                                                       cAGCCTGGTTCGTACTGGCCAGTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGGCCCGATTTGCTCGCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAg                                                                                                                                                                            <  2:154425/251‑1 (MQ=255)                                                             ggTTCGTACTGGCCACTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAATAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCACCTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCAGGTCTATATGACCGATGCGCAGGGCa                                                                                                                                                                        <  2:339522/249‑1 (MQ=255)                                                                       ggCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCg                                                                                                                                                            <  1:187775/251‑1 (MQ=255)                                                                       ggCCAGTTTTGTCAAAGAAGTTAAACCGTGCTTGCGAATAGGAACGGGGGGGACGTTGATCGACTCCTCAACGGTGCTTGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCGGGCGTTTAATCAGCTACAACGTCCCCCGGATCGCGGCAATATCGGTTGCATTAACAAAGTGCGCAATGAATATCATGTCTATAGTACCGCGGCGCAGTGGAAAGTATGGTTCg                                                                                                                                                            <  2:503267/251‑1 (MQ=255)                                                                              tttGTCAAAGAAGTTACACAGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGc                                                                                                                                                     >  1:355172/1‑251 (MQ=255)                                                                                tGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCa                                                                                                                                                    >  1:218088/1‑250 (MQ=255)                                                                                  tCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAAt                                                                                                                                                 >  1:401606/1‑251 (MQ=255)                                                                                    aaaGAAGTTACCCCGGGCGTGCGAAGAGCAACAGCGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAAAATCCCCCGGTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAGTGAATATCATGTCTATATGACCGATGCGCAGGGCAACGTATTGTTCGATTCGGCAAATa                                                                                                                                                >  2:53204/1‑250 (MQ=255)                                                                                      agaagTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGt                                                                                                                                                                                                                 <  1:288541/183‑1 (MQ=255)                                                                                      agaagTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGt                                                                                                                                                                                                                 >  2:288541/1‑183 (MQ=255)                                                                                         agTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCg                                                                                                                                          >  2:393853/1‑251 (MQ=255)                                                                                           ttAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGtt                                                                                                                                        >  1:100265/1‑251 (MQ=255)                                                                                              aaCCGGGCGTGCGAAGCGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTgg                                                                                                                                      >  2:142789/1‑250 (MQ=255)                                                                                                ccGGGGGTGCGAAGAGCAACGGTGGGGACGGTGATCGACACCGGAACGTTGCTGGCGGAGCTGGCGCGTCTCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGGCTATATGACCGATGCGCAGGGCAAGGTATTGTTCGAGTCGGCAAATAAAGCCGGTGGACa                                                                                                                                   <  2:49205/251‑1 (MQ=255)                                                                                                  gggCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAgg                                                                                                                                 >  1:195509/1‑251 (MQ=255)                                                                                                      gTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAACGCCGTTGGACAGGatt                                                                                                                              >  1:27928/1‑250 (MQ=255)                                                                                                      gTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATGCGTCAAATAAAGCCGTTGGACAGGatta                                                                                                                             <  2:468829/251‑1 (MQ=255)                                                                                                            agagCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTcgcg                                                                                                                       <  1:153293/251‑1 (MQ=255)                                                                                                              agCAACGGAGGGGACGTTGATCGACACCGAAACGTTGCTGGCGGAGCTGGGGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTAGCGCt                                                                                                                     >  2:223924/1‑251 (MQ=255)                                                                                                                              ttGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCt                                                                                                     >  1:113088/1‑251 (MQ=255)                                                                                                                                               gTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTATCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGt                                                                                    >  1:396565/1‑251 (MQ=255)                                                                                                                                                                 gcgTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGt                                                                               >  2:174241/1‑238 (MQ=255)                                                                                                                                                                 gcgTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGt                                                                               <  1:174241/238‑1 (MQ=255)                                                                                                                                                                   gTCCCGATGTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTAGTCGCGCTTGAATGACGTCTGGCTAACGTTGCGTGGTCAGTTTGGTGCGCGCAGCACGTTg                                                                <  2:366567/251‑1 (MQ=255)                                                                                                                                                                           ttGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATcc                                                        >  1:382028/1‑251 (MQ=255)                                                                                                                                                                            tGCTCTCTGGGGACCCAACGCATGGGCACCTGGCGCAGGCGTTTCATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGAATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCt                                                       >  2:420869/1‑251 (MQ=255)                                                                                                                                                                                             aCGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTc                                      <  1:47779/251‑1 (MQ=255)                                                                                                                                                                                                    ggCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAAt                                                                                                                                                 >  1:14586/1‑137 (MQ=255)                                                                                                                                                                                                    ggCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAAt                                                                                                                                                 <  2:14586/137‑1 (MQ=255)                                                                                                                                                                                                          tGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCACGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGtt                         >  1:112618/1‑251 (MQ=255)                                                                                                                                                                                                                         aTCAGATACAACATCACCCGTTTAGCGCCAATATCGGTGACATTAACAAAGCGCGCCATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTAGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATACCATCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGAAAAATCCTGCCGATCCCGACAGTTCTGTGCTGTATGTTGACGCACCGATTATg          >  2:500737/1‑251 (MQ=255)                                                                                                                                                                                                                           cAGCTACAACATCCCCCGTTTCGCGCCAATCTCGGAGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGGGCGCGCAGCACGTTGCAAAATCCTGCCGATCACGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATgg         >  1:152015/1‑250 (MQ=255)                                                                                                                                                                                                                                 caacaTCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTc  >  1:311469/1‑251 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    AATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTG‑GTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTC  >  NC_000913/4636695‑4637168

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: ‑