Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,637,151	A→G	A152A (GCA→GCG)	creC →	sensory histidine kinase CreC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,637,151	0	A	G	92.9%	109.5 / ‑0.2	42	A152A (GCA→GCG)	creC	sensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (26/13);  total (28/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.90e-01

CGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTCTCCCCGATCTCGGTAGTAGC  >  NC_000913/4636904‑4637397                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       cGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGCGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCg                                                                                                                                                                                                                                                     >  1:322413/1‑251 (MQ=255)      aaTCAGCTACAGCTACCCCCGTTTCGCGCCAAAAGCGGGGTCATTAACATAGTGCGCAAGGAATAGTAAGTCTATATGACCGTTTCGCAGGGTAAAGTATTGTGCGATTCGGCAAATAAAGCCGTTGGGCAGTATTTGTCGCGCTGGTATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGGTCTTTGCAAAATCCTGCCGATTCCGAAAGTTCGGTGATGTATGTTGCCTCGCCGATTAt                                                                                                                                                                                                                                                <  2:75929/251‑1 (MQ=255)                       cccGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAACGAATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGACTGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCtt                                                                                                                                                                                                                                >  1:179017/1‑250 (MQ=255)                          gTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTg                                                                                                                                                                                                                            >  1:310104/1‑251 (MQ=255)                                    aaTATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCGAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTg                                                                                                                                                                                                                    >  1:193692/1‑249 (MQ=255)                                               cATTAACAAAGAGCGCAATGAATATCATGTATATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGATCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAAc                                                                                                                                                                                                       >  2:68967/1‑251 (MQ=255)                                                         gTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACgcggc                                                                                                                                                                                             >  1:367142/1‑251 (MQ=255)                                                          tGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCAGCGCCGATTATGGACGGCTCGAGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACgcgggg                                                                                                                                                                                            >  2:315619/1‑249 (MQ=255)                                                                  gAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTcc                                                                                                                                                                                    <  1:315619/251‑1 (MQ=255)                                                                    atatCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGGTTTGAGCGTAGGCAAACCGACCGCGGCGATGGCTCCgg                                                                                                                                                                                  >  1:71407/1‑251 (MQ=255)                                                                             cTATATGACAGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCTCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAACGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTCTTGGCGGTTTGTACGTCGGCAAAACGAACGCGGAGATGGCTCCGGTCAttgatc                                                                                                                                                                         >  2:152/1‑247 (MQ=255)                                                                               atatGACCGATTCGCAGGGTAAAGTATTTTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTTGGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATGCCGAAAGTACGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGGTTATGGGCGTTGTGAGCGTAGGCAAACCGATCGCGGCGATGGCTCCGGTCATTAAGCGt                                                                                                                                                                       <  1:56891/251‑1 (MQ=255)                                                                                tatGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGTCAGGATTATGCGCGGTGGAATGACGTATGTCTAACTTAGCGTGGTCAGAATGGTTCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAACGTTAGGTGATGTTTGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGCCAAACAGAAGGCGGCGagatctcagttcagtagacgt                                                                                                                                                                       >  2:120509/1‑230 (MQ=255)                                                                                      ggcttcgcaggggatattATTTTTCGATTCGGAAAATAAAGCCGTTTGGCTGGATTATGCGCGTGGGAATGACTTCTGGCTAACGTTGCGGGGTCAGTTTGGTGCGCGCAGCTCGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTGTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTagcgagc                                                                                                                                                                <  2:16459/234‑1 (MQ=255)                                                                                              aGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAAtt                                                                                                                                                        <  1:47743/251‑1 (MQ=255)                                                                                                ggtaaaGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAAtttt                                                                                                                                                      >  1:34605/4‑251 (MQ=255)                                                                                                  cAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTAt                                                                                                                                                    >  1:9806/1‑251 (MQ=255)                                                                                                        aTTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGcc                                                                                                                                               >  1:157051/1‑250 (MQ=255)                                                                                                          tGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAgc                                                                                                                                            >  1:363832/1‑251 (MQ=255)                                                                                                           gTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATGCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAgc                                                                                                                                            <  1:87250/250‑1 (MQ=255)                                                                                                               gATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCAt                                                                                                                                       >  1:354692/1‑251 (MQ=255)                                                                                                                           aaaGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGAt                                                                                                                           >  1:167497/1‑251 (MQ=255)                                                                                                                             agcatttGGTCAGGAATATTCGCACTAGATTGACGTCGGGCAAGCGTTGCGGGGTCGGTATGGTGCGTGCTGTATGTTGCAAAATCCTGCCGATCCCGAAAGTCCGGTGATGTATTTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGGTTTGAGCGTAGGCAAACCGATCGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGACTTTTATGGGCCAGCGCCATTTTGCTGGGGATTg                                                                                                                         <  2:58162/246‑1 (MQ=255)                                                                                                                             aGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTg                                                                                                                         >  1:9928/1‑251 (MQ=255)                                                                                                                               cggTTGGACAGGATTTTTCGGGCTGGAATGAGGTCTGTCTAACGTTGCGTGGTCGGTATGGTGGGCGCAGGTCGTTGCAAAATCCTGCCGATCCCGAATGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGCGGGTTATTGGCGTTTTGAGCGTAGGCAAACCGAATGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGCGCCAGCGCCATTTTGTTGGGGATTGCa                                                                                                                       <  2:193692/249‑1 (MQ=255)                                                                                                                               ccGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCa                                                                                                                       <  1:249729/251‑1 (MQ=255)                                                                                                                                 gTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACt                                                                                                                     >  1:265707/1‑251 (MQ=255)                                                                                                                                   tgggcAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTgg                                                                                                                   >  1:285102/5‑251 (MQ=255)                                                                                                                                      gcAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTg                                                                                                                 >  1:242057/2‑250 (MQ=255)                                                                                                                                      gcAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTg                                                                                                                 >  1:214204/2‑250 (MQ=255)                                                                                                                                               ttCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGgcgc                                                                                                         >  1:131645/1‑249 (MQ=255)                                                                                                                                                              cGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTtggtgg                                                                                         >  1:40187/1‑250 (MQ=255)                                                                                                                                                                       aaCGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGct                                                                               >  1:20069/1‑251 (MQ=255)                                                                                                                                                                          gTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGctcc                                                                             >  1:345158/1‑249 (MQ=255)                                                                                                                                                                                    cattATGGTGCTCGCAGTACGTTGCAAAATCCTGCCTATCCCGAAAGTTCGGTTATGTATGTTGCCGCGCCGATTATTGACGGCTCGCGGCTTTTTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCt                                                                  <  1:302827/248‑1 (MQ=255)                                                                                                                                                                                                          tGCAAAATCCTGCCGATCCCGAAAGATCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTc                                            >  1:299609/1‑251 (MQ=255)                                                                                                                                                                                                             aaaaTCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACt                                         <  1:179178/251‑1 (MQ=255)                                                                                                                                                                                                               aaTCCTGCCGATCCCGAAAGGTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTTGTGGATCAACCGCTCCATTGCCCGGCCCACTCGCCATGCCGATTCCGTCACTg                                        >  2:67100/1‑250 (MQ=255)                                                                                                                                                                                                                       cGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGcc                                                                       >  1:168389/1‑211 (MQ=255)                                                                                                                                                                                                                           cccGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCg                            <  1:171011/250‑1 (MQ=255)                                                                                                                                                                                                                                         gATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCctcttcctgaa              <  1:135463/250‑8 (MQ=255)                                                                                                                                                                                                                                                    ccGCGCCGTTTATGGACGGCTTGCGGCTTATTGGCGGTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCAATAAGCGTAGCGAGCGGCGAATTTTAGGGGCCAGTGCCATTTTGCTTGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGc  <  2:40187/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       CGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTCTCCCCGATCTCGGTAGTAGC  >  NC_000913/4636904‑4637397

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