Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,132,697	G→A	E28K (GAG→AAG)	metF →	5,10‑methylenetetrahydrofolate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,132,697	0	G	A	100.0%	167.0 / NA	46	E28K (GAG→AAG)	metF	5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/19);  total (27/19)

CCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGCACAACACAACATATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGA  >  NC_000913/4132478‑4132937                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                 cccttctttttCCTCACCATCATTGTTCATTTTTCGGTTGTCGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACTGATAGATGTGCACACCACAACATATAACTACAAGCTATTGATGAGGTAAGGTATGATCTTTTTTCACGCCATCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATgg                                                                                                                                                                                                                    <  2:184009/244‑1 (MQ=255)        ttttCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGCACAACACAACATATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGAc                                                                                                                                                                                                             <  1:236893/250‑1 (MQ=255)          ttCCTCACCATCATTGGTCATTTGTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGCACAACACAACATATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGAccc                                                                                                                                                                                                           <  2:509130/250‑1 (MQ=255)             cTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGCACAACACAACATATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCtgtg                                                                                                                                                                                                       >  1:353698/1‑251 (MQ=255)               cACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGCACAACACAACATATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTCACGCCACCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGACGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTgg                                                                                                                                                                                                      >  1:113134/1‑250 (MQ=255)                         ggTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGCACAACACAACATATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCatcgc                                                                                                                                                                                           >  1:421909/1‑250 (MQ=255)                          gTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGCACAACACAACATATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCatcgat                                                                                                                                                                                          >  1:81207/1‑251 (MQ=255)                                               ttCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGCACAACACAACATATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCg                                                                                                                                                                     >  1:506976/1‑251 (MQ=255)                                                         ccTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGCACAACACAACATATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAgtatgtat                                                                                                                                                           >  1:427442/1‑251 (MQ=255)                                                         ccTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGCACAACACAACATATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTAt                                                                                                                                                           <  2:135175/251‑1 (MQ=255)                                                              aTCTTTACATCTGGACGTCTAAACGGATAGATGTGCACAACACAACATATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTAGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTg                                                                                                                                                      >  2:6356/1‑251 (MQ=255)                                                                         tGGACGTCTAAACGGATAGATGTGCACAACACAACATATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGgcgc                                                                                                                                           >  2:391099/1‑251 (MQ=255)                                                                                      ggATAGATGTGCACAACACAACATATCACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGCTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCg                                                                                                                                >  2:206181/1‑249 (MQ=255)                                                                                            aTGTGCACAACACAACATATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCatcgat                                                                                                                                                                                          >  2:109883/1‑185 (MQ=255)                                                                                            aTGTGCACAACACAACATATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCatcgat                                                                                                                                                                                          <  1:109883/185‑1 (MQ=255)                                                                                                     acacaacaTATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGcac                                                                                                               <  2:6126/251‑1 (MQ=255)                                                                                                      cacaacaTATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGcaca                                                                                                              <  1:69236/251‑1 (MQ=255)                                                                                                      cacaacaTACAACAACAAGCGACTGCTGAGGGAAGGTATGAGGTTTTTTCACGCCAGCCAGAGGGATGACCTGAATCAGAGCCTGGCAGAAGTCCGGGGGCAGATTAACGATTCGTTCAAGTTTTACGCGGCGCGGACAAGTGAAATGGAGCAGACCCTGTGGAACTCCATGGATCGCCTTAGCAGCCTGCAAAAGAAGTTTGTATCGGCGACCTCTGGCGCGAACTCCGGCGAGCGCAACAGTACGcacc                                                                                                              >  2:77878/1‑250 (MQ=255)                                                                                                        caacaTATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGc                                                                                                            >  2:389028/1‑251 (MQ=255)                                                                                                          aganatAACTACAAGCGATTGATGAGGTAAGTTATGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATGGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTTTCGTTGACCTATGGCGCGAACTCCGGCGGGCGCGACCGTACGCACAGCat                                                                                                          <  2:480416/249‑1 (MQ=255)                                                                                                                    aCAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGc                                                                                                >  1:119612/1‑251 (MQ=255)                                                                                                                          gATTGATGAGGTAAGGTATGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTaa                                                                                           >  1:361071/1‑250 (MQ=255)                                                                                                                             tgatgaGGTAAGGTATGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTaa                                                                                           >  1:456738/1‑247 (MQ=255)                                                                                                                              gatgaGGTAAGGTATGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTCTGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTaa                                                                                           >  1:479320/1‑246 (MQ=255)                                                                                                                               atgaGGTAAGGTATGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGTGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTaa                                                                                           >  1:148227/1‑245 (MQ=255)                                                                                                                               atgaGGTAAGGTATGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCACCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCg                                                                                     <  1:181782/251‑1 (MQ=255)                                                                                                                                     taaggtaTGAGCTTTTTTCACGCCTGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTACCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTgg                                                                               >  2:244384/1‑251 (MQ=255)                                                                                                                                          taTGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTgg                                                                          >  1:16392/1‑251 (MQ=255)                                                                                                                                           aTGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTgg                                                                          >  1:494474/1‑250 (MQ=255)                                                                                                                                            tGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCAtt                                                                                             <  1:325268/230‑1 (MQ=255)                                                                                                                                            tGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCAtt                                                                                             >  2:325268/1‑230 (MQ=255)                                                                                                                                              aGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGc                                                                      <  1:389028/251‑1 (MQ=255)                                                                                                                                                     ttCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCAcc                                                                <  2:353698/250‑1 (MQ=255)                                                                                                                                                                cagcGGGATGCCTTGAAGCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGAGCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTg                                                    <  2:113134/251‑1 (MQ=255)                                                                                                                                                                               aaTCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGANGcc                                     >  2:470289/1‑251 (MQ=255)                                                                                                                                                                                   agagCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCcgac                                 <  1:153242/251‑1 (MQ=255)                                                                                                                                                                                        cTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCt                            >  1:314774/1‑251 (MQ=255)                                                                                                                                                                                         tGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCTGTGAc                                                                                                                                                    <  2:214225/130‑1 (MQ=255)                                                                                                                                                                                         tGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGAc                                                                                                                                                    >  1:214225/1‑130 (MQ=255)                                                                                                                                                                                             aGAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTgcgc                        >  1:225016/1‑250 (MQ=255)                                                                                                                                                                                              gAAGTCCAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCAc                      <  1:369856/251‑1 (MQ=255)                                                                                                                                                                                                   ccAGGGGCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTg                 <  1:424228/251‑1 (MQ=255)                                                                                                                                                                                                      ggggCAGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCAc              <  2:271754/251‑1 (MQ=255)                                                                                                                                                                                                           aGATTAACGTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGac         <  2:81207/251‑1 (MQ=255)                                                                                                                                                                                                                   gTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGa  >  2:312773/1‑250 (MQ=255)                                                                                                                                                                                                                   gTTTCGTTCAAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGa  >  1:114381/1‑250 (MQ=255)                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                 CCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTCCTTCATCTTTACATCTGGACGTCTAAACGGATAGATGTGCACAACACAACATATAACTACAAGCGATTGATGAGGTAAGGTATGAGCTTTTTTCACGCCAGCCAGCGGGATGCCCTGAATCAGAGCCTGGCAGAAGTCCAGGGGCAGATTAACGTTTCGTTCGAGTTTTTCCCGCCGCGTACCAGTGAAATGGAGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACGAGCTGCGCACCATTGCACGCGACTACTGGA  >  NC_000913/4132478‑4132937

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