Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,572,928	0	C	G	12.1%	175.9 / 15.2	58	A759P (GCA→CCA)	pqqL	putative periplasmic M16 family zinc metalloendopeptidase
Reads supporting (aligned to +/- strand):  ref base C (29/22);  new base G (6/1);  total (35/23)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.59e-01
Rejected: Not supported by required number of reads on both strands.

CTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCGGCAGTAAACTGTGCAATCTGATTTTCTTGCAGTAATTTCGTGCGGTCATCAGCATAGCGCGTCTCATACATCTGCTGGGCGA  >  NC_000913/1572644‑1573227                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                              cTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGt                                                                                                                                                                                                                                                                                               <  1:292709/301‑1 (MQ=255)   cACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGAGCCATACGCGTCGGCAGATTAAAAGGTGTCCGGGAATCATAACGCTTCCAATGTGACACCTGTGCCCCAGGTTCATTTTGttt                                                                                                                                                                                                                                                                                             >  1:395702/1‑300 (MQ=255)       gATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCtttt                                                                                                                                                                                                                                                                                         >  2:428039/1‑301 (MQ=255)         tCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTAc                                                                                                                                                                                                                                                                                       >  2:333462/1‑301 (MQ=255)                   cAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGAt                                                                                                                                                                                                                                                                             <  2:396216/301‑1 (MQ=255)                              ccATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGt                                                                                                                                                                                                                                                                  >  2:426679/1‑301 (MQ=255)                                                                          ggTTGACAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAgt                                                                                                                                                                                                                     >  1:471444/1‑301 (MQ=255)                                                                                cAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTttgt                                                                                                                                                                                                               >  1:281208/1‑300 (MQ=255)                                                                                cAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTttga                                                                                                                                                                                                               >  2:80736/1‑301 (MQ=255)                                                                                cAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCCAATGGCCAATCAGAGAGTttgt                                                                                                                                                                                                               >  1:233000/1‑300 (MQ=255)                                                                                   gTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTttgattt                                                                                                                                                                                                            >  1:473905/1‑300 (MQ=255)                                                                                   gTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTttgattg                                                                                                                                                                                                            >  2:246495/1‑301 (MQ=255)                                                                                                     cTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTACTGGCGACTCAGAGTGTTTGATTGATCCTAAGTACC‑GCGTaa                                                                                                                                                                                         >  1:470673/1‑301 (MQ=255)                                                                                                       tatatCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAAtt                                                                                                                                                                                       <  2:281511/300‑1 (MQ=255)                                                                                                                         ttCAACCGAGAGGCTAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGtct                                                                                                                                                                     <  2:535591/300‑1 (MQ=255)                                                                                                                              ccGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGc                                                                                                                                                                <  2:474506/301‑1 (MQ=255)                                                                                                                              ccGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCAGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGc                                                                                                                                                                <  2:436210/301‑1 (MQ=255)                                                                                                                                     gCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTa                                                                                                                                                         <  2:295688/301‑1 (MQ=255)                                                                                                                                           aaaaaCGCGTTATGCTCCAGTAGTCTGTTCACGAATTTTAACACGTTGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACATGTGTCCGGGAATCATAACGCTTCCACTTTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACa                                                                                                                                                   <  1:333057/299‑1 (MQ=255)                                                                                                                                                  cTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACaa                                                                                                                                            <  1:426115/301‑1 (MQ=255)                                                                                                                                                         cTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGCGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCAGCGACATTACCCACAATCAAAAAAGAGat                                                                                                                                     >  1:426041/1‑301 (MQ=255)                                                                                                                                                                 gCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCAATTGTAAAACTGTTGCCAAAGTTCCATTTTTTTCTTTTACGTTACCGGATGGTTTTTCGGTCG‑GCGGGTTAATGGGTTAACTTGCGCTAAATGGCAAACCAGAGTG‑TTGAATGATTCTTAATAACGGCGTATTAAAGGCCAGAAGTTGTTTTTTTGCAACATAACCACCAATaaaaaaggaaaaacacactg                                                                                                                             >  2:467730/1‑287 (MQ=255)                                                                                                                                                                  ccTGTTAACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCTAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTgg                                                                                                                            <  1:246228/301‑1 (MQ=255)                                                                                                                                                                     gTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGgaga                                                                                                                         >  2:345550/1‑301 (MQ=255)                                                                                                                                                                                        tAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGATATAACCGCTGGAGTAGAAAACAATTgccaaacg                                                                                                      >  1:313870/1‑294 (MQ=255)                                                                                                                                                                                                 tGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATCGGcagccagc                                                                                             <  2:233243/301‑1 (MQ=255)                                                                                                                                                                                                  gCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCATAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCAGGAACTGAAAACAATTGCACATCGGcagccaggg                                                                                            >  1:371002/1‑299 (MQ=255)                                                                                                                                                                                                          gACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGCTCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTTCGACATTACCGACAATGACAAACGTGATATCAGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCGGCg                                                                                    >  1:387859/1‑301 (MQ=255)                                                                                                                                                                                                          gACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCg                                                                                    <  2:471274/301‑1 (MQ=255)                                                                                                                                                                                                               tAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCAACTGCCGCAGt                                                                               >  1:306075/1‑301 (MQ=255)                                                                                                                                                                                                                     cATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCGGCAGTAAACtg                                                                         <  2:306479/301‑1 (MQ=255)                                                                                                                                                                                                                     aaTCGAGGGCCATACGCGCCGGCAGATTAACAGGTGTCCGGGAAGCATAACGCTGTCACTGTGAAACCTGGGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCGTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCGGCAGTAAACtg                                                                         <  1:583237/300‑1 (MQ=255)                                                                                                                                                                                                                              ccATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCGGCAGTAAACTGTGCAATCTg                                                                <  1:427475/301‑1 (MQ=255)                                                                                                                                                                                                                                  aCGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACCTTACCGACAATGACAAACGTGATATCCGCTGGAGATGACAACAATTGGCGAACGGCAGCCAACGCATCTGCGGCAGTAAACTGTGCAATATGAttt                                                            >  1:556381/1‑301 (MQ=255)                                                                                                                                                                                                                                      tttgctcttccgttctttgtCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCGGCAGTAAACTGTGCAATCTGATTTTCtt                                                        <  2:608806/284‑1 (MQ=255)                                                                                                                                                                                                                                        cGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATCGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCGGCAGTAAACTGTGCAATCTGATTTTCTTGc                                                      >  2:454971/1‑301 (MQ=255)                                                                                                                                                                                                                                           cAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCGGCAGTAAACTGTGCAATCTGATTTTCTTGCAga                                                   >  2:362497/1‑300 (MQ=255)                                                                                                                                                                                                                                                  aCAGGTTTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCGGCAGTAAACTGTGCAATCTGATTTTCTTGCAGTAATTTCg                                            <  2:314274/301‑1 (MQ=255)                                                                                                                                                                                                                                                    ccccttcccgacgctcttccgttctcACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAGCGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCACTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCGGCAGTAAACTGTGCAATCTGATTTTCTTGCAGTAATTTCGTg                                          <  1:409746/276‑1 (MQ=255)                                                                                                                                                                                                                                                      gtgtCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCGGCAGTAAACTGTGCAATCTGATTTTCTTGCAGTAATTTCGTGCg                                        <  1:345087/301‑1 (MQ=255)                                                                                                                                                                                                                                                      gtgtCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCAACGAGTTTGTCTTCTGCGACATTACCGACAATGCACAACGTGATATCCGCTGGAGATGAAAACAATTGGCGCTCGGCAGCCAACGCATCTGCGGCAGTAAACTGTCCAATCTGATTTTAttagaccggcagcacggc                                        >  1:607948/1‑285 (MQ=255)                                                                                                                                                                                                                                                       cacagtggaagcgttatgattCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCGGCAGTAAACTGTGCAATCTGATTTCCTTGCAATAATTTCGTCCgg                                       >  1:281510/20‑301 (MQ=255)                                                                                                                                                                                                                                                            gggAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCGGCAGTAAACTGTGCAATCTGATTTTCTTGCAGTAATTTCGTGCGGtcata                                  >  2:130015/1‑300 (MQ=255)                                                                                                                                                                                                                                                              ggttactgtaaaagaacaaaatgaaCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCGGCAGTAAACTGTGCAATCTGATTTTCTTGCAGTAATTTCGTGCGGtcatcaa                                >  2:228857/24‑300 (MQ=255)                                                                                                                                                                                                                                                                   ctgtaaaagaacaaaatgaaCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGACACCAATTGGCGATCGGCAGCCAGCGCATCTGCGGCAGTAAACTGTGCAATCTGATTTTCTTACAGTAATTTCGTACGGTCATCAGCATat                           >  1:641902/19‑300 (MQ=255)                                                                                                                                                                                                                                                                    tAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCGGCAGTAAACTGTGCAATCAGATTTTCTTGAAGTAATTTCGTGCGGTAATCAGCATAgg                          >  1:126677/1‑300 (MQ=255)                                                                                                                                                                                                                                                                         cTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCGGCAGTACACTGTGCAATCTGATTTTCTTGCAGTAATTTAGTGCGGTCATCAGCATAGCGCGtc                     >  1:561701/1‑301 (MQ=255)                                                                                                                                                                                                                                                                           gaacaaaatgaaCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCGGCAGTAAACTGTGCAATCTGATTTTCTTGCAGTAATTTCGTGCGGTCATCAGCATAGCGCGtctc                   >  2:41649/11‑301 (MQ=255)                                                                                                                                                                                                                                                                             cACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAGCGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCACTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCGGCAGTAAACTGTGCAATCTGATTTTCTTGCAGTAATTTCGTgagatcggaagagcacaagtctgaaa                 >  2:410309/1‑276 (MQ=255)                                                                                                                                                                                                                                                                               aaaatgaaCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGCATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAACCGCGATATCCGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCGGCCGTAACCTGTGCAATCTGATTTTCTTGCAGTAATTTCGTGCGGTCATCACCATACCTCGTCTCCTAc               >  1:365968/7‑301 (MQ=255)                                                                                                                                                                                                                                                                               aaaatgaaCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCGGCAGTAAACTGTGCAATCTGATTTTCTTGCAGTAATTTCGTACGGTCATCAGCATAGCCCATCCCATAc               >  1:271625/7‑301 (MQ=255)                                                                                                                                                                                                                                                                                aaatgaaCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCGGCAGTAAACTGTGCAATCTGATTTTCTTGCAGTAATTTCGTGCGGTCATCAGCATAGCGCGTCTCATACa              >  2:395635/6‑301 (MQ=255)                                                                                                                                                                                                                                                                                   gAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGGCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCGGCAGTAAACTGTGCAATCTGATTTTCTTGCAGTAATTTCGTGCGGTCATCAGCATAGCGCGTCTCATACATct           <  2:562415/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                    gaaCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCGGCAGTAAACTGTGCAATCTGATTTTCTTGCAGTAATTTCGTGCGGTCATCAGCATAGCGCGTCTCATACATctg          <  1:41596/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                         tgtgCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCGGCAGTAAACTGTGCAATCTGATTTTCTTGCAGTAATTTCGTGCGGTCATCAGCATAGCGCGTCTCATACATCTGCTggg     >  2:460106/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                          ttgCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCGGCAGTAAACTGTGCAATCTGATTTTCTTGCAGTAATTTCGTGCGGTCATCAGCATAGCGCGTCTCATACATCTGCTGGGc    <  2:281813/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                           tgCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCGGCAGTAAACTGTGCAATCTGATTTTCTTGCAGTAATTTCGTGCGGTCATCAGCATAGCGCGTCTCATACATCTGCTGGGCg   <  2:388373/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                            gCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTTCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCGGCAGTAAACTGTGCAATCTGATTTTCTTGCAGTAATTTCGTGCGGTCATCAGCATAGCGCGTCTCATACATCTGCTGGGCGa  >  1:622574/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                              CTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGGTTGACAAGTAAAAGCCAGCAAATGACTGATATCTTTGGCCTGAGGATCAACCGAGAGGCGAGAAGAAACGCTGTATGCTCCAGATGCCTGTTCACGAATATTAACACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAAT‑GGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAAC‑GCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCGCTGGAGATGAAAACAATTGGCGATCGGCAGCCAGCGCATCTGCGGCAGTAAACTGTGCAATCTGATTTTCTTGCAGTAATTTCGTGCGGTCATCAGCATAGCGCGTCTCATACATCTGCTGGGCGA  >  NC_000913/1572644‑1573227

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑