Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,810,937	T→C	L72L (TTA→TTG)	mutM ←	DNA‑formamidopyrimidine glycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,810,937	0	T	C	100.0%	148.3 / NA	46	L72L (TTA→TTG)	mutM	DNA‑formamidopyrimidine glycosylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/26);  total (20/26)

ATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGGGTCAGCACATTATGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTGCCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGT  >  NC_000913/3810692‑3811180                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    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)                                                                                               gCCCTTCCAGCTCTTTGGGCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGAtctc                                                                                                                                                   <  1:97603/250‑1 (MQ=255)                                                                                                cccTTCCAGCCCTTTGGTCCACAGCCAGGCACAAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCAACACCAAAGCCACATGTTCATGCTTTTCAGGGGGAATTTCTTCTGGAAGGATGCTCAGGCTGCCAGACATCCCCAAATGAATGATAGTCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCtctt                                                                                                                                                 <  1:218592/251‑1 (MQ=255)                                                                                                          tctTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCAATTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCCGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTATCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTCGATCTCTTCTGAAACCgg                                                                                                                                       >  1:251263/1‑251 (MQ=255)                                                                                                                         cAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTcc                                                                                                                        <  1:378200/251‑1 (MQ=255)                                                                                                                         cAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTGTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAAGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACAACCCCAAATGAATGATAACCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTcc                                                                                                                        <  1:432915/251‑1 (MQ=255)                                                                                                                           ggccccAAAGCGGCTCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTCCAGGGGGAAGTTCTTCTGGAAGGAGGCGCAGGCGGCCAGACATCCCCAAATGAATGATAATCCAGCCCTCAGGAAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGt                                                                                                                      <  2:299775/247‑1 (MQ=255)                                                                                                                                         gcgGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGcaccaccgc                                                                                                         >  1:385339/1‑247 (MQ=255)                                                                                                                                                                 taccGTTGCTCATCACCAAATCCCCATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACgag                                                                                 <  1:139592/248‑1 (MQ=255)                                                                                                                                                                 taccGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGCCATCCCCAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCCCCACCGCATGAAGAATGGTTGCACCAACg                                                                                   >  2:180073/3‑248 (MQ=255)                                                                                                                                                                       tGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGt                                                                          >  1:46892/1‑251 (MQ=255)                                                                                                                                                                        gtccagcaCCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGtt                                                                         <  2:110930/245‑1 (MQ=255)                                                                                                                                                                          ccatcaCCAAAGCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCAAAGGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAATCCGGCCAGCGCAAGCGTCCTTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCt                                                                       <  1:425473/250‑1 (MQ=255)                                                                                                                                                                                    acccACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGCCTGCCAGACATCCCCAAATGAATGATAATCCAGCCCTCAGGCAGCTCAAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCcgcgac                                                             <  1:246262/249‑3 (MQ=255)                                                                                                                                                                                    aTCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCcgcgac                                                             >  1:255220/1‑249 (MQ=255)                                                                                                                                                                                      ccACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGTAGGCTGCCAGACATCCCCAAAGGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGcctcgacg                                                            <  1:366716/250‑7 (MQ=255)                                                                                                                                                                                          aTGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCAAATGAATGATAACCCAGCCCTCAGGCAGCGCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGACGGCTgg                                                       <  1:248893/251‑1 (MQ=255)                                                                                                                                                                                          aTGGTCAAGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGTTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTGGCACCAACGAGATGCGGTTCTATGCCGCGACGGCTgg                                                       >  2:61858/1‑251 (MQ=255)                                                                                                                                                                                               cATGCTTTTGAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCCGACATCAACACATGATTGATAATCCAGCCATCAGGCAGCTCCCGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGCTGCTCGATTAAACGGTAGATCTCTTATGACACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGCTGCGGTTCTCTGCCGCGACGGCTGGGTTc                                                   >  2:410146/1‑250 (MQ=255)                                                                                                                                                                                                      ttCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCAAAGGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGACGGCTGGTTTCAACTTCg                                            <  1:284033/250‑1 (MQ=255)                                                                                                                                                                                                             ggAAGTTCTTCTGGAAGGATGCGCAGCCTGCCAGACATCCCCAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGACGGCTGGTTTCAACTTCGGGTAATTc                                    <  1:157232/251‑1 (MQ=255)                                                                                                                                                                                                                aGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCCAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGACGGCTGGTTTCAACTTCGGGTAATTCAgg                                 <  1:52432/251‑1 (MQ=255)                                                                                                                                                                                                                         ggAAGGATGCGCAGGCTGCCAGACATCCCCAAATGACTGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGACGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCAtc                        >  2:166682/1‑251 (MQ=255)                                                                                                                                                                                                                                   gcAGGCTGCCAGACATCCCCAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATg                                                                    >  2:288200/1‑197 (MQ=255)                                                                                                                                                                                                                                   gcAGGCTGCCAGACATCCCCAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATg                                                                    <  1:288200/197‑1 (MQ=255)                                                                                                                                                                                                                                          gCCAGACATCCCCAAATGAAGGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCGTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCTCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGACGGCTGGTTTCAACTTCGGGTAATTCAGTCATAGCATCTCCAGGAATGAACAGAt       <  2:162658/251‑1 (MQ=255)                                                                                                                                                                                                                                          cccaaaCATCCCCAAATGACTGATCACCCTGCCCTGAGCCAGCTCCAGCGTCAGATATTTAGCCCGCCTCTGCACGCTAAGCACTGGTTGCTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCCCCACCGCATGAAGAATGGTTGCACCAACGAGAAGCGGGTCTATGCCTCGACGGCTGGTTTCAACTTCGGGTAAATCAGGCATAGCATCTCCAGGAATGAACAGAt       <  2:385339/246‑1 (MQ=255)                                                                                                                                                                                                                                               aCATCCCCAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGACGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCaat  >  2:267725/1‑249 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    ATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGGGTCAGCACATTATGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTGCCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGT  >  NC_000913/3810692‑3811180

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