Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,626,565	G→T	V218V (GTG→GTT)	radA →	DNA recombination protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,626,565	0	G	T	100.0%	224.7 / NA	72	V218V (GTG→GTT)	radA	DNA recombination protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (30/42);  total (30/42)

GCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGC‑TGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTC  >  NC_000913/4626323‑4626807                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    gCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGt                                                                                                                                                                                                                                             <  1:557203/251‑1 (MQ=255)   aaTGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTaa                                                                                                                                                                                                                                           <  1:363720/251‑1 (MQ=255)   aaTGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTa                                                                                                                                                                                                                                            >  1:198512/1‑250 (MQ=255)              aTCGCCTTGGCCTGCAGACTGACAATCTTAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCTTGATTGCCGAAGAAGAGCAACTGAAGCTGATGGTAATTGACTCGATCCAGGTGATTCATATGGCTGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGTTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATggc                                                                                                                                                                                                                                <  1:606464/251‑2 (MQ=255)                       gCCTTCCGACTGATAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCGGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGCTCCAGGTGATGCATATGGCGGATTTACAGTCATCGCCTGGCAGCGTGGCGGAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAATATGGCTCgctggca                                                                                                                                                                                                                       <  2:308959/251‑2 (MQ=255)                        ccTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCgctggcgg                                                                                                                                                                                                                      <  1:360040/251‑3 (MQ=255)                          tgcagACTAACAAGCTCAATATGTTGTCGGATACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGACGAGCAAACGAAGATGATGGAAATTGACTCGATCCAGGTGATGCATATGGGGGATGTACAGTCATCGCCTGGCAGCGTGGCGCTGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGATGGGG‑CGATTGT‑CCTAGTTGGACCCGTAACCAACGATGGGtggctggcgggg                                                                                                                                                                                                                    >  2:355585/5‑247 (MQ=255)                            ccGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGccc                                                                                                                                                                                                                 >  1:558923/1‑251 (MQ=255)                                  gACAATCTCAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAAACAAAGATGGCTCGCTGGC‑GGGCCCTAAAGt                                                                                                                                                                                                           >  1:476342/1‑251 (MQ=255)                                          cAATATGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTg                                                                                                                                                                                                       <  1:492655/247‑1 (MQ=255)                                              atGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACt                                                                                                                                                                                                <  1:431334/250‑1 (MQ=255)                                              atGTTGTCGGAAACCAGCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGGGCTGGAACCCTg                                                                                                                                                                                               >  1:578348/1‑251 (MQ=255)                                                gttgtCTGAAACCAGCATCGAGCAGATCTGCCGGATGGCCGAAGAAGAGCATCCGAAGCGGATGTGAATTGACTCGAGCCAGGTGATGCATATGGCGGATGTACAGGTATCGCCTGGCATCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGATGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTa                                                                                                                                                                                             <  2:328378/251‑1 (MQ=255)                                                      ggAAACCAGCATCGAGCAGATATGCATGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGGGGCGCAAGGGCGTGAAACAGCGGCATATCTGACGCGCTTCGCCAAAACGCGCGGGGGGG‑CGATTGT‑CAAGGGTGGGCACGTAAACAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGAc                                                                                                                                                                                        >  1:365655/1‑250 (MQ=255)                                                              gCATCGAGCAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTg                                                                                                                                                                               <  1:105339/251‑1 (MQ=255)                                                                   gagcAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCAGATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCtttt                                                                                                                                                                          <  1:26883/248‑1 (MQ=255)                                                                             gCCTGATTGCCGAAGAAGAGCAACCGATGCTGATAGTATTTGTCTAGATTCTGTTGCTGCATATGGCGTATGGTCAGTCTGCGCCTGGCAGCGTTGCGCCGGTGCGTGTAACGGCGTCTTATCTGACTCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTATCCAAAGATGGTTCGCTGGC‑GGGGTATGAAGTGCTGGAACACTGTATCTACTGTTCGGTGCttttttctga                                                                                                                                                                    <  2:95665/247‑7 (MQ=255)                                                                             gCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGGGATGCATATGGCGGATGTACAGGCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGc                                                                                                                                                                   <  2:371001/248‑1 (MQ=255)                                                                                tGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGAGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCCGGAACACTGTATCGCCTGTTCCGTGCTTTTGGATGGCGATGcc                                                                                                                                                             >  1:285309/1‑251 (MQ=255)                                                                                        gaagaaGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCg                                                                                                                                                     >  1:286672/1‑251 (MQ=255)                                                                                          agaagaTCAGCCGAAGTTGGTGGTAATTGACTCGTGGCAGGTGATGCATATGGCGTATGTACTGGCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAACCGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGta                                                                                                                                                   <  2:223582/251‑2 (MQ=255)                                                                                          agaagaGCAGCCGAAGCTGTTGGTAATTGACGCGTTCCAGGTTAGGCATATGGCGGAGGTCTAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGTCGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCGCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGGTCGGTGCTTTTGGATGGCGATGCCGTCTCCCGtt                                                                                                                                                   <  2:6351/251‑1 (MQ=255)                                                                                             agagTAACCGAAGCGGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATTTTCAGTCAGCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGTCGCGCGTCGCCAAAACGCGCGGTGTGG‑CTATTGG‑CATGGTTGGACACGTAACCAAATATGGCTCGCTGGC‑GGGCACTAAAGTGCGGGAACACTGTACCGACTGTGCGGTGCTGTTGGATGGCGACGCTGACTCCAGTTTTc                                                                                                                                                <  2:388471/251‑1 (MQ=255)                                                                                                  aaCCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCAct                                                                                                                                           <  1:214730/251‑2 (MQ=255)                                                                                                     cGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAATATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCActtta                                                                                                                                        <  1:394882/251‑5 (MQ=255)                                                                                                      gAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCActttac                                                                                                                                       <  1:633276/251‑6 (MQ=255)                                                                                                            gATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGc                                                                                                                                  <  1:505656/250‑1 (MQ=255)                                                                                                              gggtaatggacgtgATCCAGGTGATGCATATGGCGGAGGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGGGTGAAACGGTGGCTTATCTGACGCGCTACGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCGGTC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCAt                                                                                                                               <  2:155283/238‑1 (MQ=255)                                                                                                                  aatggatgcgagccAGGTGATGCAAATGGCGGATGGACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCGGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGTCTCCCGTTTTCGCACGTTAGGCAGCCATaaaa                                                                                                                           <  2:32440/239‑1 (MQ=255)                                                                                                                   aTTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATaaaaa                                                                                                                          <  1:509558/251‑1 (MQ=255)                                                                                                                   aTTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAAAACTGTATCGACTGTTCGGGGCTTTTGGATGGCGGTGCCGACTCCCGTTTTCGCACTTTACGCAGCCATaaaa                                                                                                                           >  1:508403/1‑250 (MQ=255)                                                                                                                     tGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAAc                                                                                                                         <  1:277807/250‑1 (MQ=255)                                                                                                                         ccGATCCGGGTGATGCAAATGGCTGATTTACAGTCATCGCCGGGCGGCGTGGCGCAGGTGCGTGAAACGGCGGTTTATCTGACGCGTGTCGTCAAAACGCGCGGTGTGG‑CGATTTG‑CATGGTTGGACACGTAACCAAAGCTGGTCGGCGGGCGTGGTC‑TAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCtt                                                                                                                    <  2:376233/250‑1 (MQ=255)                                                                                                                               cAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACCCTGTATCGACTGTTCGGGGCTTTTGGATGGCGGTGCCGACTCCCGTTTTCGCACTTTACGCCGCCCTAAAAACCGCTTCGgcgg                                                                                                              >  1:10935/1‑250 (MQ=255)                                                                                                                                  gggaggCATATGGTGTATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAACCGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGG‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGGGCGGGAACACTGTATCGACTGTTCGGGGCTTTTGGATGGCGATGCCGACTCCCGTGTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGt                                                                                                           <  2:654811/246‑1 (MQ=255)                                                                                                                                   tGATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTTTCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CCTGGTTGGACACGTAACCCAAGATGGCTCGCTTGC‑GGGCCCTAAAGTGCTGGAACACAGTATTCACTGTTCGGGGCTTTTGGATGGCGATGCCGCCTCCCGTTTTCGCACTTTTCGCCGCAATTAAAACCGCTTctgcacgggg                                                                                                          >  2:99834/1‑242 (MQ=255)                                                                                                                                    gATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTg                                                                                                          <  1:202487/250‑1 (MQ=255)                                                                                                                                    gATGCATATGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGa                                                                                                         <  1:6086/251‑1 (MQ=255)                                                                                                                                         atatGGCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGAAACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGTCTGTTCGGTGCTTTTGGATGGCGATGCCGAATCCCGTTTTTGCACTTTAAGCAGCCCTAAAAAACGCTTCGGCGAGGTGAATg                                                                                                      >  1:300109/1‑249 (MQ=255)                                                                                                                                              gCGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGCTGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGCAAACTGTATCGGCTGTTCGGTGCTTTTGGATGGCGCTGCCGACTCCCGTTTTCGCACTTTACGCAGGCATAAAAAACGCTTCGGCGCGGTGAATGAGCTggg                                                                                               >  1:410471/1‑251 (MQ=255)                                                                                                                                               cGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACTGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATTGCTCgctggc                                                                                                                                                                                                                        <  2:446261/130‑1 (MQ=255)                                                                                                                                               cGGATGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCCAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCgctggc                                                                                                                                                                                                                        >  1:446261/1‑130 (MQ=255)                                                                                                                                                  aTGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTc                                                                                           >  1:72568/1‑251 (MQ=255)                                                                                                                                                  aTGTACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGGGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTc                                                                                           >  1:527580/1‑251 (MQ=255)                                                                                                                                                     tACAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTc                                                                                        >  1:492594/1‑251 (MQ=255)                                                                                                                                                      aCAGTCATCGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTc                                                                                        <  1:119060/250‑1 (MQ=255)                                                                                                                                                              cGCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGCTGCCGACTCCCGTTTTCGCACTTTACGCCGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGAcc                                                                               >  1:221082/1‑251 (MQ=255)                                                                                                                                                               gCCTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTTGACACGTAACCAAAGATAGCTCGCTGGC‑GGGACCTAAAGTGCTGGAACACTGTATCGGATGTTCGGTGCTTTTGGAAGGCGCGGCCGCCTCCCGTTTTCGCCCTTTAAGCCGCCATAAAAACCGCTTCGGCGCGGGGCAAGGGCTGGGCGTCTTcactagtaccc                                                                              >  2:458050/1‑241 (MQ=255)                                                                                                                                                               gCCTGGCAGCGGGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGA‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGGACACTGGATCGACTGTTCGGTGGTTTTGGATGGCGATGCCGACCCCCGTTTTCGCCCTTTACCCCGCAAAAAAAACACCTTCGGCGCGGTGAATGAGCAGGGCGTCTTccgcatggccc                                                                              >  2:348741/1‑241 (MQ=255)                                                                                                                                                                ccTGGCAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGGGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGGGCTGGGACACTGTATCGACTGGTCGGGGCTTTTGGATGGCAATGCCGACTCCCGTTTTCGCACTTTACGCAACCATAAAAACCGCTTCGGCGCGGGGGATGAGGTGGGCGTCTTCGCGAttgccg                                                                              >  2:545849/1‑245 (MQ=255)                                                                                                                                                                   ggcagggggtggagggggCGTGAAACGGCGGCTTATCTGACGAGATTCCGCAAAACGCGTGGTGTGG‑GGAT‑GTGAATGGTTGGACACGTAACCAAAGATGTCTCGCTGGT‑GGGACCTAAAGTGCTGGAAAAATGTATCGACTGTTCTGTGCTTTTTGAAGGCGATGGCGACTCCCGTTTGCGTACTTTACGCAGCAATAAAACCCGCTGCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGCCCGaaaa                                                                          <  2:265528/234‑3 (MQ=255)                                                                                                                                                                    gcaGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAg                                                                         <  1:120830/249‑1 (MQ=255)                                                                                                                                                                                  gTGCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTGGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGt                                                           <  1:296828/251‑1 (MQ=255)                                                                                                                                                                                    gCGTGAAACGGCGGCTTATCTGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCa                                                         <  1:545849/251‑1 (MQ=255)                                                                                                                                                                                                 gctattctGACGCGCTTCGCCAAAACGCGCGGGGTGG‑CGATAGT‑CATGGGTGGACAAGTAACCAAAGATGGCTCGCTGTC‑GGGCCCTAAGGTGCTGGAACACTGTATCGACTGGTCGGTGCTTTTGGAGGGCGATGCCGACTCCCGTTTTCGCACTGTACGCAGCCATAAAAACCGCTGCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGc                                            <  2:6321/244‑1 (MQ=255)                                                                                                                                                                                                    tatctGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCa                                           <  1:99834/245‑1 (MQ=255)                                                                                                                                                                                                     atctGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCa                                                                                                                                <  2:371176/160‑1 (MQ=255)                                                                                                                                                                                                     atctGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCa                                                                                                                                >  1:371176/4‑163 (MQ=255)                                                                                                                                                                                                     atctGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAAtt                                        >  1:338213/4‑251 (MQ=255)                                                                                                                                                                                                      tctGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAAttt                                       >  1:162822/3‑251 (MQ=255)                                                                                                                                                                                                        tGACGCGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTc                                     >  1:648762/1‑251 (MQ=255)                                                                                                                                                                                                              gCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGt                               <  1:257533/251‑1 (MQ=255)                                                                                                                                                                                                                gtggCCAAAACGCGCGG‑GTGGGCGATTGG‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCGTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTGCTTAAGTcg                             <  2:545962/248‑1 (MQ=255)                                                                                                                                                                                                                    ccAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCgcgga                         <  1:498873/251‑2 (MQ=255)                                                                                                                                                                                                                    ccAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGGGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCgcggg                         >  1:453332/1‑250 (MQ=255)                                                                                                                                                                                                                        aactcgtggTGTGG‑AGATTGG‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCATAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGAGGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATg                     <  2:241219/244‑1 (MQ=255)                                                                                                                                                                                                                                  gtgG‑CGATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCt            <  1:140892/250‑1 (MQ=255)                                                                                                                                                                                                                                  gtgG‑CGATTGG‑CATGGTTGGACACGTAACCAAAGATGGCTCGCGGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCAACTGTTCGGGTCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTTCGCAGCCATAAAAACCGCTTCGGCGCGGGGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGGGAAGACAGCAAACCTTCGGCAATTTTCTTAAGTGGCGGGGCTGAGGTGACCTc           >  2:87298/1‑251 (MQ=255)                                                                                                                                                                                                                                     ggtgATTTT‑CATGGTTGGATACGTACCCAAAGACGCATCGCTGGC‑GGGCCCTAATGTGCTGGAACACTGTAGCGACGGTTCGGTGATTGTGGATGGCGATGCAGTCTCCCTTCTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGTGCTGTGAATGAGCTTGGCTTCTTCGCGATGACCGAACCGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAAGTTTCTTAAGGCGCGTAGATTAAGTGACCTcc          <  2:406871/247‑1 (MQ=255)                                                                                                                                                                                                                                        gATTGT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGGCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCggcc      >  1:363401/1‑249 (MQ=255)                                                                                                                                                                                                                                          ttGT‑CATGGTTGGTCACGTAACCAAAGCTGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCGTTTGGATGGCGATGCCGACTCCCGTATTCGCACTTTACGCAGCCATACAAACCACTTCGGCGCGGTGAATGAGCTCGGCTTATTCGCGTTGACAGCAAAGGGGCTGCGGGAGGGCAGCGACCATCAGGCAGTTTTCTGAAGTGGGGGAGATGAAGAGAACTCCggagac    >  2:50588/1‑247 (MQ=255)                                                                                                                                                                                                                                            gT‑CATGGTTGGACACGTAACCAAAGATGGCTCGCTGGC‑GGGCCCTAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACTTTACGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGAGATGAAGTGACCTCCGGCAGCTc  >  1:265148/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    GCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGG‑CGATTGT‑CATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGC‑TGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATGCCGACTCCCGTTTTCGCACCTTGCGCAGCCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTC  >  NC_000913/4626323‑4626807

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