Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	154,818:1	+G	10.0%	coding (609/2598 nt)	htrE ←	putative outer membrane usher protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	154,817	1	.	G	10.0%	146.2 / 13.4	40	coding (610/2598 nt)	htrE	putative outer membrane usher protein
Reads supporting (aligned to +/- strand):  ref base . (17/19);  new base G (2/2);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

TAGTGACTTTGGCGTTGGTATTGGCAACGCCATGAATGATAGGCGCAAAGCTGGCTAAAGTCGGAGGCAACATGCGGCTGTCACTGTATAAACGAATGCCTCGGATACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGTTGTATGACAACATGGCCGCATTAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTCATAACGAACAGAAGCCTGAGGGATAATTTCCGTCAAATTGAGGCAATTGCCGAGCGTTTCATCCCTGGCAAGCAGAACGGCTTTTT  >  NC_000913/154532‑155089                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                 tAGTGACTTTGGCGTTGGTATTGGCAACGCCATGAATGATAGGCGCAAAGCTGGCTAAAGTCGGAGGCAACATGCGGCTGTCACTGTATAAACGAATGCCTCGGATACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑Gtt                                                                                                                                                                                                                                                                               <  1:91056/289‑1 (MQ=255)      aCTTTGGCGTTGGTATTGGCAACGCCATGAATGATAGGCGCAAAGCTGGCTAAAGTCGGAGGCAACATGCGGCTGTCACTGTATAAACGAATGCCTCGGATACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAAt                                                                                                                                                                                                                                                                           >  2:580244/1‑288 (MQ=255)      aCTTTGGCGTTGGTATTGGCAACGCCATGAATGATAGGCGCAAAGCTGGCTAAAGTCGGAGGCAACATGCGGCTGTCACTGTATAAACGAATGCCTCGGATACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAAt                                                                                                                                                                                                                                                                           >  1:82199/1‑288 (MQ=255)       cTTTGGCGTTGGTAGTGGCAACGCCATGAATGATAGGGTCAAAGCTGGCTAAAGTCGGAGGCAACATGCGGCTGTCACTGTATAAACGAATGCCTCGGATACTGACGGAATCAAAGGGTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATg                                                                                                                                                                                                                                                                          <  1:436679/288‑1 (MQ=255)           ggCGTTGGTATTGGCAACGCCATGAATGATAGGCGCAAAGCTGGCTAAAGTCGGAGGCAACATGCGGCTGTCACTGTATAAACGAATGCCTCGGATACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATgcagca                                                                                                                                                                                                                                                                     >  2:248456/1‑289 (MQ=255)                             gCCATGAATGATAGGCGCAAAGCTGGCTAAAGTCGGAGGCAACATGCGGCTGTCACTGTATAAACGAATGCCTCGGATACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCAtttttt                                                                                                                                                                                                                                                   <  1:248456/289‑1 (MQ=255)                                  gaatgaTAGGCGCAAAGCTGGCTAAAGTCGGAGGCAACATGCGGCTGTCACTGTATAAACGAATGCCTCGGATACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGAcc                                                                                                                                                                                                                                              <  2:553912/289‑1 (MQ=255)                                            cgcAAAGCTGGCTAAAGTCGAGGGCAACATGCGGCTGTCACTGTTTAAACGAATGCCTCGGATACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTc                                                                                                                                                                                                                                     <  1:435624/288‑1 (MQ=255)                                                       cTAAAGTCGGAGGCAACATGCGGCTGTCACTGTATAAACGAATGCCTCGGATACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTTATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATc                                                                                                                                                                                                                          <  2:11744/288‑1 (MQ=255)                                                              cGGAGGCAACATGCGGCTGTCACTGTATAAACGAATGCCTCGGATACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAgg                                                                                                                                                                                                                  <  2:224885/289‑1 (MQ=255)                                                                gAGGCAACATGCGGCTGTCACTGTATAAACGAATGCCTCGGATACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGtt                                                                                                                                                                                                                >  2:38035/1‑289 (MQ=255)                                                                                                 tGCCTCGGATACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACCGGTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGTTGTATGACAACATGGCCGCATTAATGCCGtttt                                                                                                                                                                                >  1:126272/1‑289 (MQ=255)                                                                                                 tGCCTCGGATACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACCGGTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGTTGTATGACAACATGGCCGCATTAATGCCGttt                                                                                                                                                                                 >  1:381200/1‑288 (MQ=255)                                                                                                 tGCCTCGGATACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCa                                                                                                                                                                                                                                    <  2:416178/236‑1 (MQ=255)                                                                                                 tGCCTCGGATACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCa                                                                                                                                                                                                                                    >  1:416178/1‑236 (MQ=255)                                                                                                 tGCCTCGGATACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGTTGTATGACAACATGGCCGCATTAATGCCGtttt                                                                                                                                                                                >  1:214759/1‑288 (MQ=255)                                                                                                          tACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGTTGTATGACAACATGGCCGCATTAATGCCGTTTTCCCATAACGa                                                                                                                                                                      <  1:580244/289‑1 (MQ=255)                                                                                                                cGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGTTGTATGACAACATGGCCGCATTAATGCCGTTTTCCCATAACGATGGAt                                                                                                                                                                 <  2:81418/288‑1 (MQ=255)                                                                                                                     tCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGTTGTATGACAACATGGCCGCATTAATGCCGTTTTCCCATAACGATGGATCAACAt                                                                                                                                                           <  1:38035/289‑1 (MQ=255)                                                                                                                                                tCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGtt                                                                                                                                                                                                                >  1:61595/1‑209 (MQ=255)                                                                                                                                                tCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGtt                                                                                                                                                                                                                <  2:61595/209‑1 (MQ=255)                                                                                                                                                tCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGTTGTATGACAACATGGCCGCATTAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAg                                                                                                                                 >  2:219039/1‑288 (MQ=255)                                                                                                                                                  aCCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTAt                                                                                                                                                                                                                                <  1:277788/191‑1 (MQ=255)                                                                                                                                                  aCCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTAt                                                                                                                                                                                                                                >  2:277788/1‑191 (MQ=255)                                                                                                                                                          tGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACCGGTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGTTGTATGACAACATGGCCGCATTAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGaa                                                                                                                       <  2:126272/289‑1 (MQ=255)                                                                                                                                                          tGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACCGGTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGTTGTATGACAACATGGCCGCATTAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGaa                                                                                                                       <  2:381200/289‑1 (MQ=255)                                                                                                                                                                  gaACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGTTGTATGACAACATGGCCGCATTAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCtata                                                                                                               >  2:440918/1‑288 (MQ=255)                                                                                                                                                                                   atCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGTTGTATGACAACATGGCCGCATTAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGTTa                                                                                             <  2:82199/289‑1 (MQ=255)                                                                                                                                                                                      aCGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGTTGTATGACAACATGGCCGCATTAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACg                                                                                          <  2:214759/289‑1 (MQ=255)                                                                                                                                                                                       cGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGTTGTATGACAACATGGCCGCATTAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                                                                                                        <  2:245771/258‑1 (MQ=255)                                                                                                                                                                                       cGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGTTGTATGACAACATGGCCGCATTAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                                                                                                        >  1:245771/1‑258 (MQ=255)                                                                                                                                                                                            aaCATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGTTGTATGACAACATGGCCGCATTAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGaa                                                                                                                       >  2:560362/1‑254 (MQ=255)                                                                                                                                                                                            aaCATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGTTGTATGACAACATGGCCGCATTAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGaa                                                                                                                       <  1:560362/254‑1 (MQ=255)                                                                                                                                                                                                    gATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGTTGTATGACAACATGGCCGCATTAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTGGATCGTTAACGTCATAACGAACAGa                                                                            >  2:77756/1‑289 (MQ=255)                                                                                                                                                                                                      ttCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGTTGTATGACAACATGGCCGCATTAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTc                                                                                        >  2:550967/1‑275 (MQ=255)                                                                                                                                                                                                      ttCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGTTGTATGACAACATGGCCGCATTAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTc                                                                                        <  1:550967/275‑1 (MQ=255)                                                                                                                                                                                                                                               gTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGTTGTATGACAACATGGCCGCATTAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGaa                                                                                                                       <  2:547186/203‑1 (MQ=255)                                                                                                                                                                                                                                               gTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGTTGTATGACAACATGGCCGCATTAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGaa                                                                                                                       >  1:547186/1‑203 (MQ=255)                                                                                                                                                                                                                                               gTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGTTGTATGACAACATGGCCGCATTAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTCATAACGAACAGAAGCCTGAGGGATAATTTCCGTCAAATTGAGGCAATTGCCGAGc                                 >  2:115332/1‑289 (MQ=255)                                                                                                                                                                                                                                                                              ccTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGTTGTATGACAACATGGCCGCATTAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTCATAACGAACAGAAGCCTGAGGGATAATTTCCGTCAAATTGAGGCAATTGCCGAGCGTTTCATCCCTGGCAAGCAGAACGGCttttt  >  2:322443/1‑289 (MQ=255)                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                 TAGTGACTTTGGCGTTGGTATTGGCAACGCCATGAATGATAGGCGCAAAGCTGGCTAAAGTCGGAGGCAACATGCGGCTGTCACTGTATAAACGAATGCCTCGGATACTGACGGAATCAAAGGTTTCGCCCGTCGTATAAGACTCACCAAGAATGAGTTGAGAACGCAGCGAGGCGATATCACGCTGAACATACCGATTCTTAAAATCATAATTACTGCCAGAATCGGTCATCCAGTTGTAGTTGCCCGAGGCACGCAGTCGCCATGCACCTAAATTCATCCCACC‑GTTAAATGCAGCATAAATGCTTTCATTTTTTCGACCAGGGGTTTCACTATGATATCCGTTGAGGTTGTATGACAACATGGCCGCATTAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTCATAACGAACAGAAGCCTGAGGGATAATTTCCGTCAAATTGAGGCAATTGCCGAGCGTTTCATCCCTGGCAAGCAGAACGGCTTTTT  >  NC_000913/154532‑155089

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