Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,398,292	T→C	I294T (ATC→ACC)	mutL →	DNA mismatch repair protein MutL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,398,292	0	T	C	100.0%	190.6 / NA	56	I294T (ATC→ACC)	mutL	DNA mismatch repair protein MutL
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/27);  total (29/27)

GCGCGATTTGCGGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCC  >  NC_000913/4398049‑4398532                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 gcgcGATTTGCGGCACCGCTTTTCTTGAACAAGCGCTGGCGATGGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGAGGCGCGATCGCCTGATCAATCACGCGATCCGCTAGGCCTGCTAAGACAAACGGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCca                                                                                                                                                                                                                                           <  2:363632/251‑1 (MQ=255)  cgcgATTTGCGGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCcac                                                                                                                                                                                                                                          <  1:322703/251‑1 (MQ=255)           cGGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTgg                                                                                                                                                                                                                                 >  1:288955/1‑251 (MQ=255)             gCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGAGGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGAc                                                                                                                                                                                                                               <  1:362199/251‑1 (MQ=255)                 cGCTTTTCTTGAACAAGCGCTGGCGATTGAAGGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCa                                                                                                                                                                                                                           <  1:19226/251‑1 (MQ=255)                   cTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAAc                                                                                                                                                                                                                         >  2:116109/1‑251 (MQ=255)                        cTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGc                                                                                                                                                                                                                     >  2:90259/1‑250 (MQ=255)                            aacaatgtcTGGCGATTGAATGGCAAAACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACGCCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCTGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCAcccc                                                                                                                                                                                                                <  2:275098/243‑1 (MQ=255)                               aaGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGcc                                                                                                                                                                                                             <  1:82334/251‑1 (MQ=255)                                         gATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACgc                                                                                                                                                                                                     >  1:310106/1‑248 (MQ=255)                                         gATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGa                                                                                                                                                                                                     >  1:373726/1‑249 (MQ=255)                                            tGAAGGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTGGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTg                                                                                                                                                                                                 <  2:178468/250‑1 (MQ=255)                                             gAATGGCAAAACGGCGATCTCACGCTACGCGGCTGCGCGAGCGATACAAATGACATCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGCTCGCATGATGCGCGAACGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACAGGGGGCCGATCCGCAACCGGCATTTGTGTTGTATCTGGAGACCGACACACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGcgcg                                                                                                                                                                                               >  1:18342/1‑251 (MQ=255)                                                  gCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTcc                                                                                                                                                                                          >  2:54704/1‑251 (MQ=255)                                                       acGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATc                                                                                                                                                                                       >  1:68288/1‑249 (MQ=255)                                                            gATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTcgcg                                                                                                                                                                                >  2:127090/1‑251 (MQ=255)                                                                tcACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGCTCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTg                                                                                                                                                                            >  1:359514/1‑251 (MQ=255)                                                                 cACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTcgcg                                                                                                                                                                                <  2:358236/246‑1 (MQ=255)                                                                              tGGGTGGCCGATCAAAATCACACCACGCCCGCACTGGCAGACATTCAGTCTTGCTACGTGAACGGTCGCATGGTGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCCGGCGGACGTCAACGTGCACCCCGCCAAACACGACGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTtat                                                                                                                                                              >  2:110953/1‑251 (MQ=255)                                                                               gggTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGAGGCGCGATCGCCTGATCAATCACGCGATCCGACAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTGCCATCAGTCGCGTCTGGTGCATGATTTta                                                                                                                                                               <  1:288708/249‑1 (MQ=255)                                                                                    gCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTtatctatc                                                                                                                                                         >  1:161542/1‑250 (MQ=255)                                                                                       gATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAggg                                                                                                                                                     <  1:252150/251‑1 (MQ=255)                                                                                             aaTCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCt                                                                                                                                               <  2:76295/251‑1 (MQ=255)                                                                                                        gCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTac                                                                                                                                    <  1:137508/251‑1 (MQ=255)                                                                                                           cGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTacaac                                                                                                                                 >  1:316205/1‑251 (MQ=255)                                                                                                                tGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGATACAACAGCa                                                                                                                             >  1:253220/1‑250 (MQ=255)                                                                                                                          tCAGTATTGCTACGTGAACGGTCGCATGATGCGCGAGCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAAcg                                                                                                                   <  2:145258/250‑1 (MQ=255)                                                                                                                               aTTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTGGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGAGCTGAGCTAGCTACAACAGCAACTGGAAACGCCGCt                                                                                                              <  1:110953/250‑1 (MQ=255)                                                                                                                                       gTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTgg                                                                                                      >  2:73679/1‑250 (MQ=255)                                                                                                                                          aaCGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACg                                                                                                   >  1:264581/1‑250 (MQ=255)                                                                                                                                                           cgATCGCCTGATCAATCACGCGTTCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCAc                                                                                 <  2:264581/251‑1 (MQ=255)                                                                                                                                                             aTCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGCCATTGCTCCCCAACCTGCACcg                                                                               <  1:75232/251‑1 (MQ=255)                                                                                                                                                                gCCTGATCAATCACGCAATCCGCACGGCCTGCGCAGACAAACTGGGGGCCGATCAGCAACCGGAATTTGTGTTGTATCTGGAGACCGACCCACATCAGGAGGACGTCCACGTGCACCCCGCCAAACACGCAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGAGTGCTACCACAGCAACTGGAACCGCCGCTACCGCTGGACGATGAACCCCACCCTGccacgcgt                                                                            >  2:365591/1‑244 (MQ=255)                                                                                                                                                                 ccTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGtt                                                                           <  2:99955/251‑1 (MQ=255)                                                                                                                                                                 ccTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGtt                                                                           >  1:158368/1‑251 (MQ=255)                                                                                                                                                                   tGATCAATCACGCGATCCGCCAGGCCTGCGAAGACCAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCAGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCGGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAACCGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGTTcc                                                                         >  2:330362/1‑251 (MQ=255)                                                                                                                                                                                 aTCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACcg                                                           >  2:288316/1‑251 (MQ=255)                                                                                                                                                                                   ccGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACcgc                                                          <  2:310167/250‑1 (MQ=255)                                                                                                                                                                                    cGCCAGGCCTGCGAAGACACACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACACCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAACCcgcgg                                                        >  2:330364/1‑250 (MQ=255)                                                                                                                                                                                     gCCAGGCCTGCGAAGACAAACGGGGGGCCGATCAGCAACCGGCATTTGTGTGGTATCTGGAGACCGACCCACATCAGGGGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTg                                                       <  2:161542/251‑1 (MQ=255)                                                                                                                                                                                           ccTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGgcggcg                                                 >  2:286312/1‑251 (MQ=255)                                                                                                                                                                                                    aCAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTACACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCaa                                         >  2:185387/1‑250 (MQ=255)                                                                                                                                                                                                      aaaCTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCa                                      <  2:328969/251‑1 (MQ=255)                                                                                                                                                                                                            ggggCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGc                                >  1:181781/1‑251 (MQ=255)                                                                                                                                                                                                              ggCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAg                              <  2:26240/251‑1 (MQ=255)                                                                                                                                                                                                                 cGATCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAAc                           >  1:116848/1‑251 (MQ=255)                                                                                                                                                                                                                   aTCAGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCCCCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCg                         <  2:75772/251‑1 (MQ=255)                                                                                                                                                                                                                      ccgaaCCGGCATTTGTGTTGTATCTGGAGACCGACCCCCATCCGGTGGACGTCAACGTGCCCCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCa                      <  1:326618/248‑1 (MQ=255)                                                                                                                                                                                                                      aGCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCa                      <  1:38870/251‑1 (MQ=255)                                                                                                                                                                                                                       gCAACCGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGTGCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAg                     <  2:115837/251‑1 (MQ=255)                                                                                                                                                                                                                           ccGGCATTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCg                 <  1:127090/251‑1 (MQ=255)                                                                                                                                                                                                                                aTTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGc            >  1:114824/1‑251 (MQ=255)                                                                                                                                                                                                                                aTTTGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGc            >  1:100460/1‑251 (MQ=255)                                                                                                                                                                                                                                 tttGTGTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGcc           >  1:113987/1‑251 (MQ=255)                                                                                                                                                                                                                                   tgtgTTGTATCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCgg         >  1:51115/1‑251 (MQ=255)                                                                                                                                                                                                                                           aTCTGGAGACCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTcc  <  1:288316/250‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 GCGCGATTTGCGGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGCCTGATCAATCACGCGATCCGCCAGGCCTGCGAAGACAAACTGGGGGCCGATCAGCAACCGGCATTTGTGTTGTATCTGGAGATCGACCCACATCAGGTGGACGTCAACGTGCACCCCGCCAAACACGAAGTGCGTTTCCATCAGTCGCGTCTGGTGCATGATTTTATCTATCAGGGCGTGCTGAGCGTGCTACAACAGCAACTGGAAACGCCGCTACCGCTGGACGATGAACCCCAACCTGCACCGCGTTCCATTCCGGAAAACCGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCC  >  NC_000913/4398049‑4398532

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