Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,810,853	G→A	G100G (GGC→GGT)	mutM ←	DNA‑formamidopyrimidine glycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,810,853	0	G	A	96.0%	155.7 / NA	50	G100G (GGC→GGT)	mutM	DNA‑formamidopyrimidine glycosylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (24/24);  minor base C (1/1);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ATGGATCCCCGCCGCAAACAGTGATTCGCTGGCATAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGGGTCAGCACATTATGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTGCCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGC  >  NC_000913/3810607‑3811077                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                 atggatCCCCGACGCAAACAGAGATTCGCAGGAATAGATATTCCCTACCCCTCCCACCAGCTTGTTATCCATCAGCCACGGTTTACTCGCCGTCTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTACGCGGCTCCGGTCCAAGATGAGCCACAACGTTGTGCCCTTCCAGCTCTTTGGCCCCAAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCCCTTTACCGt                                                                                                                                                                                                                              >  2:257528/1‑251 (MQ=255)      cccccGCCGCAAACAGTGATTCGCTGGCATAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGGCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCtc                                                                                                                                                                                                                         <  2:294234/250‑1 (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                                                                                                                                                                          >  2:35756/1‑230 (MQ=255)                                                                          aGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGAGCCAGCACGTTGTGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTACCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTgg                                                                                                                                                                          <  1:35756/230‑1 (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)                                                                                                                                                                                    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)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                 ATGGATCCCCGCCGCAAACAGTGATTCGCTGGCATAGATATTCCCTACCCCTACCACCAGCTTGTTATCCATCAGCCACGGTTTAATCGCCGTTTTTTTCTTCGCGCACTTCTGATGCAGATACTCACCATTGAAATCGTCGCTAAGCGGCTCCGGTCCAAGATGGGTCAGCACATTATGCCCTTCCAGCTCTTTGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTGCCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGC  >  NC_000913/3810607‑3811077

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: ‑