Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	961,772	C→T	A191V (GCG→GTG)	cmk →	cytidylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	961,772	0	C	T	88.7%	136.6 / 3.6	53	A191V (GCG→GTG)	cmk	cytidylate kinase
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base T (13/34);  total (16/37)
Fisher's exact test for biased strand distribution p-value = 3.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

GAAGCATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCAATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGATTAAACATGACTGAATCTTTTGCTCAACT  >  NC_000913/961531‑962017                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       gAAGCATTAGTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGt                                                                                                                                                                                                                                               <  2:219073/251‑1 (MQ=255)   aGCATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGACGATGGCCGCGACATGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTACTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGCGAAGGGATTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGAGGTAGCGCCACTGGTTc                                                                                                                                                                                                                                             >  1:314275/1‑251 (MQ=255)      attattGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCgg                                                                                                                                                                                                                                          >  2:103489/1‑251 (MQ=255)         attGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCa                                                                                                                                                                                                                                        <  2:326324/250‑1 (MQ=255)          ttGCTTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAg                                                                                                                                                                                                                                       <  2:94843/250‑1 (MQ=255)                cGCCAACGCGCGTTTCGCGAATTTCCAGGTCTGATTTCCGATGGCCGCGACATGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCAGCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTTAGCGCCGTTTGGCCGAGATCAAAGAACGCGACGGCCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGc                                                                                                                                                                                                                                <  1:272569/251‑1 (MQ=255)                cGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGc                                                                                                                                                                                                                                <  1:284821/251‑1 (MQ=255)                 gCCAACGCGCGTTTCGCGCATTACCAGGTCTGATTGCAGATGGCCGCGACATGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCgctgct                                                                                                                                                                                                                               >  2:381311/1‑251 (MQ=255)                 gCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGCt                                                                                                                                                                                                                               >  1:63583/1‑251 (MQ=255)                  ccAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGCtt                                                                                                                                                                                                                              <  2:44732/251‑1 (MQ=255)                         gcgTTTCGCGAATTACAAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGGGGTATTCCATGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGACCGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGACCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGATTTAGTGtt                                                                                                                                                                                                                       >  1:79012/1‑251 (MQ=255)                                   aaTTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATTGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTccacc                                                                                                                                                                                                             >  1:165176/1‑251 (MQ=255)                                   aaTTACCAGGTCTGATCGCCGATGGTCGGGACATGGGAACGGTGGTATTCTCTGATGCCCCAGTGAAAATTTTCCTTGACCCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTGAGTGTTAACGTTGAGCGCCTTTTGGCCGAGATCAAATAACGCGACGACCGCGATCGTAACCTAGCGGTAGTTCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTTGATTccacc                                                                                                                                                                                                             <  1:298616/251‑1 (MQ=255)                                        ccAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAAGTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTGCAGTTGCAGGAGAAGGGCTTTAGTGTTAACTGTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGACACTGTTTCCGGCAGCCTATGCTTTAGTTTTGGATTCCACCACCtt                                                                                                                                                                                                        <  1:136493/251‑1 (MQ=255)                                                                gACATGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGaaaa                                                                                                                                                                               <  2:341868/251‑1 (MQ=255)                                                                  cATGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAg                                                                                                                                                                             <  1:391864/251‑1 (MQ=255)                                                                   aTGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGGGTGGGATTCCACCACCTTAAGCAGTGAGCAAG‑TGATTGAAAAAgc                                                                                                                                                                            <  2:82452/251‑1 (MQ=255)                                                                       gAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTa                                                                                                                                                                        <  2:192940/251‑1 (MQ=255)                                                                        aaCGGGGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCGTTAGTGTTAACTGTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTTCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TTATTGAAAAAGCGCTAc                                                                                                                                                                       <  2:424903/251‑1 (MQ=255)                                                                               gtATTCCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATAcgc                                                                                                                                                                <  2:63583/251‑1 (MQ=255)                                                                               ggattcgcggATGCACCCGTGAAAAGTTTACTTGACGCCTCCTCTGAAGAACGTGCGCATCGTCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCGTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTACGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATAcgc                                                                                                                                                                <  2:16104/242‑1 (MQ=255)                                                                                 aTTCCCTGATGAACCAGTGAAAATTTTCCTTGTCGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCGCCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATAcgcgc                                                                                                                                                              <  2:139543/251‑1 (MQ=255)                                                                                           gCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCGCGCCAGAAAtt                                                                                                                                                    <  2:361107/251‑1 (MQ=255)                                                                                           gCACCAGGGAAAATTTTCCGTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCGCGCCAGAAAtt                                                                                                                                                    <  2:354848/251‑1 (MQ=255)                                                                                                            ccTTGACGCCGCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCg                                                                                                                                    <  2:269811/250‑1 (MQ=255)                                                                                                               tGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGAcc                                                                                                                                 <  1:333045/250‑1 (MQ=255)                                                                                                                gACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCCACAAAAGCAGGAGAAGGGCTTTAGAGTTACCTTTGAGAGAATTTTGGCCGAGATCACAGAACGCGACGTCCGCGATCGTAACCGAGCGGGAGCGCAACTGGTTCAGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGAATTCAGCCAG‑CGATAGAAAAAGCGCTACAATACGCGCGCCCGAAATTGGCTCTCGCATAAGCGACCg                                                                                                                                >  1:384901/1‑250 (MQ=255)                                                                                                                   gCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCGCGcc                                                                                                                                                           >  1:80082/1‑220 (MQ=255)                                                                                                                   gCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCGCGcc                                                                                                                                                           <  2:80082/220‑1 (MQ=255)                                                                                                                    gcgcctgggAAGAGCGTGCGCATCGCCGGATGCTACAGTTGCAGTAGAAGGGCTTTATTGGTAACTTTGAGCGTCTTTTGGCCGAGATCAAAGAACGCGACGATCGCGATCGTGACGTAGCGGTAGTGCGAGTGGTGCCGGCAGCCGAGTCTTTATTGTGGGATGCCTCCATCTTGATCATTGTGTAAG‑TGATGGAAAAGGTGCTACTAGACTGGCGCCATGAATTGTCTCTCGcgctatcgactgaattt                                                                                                                           <  2:380157/244‑17 (MQ=255)                                                                                                                       cctcGGTAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCGGGCCAGAAATTGGCTCGCGGATAAGCGACCGAAGTTGCa                                                                                                                        <  2:288765/251‑1 (MQ=255)                                                                                                                          cGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGt                                                                                                                      >  1:152778/1‑250 (MQ=255)                                                                                                                               gaaCGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAACGGCTTTAGTGTTCACTTTGAGCGCCTTTTGGCAGAGATCAAAGAACGAGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCAGATGCTTTAGTGTTGGATTACACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTAccccc                                                                                                                >  1:174824/1‑251 (MQ=255)                                                                                                                                  cGTGCGCATTGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGtt                                                                                                             <  2:165176/251‑1 (MQ=255)                                                                                                                                        cgtagcgggctGCTACAGGTGGAGGAGAGGGGGTATGGTGTTAACTTTGAGCGGTTTTTGGCCGAGATCAAAGAACGCGACGACCGCGCACGTAACCGAGAGGTAGTGCCACTGGGTCCGGCAGCCGATGCTTTAGGGTGGGAATCAGCCTGCCTAAGCATTGAGAAAGGTG‑TTGTAGAAGCGCTCCAATACGCGCGCCAGAAATCGGCTCTCGCATAAGCGACCGAATATGCAGCACCGCAGTGGCAATg                                                                                                       >  2:131295/11‑251 (MQ=255)                                                                                                                                            gccgcATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCAATGGAAt                                                                                                   >  1:87004/1‑251 (MQ=255)                                                                                                                                              cgcctGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCGCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCAATGGAATGa                                                                                                 <  2:310740/247‑1 (MQ=255)                                                                                                                                                              aGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTCGCCGAGATCACAGAACGCGACGACAGCGCTCGTACCCGAGCGGTAGTCCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGCAAAAGCGCTACAATACGCGCGCCAGAAATTGGATCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCAATGGAATGACAGCGGGTATGTTaa                                                                                  >  1:202521/1‑250 (MQ=255)                                                                                                                                                                          ttAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTGGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATGGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTGCCCCCGTTGCAATGGAATGACAGCGGGTATGTTAAACAACCCCATcc                                                                      <  1:272123/250‑1 (MQ=255)                                                                                                                                                                                  aaCTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCTCCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCAATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAg                                                             <  1:5173/251‑1 (MQ=255)                                                                                                                                                                                        gAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTTCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGGATAAGCGACCGAATTTGCGGTACCCCCGTTGCAAGGGAATGGCAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAgg                                                        <  1:254744/250‑1 (MQ=255)                                                                                                                                                                                        gAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCAATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGt                                                       <  2:193541/251‑1 (MQ=255)                                                                                                                                                                                          gcgcCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATAGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCAATGGAATGAAAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTgg                                                     >  1:223624/1‑251 (MQ=255)                                                                                                                                                                                            cccTTTTGGCTGAGAGCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCCCTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCTCGCCAGAAAGTGGCTCTCGCATAAGCGACCGAAGTTGCAGTGCCCCCGTTGCAATGGAAGGACAGCGTGTATGGTAAACAACCCCATCCGGCATGGAGCCAGGTGGAc                                                   <  2:268114/250‑1 (MQ=255)                                                                                                                                                                                                   ggCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCAATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTaaa                                             <  1:381311/250‑1 (MQ=255)                                                                                                                                                                                                             aaaGAACGCGACGACCGCGACCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATGGGCTCGCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCAATGGAATGACAGCGGGTAGGTGAAACAACCCCAGCCGGCATGGAGCCAGGTGGTGGTTAAAGATAAACCTGa                                  <  2:206533/251‑1 (MQ=255)                                                                                                                                                                                                                 aaCGCGACGACCGCGATCGTAACCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCAATGGAATGACAGCGGGTATGTTAAACATCCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGAt                              <  1:423188/251‑1 (MQ=255)                                                                                                                                                                                                                           ccGCGATCGTAACCGAGCGATAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCGCGCCCGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCAATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAACTATAAACCTGAAGCTTAAACATGAc                    >  2:142037/1‑251 (MQ=255)                                                                                                                                                                                                                                     aaCCGAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCAATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGATTAAACATGACTGAATCtttt          <  1:253339/251‑1 (MQ=255)                                                                                                                                                                                                                                         gAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCAATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGATTAAACATGACTGAATCttt           <  1:204579/246‑1 (MQ=255)                                                                                                                                                                                                                                         gAGCGGTAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCAATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGATTAAACATGACTGAATCttt           >  2:204579/1‑246 (MQ=255)                                                                                                                                                                                                                                            cGGTAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCAATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGATTAAACATGACTGAATCTTTTGCTCAAc   <  2:51290/251‑1 (MQ=255)                                                                                                                                                                                                                                             ggTAGTGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCAATGGAATGACAGCGGGTATGTTAAACAACCACAACCTGCATGGAGCCAGGTGGGCGTTAAACATAAACCTGAAGATTAAACATGACTGAATCTTTTGCTCAAct  <  1:394174/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       GAAGCATTATTGCGTCGCCAACGCGCGTTTCGCGAATTACCAGGTCTGATTGCCGATGGCCGCGACATGGGAACGGTGGTATTCCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAG‑TGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAGTACCCCCGTTGCAATGGAATGACAGCGGGTATGTTAAACAACCCCATCCGGCATGGAGCCAGGTGGACGTTAAATATAAACCTGAAGATTAAACATGACTGAATCTTTTGCTCAACT  >  NC_000913/961531‑962017

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