Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,349,573	G→A	R154R (CGC→CGT)	dcuR ←	DNA‑binding transcriptional activator DcuR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,349,573	0	G	A	100.0%	201.0 / NA	59	R154R (CGC→CGT)	dcuR	DNA‑binding transcriptional activator DcuR
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/28);  total (31/28)

GAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTATTTACGGCAGGAAACACGCGAAATGTTAACCTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAGCGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTTGATCAGGTAATCCACGACACCGTAATGCAGCGAATCTTTAATGGTTGCCGCATCGGCTGCGGAGGAGATGACAATCACATCACTTTTGCAACGCGC  >  NC_000913/4349340‑4349815                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                   gagtagtgTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTATTTACGGCAGGAAACACGCGAAATGTTAACCTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTc                                                                                                                                                                                                                                   <  2:198929/251‑1 (MQ=255)  agtagtgTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTATTTACGGCAGGAAACACGCGAAATGTTAACCTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCg                                                                                                                                                                                                                                  >  1:491540/1‑251 (MQ=255)        gtTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTATTTACGGCAGGAAACACGCGAAATGTTAACCTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTgg                                                                                                                                                                                                                              >  1:302366/1‑249 (MQ=255)             gCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTATTTACGGCAGGAAACACGCGAAATGTTAACCTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAgctgct                                                                                                                                                                                                                       >  1:596147/1‑251 (MQ=255)                      cggtaacggtaaACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTATTTACGGCAGGAAACACGCGAAATGTTAACCTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCCAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTTGAt                                                                                                                                                                                                              >  2:297339/1‑251 (MQ=255)                          aacggtaaACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTATTTACGGCAGGAAAAAAGCGAAATGTTAACCTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGAAGAGGATCCTGTTCGTTGGCGCTGCTGCCGTGGATTAGc                                                                                                                                                                                                          >  2:45024/1‑251 (MQ=255)                                         gCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTATTTACGGAAGGAAAAACGCGAAATGTTAACCTTGTTGGCTAGCTCGGCGGTTGAAAATTAATAGCCCTGATTCGCGTCCATCCACTGGCACAGTGAGC‑GTAAAGCCTGCGGCGTTCAGCCTTTTGGCAAACCCAGAGGATTCTTTTCGGTGGAGCTGCTGCCGTGGGATAGCTGATCAAgatccgcc                                                                                                                                                                                           >  2:166756/1‑244 (MQ=255)                                           cccGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTATTTACGGCAGGAAACACGCGAAATGTTAACCTCGTTGGCTAGCTCGTCGGTTTAAAATTCATAGTCCTTATTCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTTGCAAACGCCGAGGATCCTGTTCGTTGGAGCTTCTGCCGt                                                                                                                                                                                                                  >  2:496807/1‑226 (MQ=255)                                           cccGTGACGCCATAAGGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTATTTACGGCAGGAAACACGCGAAATGTTAACCTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGt                                                                                                                                                                                                                  <  1:496807/226‑1 (MQ=255)                                             cGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTATTTACGGCAGGAAACACGCGAAATGTTAACCTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGt                                                                                                                                                                                       >  1:528351/1‑251 (MQ=255)                                                       tAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTATTTACGGCAGGAAACACGCGAAATGTTAACCTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTg                                                                                                                                                                             >  1:621566/1‑251 (MQ=255)                                                               aCTGGTGAACAAGATGTGGCAGGTGACCAGCCAGATGAGGTATTTACGGCAGGAAACACGCGAAATGTTAACCTCGTTGGCTAGCGCGTTGGTTGAAAATTCATAGTCCTGATGCGCGTCAAGCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGTAGCTGTTGCCGTGGATTAGCTGATCATGCTCGGCCTGGTCATAATACTGATGtttt                                                                                                                                                                      <  2:115585/250‑1 (MQ=255)                                                                cTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTATTTACGGAAGGAAACACGGGAAATGTTAACCTCGTTGACTAGCTCGACGGCTGCAAATTCATAGTCCTGATGCGAGTCAATCCCCTGGCCAAGAG‑GCAGGAAAGTCTGCGGCGTTAAGACTTTTGGCAAAAGACGAGGATCCTGTTAGTTGGAGATGTTGCCGTGGATTAGCTGATCAAGCCCGGCCTGGTCATACTAATGGTGtttttt                                                                                                                                                                    >  1:51777/1‑250 (MQ=255)                                                                  ggTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTATTTACGGCAGGAAACACGCGAAATGTTAACCTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCg                                                                                                                                                                                                                                  <  2:51864/187‑1 (MQ=255)                                                                  ggTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTATTTACGGCAGGAAACACGCGAAATGTTAACCTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCg                                                                                                                                                                                                                                  >  1:51864/1‑187 (MQ=255)                                                                        cAAGATGTGGCAGTTGACCAGCCAGATGAGGTATTTACGGCAGGAAACACGCGAAATGTTAACCTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCa                                                                                                                                                            <  1:78313/251‑1 (MQ=255)                                                                                   aGTTGACCAGCCAGATGAGGTATTTACGGCAGGAAACACGCGAAATGTTAACCTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGctgctg                                                                                                                                                                                                                      <  2:13940/182‑1 (MQ=255)                                                                                   aGTTGACCAGCCAGATGAGGTATTTACGGCAGGAAACACGCGAAATGTTAACCTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGctgctg                                                                                                                                                                                                                      >  1:13940/1‑182 (MQ=255)                                                                                       gACCAGCCAGATGAGGTATTTACGGCAGGAAACACGCGAAATGTTAACCTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTgccgcc                                                                                                                                             <  2:302270/251‑1 (MQ=255)                                                                                             ccagATGAGGTATTTACGGCAGGAAACACGCGAAATGTTAACCTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCgg                                                                                                                                       >  1:414723/1‑251 (MQ=255)                                                                                                gATGAGGTATTTACGGCAGGAAACACGCGAAATGTTAACCTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCCGCGCCCTTTTCTTTTGCCGCCAGCCGGTgg                                                                                                                                    >  2:92651/1‑250 (MQ=255)                                                                                                  tGAGGAATTTACGGCAGGAAACACGCGAAATGTTAACCTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAgc                                                                                                                                  >  2:37109/1‑251 (MQ=255)                                                                                                        aTTTACGGCAGGAAACACGCGAAATGTTAACCTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCCGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCctct                                                                                                                            >  2:468666/1‑251 (MQ=255)                                                                                                         tttACGGCAGGAAACACGCGAAATGTTAACCTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCCGCGCCATTTTCTTTTGCCGCCAACCGGTGAGCGCCNCtt                                                                                                                           >  2:595765/1‑251 (MQ=255)                                                                                                                caggacacacgtgcgattttAACCTCGTTGGCTAGCTCCACGGTTGCAAATTCATAGTCTTGATGCGCGCCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCACACGACGAGGATCCTGTCCGTTGGACCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATCATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCg                                                                                                                    >  2:87066/19‑251 (MQ=255)                                                                                                                 aGGAAACACGCGAAATGTTAACCTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCgg                                                                                                                   >  1:248823/1‑251 (MQ=255)                                                                                                                  ggAAACACGCAATATGTTAACATCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGGGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCggg                                                                                                                  <  2:629368/251‑1 (MQ=255)                                                                                                                     aaCACGCGAAATGTTAACCTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAg                                                                                                               <  2:171452/251‑1 (MQ=255)                                                                                                                        acGCGAAATGTTAACCTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCt                                                                                                            <  1:279584/251‑1 (MQ=255)                                                                                                                             aaaTGTTAACCTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAg                                                                                                       <  2:621566/251‑1 (MQ=255)                                                                                                                                     aCCTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTTGAt                                                                                               <  1:549700/251‑1 (MQ=255)                                                                                                                                      ccTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTTGATc                                                                                              <  1:144762/251‑1 (MQ=255)                                                                                                                                         cGTTGGCTAGCTAGTCGGTTGAAAATTCATAGTCCAGATGCGCGTCAATCAACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGCTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTGTTGCAGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTGGATCAgg                                                                                           >  2:484552/1‑251 (MQ=255)                                                                                                                                              gCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTTGATCAGGTAATc                                                                                      <  1:262760/251‑1 (MQ=255)                                                                                                                                                tAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTTGATCAGGTAATcc                                                                                     >  1:587293/1‑250 (MQ=255)                                                                                                                                                                aaTTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTTGATCAGGTAATCCACGACACCGTAATGc                                                                      >  2:397601/1‑249 (MQ=255)                                                                                                                                                                 aTTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTTGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCTGGAAGCCTGGAAGGGTTTGATCAGGTAATCCACGTCACCGTAATGCAGc                                                                   <  1:206699/251‑1 (MQ=255)                                                                                                                                                                 aTTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTTGATCAGGTAATCCACGACACCGTAATGCAGc                                                                   >  1:570680/1‑251 (MQ=255)                                                                                                                                                                       aGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTTGATCAGGTAATCCACGACACCGTAATGCAGCGAATCt                                                             >  1:172390/1‑251 (MQ=255)                                                                                                                                                                             gATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTTGATCAGGTAATCCACGACACCGTAATGCAGCGAATCTTTAATg                                                       >  1:476394/1‑251 (MQ=255)                                                                                                                                                                                     tCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTTGATCAGGTAATCCACGACACCGTAATGCAGCGAATCTTTAATGGTTGCCGc                                               >  2:266791/1‑251 (MQ=255)                                                                                                                                                                                          ccACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTTGATCAGGTAATCCACGACACCGTAATGCAGCGAATCTTTAATGGTTGCCTCATCg                                           <  1:516846/250‑1 (MQ=255)                                                                                                                                                                                           cACNGGCACAGTGTGC‑GTAACGGCTGCGTCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTTGATCAGGTAATCCACGACACCGTAATGCAGCGAATCTTTAATGGTTGTCGCATCGGc                                         <  2:596147/251‑1 (MQ=255)                                                                                                                                                                                             cTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCCAACGACGCGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTTg                                                                                                 >  2:88043/1‑193 (MQ=255)                                                                                                                                                                                             cTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTGCTTTGGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTTg                                                                                                 <  1:88043/193‑1 (MQ=255)                                                                                                                                                                                              tGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTTGATCAGGTAATCCACGACACCGTAATGCAGCGAATc                                                              >  2:129561/1‑227 (MQ=255)                                                                                                                                                                                              tGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTTGATCAGGTAATCCACGACACCGTAATGCAGCGAATc                                                              <  1:129561/227‑1 (MQ=255)                                                                                                                                                                                                  acaGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCgg                                                                                                                   >  1:300912/1‑170 (MQ=255)                                                                                                                                                                                                  acaGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCgg                                                                                                                   <  2:300912/170‑1 (MQ=255)                                                                                                                                                                                                      tgtgC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTTGATCAGGTAATCCACGACACCGTAATGCAGCGAATCTTTAATGGTTGCCGCATCGGCTGCggagg                                 >  1:319161/1‑248 (MQ=255)                                                                                                                                                                                                              aaCGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTTGATCAGGTAATCCACGACACCGTAATGCAGCGAATCTTTAATGGTTGCCGCATCGGCTGCGGAGGAGATGACaa                        <  2:248823/250‑1 (MQ=255)                                                                                                                                                                                                                cGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTTGATCAGGTAATCCACGACACCGTAATGCAGCGAATCTTTAATGGTTGCCGCATCGGCTGCGGAGGAGATGACAat                       <  1:561093/249‑1 (MQ=255)                                                                                                                                                                                                                cGTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTTGATCAGGTAATCCACGACACCGTAATGCAGCGAATCTTTAATGGTTGCCGCATCGGCTGCGGAGGAGATGACAat                       >  2:561093/1‑249 (MQ=255)                                                                                                                                                                                                                 gTCTGCGGCGTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTTGATCAGGTAATCCACGACACCGTAATGCAGCGAATCTTTAATGGTTGCCGCATCGGCTGCGGAGGAGATGACAatcac                    <  1:92651/251‑1 (MQ=255)                                                                                                                                                                                                                     gcggcgTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTTGATCAGGTAATCCACGACACCGTAATGCAGCGAATCTTTAATGGTTGCCGCATCGGCTGCGGAGGAGATGACAatcacatc                 <  1:193606/250‑1 (MQ=255)                                                                                                                                                                                                                        gcgTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTTGATCAGGTAATCCACGACACCGTAATGCAGCGAATCTTTAATGGTTGCCGCATCGGCTGCGGAGGAGATGACAATCACATCACtt             <  2:178354/251‑1 (MQ=255)                                                                                                                                                                                                                          gTTAAGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTTGATCAGGTAATCCACGACACCGTAATGCAGCGAATCTTTAATGGTTGCCGCATCGGCTGCGGAGGAGATGACAATCACATCACtttt           <  2:172390/251‑1 (MQ=255)                                                                                                                                                                                                                             aaGCCTTTTGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTTGATCAGGTAATCCACGACACCGTAATGCAGCGAATCTTTAATGGTTGCCGCATCGGCTGCGGAGGAGATGACAATCACATCACTTTTGc         <  2:491/250‑1 (MQ=255)                                                                                                                                                                                                                                    ttGGCAAACGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTTGATCAGGTAATCCACGACACCGTAATGCAGCGAATCTTTAATGGTTGCCGCATCGGCTGCGGAGGAGATGACAATCACATCACTTTTGCAAcgcgc  <  2:491037/250‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                   GAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTATTTACGGCAGGAAACACGCGAAATGTTAACCTCGTTGGCTAGCTCGTCGGTTGAAAATTCATAGTCCTGATGCGCGTCAATCCACTGGCACAGTGTGC‑GTAACGTCTGCGGCGTTAAGCCTTTTGGCAAGCGACGAGGATCCTGTTCGTTGGAGCTGCTGCCGTGGATTAGCTGATCAAGCTCGGCCTGGTCATAATACTGATGTTTTTCCAGCGCCATTTTCTTTTGCCGCCAGCCGGTGAGCGCCTCTTCAAAGCGGGAAGCCTGGAAGGGTTTGATCAGGTAATCCACGACACCGTAATGCAGCGAATCTTTAATGGTTGCCGCATCGGCTGCGGAGGAGATGACAATCACATCACTTTTGCAACGCGC  >  NC_000913/4349340‑4349815

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