Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,054,855	G→A	Q442* (CAG→TAG)	torS ←	sensor histidine kinase TorS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,054,855	0	G	A	98.4%	191.7 / NA	63	Q442* (CAG→TAG)	torS	sensor histidine kinase TorS
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (31/31);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTGGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAGGGGAGTCAATATCACCGTCCAGCAGCCGTTGCAGCGCCTGCGTTTGTTCG  >  NC_000913/1054608‑1055096                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  agcGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCctagct                                                                                                                                                                                                                                                >  2:462981/1‑251 (MQ=255)           gTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCt                                                                                                                                                                                                                                      <  1:392217/251‑1 (MQ=255)              cAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCACATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCt                                                                                                                                                                                                                                   >  1:310042/1‑251 (MQ=255)               aaaGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCtt                                                                                                                                                                                                                                  >  1:318835/1‑251 (MQ=255)               aaaGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCtt                                                                                                                                                                                                                                  <  1:275884/251‑1 (MQ=255)                acgtgCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCA‑GGCATCATCGGCGCCAGAAACCCCGACTTTGCCTAGATGGCTTTTTCTgtatc                                                                                                                                                                                                                                 >  2:341145/5‑247 (MQ=255)                  gggCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCg                                                                                                                                                                                                                               <  2:143221/251‑1 (MQ=255)                       cATCGCTGACCGAAACATTGTTGCCACCTGCATCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCTGAGTCAGTAAGTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGGGCCGAGAATACCGTGCAGCGGTGTGCGGATTTCATGGC‑TCTTCGCCGCCAGAAACGCCGGTTTGTCCGAGCTGGCGATGTCTGCTTCCGCCCGt                                                                                                                                                                                                                          <  2:84115/251‑1 (MQ=255)                       cATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCg                                                                                                                                                                                                                           <  1:185901/250‑1 (MQ=255)                               aCCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCg                                                                                                                                                                                                                  <  1:598441/251‑1 (MQ=255)                                 cGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGgt                                                                                                                                                                                                                <  1:195978/251‑1 (MQ=255)                                  gAAATATTCTTGCCACCGGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCAGCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCATCAGTTGAGCAGTGCCGAGAATTCCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGgtg                                                                                                                                                                                                               <  1:576293/251‑1 (MQ=255)                                         tCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTAt                                                                                                                                                                                                         <  1:530348/250‑1 (MQ=255)                                               cACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGt                                                                                                                                                                                                  <  2:2980/251‑1 (MQ=255)                                                 ccTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTGCTATCACCAGTTc                                                                                                                                                                                                <  1:104432/251‑1 (MQ=255)                                                  cTGCTTCGATAACGGAATAATAGCGAATATCGTTGAGGTTGGTCAGCCACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAGGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTAACGCGGTGTGCGGCTCTCCTGGC‑TCCTCGCCGCCAGAAACGCAGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTcc                                                                                                                                                                                               >  2:556845/1‑251 (MQ=255)                                                    gCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTcc                                                                                                                                                                                               <  2:337728/249‑1 (MQ=255)                                                           tAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTc                                                                                                                                                                                       >  1:162783/1‑250 (MQ=255)                                                             gCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGc                                                                                                                                                                                    >  1:10993/1‑251 (MQ=255)                                                               ggAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGTCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTg                                                                                                                                                                                  <  1:364478/251‑1 (MQ=255)                                                                       cGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCt                                                                                                                                                                          >  2:480732/1‑251 (MQ=255)                                                                         agaATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCttt                                                                                                                                                                        <  1:566218/251‑1 (MQ=255)                                                                             tatCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGAcc                                                                                                                                                                    >  1:508576/1‑251 (MQ=255)                                                                             tatCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGAcc                                                                                                                                                                    >  1:32505/1‑251 (MQ=255)                                                                                agTTGAGGATGGGCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTgc                                                                                                                                                                 >  2:193743/2‑251 (MQ=255)                                                                                        aTGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCAGGGCGTTGAGGGGGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTGATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAg                                                                                                                                                         <  2:180209/251‑1 (MQ=255)                                                                                         tGGTCAGCAACGATTCAACAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGc                                                                                                                                                        >  1:189730/1‑251 (MQ=255)                                                                                                cAACGATTCACCAGAGTCAGTAATTGCACGCAAATCATCACGCTGGGCGTTACGTGCGGGGTTATCTGACAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTcgc                                                                                                                                                 >  1:365079/1‑251 (MQ=255)                                                                                                   cGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGt                                                                                                                                              >  1:245460/1‑251 (MQ=255)                                                                                                     aTTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCt                                                                                                                                                       >  1:83196/1‑240 (MQ=255)                                                                                                     aTTCACCAGAGTCAGTAATTGCCCGCAAAGCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCt                                                                                                                                                       <  2:83196/240‑1 (MQ=255)                                                                                                           cAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTc                                                                                                                                      <  1:429010/251‑1 (MQ=255)                                                                                                             gagTCAGTAATTGCCCGCAAATCATCACGCGGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCat                                                                                                                                    <  2:4604/251‑2 (MQ=255)                                                                                                                     aaTTTCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCGCATGGC‑GCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTTGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCTTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTCCGCGGTGGCGATTCAGTGCATGa                                                                                                                             <  1:600615/250‑1 (MQ=255)                                                                                                                                            gggCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCCGTGCATGAACATTGCTGCGCAACGcatccatc                                                                                                     >  2:100645/1‑251 (MQ=255)                                                                                                                                                    aGTGCGGGTTTATCTGCCAGCAGTTGAGCAGTGCCTAGAATACCGTACAGCGGTGTGCGGATCGCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACTCGCTTTTACCTTCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGccc                                                                                             <  2:629907/251‑1 (MQ=255)                                                                                                                                                     gTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCg                                                                                            <  1:285025/251‑1 (MQ=255)                                                                                                                                                      tGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCg                                                                                            >  1:61310/1‑250 (MQ=255)                                                                                                                                                      tGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCg                                                                                            >  1:516957/1‑250 (MQ=255)                                                                                                                                                        cGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTTCCTAGCTGGCTTTTTCTGCGTCCGCCCGTGCCTGGCGGTGTGCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGAGTCAGTGCATGAACATGGCTGCGAAACGCATCCATCAGCCGCCCGATg                                                                                         <  2:549969/251‑1 (MQ=255)                                                                                                                                                            gTTATCTGCCAGCAGTTGAGAAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTAt                                                                                     >  1:532802/1‑251 (MQ=255)                                                                                                                                                             ttATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGTCCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATc                                                                                    >  1:554836/1‑251 (MQ=255)                                                                                                                                                               agcTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCTGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCa                                                                                  <  2:512590/249‑1 (MQ=255)                                                                                                                                                                tCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGTTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATcc                                                                                   >  2:367298/1‑249 (MQ=255)                                                                                                                                                                 cTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAg                                                                                                                                    >  1:187828/1‑199 (MQ=255)                                                                                                                                                                 cTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCAGCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAg                                                                                                                                    <  2:187828/199‑1 (MQ=255)                                                                                                                                                                      agcagTTGAGCAGTGCCGAGACTACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCg                                                                           >  2:106921/1‑251 (MQ=255)                                                                                                                                                                         agTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTAc                                                                        >  1:568752/1‑251 (MQ=255)                                                                                                                                                                                    gCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGttt                                                             <  2:47626/251‑1 (MQ=255)                                                                                                                                                                                           aTACCGTACAGCGGTGTGCGGATCTCATTGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGa                                                       >  1:241700/1‑250 (MQ=255)                                                                                                                                                                                           aTACCGGACAGCGGAGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCAGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCGGCGAAACGCATCCATCAGCCGCCCGAGGGTTTCCAGCTCCCGTACGCCCGCGGTTTCCGGGa                                                       >  2:279937/1‑250 (MQ=255)                                                                                                                                                                                             aCCGGACAGCGGTGTGCTGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCGGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGTTGGCGATTCAGTGCATGAACATTGCGGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGaaa                                                     <  2:541828/250‑1 (MQ=255)                                                                                                                                                                                                    aGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAGGGGAGTc                                             >  1:14812/1‑251 (MQ=255)                                                                                                                                                                                                      cGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAGGGGAggca                                            >  1:74528/1‑247 (MQ=255)                                                                                                                                                                                                        gctttcggcgtacttGGC‑TCATCGCGGCCAAAAAGGCCTATTGTGCCTGGCTGGCTTTTTCTGCTCCCTGCCGTGCCTGTCGCGGGTCTATCACCATTTCCCGCAATTCCGCTGTACGCGCTTTGACCGGCGCCGCCATCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATACATCAGCCGCCCGATGGTATCCAGCTCCCGTACGACCGCGGTTTCCGGGAAAGTGGAGTCAata                                         <  2:338091/237‑1 (MQ=255)                                                                                                                                                                                                           ggCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTGTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATGCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGAGGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAGGTGAGTCACTATCa                                      <  2:445067/250‑1 (MQ=255)                                                                                                                                                                                                              ggATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCCGCGGTTTCCGGGAACGGGGAGTCAATATCACCg                                   >  1:164521/1‑251 (MQ=255)                                                                                                                                                                                                                        ggcgcATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCAGCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAGGGGAGTCAATATCACCGTCcagcag                           <  2:621473/246‑1 (MQ=255)                                                                                                                                                                                                                        gC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAGGGGAGTCAATATCACCGTCCAGCAGcc                         >  1:283591/1‑251 (MQ=255)                                                                                                                                                                                                                            cATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAGGGGAGTCAATATCACCGTCCAGCAGCCGtt                      <  1:167373/251‑1 (MQ=255)                                                                                                                                                                                                                                  cgccAGAAACGCCGCTTTTGCCCAGCTGGCTTTTTCTGCTTCCGCACGTGCCTGTCGGTGTTCTATCACCATTTACTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAACGGGGAGGCACTCTCACCGTCAAGCAGCCGTTGCAgcg                >  2:262477/1‑251 (MQ=255)                                                                                                                                                                                                                                      aGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAGGGGAGTCAATATCACCGTCCAGCAGCCGTTGCAGCGCCTg            <  2:318835/251‑1 (MQ=255)                                                                                                                                                                                                                                               gATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATGCAGCTGTACGCGCGTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAGGGGAGTCAATATCACCGTCCAGCAGCCGTTGCAGCGCCTGCGTTTGtgc   <  2:90093/251‑3 (MQ=255)                                                                                                                                                                                                                                                aTTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTCCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGCGGCGATTCAGTGCATGCACATTGCTGCGAAACGCCTCCCTCAGCCGCCCGATGGGATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAACGGGGAGCCAATATCACCGACCAGCAGCCCTTGAAGCGCCAGCGTTTGTTCg  >  2:101549/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  AGCGGGCGCGGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGC‑TCATCGCCGCCAGAAACGCCGATTTTGCCTGGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAGGGGAGTCAATATCACCGTCCAGCAGCCGTTGCAGCGCCTGCGTTTGTTCG  >  NC_000913/1054608‑1055096

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