Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	536,173	T→C	L155P (CTC→CCC)	hyi →	hydroxypyruvate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	536,173	0	T	C	93.5%	150.0 / NA	46	L155P (CTC→CCC)	hyi	hydroxypyruvate isomerase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (21/22);  minor base A (3/0);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAA  >  NC_000913/535929‑536410                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              gCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGCTGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCCTTTTGATATTCCTGGTTTCCATCCCACCgg                                                                                                                                                                                                                                         >  2:266328/1‑251 (MQ=255)       aTTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTGATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCg                                                                                                                                                                                                                                   <  2:538509/251‑1 (MQ=255)                 gTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAAAGCTTGTAGAAAAACTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCATGCGCTg                                                                                                                                                                                                                         >  2:294512/1‑251 (MQ=255)                   gaagaaGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGAGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACCGATTCACGCAACGCTTGTAGAAAACCTGCGTTATTCCGCGAATATGCTGATGAAAGGAGATATTTTATTACTGCTTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCACACCGGAACTCGGCAGGCGCTgc                                                                                                                                                                                                                        >  2:640385/1‑249 (MQ=255)                              tCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTTATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATtgattgat                                                                                                                                                                                                             <  1:669119/250‑1 (MQ=255)                               cGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATtgatg                                                                                                                                                                                                            <  2:703850/250‑1 (MQ=255)                                  gATGGCGTAGCATCAGCGATTCGTTATGCCCGTGCGCTTGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTACCCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATgt                                                                                                                                                                                                        <  2:550863/251‑1 (MQ=255)                                        gTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATgttggttg                                                                                                                                                                                                  <  1:236612/251‑1 (MQ=255)                                        gTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATgttggttg                                                                                                                                                                                                  >  1:721681/1‑251 (MQ=255)                                          agcagcagcGATTCGATATGCACGTGCGCTTGGTAATAAAAAAATTAACTGTCTGGTCGGTACCACGAAGGCTGGTTTCAGTAGTGAACAGATTCACGCAACACTTGTAGAAAACCTGCGTTATGCCGCGACTCTGATGATGAAAGAAGATGTTTTATTACTGGTTGAACATATTAACCATTTTAGTATTCCTGGTTTCCATCACCCCGGAACTGAGGAGGCGCTGATATTGATTGATTTTGTTGGAtgct                                                                                                                                                                                                >  2:24966/1‑251 (MQ=255)                                            cagcagcTATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTtgctgc                                                                                                                                                                                              <  1:365440/251‑1 (MQ=255)                                               cagcGATGCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAAt                                                                                                                                                                                           <  1:705086/251‑1 (MQ=255)                                                   gATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTg                                                                                                                                                                                        >  1:516109/1‑250 (MQ=255)                                                             gCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTTATNATGTTGGTTGCTGCAATTTGAAAATTCAGt                                                                                                                                                                              >  2:362669/1‑250 (MQ=255)                                                                    cgcTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACAtt                                                                                                                                                                      >  1:665744/1‑251 (MQ=255)                                                                     gcTGGGTAATAAAAAACTTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCAACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACAttt                                                                                                                                                                     >  1:252277/1‑251 (MQ=255)                                                                        gggTAATAAAAAAATTAACTGCCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGAGTGAACCTATTAACCATATTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCACTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATc                                                                                                                                                                  <  2:437072/251‑1 (MQ=255)                                                                                ctaaaaTTAACTGTCTGGGCGGTATAACGCCGGCTGGTTTCAGAAGTGAACAGCTTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCAATCCAACAGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTGGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGc                                                                                                                                                            >  2:608144/3‑249 (MQ=255)                                                                                                 gTCGTTAAAACGCCGGCTAGTTTCAGCAGGGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCTAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGTTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGTTGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAATGTGAAtt                                                                                                                                         <  2:622654/251‑1 (MQ=255)                                                                                                   cGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGCTATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGGTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGACGGAAGGTGAATTa                                                                                                                                        >  2:79497/1‑250 (MQ=255)                                                                                                          aCGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAacac                                                                                                                                <  1:640726/251‑1 (MQ=255)                                                                                                                    gTTTCAGCAGTGAACAGAGTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGTAAGGTGAATTAACCAACACCATGACTc                                                                                                                        >  2:600137/1‑249 (MQ=255)                                                                                                                         agcagTGAACAGATTCACGCAACGCTTGTAGNAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTAATAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGATGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAAGTCAGTATGACATTTATCATAGGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGc                                                                                                                 <  2:351282/251‑1 (MQ=255)                                                                                                                             gTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGAt                                                                                                             >  1:350713/1‑251 (MQ=255)                                                                                                                               gAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATa                                                                                                            <  2:140647/250‑1 (MQ=255)                                                                                                                                    gATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTg                                                                                                      >  1:159009/1‑251 (MQ=255)                                                                                                                                                 attGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAAt                                                                                         <  2:516109/250‑1 (MQ=255)                                                                                                                                                   tGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAAt                                                                                         <  1:630490/249‑1 (MQ=255)                                                                                                                                                                        cgcgAATATGCTGCTGAAAGAAGATATTTTCTTCATGATTGACCCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCgcggcg                                                                  >  1:559646/1‑251 (MQ=255)                                                                                                                                                                            aaTATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTCACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCgaacc                                                              >  2:18700/1‑251 (MQ=255)                                                                                                                                                                            aaTATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCgaacc                                                              <  2:721681/251‑1 (MQ=255)                                                                                                                                                                               gtGCTGATGAAAGAAGATATTTTATGACGGATTGAGCCGAGTAACCATGTTGATATTCCTGGTTTCCATCCCGCCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACAGTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACGCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCgaaccggg                                                           <  2:270648/250‑1 (MQ=255)                                                                                                                                                                                       ggaagaagAGATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGACAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGACAATTCAGTATGGCATTTATCATATGCAGGGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGCTACTCCGCATCGCGGCggaccgggcaccgg                                                     >  2:132397/3‑249 (MQ=255)                                                                                                                                                                                        aaaGAAGAGATTTTATTACTGATGGAACCTATGAACCATATTGATATTCCTGGTGTACATCCCACCGGAACTAGGCAGGCGCTGAGATTGATTGATGATGTTGGTTGCTGCAATTAGAAAATTCAGTATGACATTGATCATATGCAGCGGATGGAAGGTGAATTAACCGACACCATGACTCAGTGGGCAGATAAAATTGGTCACCTGCAAAGGGCCGATAATCCGCATCGCGGCGAACCGGGAACCGgaga                                                  <  2:163028/251‑1 (MQ=255)                                                                                                                                                                                         aagaagATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGgaga                                                  >  1:14936/1‑250 (MQ=255)                                                                                                                                                                                          agaagaTATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGaaa                                                <  1:192082/251‑1 (MQ=255)                                                                                                                                                                                            aagaTATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAaatt                                              <  1:192218/251‑1 (MQ=255)                                                                                                                                                                                                    ttattaCTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATg                                       >  1:726528/1‑250 (MQ=255)                                                                                                                                                                                                    ttattaCTGATTGAACCTAGTAACCATTTTNATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATg                                       <  2:357251/250‑1 (MQ=255)                                                                                                                                                                                                            gATTGAACCTATTACCCATTTTGATCTTCATGGTTTCCATCACACCGGAACTCGGCAGGCGCTGAAATTCATTGATGATGTTGGATGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGAGAATTAACCAACACCATGACTCTGTGGGCTGATAAACTTGGTCACCTGCAACTTGCCGTTAATCCGCATCGCGGCGAACCGGGAACCGGCGAAATTAATTATGCTAATctct                              >  2:675128/1‑251 (MQ=255)                                                                                                                                                                                                             aTTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCtt                             >  1:244315/1‑251 (MQ=255)                                                                                                                                                                                                                         aaCCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCg                  >  1:142057/1‑250 (MQ=255)                                                                                                                                                                                                                         aaCCATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCg                  <  1:267823/250‑1 (MQ=255)                                                                                                                                                                                                                            cATTTTGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGaaa               >  1:379180/1‑250 (MQ=255)                                                                                                                                                                                                                              ttttGATATTCCTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGCTGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAAtt            >  2:376423/1‑251 (MQ=255)                                                                                                                                                                                                                                        ccTGGTTTCCATCCCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACaa  <  2:726528/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              GCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCACGCAACGCTTGTAGAAAACCTGCGTTATGCCGCGAATATGCTGATGAAAGAAGATATTTTATTACTGATTGAACCTATTAACCATTTTGATATTCCTGGTTTCCATCTCACCGGAACTCGGCAGGCGCTGAAATTGATTGATGATGTTGGTTGCTGCAATTTGAAAATTCAGTATGACATTTATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAA  >  NC_000913/535929‑536410

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