Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	92,687	A→G	P425P (CCA→CCG)	ftsI →	peptidoglycan DD‑transpeptidase FtsI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	92,687	0	A	G	95.1%	194.6 / 1.1	62	P425P (CCA→CCG)	ftsI	peptidoglycan DD‑transpeptidase FtsI
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (33/25);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01

GTGGCACGCTACAGCGAATTAACCCTGACCGGGGTATTACAGAAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCATTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACTGTGGTGCATATGATGGAAAGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTA  >  NC_000913/92442‑92932                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      gTGGCACGCTACAGCGAATTAACCCTGACCGGGGTATTACAGAAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTAc                                                                                                                                                                                                                                                   >  1:77770/1‑250 (MQ=255) gTGGCACGCTACAGCGAATTAACCCTGACCGGGGTATTACAGAAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACa                                                                                                                                                                                                                                                  <  2:198773/251‑1 (MQ=255)       cGCTACAGCGAATTAACCCTGACCGGGGTATTACCGAAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTTGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATGTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAgc                                                                                                                                                                                                                                            <  2:42946/251‑1 (MQ=255)         cTACAGCGAATTAACCCTGACCGGGGTATTACAGAAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAgcgc                                                                                                                                                                                                                                          >  1:44753/1‑251 (MQ=255)            cAGCGAATTAACCCTGACCGGGGTATTACAGAAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAg                                                                                                                                                                                                                                       <  1:170794/251‑1 (MQ=255)             aGCGAATTAACCCTGACCGGGGTATTACAGAAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCAGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGGCTGACATAGAGAGGGCCACCTTCTCTTTAGGCTACGGGCTAATGGTAACACCGCTACAGTTAGCGCGAGt                                                                                                                                                                                                                                      >  2:179043/1‑251 (MQ=255)               cGAATTAACCCTGACCGGGGTATTACAGAAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCt                                                                                                                                                                                                                                    >  1:447807/1‑251 (MQ=255)                 aaTTAACCCTGACCGGGGTATTACAGAAGTCGAGTAACGTCGGTGGTTCCAAGCTGGCGTTAGCGACGCCGTCCTCAGCGTTAGTAGAGACTTACCCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGGTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACGCCGTTACAGTTAGCGCGAGTCttc                                                                                                                                                                                                                                  <  2:232497/251‑3 (MQ=255)                           gACCGGGGTATTACAGAAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCATTACAGTTAGCGCGAGTCTACGCAACTATCg                                                                                                                                                                                                                        >  1:152869/1‑251 (MQ=255)                           gACCGGGGTATTACAGAAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAAAAACGGTGGTCTGACATAGAGAGGGCGACCTTCTCTTTCGGCTACGGGCTAATGGTAACCCCGTTACAGTTAGCGCGGGTCTACGCAACTATCg                                                                                                                                                                                                                        >  2:431094/1‑251 (MQ=255)                               ggggTATTACAGAAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGc                                                                                                                                                                                                                      >  1:14971/1‑249 (MQ=255)                                  gTATTACAGAAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGACAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCt                                                                                                                                                                                                                  >  1:265862/1‑250 (MQ=255)                                      tACAGAAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGc                                                                                                                                                                                                             >  1:99061/1‑251 (MQ=255)                                         aGAAGTCGAGTAACGTCGGTGTTTCCAAGCTTGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCATTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCAtt                                                                                                                                                                                                          <  2:152869/251‑1 (MQ=255)                                               cGAGTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCGCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGc                                                                                                                                                                                                    <  1:224270/251‑1 (MQ=255)                                                 aGTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGccc                                                                                                                                                                                                  >  1:533427/1‑251 (MQ=255)                                                  gTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCa                                                                                                                                                                                                 <  1:182410/251‑1 (MQ=255)                                                                       cggCGTTAGCGATGCCTCCCTCAGCGTTAGTAGATTCTTACTCACGCTTTGGACTGGGAAAAGCGACCAATTTGGTGTTGGTCGGAGAACGCAGTGGTTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCGCCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTGCAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGAc                                                                                                                                                                            <  1:103182/250‑1 (MQ=255)                                                                          cGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGAcccc                                                                                                                                                                         >  2:191696/1‑251 (MQ=255)                                                                                      cGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTAAGGCCTTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCggg                                                                                                                                                             >  2:51523/1‑250 (MQ=255)                                                                                           tCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACg                                                                                                                                                        >  1:497978/1‑251 (MQ=255)                                                                                                gTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGCGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCACCTATCGGCCGCTACGGCATTTATCGCCCCCTGTCGATTACCCAAGTTGACCCCCCGGTTCCCGGTGAACGGGTCt                                                                                                                                                   >  2:254219/1‑251 (MQ=255)                                                                                                   agtagATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTACTGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTTACCCCCCGGTTCCCAGTGAACGTGTCTTcc                                                                                                                                                >  1:429220/1‑251 (MQ=255)                                                                                                   agtagATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTcc                                                                                                                                                >  1:538033/1‑251 (MQ=255)                                                                                                             tactCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATcc                                                                                                                                        <  1:523723/249‑1 (MQ=255)                                                                                                                 cACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTg                                                                                                                                    <  1:51523/249‑1 (MQ=255)                                                                                                                          gACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACcgtg                                                                                                                         <  2:54372/251‑1 (MQ=255)                                                                                                                          gACTGGGAAAAGCGACCAAGTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACcgtg                                                                                                                         <  1:62382/251‑1 (MQ=255)                                                                                                                           aCTGGGAAAAGCGACCAATGTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCgtgg                                                                                                                        <  2:14971/251‑1 (MQ=255)                                                                                                                                  aaaGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGATTCCATTGTCCGCCCCGTGGTGCatat                                                                                                                 >  2:412990/1‑251 (MQ=255)                                                                                                                                                ggggTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGCCATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTTGATTACCAAAGTTTACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCAATGTCCGCCCCGTGGTGCATATGATGGAAAGCgggg                                                                                                   >  2:429046/1‑248 (MQ=255)                                                                                                                                                      ggTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCCTTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTATTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTAc                                                                                             >  1:245520/1‑251 (MQ=255)                                                                                                                                                             gaACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGggcgg                                                                                      >  1:168599/1‑251 (MQ=255)                                                                                                                                                               aCGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCa                                                                                                                                                                                                            >  1:235814/1‑131 (MQ=255)                                                                                                                                                               aCGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCa                                                                                                                                                                                                            <  2:235814/131‑1 (MQ=255)                                                                                                                                                                  cAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGggcggcggcg                                                                                 >  2:268409/1‑251 (MQ=255)                                                                                                                                                                    gTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGggcggcggcggc                                                                               <  1:397360/251‑1 (MQ=255)                                                                                                                                                                           atatCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAg                                                                         <  2:67644/250‑1 (MQ=255)                                                                                                                                                                            tatCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAggc                                                                       >  1:436464/1‑251 (MQ=255)                                                                                                                                                                               ccTCAAAAACAACGGTGGTCGGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCTGTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAggcggc                                                                    <  2:538033/251‑1 (MQ=255)                                                                                                                                                                                 tCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGgcggcg                                                                   >  1:354711/1‑250 (MQ=255)                                                                                                                                                                                   aaaaaCAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGAAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGTCGAtt                                                                >  2:338339/1‑251 (MQ=255)                                                                                                                                                                                    aaaaCAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTa                                                               <  2:44753/251‑1 (MQ=255)                                                                                                                                                                                     aaaCAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTa                                                               >  1:175918/1‑250 (MQ=255)                                                                                                                                                                                         aacGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGc                                                          <  1:431094/251‑1 (MQ=255)                                                                                                                                                                                            ggtggtCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCtat                                                       <  2:31590/251‑1 (MQ=255)                                                                                                                                                                                            ggtggtCTGACATAGAGAGGGCCACCTTCTCTTTCAGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCACAGTTGACCCCCCGGTTCCCGGAGAACGTGTCTTACCGGAATCCATTGTCCGCACCGTGTTGCATATGATGGAAAGCGTGACGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAATGCtat                                                       >  1:419351/1‑251 (MQ=255)                                                                                                                                                                                             gtggtCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCtatc                                                      >  1:557880/1‑251 (MQ=255)                                                                                                                                                                                              agacggggtctcttccgatcttccTTCTCCTTCGGCTTCGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCtatcg                                                     <  1:170257/229‑1 (MQ=255)                                                                                                                                                                                               ggtCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCtatcg                                                     <  1:445806/250‑1 (MQ=255)                                                                                                                                                                                                      cATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCtat                                                       >  1:237484/1‑241 (MQ=255)                                                                                                                                                                                                      cATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCtat                                                       <  2:237484/241‑1 (MQ=255)                                                                                                                                                                                                                  cACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAACCCGGt                                  >  2:340400/1‑250 (MQ=255)                                                                                                                                                                                                                   aCCTTCTCTTTCGGCTACGGCCTAATGGAACCACCGTTACAGTTAGCCCGGTTCTCCCCAACTATCGGCAGCTCCGGCATTTATCGCCCCCTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGACCGTGTCTTCCCGGAATCCATTGTCCGCCCCGTGGTCCAAATGTTGGAAACCGTGGCGCTCCCGGCCGGCGGCGGCGTGAGGGCGGCGATTAAAGGCtatagagaccggaagcggttcggttg                                >  2:170257/1‑228 (MQ=255)                                                                                                                                                                                                                        ctctTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTg                                                                                                    >  1:434815/1‑178 (MQ=255)                                                                                                                                                                                                                        ctctTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTg                                                                                                    <  2:434815/178‑1 (MQ=255)                                                                                                                                                                                                                                gCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGGCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTc                   <  2:286348/251‑1 (MQ=255)                                                                                                                                                                                                                                 cTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCg                  >  2:445089/1‑251 (MQ=255)                                                                                                                                                                                                                                   aCGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTAACGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCgg                 <  2:427872/250‑1 (MQ=255)                                                                                                                                                                                                                                                aaCACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATAGCCCACTGACGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCCTTGTCCGCACCGTGGTGCCTATGATGGAAAGCGTGGCGCTCCCGGGCGGCGGCGGCGTGAAGGCGGCGCTTAAAAGCTATCGTACCGCCATTAAAACCGGTACCGCGAACAAGGTCGGGCCGGCCGGgcgcc   >  2:147466/1‑250 (MQ=255)                                                                                                                                                                                                                                                 acacCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTa  <  1:491118/251‑1 (MQ=255)                                                                                                                                                                                                                                                     cATTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGc                                                                                                                               <  2:188859/122‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      GTGGCACGCTACAGCGAATTAACCCTGACCGGGGTATTACAGAAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCATTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACTGTGGTGCATATGATGGAAAGCGTGGCGCTACCAGGCGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACCGCGAAAAAGGTCGGGCCGGACGGTCGCTA  >  NC_000913/92442‑92932

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