Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	292,019	C→T	L174L (CTA→TTA)	yagJ →	protein YagJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	292,019	0	C	T	96.3%	156.4 / ‑2.2	54	L174L (CTA→TTA)	yagJ	protein YagJ
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (29/23);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

CCATTTGTTGTGTGACGAAAGCGGGAACTATCCGGCTCTGTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGC  >  NC_000913/291772‑292259                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 ccATTTGTTGTGTGACGAAAGCGGGAACTATCCGGCTCTGTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAc                                                                                                                                                                                                                                               >  1:354507/1‑251 (MQ=255)      tgttgtGTGACGAAAGCGGGAACTATCCGGCTCTGTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCt                                                                                                                                                                                                                                           <  1:227143/250‑1 (MQ=255)       gttgtGTGACGAAAGCGGGAACTATCCGGCTCTGTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTg                                                                                                                                                                                                                                          <  1:308625/250‑1 (MQ=255)                 gAAAGCGGGAACTATCCGGCTCTGTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAGACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCttt                                                                                                                                                                                                                               >  1:149204/1‑251 (MQ=255)                      cGGGAACTATCCGGCTCTGTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGccc                                                                                                                                                                                                                           >  1:45273/1‑250 (MQ=255)                       gggAACTATCCGGCTCTGTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCa                                                                                                                                                                                                                          >  1:632839/1‑250 (MQ=255)                       gggAACTATCCGGCTCTGTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGAGTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCttt                                                                                                                                                                                                                               <  1:239439/245‑1 (MQ=255)                       gggAACTATCCGGCTCTGTTGAACCATTGGGAAACCAAGGCTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGAATCGTAATACACAGTACCCAGGGCAGTCGTAACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGCTTTCCTGTTCTTTGCCGAACAGGATGGCCAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCttt                                                                                                                                                                                                                               >  2:239439/1‑245 (MQ=255)                          aaCTATCCGGCTCTGTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACt                                                                                                                                                                                                                      >  2:420984/1‑251 (MQ=255)                                     tctGTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGACGGCTTTGCTTTCTGGTATCGTAATCCACCGTACACAGGGCAGTAGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCAAGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGCTCCACATAGCTTACATCAGGCTGATGCTTTGCCCAAACTGGAAGGACTTg                                                                                                                                                                                                           >  2:591900/1‑251 (MQ=255)                                        gTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAGACCAAACGGGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTg                                                                                                                                                                                                           <  2:149204/248‑1 (MQ=255)                                                      aaaCCAAGTTTTGTGAGGTCTAAACCAAACGTGATGGATTTGCTGTCTGGTATCGTAAGCCACAGGACACATGGCAGTCTTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTACCGATTTCTTGTTCTTTGCCGATCAGGATGTTAATATGGTTGTGGACCTGGTAGATCCACATAGTTTACATCTGGCTGATGCTTTGCCTAAACGTGATGGACTTGCGCTATATGTCGaatat                                                                                                                                                                                          <  2:493527/251‑4 (MQ=255)                                                         ccAAGGATTTTGCGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTAGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGCCGGCCTTGCGCTATATGCCGAAcatcat                                                                                                                                                                                       >  2:399678/1‑251 (MQ=255)                                                          cAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTATGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCAtt                                                                                                                                                                                      >  2:592527/1‑251 (MQ=255)                                                           aaGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACGGGAAGGACTTGCGCTATATGCCGAACATCATTc                                                                                                                                                                                     <  1:362373/251‑1 (MQ=255)                                                            aGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCt                                                                                                                                                                                    <  2:371555/251‑1 (MQ=255)                                                            aGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCAACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCt                                                                                                                                                                                    >  1:323721/1‑251 (MQ=255)                                                                           aaaCCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATAGCAAAATGGTTGTGGACATGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGcg                                                                                                                                                                      >  1:623886/1‑250 (MQ=255)                                                                                        aGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGACGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACCGGCGAATCGAATCTGTCg                                                                                                                                                        >  2:235935/1‑251 (MQ=255)                                                                                                      ggTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAg                                                                                                                                             <  1:158705/248‑1 (MQ=255)                                                                                                           cGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGAt                                                                                                                                                                                 <  2:133623/207‑1 (MQ=255)                                                                                                           cGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTCACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGAt                                                                                                                                                                                 >  1:133623/1‑207 (MQ=255)                                                                                                                            aGGGCAGGCGTCACTGGGACTAGCTTATGTTGCAGCTGAAAAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCCCATCGCTTACATCTGGCTGATGCTTTGCCAAAACTGGAAGGACTTGCGCTATATGCCGAACATAATTCTGATGCTTACAGGCGAATCGAATATGTCGCCGAAGTAAAGGGTAATTTAAGAGTGTTAAATTTg                                                                                                                     >  2:663866/1‑250 (MQ=255)                                                                                                                             gggCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTTGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAGTCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGaa                                                                                                                   <  1:445668/251‑1 (MQ=255)                                                                                                                                    cgtcACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTTGATTTGAAACGGc                                                                                                              >  2:174707/1‑249 (MQ=255)                                                                                                                                       cACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGAt                                                                                                         >  1:307873/1‑251 (MQ=255)                                                                                                                                        aCTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGAtg                                                                                                        >  1:15148/1‑251 (MQ=255)                                                                                                                                             gAATTGCTTATGTTGAAGCTGAATAGGACAAGATTGTTCGTCCCGATGTCCTTTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGTAAGGACTTGCGCTATATTCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACTAGTGTTAGATTTGAAACGGCATGATGTGCAg                                                                                                   <  2:446977/251‑1 (MQ=255)                                                                                                                                                 cGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGTATGTGCAGGAt                                                                                                <  2:632839/250‑1 (MQ=255)                                                                                                                                                 cGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGAAGGATg                                                                                               >  1:612816/1‑251 (MQ=255)                                                                                                                                                  gCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGc                                                                                              >  2:243416/1‑251 (MQ=255)                                                                                                                                                      aTGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGtt                                                                                          <  1:350834/251‑1 (MQ=255)                                                                                                                                                      aTGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGt                                                                                           <  1:243416/250‑1 (MQ=255)                                                                                                                                                       tGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTAAAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGGGTGTTAGATTTGAAACGGCAGGCTGTGCAGGATGCTGttt                                                                                         >  2:350599/1‑250 (MQ=255)                                                                                                                                                        gTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTGAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGc                                                                                        <  2:344645/251‑1 (MQ=255)                                                                                                                                                           gAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTa                                                                                      <  1:165483/250‑1 (MQ=255)                                                                                                                                                                 gAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTg                                                                                <  1:399678/250‑1 (MQ=255)                                                                                                                                                                    cAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGaaa                                                                             >  1:166419/1‑250 (MQ=255)                                                                                                                                                                      gTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAAATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGCGTGTTAGATTTGCAACGGCAGGATGTGCAGGATGCGGTTGCTACCGCTGAACATg                                                                          >  2:509162/1‑251 (MQ=255)                                                                                                                                                                       tACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGAtt                                                                                                                       <  2:237556/205‑1 (MQ=255)                                                                                                                                                                       tACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGAtt                                                                                                                       >  1:237556/1‑205 (MQ=255)                                                                                                                                                                                gTTCGTCCCGATTGCCTGTTCTTTGCCGAACAGGATGGCAAAATGGATGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGCGTGTTAGATTTGAAACGGCAGGATGAGCAGGATGCTGTTGCTACCGCTGCAAATGCAGAAACGtt                                                                >  1:610069/1‑251 (MQ=255)                                                                                                                                                                                    gTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGCAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGATAttt                                                            >  2:249617/1‑251 (MQ=255)                                                                                                                                                                                              ccTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGAc                                                  >  1:517023/1‑251 (MQ=255)                                                                                                                                                                                              ccTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGAc                                                  <  1:70264/251‑1 (MQ=255)                                                                                                                                                                                                       ggCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCtgatg                                         <  2:37913/250‑1 (MQ=255)                                                                                                                                                                                                                 gATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAgg                                >  2:312736/1‑249 (MQ=255)                                                                                                                                                                                                                     gCAAAATGGTTGTGGACCTGGTAGAGCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGAGTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCt                           <  1:591900/251‑1 (MQ=255)                                                                                                                                                                                                                       aaaaTGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCtat                         >  1:641387/1‑251 (MQ=255)                                                                                                                                                                                                                         aaTGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCtt                                                                                                                                                                             <  1:641510/101‑1 (MQ=255)                                                                                                                                                                                                                         aaTGGTTGTGGACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCtt                                                                                                                                                                             >  2:641510/1‑101 (MQ=255)                                                                                                                                                                                                                                   gACCTGGTAGGACCGCATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGt             <  1:663866/251‑1 (MQ=255)                                                                                                                                                                                                                                   gACCTGGTAGATCCACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGCCTTGTTGATGACTATCAGTAATCTATAGAAAATTgagt             >  2:649998/1‑248 (MQ=255)                                                                                                                                                                                                                                               ccACATAGCTTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGc  <  1:617914/250‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 CCATTTGTTGTGTGACGAAAGCGGGAACTATCCGGCTCTGTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTACATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGC  >  NC_000913/291772‑292259

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