Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,215,670	G→T	D65Y (GAT→TAT)	aceB →	malate synthase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,215,670	0	G	T	94.9%	113.0 / 0.3	39	D65Y (GAT→TAT)	aceB	malate synthase A
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (21/16);  total (22/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.77e-01

TCGTTCACAGTGGGGAAGTTTTCGGATCCATGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTGATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTGATGCGGTTAACGGCACCATC  >  NC_000913/4215431‑4215909                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                tCGTTCACAGTGGGGAAGTTTTCGGATCCATGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTtattttattt                                                                                                                                                                                                                                        <  2:147413/249‑1 (MQ=255)      cacaGTGGGGAAGTTTTCGGATCCATGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAAc                                                                                                                                                                                                                                 >  1:118775/1‑251 (MQ=255)                aaGTTTTCGGATCCATGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCAtt                                                                                                                                                                                                                       <  1:826/251‑1 (MQ=255)                 aGTTTTCGGATCCATGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTc                                                                                                                                                                                                                      <  1:204947/251‑1 (MQ=255)                      ccGGATCCATGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTGCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCAGTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATGGATAACGGAACGTTGCCTTATTATATTTCGGAAACAGCTTCCATTCTCGAt                                                                                                                                                                                                                 <  2:244684/250‑1 (MQ=255)                           tCCATGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTgt                                                                                                                                                                                                            >  2:388426/1‑250 (MQ=255)                              aTGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTg                                                                                                                                                                                                         >  1:29959/1‑251 (MQ=255)                               tGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTgg                                                                                                                                                                                                        >  1:13114/1‑251 (MQ=255)                                       agCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGGATGGCGAGAAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCAACAACGCAATAAACTACTGGCAGCGCGCATTCAGCAGCAGCAAGATATTATTAACGGAACGTTGCATTATTTTATTTCGGAACCAGCTTCCATTCGCGATGCTACTTGGaaattt                                                                                                                                                                                                  >  2:405175/1‑246 (MQ=255)                                                aTGACTGAACAGGCAACATCAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTcc                                                                                                                                                                                       <  1:68663/251‑1 (MQ=255)                                                                            aaCTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCCTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTagag                                                                                                                                                           <  1:19284/251‑1 (MQ=255)                                                                                         aaGGCCGTATGGCGAGCAGGAGGAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTGATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCACGGGAt                                                                                                                                                                                          <  1:238392/207‑1 (MQ=255)                                                                                         aaGGCCGTATGGCGAGCAGGAGGAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGAGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTGATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCACGGGAt                                                                                                                                                                                          >  2:238392/1‑207 (MQ=255)                                                                                                       aGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGGGACGCATTTTACGCCACAAAGCAATAAACTTATGGCCGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGCACGTTGCCTTATTTTATTTCGTAAAAAGCTTCCATTCTCGATACTGATTGTTAAATTCGCGGGTTACATGCAAACGAAGATGACCGCCGAGTAGATATAACTGAACCGGTAGAGCGAAAggtg                                                                                                                                >  2:381146/1‑248 (MQ=255)                                                                                                              gcaGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATc                                                                                                                          >  2:103155/3‑250 (MQ=255)                                                                                                                               ccGAAGCAGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCa                                                                                                          <  2:113656/249‑1 (MQ=255)                                                                                                                                     cGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGaa                                                                                                   >  1:349417/1‑250 (MQ=255)                                                                                                                                          gAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCtt                                                                                             >  2:229579/1‑251 (MQ=255)                                                                                                                                          gAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCAATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGCTGGTGATCAAAGCGCTCAACGCCAATGTGAAAGtatt                                                                                             >  2:336795/1‑248 (MQ=255)                                                                                                                                                  gactgaGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCg                                                                                     <  2:13114/251‑1 (MQ=255)                                                                                                                                                           ggTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGaa                                                                             <  1:240174/250‑1 (MQ=255)                                                                                                                                                              gACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGAtt                                                                         >  1:125552/1‑251 (MQ=255)                                                                                                                                                               aCGCATTTTACGCCACAACGCAATAAACGTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAATACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTc                                                                        <  1:493470/251‑1 (MQ=255)                                                                                                                                                                 gCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGAt                                                                                   <  1:386475/238‑1 (MQ=255)                                                                                                                                                                 gCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGAt                                                                                   >  2:386475/1‑238 (MQ=255)                                                                                                                                                                    ttttACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTag                                                                                                                                           <  1:427559/179‑1 (MQ=255)                                                                                                                                                                    ttttACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTag                                                                                                                                           >  2:427559/1‑179 (MQ=255)                                                                                                                                                                            cacaACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGAc                                                           <  1:388426/251‑1 (MQ=255)                                                                                                                                                                                            tCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATc                                            <  1:229579/250‑1 (MQ=255)                                                                                                                                                                                                    cgcgcATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCa                                    >  1:137520/1‑250 (MQ=255)                                                                                                                                                                                                      cgcATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAAt                                 >  2:277561/1‑251 (MQ=255)                                                                                                                                                                                                            cagcagcagCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCt                           <  1:352631/251‑1 (MQ=255)                                                                                                                                                                                                                agcagcaAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGt                       >  1:426434/1‑251 (MQ=255)                                                                                                                                                                                                                agcagcaAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGt                       >  1:198424/1‑251 (MQ=255)                                                                                                                                                                                                                    gcaAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTGATg                   <  1:215900/251‑1 (MQ=255)                                                                                                                                                                                                                     caAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTGATGc                  >  2:106245/1‑251 (MQ=255)                                                                                                                                                                                                                     caAGATATTGATAACGGAACGTTGCCTTATTTTAGTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTGATGc                  <  2:7493/251‑1 (MQ=255)                                                                                                                                                                                                                     caAGATATTGATAACGGAACGTTGACTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGTTTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGAAAATTAACCTGAGTGGTGc                  >  1:438767/1‑251 (MQ=255)                                                                                                                                                                                                                               aTAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGCAGATTCACTGGCACCAGACTGGAACAAAGAGATCGCCGGGCAAATTAACCTGCGTGATGCGGTTAACGGc        >  1:338797/1‑251 (MQ=255)                                                                                                                                                                                                                                     gAACGTTGCCTTATTTTATTTCGGCAACAGCTTCCATTCGCGATGATGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGTCTATTTTACCTGCGTGATGCGGTTAACGGCaaaatc  >  2:148776/1‑246 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                TCGTTCACAGTGGGGAAGTTTTCGGATCCATGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTGATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTGATGCGGTTAACGGCACCATC  >  NC_000913/4215431‑4215909

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