Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,810,784	A→G	H123H (CAT→CAC)	mutM ←	DNA‑formamidopyrimidine glycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,810,784	0	A	G	100.0%	198.6 / NA	58	H123H (CAT→CAC)	mutM	DNA‑formamidopyrimidine glycosylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/24);  total (34/24)

AGCTAACAATTCACACTCTGCCAGCGACAGTGATGACGCCAGCCGATCCGGATGGATCCCCGCCGCAAACAGTGATTCGCTGGCATAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGGGTCAGCACATTATGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTGCCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCT  >  NC_000913/3810556‑3811020                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             aGCTAACAATTAACACTCTGCCAGCGACAGTGATGACGCCAGCCGATCCGGATGGATCCCCGCCGCAAACAGTGATTCGCTGGCATAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATAAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTcc                                                                                                                                                                                                                        >  2:184659/1‑251 (MQ=255)     aaCAATTCACACTCTGCCAGCGACAGTGATGACGCCAGCCGATCCGGATGGCTCCCCGCCGCAAACAGCGATTCGCTGGCATAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGCAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCcac                                                                                                                                                                                                                      >  1:30738/1‑249 (MQ=255)     aaCAATTCAAACTCTGCCAGCGACAGAGATGACGCCAGCCGCTCCGGATGGATCCCCGCCGCAAACCGTGATTCGCTGGCATAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGCTGAAGATACTCACCATTGAAATCGTCGCTACGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCcacac                                                                                                                                                                                                                    >  2:205688/1‑250 (MQ=255)       cAATTCACACTATGCCAGCGACAGTGATGACGCCAGCCGATCCGGATGGATCCCCGCCGCAAACAGTGATTCGCTGGCATAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGcc                                                                                                                                                                                                                  >  1:421395/1‑251 (MQ=255)         aTTCACACTCTGCCAGCGACAGTGATGACGCCAGCCGATCCGGATGGATCCCCGCCGCAAACAGTGATTCGCTGGCATAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGcc                                                                                                                                                                                                                  >  1:149599/1‑249 (MQ=255)              cacTCTGCCAGCGACAGTGATGACGCCAGCCGATCCGGATGGATCCCCGCCGCAAACAGTGATTCGCTGGCATAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCGCCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGTTCCACAGCCAGGCAcc                                                                                                                                                                                                           >  1:149560/1‑251 (MQ=255)                           aCAGTGATGACGCCAGCCGATCCGGATGGATCCACGCCGCAAACAGTGATTCGCTGGCATAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCATTCCAGCTCTTTGGTCCACAGCCAGGCCCCAAAGCGGCGCgg                                                                                                                                                                                               >  2:246262/1‑250 (MQ=255)                             gttgaggaccccggccgaccggtatggacCCCCGCCGCAAACAGTGATTCGCTGGCATAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATc                                                                                                                                                                                            <  2:263571/228‑1 (MQ=255)                             aGTGATGCCGCCAGCCGATCCGGATGGATCCCCGCCGCAGACAGTGATTCGCTGGCATAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCCCGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATc                                                                                                                                                                                            >  1:274822/1‑251 (MQ=255)                               tgatgaCGCCAGCCGATCCGGATGGATCCCCGCCGCAAACAGTGATTCGCTGGCATAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCgg                                                                                                                                                                                          >  1:162658/1‑251 (MQ=255)                                 atgaCGCCAGCCGATCCGGATGGATCCCCGCCGCAAACAGTGATTCGCTGGCATAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGgtg                                                                                                                                                                                        >  1:294234/1‑251 (MQ=255)                                        cAGCCGATCCGGATGGACCCCCGCCGCAAACAGTGATTCGCTGGCATAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCgca                                                                                                                                                                                 <  1:53144/251‑1 (MQ=255)                                                 gggatggacCCCCGCCGCAAACAGTGATTCGCTGGCATAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCGCCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACtttac                                                                                                                                                                        <  2:149560/250‑3 (MQ=255)                                                    atggatCCCCGACGCAAACAGAGATTCGCAGGAATAGATATTCCCTACCCCTCCCACCAGCTTGTTATCCATCAGCCACGGTTTACTCGCCGTCTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTACGCGGCTCCGGTCCAAGATGAGCCACAACGTTGTGCCCTTCCAGCTCTTTGGCCCCAAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCCCTTTACCGt                                                                                                                                                                     >  2:257528/1‑251 (MQ=255)                                                         cccccGCCGCAAACAGTGATTCGCTGGCATAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGGCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCtc                                                                                                                                                                <  2:294234/250‑1 (MQ=255)                                                              gccgcAAACAGTGATTCGCTGGCATAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCATGCACCAAAGCGTCGCGGAt                                                                                                                                                                                             <  2:270787/217‑1 (MQ=255)                                                              gccgcAAACAGTGATTCGCTGGCATAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGAt                                                                                                                                                                                             >  1:270787/1‑217 (MQ=255)                                                                 gcAAACAGTGATTCGCTGGCATAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCaa                                                                                                                                                        <  2:176154/251‑1 (MQ=255)                                                                   aaaCAGTGATTCGCTGGCATAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAAt                                                                                                                                                      <  1:121093/251‑1 (MQ=255)                                                                    aaCAGTGATTCGCTGGCATAGATATTCCCTACCCCTACCTCCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCTTTTTTTTCTTCGCGCACTTCTGATGCAGATACGCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGGCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATc                                                                                                                                                     <  2:399042/251‑1 (MQ=255)                                                                      aaGTGATTCGCAGGCATAGATATTCCCTACACCTACCACCCGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAACTCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCCGCTCGGTGTAGCGCAGCACTTTCCCGTTGCTCATCACCACATcc                                                                                                                                                    >  2:233406/2‑250 (MQ=255)                                                                        gTGATTCGCTGGCATAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCcac                                                                                                                                                  >  1:140306/1‑250 (MQ=255)                                                                           aTTCGCTGGCATAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATgg                                                                                                                                              >  1:268060/1‑251 (MQ=255)                                                                             tCGCTGGCATAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGt                                                                                                                                             <  1:313566/250‑1 (MQ=255)                                                                                cTGGCATAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATg                                                                                                                                         >  1:299775/1‑251 (MQ=255)                                                                                     aTAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCtttt                                                                                                                                    >  2:426171/1‑251 (MQ=255)                                                                                          tatTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAggg                                                                                                                               >  2:343411/1‑251 (MQ=255)                                                                                            ttCCCTACCCCTACCACCAGCTTGTTATCCCTCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGCGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCCTGCTTTTCAggggg                                                                                                                             >  2:401279/1‑251 (MQ=255)                                                                                            ttCCCTACACCTACCACCAGATTGTTATACATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTACGGTCCAAGATGAGACAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAggggg                                                                                                                             >  2:422956/1‑251 (MQ=255)                                                                                                ctacccctaccACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGCAGt                                                                                                                         >  2:52432/1‑251 (MQ=255)                                                                                                              aGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGAGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGAt                                                                                                           >  1:272018/1‑251 (MQ=255)                                                                                                              aGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTGCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGAt                                                                                                           <  1:436863/251‑1 (MQ=255)                                                                                                               gcctggtggccATCAGTCACTGTTTATTCGCCGTTTTTTTCTTCGCGCGCTTCTGATTCATCTACTAACCATTGAATGCGTCGCGAAGCGGCTCCGGTCCAAGATGAGACAGCACGTTGTGCCCTTCCAGCTTTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTCGCGGAGCACTTTACCGTTACTCATCACCAAACCCACATGGTCATGCTTTGCAGGGGGAAGTTCTGCTGTAAGGAt                                                                                                           <  2:245958/241‑1 (MQ=255)                                                                                                                cTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAACTCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATgc                                                                                                         >  2:97603/1‑251 (MQ=255)                                                                                                                  tGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATgc                                                                                                         <  1:177346/249‑1 (MQ=255)                                                                                                                  tGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATgc                                                                                                         <  1:274888/249‑1 (MQ=255)                                                                                                                             aGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTgg                                                                                                                 <  1:35756/230‑1 (MQ=255)                                                                                                                             aGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTgg                                                                                                                 >  2:35756/1‑230 (MQ=255)                                                                                                                                 aCGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCAGCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCGGGCTGCCAGACATc                                                                                        <  1:326027/251‑1 (MQ=255)                                                                                                                                         aTCGCCTTTTTTTTCTTCGCGCACTTCTGATGCAGAGACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTGCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCAaatg                                                                                <  1:278409/251‑1 (MQ=255)                                                                                                                                               gTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCAAATGAATGAt                                                                          >  1:384101/1‑251 (MQ=255)                                                                                                                                                   ttttCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCAAATGAATGATAATc                                                                      >  2:112050/1‑251 (MQ=255)                                                                                                                                                          cgcgCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCAAATGAATGATAATCCAGCCCt                                                               >  1:431957/1‑251 (MQ=255)                                                                                                                                                                  tCTGATGCATATACTCACCGTTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTTGTCCACAGCCATGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCTCTGTACCGTTGCTCATCACCAAATCCACATGTTCATGCTTTGCAGGTGGAAGTTCTTCTGGGAGGATGCGCAGGCTGGCAGATATTGCCAAATGAATGATAATCCAGCCCgcaggcag                                                        <  2:380800/250‑1 (MQ=255)                                                                                                                                                                     gATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCAAATGAATGATAATCCAGCCCTCAGGCAGCTc                                                     <  1:426171/250‑1 (MQ=255)                                                                                                                                                                     gATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCAAATGAATGATAATCCAGCCCTCAGGCAGCTc                                                     <  1:343411/250‑1 (MQ=255)                                                                                                                                                                      aTGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCAAATGAATGATAATCCAGCCCTCAGGCAGCTCca                                                   <  1:410636/251‑1 (MQ=255)                                                                                                                                                                                         aaaTCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCccg                                >  2:350193/1‑251 (MQ=255)                                                                                                                                                                                             cgtcgCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCcgccgc                            >  1:65647/1‑251 (MQ=255)                                                                                                                                                                                                    aaGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGc                     >  2:90005/1‑251 (MQ=255)                                                                                                                                                                                                         gcgccGGTCCAAGATGAGCCAGCACGTTGTGCCCTGCCAGCTCTGGGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGGAGCGCAGCACTTGACCGTTGCTCATCACCAAATCCACATGGTCAGGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCGGCCAGACATCCCCAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGAGATTTATCCCGCCGCTGCACGCTAAGc                <  2:268060/248‑1 (MQ=255)                                                                                                                                                                                                          cTCCGGTCCAAGATAAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGc                >  1:74840/1‑250 (MQ=255)                                                                                                                                                                                                           tCCGGTCCACGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGACGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCACATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCCGATATTTAGCCCGCCGCTGCACGCTAAGCAc              >  2:284033/1‑251 (MQ=255)                                                                                                                                                                                                            ccGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCACTTTGGTCCACATCCAGGCACCAAAGCGGCGCGGATCGGTGGAGCGCAGCACTTTACCGTTGCTCAGCACCAAATCCACATGTTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGATATCCCCAATTGAATGATAATCCAGCCCTCATGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCa               <  1:35773/249‑1 (MQ=255)                                                                                                                                                                                                               gTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGATAAGCACtggt          >  2:157232/1‑251 (MQ=255)                                                                                                                                                                                                                tCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTCGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACtggtt         >  2:366716/1‑251 (MQ=255)                                                                                                                                                                                                                 ccAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTGGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCAAAGGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACtggttg        <  2:149599/251‑1 (MQ=255)                                                                                                                                                                                                                        gagccAGCACGTTGTGCCCTCCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGTTCGGTGTAGCGCAGCCCTTTCCCGTTGCTCATCCCCAAAGCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCGGGAAGGAGGCGCAGGCTGCCAGACATCCCCAAATGAAGGAGAATCCAGCCCGCAGGCAGCTCCAGCAGCAGATATGTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTctct  <  2:431957/246‑4 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             AGCTAACAATTCACACTCTGCCAGCGACAGTGATGACGCCAGCCGATCCGGATGGATCCCCGCCGCAAACAGTGATTCGCTGGCATAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGGGTCAGCACATTATGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTGCCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCT  >  NC_000913/3810556‑3811020

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