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	98.5%	220.4 / NA	69	D65Y (GAT→TAT)	aceB	malate synthase A
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base T (33/34);  minor base C (1/0);  total (34/35)
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.98e-01

GATCCTTCGTTCACAGTGGGGAAGTTTTCGGATCCATGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTGATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTGATGCGGTTAACGGCACCATCAG  >  NC_000913/4215425‑4215911                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  gATCCTTCGTTCACAGTGGGGAAGTTTTCGGATCCATGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTAACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTAtttt                                                                                                                                                                                                                                              >  1:163337/1‑251 (MQ=255)            cacaGTGGGGAAGTTTTCGGATCCATGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAAc                                                                                                                                                                                                                                   >  1:99175/1‑251 (MQ=255)              caGTGGGGAAGTTTTCGGATCCATGACGAGGAGCTGCACGATGACTGAAAAGGCAACAACAACCGTTGAACTGGCATTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCAGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGAGCATTCAGCAACAGAAAGATATTGATAACGTAACGTTGCCTTATTTTATTTCGGCAActg                                                                                                                                                                                                                                 >  2:495607/1‑249 (MQ=255)               aGTGGGGAAGTTTTCGGATCCATGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGc                                                                                                                                                                                                                                <  1:187467/251‑1 (MQ=255)                gTGGGGAAGTTTTCGGATCCATGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAAAAGCt                                                                                                                                                                                                                               >  1:352264/1‑251 (MQ=255)                      aaGTTTTCGGATCCATGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCAtt                                                                                                                                                                                                                         <  1:620328/251‑1 (MQ=255)                              ggATCCATGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGc                                                                                                                                                                                                                 >  2:614119/1‑251 (MQ=255)                              ggATCCATGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGACACAACGCAATACACTTCTGGAAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGAttc                                                                                                                                                                                                                 >  1:383344/1‑249 (MQ=255)                               gATCCATGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCt                                                                                                                                                                                                                <  1:548783/251‑1 (MQ=255)                                        cGAGGAGCTGCACGATGACTGACCAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCCAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCCAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGaa                                                                                                                                                                                                        >  2:292925/1‑250 (MQ=255)                                            gagCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCAGTCGCGATGCGGATTGGAAAATTc                                                                                                                                                                                                   <  2:209250/251‑1 (MQ=255)                                                    cGATGCCTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGACGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCAGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTAGCGATGCTGATTGGAAAATTCGCGGGAtt                                                                                                                                                                                           >  2:94315/1‑251 (MQ=255)                                                            gAACAGTCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATGTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGGAGCAAGATATTTAGAACGGAACGTTGCCTTATTTGATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAACTCGCGTGATGCCTGCg                                                                                                                                                                                     <  2:603675/249‑1 (MQ=255)                                                            gAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATAGTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCGGATTGGAAAATTCGCGGGATTCCTGCgg                                                                                                                                                                                    <  2:615161/250‑1 (MQ=255)                                                             aaCAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGAc                                                                                                                                                                                  >  2:104917/1‑251 (MQ=255)                                                                  gCAACAACAACCGATGAATTGTCTTTCACAAGGTCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTGTACGCCGCAACGCAATAAACATCTGGTAGCGCGCTTTCAGCAGCAGCAAGATATTGATAACGGAACGTGGCCTTATTTTATTTCGGACACAGCTTCTATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTaga                                                                                                                                                                             <  1:611268/251‑1 (MQ=255)                                                                      caacaacCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGAc                                                                                                                                                                         >  1:495439/1‑251 (MQ=255)                                                                             cGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCcgcg                                                                                                                                                                  <  1:52708/251‑1 (MQ=255)                                                                              gATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTGCTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCtt                                                                                                                                                                                                                              <  1:381926/190‑1 (MQ=255)                                                                              gATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGCCTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCtt                                                                                                                                                                                                                              >  2:381926/1‑190 (MQ=255)                                                                                 gAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGt                                                                                                                                                                 >  1:473873/1‑248 (MQ=255)                                                                                  aaCTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCTAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATGCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTagag                                                                                                                                                             <  1:578425/251‑1 (MQ=255)                                                                                   aCTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAgaga                                                                                                                                                            <  1:104917/251‑1 (MQ=255)                                                                                     tGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGAt                                                                                                                                                           >  2:161196/1‑250 (MQ=255)                                                                                                 ggCCGTATGGCGAGCATGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTa                                                                                                                                              <  2:245014/251‑1 (MQ=255)                                                                                                       aTGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAgcgc                                                                                                                                        >  1:574836/1‑251 (MQ=255)                                                                                                            gAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCa                                                                                                                                       >  1:159950/1‑247 (MQ=255)                                                                                                                aGGAGAAGCAAATTCTTTCTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCATCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTg                                                                                                                               <  2:352264/251‑1 (MQ=255)                                                                                                                   agaAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATAGTGATAACGGAACGTTGCCTTATTTTATGTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGGCTTAGAAGACCGCCGCGTAGAGATAACGGGCCCTGTAGAGCGCAAGACGGTGATc                                                                                                                            <  2:608789/251‑1 (MQ=255)                                                                                                                      aGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAACATGGTGATCAAc                                                                                                                         <  1:162083/251‑1 (MQ=255)                                                                                                                       gCAAATTATTACGGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTCATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGATGACCGCCGCGTAGAGATAACTGACCCGGTAGAGCGCAAGATGGTGATCAAcg                                                                                                                        >  2:133926/1‑251 (MQ=255)                                                                                                                         aaaTTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAAcgcg                                                                                                                      <  1:367260/251‑1 (MQ=255)                                                                                                                            ttCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGATGGTGACGCATTTTACGCCACACCGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATACCGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGAAGAGCGCAAGATGGTGATCAACGCGCTc                                                                                                                   >  2:238358/1‑251 (MQ=255)                                                                                                                             tCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTc                                                                                                                   <  2:76255/250‑1 (MQ=255)                                                                                                                                tACTGCCGAAGCAGTAGAATTTCTGGCTTAGCTGGTGACGCATTTTACACCACAACGCAATAAACTACTGGCAGCGCGTAGTGAGCAGCATCAAGATATTGATAACGGAACGTTGCCCTAGTTTAGTTCGGAAACAGCTTCTATGCGCGATGCTTATTGGAGAATTCGCGCGAGCCCTGCGTACTTAGAAGACCGCCGCGTAGAGTTAACTGGCCCTGTAGAGCGCAAGATGGTGATCAACGCGCTCAAc                                                                                                                <  1:43455/250‑1 (MQ=255)                                                                                                                                    gCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTGTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGcc                                                                                                             <  1:249906/249‑1 (MQ=255)                                                                                                                                     ccGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAAt                                                                                                          <  1:8707/251‑1 (MQ=255)                                                                                                                                         aGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAAtgtg                                                                                                       >  2:369531/1‑250 (MQ=255)                                                                                                                                          gCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGa                                                                                                      >  2:120416/1‑250 (MQ=255)                                                                                                                                           cGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTGCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGaaa                                                                                                    <  1:292925/251‑1 (MQ=255)                                                                                                                                            ggTAGAATTTCTGACTGAGCTGGCGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGaa                                                                                                     >  1:183984/1‑249 (MQ=255)                                                                                                                                             gTAGAATATCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAAATTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGATTCCATTCGCGATGATGATTGAAAAATTCGCGGGATTCCTGCGGCCTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGAGATCAACGCGCTCAACGCAAATGTGAAAGt                                                                                                  >  1:76073/1‑251 (MQ=255)                                                                                                                                                 aaTTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCttt                                                                                              <  1:57051/251‑1 (MQ=255)                                                                                                                                                    ttCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGCTTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATg                                                                                           >  2:128303/1‑251 (MQ=255)                                                                                                                                                    ttCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTAt                                                                                            <  1:151562/250‑1 (MQ=255)                                                                                                                                                             aGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTc                                                                                  <  2:397743/251‑1 (MQ=255)                                                                                                                                                                         atattACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTCGCGAAAACTGACACGGAAGAGCGAAAGATGGTGATAAACGCGCTTAACGCCAATGTGAAAGTCTTTATGGACGATTTCGAAGATTCACTg                                                                      >  2:178618/4‑251 (MQ=255)                                                                                                                                                                          ttttACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTgg                                                                     >  2:11017/1‑251 (MQ=255)                                                                                                                                                                                   aaaaCGAAATCAACTTCTCGCAGCCCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGAGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCACGATGTTGATCAACGCGCTCAACGCCAATGCGAATGTATTAATGGCTGATTTCGCAGATTCACTGGCCCCaaac                                                             >  2:72215/3‑247 (MQ=255)                                                                                                                                                                                       cGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGa                                                         >  1:173911/1‑250 (MQ=255)                                                                                                                                                                                        gCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAAc                                                       <  2:183984/251‑1 (MQ=255)                                                                                                                                                                                         caagaaaCTGCTGGCAGCGCGCAGTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGTCTGGAACa                                                      <  2:495439/247‑1 (MQ=255)                                                                                                                                                                                           aTAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACaaa                                                    <  2:306597/251‑1 (MQ=255)                                                                                                                                                                                               aCTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATGCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGa                                                <  1:543665/251‑1 (MQ=255)                                                                                                                                                                                                  tCTGGCAGAGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGCTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGAt                                               <  2:412899/249‑1 (MQ=255)                                                                                                                                                                                                   cTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCg                                             >  2:107789/1‑250 (MQ=255)                                                                                                                                                                                                     ggCAGCGCGCATGCAGCAGCAGCAATATATTGATAACGGAACGTTGCCTGAATTTATTTCTGAAACAGCTTCCATTCGCGATGCTGATTGTAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGGGCTCAACGCCAATGAGAAAGTCTTTATTGCCGATTTCGAAGATTCACTGGCACCAGACTGTAACAAAGTGATCg                                             <  1:94315/248‑1 (MQ=255)                                                                                                                                                                                                            cgcATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAAt                                   <  1:69878/251‑1 (MQ=255)                                                                                                                                                                                                              cATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTa                                 <  1:151815/251‑1 (MQ=255)                                                                                                                                                                                                              cATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCGTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCGTTATGGCCGACTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGTTCGACGGGCAAATTa                                 <  2:542191/251‑1 (MQ=255)                                                                                                                                                                                                               aTTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTa                                 >  2:319877/1‑250 (MQ=255)                                                                                                                                                                                                                      agcagcaAGATATTGATAACGGAACGTTGCCTTATTTTAGTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGAGCGACGGGCAAATTAACCGGCGt                         <  2:99175/251‑1 (MQ=255)                                                                                                                                                                                                                                  ttGATAACGGAACGTTGCCTTATTTTATTTCGGCAACAGCTTCCATTCGCGATGCTGATTGGAATATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTGAGGCGGTTAAc             >  2:296560/1‑251 (MQ=255)                                                                                                                                                                                                                                   tGATAACGGAACGTTGCATTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGAGATGCGGTTAACg            >  2:303022/1‑251 (MQ=255)                                                                                                                                                                                                                                     ataccGGAACGTTGCCTTCTTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGCTTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTGATGCGGTTAACGGc          >  2:503476/5‑251 (MQ=255)                                                                                                                                                                                                                                       aaCGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTGATGCGGTTAACGGCAc        >  1:398467/1‑251 (MQ=255)                                                                                                                                                                                                                                          ggAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGCCTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGCTTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTGATGCGGTTAACGGCACCAt     >  2:238807/1‑251 (MQ=255)                                                                                                                                                                                                                                             aCGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTGATGCGGTTAACGGCACCATc    >  2:578820/1‑249 (MQ=255)                                                                                                                                                                                                                                             aCGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTGATGCGGTTAACGGCACCATCAg  <  2:264813/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  GATCCTTCGTTCACAGTGGGGAAGTTTTCGGATCCATGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTGATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTGATGCGGTTAACGGCACCATCAG  >  NC_000913/4215425‑4215911

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