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	95.3%	186.4 / ‑3.5	65	Q442* (CAG→TAG)	torS	sensor histidine kinase TorS
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (36/25);  total (39/26)
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

GGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTGGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAGGGGAGTCAATATCACCGTCCAGCAGCCGTTGCAGCGCCTGCGTTTG  >  NC_000913/1054617‑1055092                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              ggTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCtttttt                                                                                                                                                                                                                                   >  1:450287/1‑250 (MQ=255) ggTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTc                                                                                                                                                                                                                                   >  1:86857/1‑251 (MQ=255)     cAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAGTTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGACTCCTCGACGCCAGAAACGCACATTTAGCCTGGCTGGCTTTATCTgtt                                                                                                                                                                                                                               >  2:106720/1‑249 (MQ=255)         gggCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAACTCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCg                                                                                                                                                                                                                           >  1:239940/1‑251 (MQ=255)             tCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCg                                                                                                                                                                                                                       >  2:166406/1‑251 (MQ=255)              cATCGCTGACCGAAACATTCTTGCCCCCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCGGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATTGCTCACCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGccc                                                                                                                                                                                                                        <  2:368492/249‑1 (MQ=255)              cATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGccc                                                                                                                                                                                                                        <  1:81226/249‑1 (MQ=255)                                    gcccccTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCa                                                                                                                                                                                                <  2:456437/247‑1 (MQ=255)                                    gCCACCTGCTTCGATAGCGGAATAATCGAGAATATCTTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGTGCTGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCTCCAGAAACGCCGATTTTGGCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCAcc                                                                                                                                                                                                 <  1:26264/250‑1 (MQ=255)                                      cacatGCTTCGATAGCGGCATAATCGAGAATATCGTTGAGGATGGTCAGCACCGGTTCACCCGCGTCAGTAATTGCACGCAACTCATCCAGCTGGGCGTTACGTGCGGGGTTATCTGCCAGCAGTTGAGAAGTGCCGCGACTACCGTACAGCGGTGTGCGGATCTAATGGCTCATCGCCGCCAGAAACGCCGGTTTTGCCTAGCTTGCTTTTTCTGCTTCCGCCCATGCCTGGCGTTGTTCttgccccagt                                                                                                                                                                                              >  2:397980/5‑242 (MQ=255)                                         cTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCttttt                                                                                                                                                                                                                                    <  2:154590/210‑1 (MQ=255)                                         cTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCttttt                                                                                                                                                                                                                                    >  1:154590/1‑210 (MQ=255)                                           gCTTCGCTAGCGGAATAATCGAGAATATCGTTGAGCATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGAAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCCGTGCCGAGAATACCGTACAGCGGTGTTCGGATCTCATGGCTCATCGCCGCCAGACACCCCGCTCTTGCCTCGCTGGCTTTTTCTGCGTCCCCCAGTGCCTGGCGGGGTTTCATCACAAGTTCCTg                                                                                                                                                                                         >  2:128728/1‑251 (MQ=255)                                                        aataatCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACAAGTTCCTGCAATTCAGCTgac                                                                                                                                                                            >  2:467397/1‑249 (MQ=255)                                                         ataatCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTAc                                                                                                                                                                           >  1:440076/1‑251 (MQ=255)                                                           aatCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTAcg                                                                                                                                                                          <  1:467397/250‑1 (MQ=255)                                                                  aatttCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGa                                                                                                                                                                  <  2:440076/247‑1 (MQ=255)                                                                         ttGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTGGCTGGCTTTTTCTGCTTCCGCCCGGGCCTGTCGGTGTTCTATCACCAGTTCCTGCATTTCAGCTGTACGCGCTTTGACCTgccc                                                                                                                                                           >  2:329689/1‑249 (MQ=255)                                                                           gAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGcgccg                                                                                                                                                         <  1:223949/251‑1 (MQ=255)                                                                            aGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGTTCTCATGGCTTATCGCCGCCAGAAACGCCGCTTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGcgccgc                                                                                                                                                        >  1:196051/1‑251 (MQ=255)                                                                            aGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGcgccgc                                                                                                                                                        >  2:259978/1‑251 (MQ=255)                                                                             ggATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGCATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCAGCCAGAAACGCCGATTATGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTCCGCGCTTTGCCCTGcgccacc                                                                                                                                                       >  2:381606/1‑248 (MQ=255)                                                                                     aGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCAGCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTc                                                                                                                                               <  2:101239/251‑1 (MQ=255)                                                                                            aTTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTgg                                                                                                                                        <  2:5248/251‑1 (MQ=255)                                                                                            aTTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCACCAGCTGTTCGCGGTgg                                                                                                                                        >  2:11052/1‑251 (MQ=255)                                                                                                 ccAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGAtt                                                                                                                                   >  1:388555/1‑251 (MQ=255)                                                                                                   agagTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGTCCTGCGCCGCCAGCTGTTCGCGGTGGCGTTTc                                                                                                                                  >  2:156009/1‑250 (MQ=255)                                                                                                         aGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCaa                                                                                                                                                                                      >  1:369715/1‑192 (MQ=255)                                                                                                         aGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCaa                                                                                                                                                                                      <  2:369715/192‑1 (MQ=255)                                                                                                         aGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCAt                                                                                                                           >  1:471546/1‑251 (MQ=255)                                                                                                                 cccGCAAATCATCCCGCTGGGCGTTAAGTGCGGGGTTATCGGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTGCGCGGTGGCGATTCAGTGCATGAACATtg                                                                                                                   <  1:288136/251‑1 (MQ=255)                                                                                                                   cGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATtgct                                                                                                                 <  1:23813/251‑1 (MQ=255)                                                                                                                    ccAAATCACCACGCTGGTCGCTAAGTGAGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGGACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTCCCTAGCTGGCTGTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATtgctg                                                                                                                <  2:334224/250‑1 (MQ=255)                                                                                                                        atcatcACGCTGGGCGTTACGTGCGGGGTTCTCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGa                                                                                                             >  2:58962/1‑250 (MQ=255)                                                                                                                        atcatcACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGATGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGa                                                                                                             >  1:54904/1‑250 (MQ=255)                                                                                                                             caCGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCTCCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGca                                                                                                       <  1:222420/251‑1 (MQ=255)                                                                                                                             caCGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGca                                                                                                       <  1:470636/251‑1 (MQ=255)                                                                                                                                    ggCGTTAAGGGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTGGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCa                                                                                                <  2:141548/251‑1 (MQ=255)                                                                                                                                      cGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAgc                                                                                              >  1:350002/1‑251 (MQ=255)                                                                                                                                          aagcgcggtgtttacTGCCAGCAGTTGATGTGTGCCGTGAAAACCGTACAGCGGTTGGCGGATCTCATGGCTAAGCTACGCCAGATACGCCGAGTTTGCCTAGCTAGCTTTTGCTGCTTTCGCCCGTGCCTGTCGGTGTTCTATCTCCAGTTCCTGCAATTCAGCGGTACGCGCTTTGACCTGCGCCGCCAGCTGTGCGCGGTGCCGATTCAGTGCATGAACATTGCTTCGAAACGCATCCATCAtccgcc                                                                                          <  2:196051/237‑1 (MQ=255)                                                                                                                                           agagCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCAAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAAAATTGCTGCGAAACGCATCCATCAGCCGccc                                                                                         >  2:352959/4‑251 (MQ=255)                                                                                                                                              gcgggtgtATCTGCCAGCAGGTTATCAGTGCCGAGAATACCGTACACCTGTGTTCGTATCTCATTACTCATCGCCACCAGAACCGCCGATTTTGCCTATCTGGCTTTTTCTTCTGCCGCCCTTGCCTGTCTGTGTTCTATCCTCAGTTCCTTCAATGCAGCTGTACGCTCTTTTACCTGCGCCGCCAGCTGTTCGCGGTGGCGAGTCAGTTCACGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGAt                                                                                      <  1:133018/244‑1 (MQ=255)                                                                                                                                              gCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGAt                                                                                      >  2:278598/1‑251 (MQ=255)                                                                                                                                                   gTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACa                                                                                                                      >  1:163316/1‑214 (MQ=255)                                                                                                                                                   gTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGGGTGCGGATCGCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACa                                                                                                                      <  2:163316/214‑1 (MQ=255)                                                                                                                                                          gCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTCTCCAGCTc                                                                          >  1:428963/1‑251 (MQ=255)                                                                                                                                                           ccAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTcc                                                                         >  1:438150/1‑251 (MQ=255)                                                                                                                                                                agTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTAc                                                                    <  1:278598/251‑1 (MQ=255)                                                                                                                                                                 gTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACg                                                                   >  1:32887/1‑251 (MQ=255)                                                                                                                                                                      gCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGc                                                              >  1:439567/1‑251 (MQ=255)                                                                                                                                                                           gCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGtt                                                          <  1:146058/250‑1 (MQ=255)                                                                                                                                                                                  aTACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGa                                                   >  1:115254/1‑250 (MQ=255)                                                                                                                                                                                           aGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAGGGGAGTc                                         >  1:175943/1‑251 (MQ=255)                                                                                                                                                                                               gtgtGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAGGGGAGTCAat                                      >  1:306933/1‑250 (MQ=255)                                                                                                                                                                                                        tctcATGGATCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCAAGCGGTTTCCGGGAAAGGGGAGTCAATATCACCGTcc                            >  1:66025/1‑251 (MQ=255)                                                                                                                                                                                                           cATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAGGGGAGTCAATATCACCGTCca                           <  1:421530/249‑1 (MQ=255)                                                                                                                                                                                                            aTGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGCTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAGGGGAGTCAATATCACCGTCcagc                         >  1:204088/1‑250 (MQ=255)                                                                                                                                                                                                               gCTCATCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAGGGGAGTCAATATCACCGTCCAGCAGcc                     >  1:305028/1‑251 (MQ=255)                                                                                                                                                                                                                   aTCGCCGCCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAGGGGAGTCAATATCACCGTCCAGCAGCCGTTg                 <  1:156009/251‑1 (MQ=255)                                                                                                                                                                                                                    ccgctgccAGAAACGCCGATTTTGCCTAGCTGGCTGTTTTTGATTCCGCCCGTGCTGGTGGGTGTCATATCACCAGTTTCTGCAATTACGCTGCACGCGCTTTGACCTGCGCCGCCAGCGGCTCGCGGTGGAGATTCAGTGCATGAACATCGCTGCGAAACGCCCCCATCAGCCGCCCGATGGTATCCTGCTCCCGGACGCTAGCGGTGGCCGGGAAAGGGGAGTCAATATCACCGTCCAGCAGTCGTTGc                <  2:371888/246‑1 (MQ=255)                                                                                                                                                                                                                     cgccgcCAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCATCTCCCGTACGCCAGCGGTTTCCGGGAAAGGGGAGTCAATATCACCGTCCAGCAGCCGTTg                 <  1:381362/249‑1 (MQ=255)                                                                                                                                                                                                                      gccgccAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAGGGGCGTCAATATCACCGTCCAGCAGCCGTTGc                >  2:99454/1‑249 (MQ=255)                                                                                                                                                                                                                        cgccAGAAACGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGACAGGGGAGTCAATATCACCGTCCAGCAGCCGTTGCAgcg            >  1:268715/1‑251 (MQ=255)                                                                                                                                                                                                                               aaCGCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTAACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAGGGGAGTCACTATCACCGTCCAGCAGCAGTTGCAGCGCCTGCGt     >  2:289257/1‑251 (MQ=255)                                                                                                                                                                                                                                  gCCGATTTTGCCTAGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAGGGGAGTCAATATCACCGTCCAGCAGCCGTTGCAGCGCCTGCGTTTg  >  1:103416/1‑251 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              GGTTCAAAGGGCTCATCGCTGACCGAAACATTCTTGCCACCTGCTTCGATAGCGGAATAATCGAGAATATCGTTGAGGATGGTCAGCAACGATTCACCAGAGTCAGTAATTGCCCGCAAATCATCACGCTGGGCGTTAAGTGCGGGGTTATCTGCCAGCAGTTGAGCAGTGCCGAGAATACCGTACAGCGGTGTGCGGATCTCATGGCTCATCGCCGCCAGAAACGCCGATTTTGCCTGGCTGGCTTTTTCTGCTTCCGCCCGTGCCTGTCGGTGTTCTATCACCAGTTCCTGCAATTCAGCTGTACGCGCTTTGACCTGCGCCGCCAGCTGTTCGCGGTGGCGATTCAGTGCATGAACATTGCTGCGAAACGCATCCATCAGCCGCCCGATGGTATCCAGCTCCCGTACGCCAGCGGTTTCCGGGAAAGGGGAGTCAATATCACCGTCCAGCAGCCGTTGCAGCGCCTGCGTTTG  >  NC_000913/1054617‑1055092

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