Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,204,617	G→A	intergenic (‑28/‑77)	folB ← / → plsY	dihydroneopterin aldolase/putative glycerol‑3‑phosphate acyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,204,617	0	G	A	100.0%	219.3 / NA	66	intergenic (‑28/‑77)	folB/plsY	dihydroneopterin aldolase/putative glycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/39);  total (27/39)

TCCACCAGCGCAAAACGCGCCCCCTCGACGTGGCTGACCACCGTTTCTGCAATGTCAGCGTAACTGAGGCAATCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCGGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCCGCGAACCAGCGGCTCCGGCAATCCAGGCGCAACCAATGTGTTACGTATCGGTGGCAAGGGAGCAGCC  >  NC_000913/3204372‑3204858                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  tCCACCAGCGCAAAACGCGCCCCCTCGACGTGGCTGACCACCGTTTCTGCAATGTCAGCGTAACTGAGGCAATCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATAc                                                                                                                                                                                                                                              >  2:103515/1‑251 (MQ=255)  ccaccaGCGCAAAACGCGCCCCCTCGACGTGGCTGACCACCGTTTCTGCAATGTCAGCGTAACTGAGGCAATCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATAcc                                                                                                                                                                                                                                             >  1:497314/1‑251 (MQ=255)      cagcgaaacacgagaccccTCGACGTGGCTGACCACCGTTTCTGCAAGGTCAGCGTAACTGAGGCAATCCGCAACATCATCACTTTTCGCCGCTTTACGGTTAACCCACGCCATTTCGATATCGACCACTAACTTCTGTTCGGTGGTCTGTTCCCAGTCGTAAACAACAATAGTGGTGATTACCGCAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACtt                                                                                                                                                                                                                                         >  2:513975/16‑251 (MQ=255)        gctcAAAACGCGCCCCCACGACGTGGCTGACCACCGTTTCTCCAATGTCAGCGTAACTGAGGCAACCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTGTTGGCTTACCAGATACCACATcc                                                                                                                                                                                                                                       <  2:20911/248‑1 (MQ=255)           cAAAACGCGCCCCCTCGACGTGGCTGACCACCGTTTCTGCAATGTCAGCGTAACTGAGGCAATCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGc                                                                                                                                                                                                                                    >  2:153805/1‑251 (MQ=255)           cAAAACGCGCCCCCTCGACGGGGCTGACCACCGTCTCTGCAATAACAGCGTACCTGAGGCAATCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACCCCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACACTATCCATCACGTCCTTCCTGCTTTTTGGCTAACCAGATACCACTTCCGGc                                                                                                                                                                                                                                    >  2:362047/1‑251 (MQ=255)              aaCGCGACCACTCGACGTGGGTGACCACCGTTTCTGCAATATCAGCGTAACTGAGGCAATCCGCCACCTCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGGACACTAACTTCTGTTATCTGGTCGGTTCCCAGTCGTAAACACCAATAGGGTTGGTTACCGACAGTTGCTCTATAAATACAATATCCATCACGACCTGCCTTCTTTTTGGCTAACCAGATACCACTTCAGGCGa                                                                                                                                                                                                                                  >  2:117320/1‑250 (MQ=255)                   gcCCCCTCGACGTGGCTGACCACCGTTTCTGCAATGTCAGCGTAACTGAGGCAATCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAAtgtg                                                                                                                                                                                                                            <  1:454266/251‑1 (MQ=255)                     ccacTCGACGGCGCTGACCACCGTTTCTGCAATGGCAGCTAAACTGAGGCAATCCGCCACATCATCACTTTTCTCCGCTTTACGGTTAGCCCACGCCATTTCGATATCGAACACTAACTGCTTTTCGATGGTATGTTCCCAGTCGTAAACTCCAATAGTGGTGATTACCGAAAGTTGCTCTAGAAATACAATATCCATCACGTCATGGCTGCTTGTTGGCTAATCAGATTCCCCTTCCGGCGAAACGTGCg                                                                                                                                                                                                                          <  2:48151/248‑1 (MQ=255)                      cccTCGACGTGGCTGACCACCGTTTCTGCAATGTCAGCGTAACTGAGGCAATCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGt                                                                                                                                                                                                                         >  1:532717/1‑251 (MQ=255)                           gACGTGGCTGCCCACCGTTTCTGCAATGTCAGCGTACCTGAGGCAATCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAAGATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGgattat                                                                                                                                                                                                                    <  2:555237/251‑1 (MQ=255)                                ggCTGACCCCCGTTTCTGCAATGTCAGCGTAACTGAGGAAATCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCcac                                                                                                                                                                                                                <  1:392120/250‑1 (MQ=255)                                   ggccccccGTTTCTGCAATGTCAGCGTACCTGAGGAAATCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGAt                                                                                                                                                                                                            <  2:38488/245‑1 (MQ=255)                                    gACCACCGTTTCTGCAATGTCAGCGTAACTGAGGCAATCCGCCACATCATCACTTTTCGTCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGAt                                                                                                                                                                                                            <  1:54259/250‑1 (MQ=255)                                        acgggttgtgggaTGTCAGCGTAACTGGGGCGGTCCGCGACGGCAGGACTTTTCGCCGCGTTGCGTGTAGGCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGTAAGTTGCTCTATAAATACAATTTCCATCACGTCCTGCCTGCTTTTTTGTTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATc                                                                                                                                                                                                       <  1:394513/239‑1 (MQ=255)                                        accGTTTCTGCAATGTCAGCGTAACTGAGGCAATCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATc                                                                                                                                                                                                       <  1:37016/251‑1 (MQ=255)                                           gttccTGCAATGTCAGCGTAACTGAGGCAATCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGtt                                                                                                                                                                                                    <  2:167317/247‑1 (MQ=255)                                                 gCAATGTCAGCGTAACTGAGGCAATCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGaa                                                                                                                                                                                                 <  1:8886/248‑1 (MQ=255)                                                        tagctacacgtgggaaaTCCGCCAGATCATCACTGTTCGCCTCTTTACGTGTATCCCTTGCCAGTGCGATATCGAACTCTAACTTCTGGTCGAAGGTCTGTTCCCAGTTGTAACCACCAATAGTGGTGATTACCGAAAGTTGCTCTAGAAATACAATCTCCATCACGGCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCTGCGAAATGTGCGTATTATCCACGGATTCATCTTTGAACACGaagtg                                                                                                                                                                                        <  2:131472/235‑4 (MQ=255)                                                          gCGTAACTGAGCCAATCCGCCACATCATCACGTTTCGTCGCTTTACGGTTAGCCCCCGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATTTGCGTATTATCCACAGATTCATCGTTGAACACGAAtttt                                                                                                                                                                                       <  2:497314/249‑1 (MQ=255)                                                           cGTAACTGAGGCAATCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCCCAGATTCATCGTTGAACACGAATTTTCa                                                                                                                                                                                     >  2:347454/1‑250 (MQ=255)                                                                   aGGCAATCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAAc                                                                                                                                                                            >  2:183368/1‑251 (MQ=255)                                                                       aaTCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCt                                                                                                                                                                        >  1:291290/1‑251 (MQ=255)                                                                        aTCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCtt                                                                                                                                                                       <  1:487237/251‑1 (MQ=255)                                                                          ccGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGCAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTAt                                                                                                                                                                     >  1:555355/1‑251 (MQ=255)                                                                              cacaTCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGt                                                                                                                                                                 >  2:154412/1‑251 (MQ=255)                                                                                      caCTATTCGCCGCTTTACGGTTATCCCACGCCAATTCGATATCGAACACTTAGTTCTGTTCGATGGTCTGTTCACAGTCGTAAACGCCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAAGATCCGTCACTTGCTGCCTGTTTTTTTGCTAATCAGATACCAGTCGCGGCGAAATGTGCGTACTATCCACAGGGTCATCGTTGAAGGCGAGTTTTCAAAACGGATCAGTTTATGAGTGCaagcg                                                                                                                                                          <  2:88746/250‑4 (MQ=255)                                                                                          tttCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGcc                                                                                                                                                      <  1:235978/250‑1 (MQ=255)                                                                                              cccgcTTTACGGTTATCCCGCGCCATGTCGATATCGAACACTAACTTCCGTTCGATGGTCTGTTCCCAGTCGTAAACTCCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGATTGCAATCGCGCCTGGaa                                                                                                                                                 <  2:280843/250‑1 (MQ=255)                                                                                               ccgcTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAAt                                                                                                                                                <  1:10832/251‑1 (MQ=255)                                                                                                cgcTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATg                                                                                                                                               <  1:281431/251‑1 (MQ=255)                                                                                                  cTTTACGGATATCCCACGCCATTTCGCTATCGACCACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGCAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCAACAGATTCATCGTTGAACACGACTTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGCATGAt                                                                                                                                             >  2:561583/1‑251 (MQ=255)                                                                                                         gTTATGCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGACCTATAACTACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAATCAGAAACCACTTCCGGCGAAATGTGCGTATTATCCACAGCTTCATCGTTCAACACGAACTGTCAACACGGAACAGCTTATGAGTGCAATCGCACCTGGAATGATCCTCaac                                                                                                                                      >  1:521200/1‑249 (MQ=255)                                                                                                               ccACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCCGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAAGATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGttc                                                                                                                                <  2:154807/251‑3 (MQ=255)                                                                                                                  cGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACctc                                                                                                                             <  1:9219/251‑1 (MQ=255)                                                                                                                           gATATCGAACACTAACTTCTGGTCGATGGTCTGTTCCCAGGCGGAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTcc                                                                                                                    <  2:30420/251‑1 (MQ=255)                                                                                                                           gATATCGAACACTAACTGCTGTTCGATGGTCTGTTCCCAGTCGGAAACACCAATAGTGGTGATGACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAAGTTTCAAAACGGAACAGCTTATGAGGGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTcc                                                                                                                    <  2:235254/251‑1 (MQ=255)                                                                                                                                   acacTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTcc                                                                                                              >  1:7521/1‑249 (MQ=255)                                                                                                                                        aaCTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACATCAATAGTGGTGCTTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCGCCATTTCCAgggc                                                                                                         <  2:25070/249‑4 (MQ=255)                                                                                                                                                tgcGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGt                                                                                                <  1:542741/248‑1 (MQ=255)                                                                                                                                                            ttCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACCGCTTATGATTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTtgtgg                                                                                   >  2:12071/1‑250 (MQ=255)                                                                                                                                                                    cGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTcc                                                                                                                    <  1:384713/210‑1 (MQ=255)                                                                                                                                                                    cGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTcc                                                                                                                    >  2:384713/1‑210 (MQ=255)                                                                                                                                                                     gTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGccc                                                                          <  1:455659/251‑1 (MQ=255)                                                                                                                                                                      tAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCg                                                                         >  1:420438/1‑251 (MQ=255)                                                                                                                                                                       aaaCACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCCTCCCGTGCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCGCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCg                                                                         <  1:117320/250‑1 (MQ=255)                                                                                                                                                                            ccAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTtg                                                                                      >  1:568632/1‑232 (MQ=255)                                                                                                                                                                            ccAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGGGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTtg                                                                                      <  2:568632/232‑1 (MQ=255)                                                                                                                                                                             cAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCcgcg                                                                  >  2:487935/1‑251 (MQ=255)                                                                                                                                                                                    ggtgATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCCGCGAACCAg                                                            <  1:153805/250‑1 (MQ=255)                                                                                                                                                                                       gATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCg                                                                                                                         >  1:442889/1‑186 (MQ=255)                                                                                                                                                                                       gATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCg                                                                                                                         <  2:442889/186‑1 (MQ=255)                                                                                                                                                                                        aTTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCCGCGAACCAGCGGCt                                                       <  1:280245/251‑1 (MQ=255)                                                                                                                                                                                          tACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCCGCGAACCAGCGGCTc                                                      <  1:320528/250‑1 (MQ=255)                                                                                                                                                                                                aaGTTGCTCTATAACTACCATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAAAAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCCGCGAACCAGCGGCTCCGGCatt                                               >  2:83768/1‑249 (MQ=255)                                                                                                                                                                                                        ctATAAATACAATATCCATCTCGTCCTGCCTGCTTTTTGGCTAACCAGAGACCACTTCCGGCGAAATGTGCGGATTATCCACAGATGCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTGCCAGCGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCCGCGAACCAGCGGCTCCGGCAATCCAgg                                          <  2:205517/248‑1 (MQ=255)                                                                                                                                                                                                              aTACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCCGCGAACCAGCGGCTCCGGCAATCCAGGCGCAACCa                                  >  1:94539/1‑250 (MQ=255)                                                                                                                                                                                                                    tatCCATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTAGTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGTAACAGTTTATGAGTGCAATCGCGCCTGTAATGATCCTCATcg                                                                                                                                     <  2:221813/145‑1 (MQ=255)                                                                                                                                                                                                                    tatCCATCACGTCCTGCCTGATTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATcg                                                                                                                                     >  1:221813/1‑145 (MQ=255)                                                                                                                                                                                                                        cATCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCCGCGAACCAGCGGCTCCGGCAATCCAGGCGCAACCAATGTGTTACGt                       >  1:124380/1‑251 (MQ=255)                                                                                                                                                                                                                          tCACGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCCGCGAACCAGCGGCTCCGGCAATCCAGGCGCAACCAATGTGTTACGTAt                     >  1:237169/1‑251 (MQ=255)                                                                                                                                                                                                                            aCGTCCTGCCTGCTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCCGCGAACCAGCGGCTCCGGCAATCCAGGCGCAACCAATGTGTTACGTATCg                   <  1:206366/251‑1 (MQ=255)                                                                                                                                                                                                                                   gcctgcTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCCGCGAACCAGCGGCTCCGGCAATCCAGgc                                         <  2:497035/222‑1 (MQ=255)                                                                                                                                                                                                                                   gcctgcTTTTTGGCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCCGCGAACCAGCGGCTCCGGCAATCCAGgc                                         >  1:497035/1‑222 (MQ=255)                                                                                                                                                                                                                                       gattgttGGATAACCAGATACCACTTCCTGCGAAATGTGCGTATTAGCCTCCGATTCATCGGTGAACACTAATTTTCAAAACGCAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCGCATCGCGTACCTCTGCGGCTCCAGTCCCAGTGCCAGGCTGGTGTGCCTCTTGTGTGGGCTGCCCGATCCGCGAACTAGCGGCTCTGGCAATCCTGGCGCAACCACTGATTTACGTATCGGtcgaaagtga        <  2:398988/246‑10 (MQ=255)                                                                                                                                                                                                                                              ggCTAACCAGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCCGCGAACCAGCGGCTCCGGCAATCCAGGCGCAACCAATGGGTGACGTATCGGTGTCAAGGGATCAGcc  <  2:397269/250‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  TCCACCAGCGCAAAACGCGCCCCCTCGACGTGGCTGACCACCGTTTCTGCAATGTCAGCGTAACTGAGGCAATCCGCCACATCATCACTTTTCGCCGCTTTACGGTTATCCCACGCCATTTCGATATCGAACACTAACTTCTGTTCGATGGTCTGTTCCCAGTCGTAAACACCAATAGTGGTGATTACCGAAAGTTGCTCTATAAATACAATATCCATCACGTCCTGCCTGCTTTTTGGCTAACCGGATACCACTTCCGGCGAAATGTGCGTATTATCCACAGATTCATCGTTGAACACGAATTTTCAAAACGGAACAGCTTATGAGTGCAATCGCGCCTGGAATGATCCTCATCGCGTACCTCTGCGGCTCCATTTCCAGTGCCATTCTGGTTTGCCGCTTGTGTGGGCTGCCCGATCCGCGAACCAGCGGCTCCGGCAATCCAGGCGCAACCAATGTGTTACGTATCGGTGGCAAGGGAGCAGCC  >  NC_000913/3204372‑3204858

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