Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,518,339	G→A	L50F (CTC→TTC)	aroB ←	3‑dehydroquinate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,518,339	0	G	A	100.0%	127.4 / NA	37	L50F (CTC→TTC)	aroB	3‑dehydroquinate synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/23);  total (14/23)

CGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAGGCTTTT‑ATACTGCTCGCCGTCAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACCTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACT  >  NC_000913/3518100‑3518585                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       cgtcgGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAGGCTTTT‑ATACTGCTCGCCGTCAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAgagg                                                                                                                                                                                                                                               >  2:340507/1‑249 (MQ=255) cgtcgGGACTTGAAGGAACCGGGCCCCGCGCTGATAACCCGCCGCCGCAAAGCCGGTCAAACCGCCCACTACGCCCCCGCCAAGCTCCCCCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACACCCAGGCTTTT‑ATCCGGCTCGCCGTCAGGGAGGATAACTCTATTGTCGTGAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAgaggag                                                                                                                                                                                                                                             <  1:49950/251‑1 (MQ=255)     gggACTTAAATGAAACGGCCCCCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGCCCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAGGCTTTT‑ATACTGCTCGCCGTCAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAg                                                                                                                                                                                                                                         <  1:160598/251‑1 (MQ=255)      ggACTTGAATGAACCGGACCCCGCGTGGATAACCCGCCGCCGCGAAGCCGGTCGGCTCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAGGCTTTT‑ATACTGCTCGCCGTCAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAgg                                                                                                                                                                                                                                        <  2:319955/251‑1 (MQ=255)            gaagaaaaCGGCCACCGCGCTGATACCTCGCCGCCGCGAAGCCGGTCAGATCGCCCCCAACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACACCGCCGTAAAGACGGTATCGATTACAGCCAGGCTTTT‑ATACTGCTCGCCGTCAGGGAGGATACCGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTc                                                                                                                                                                                                                                   <  1:280630/245‑1 (MQ=255)                   acggcccccGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGCGCGCCCACTACGCCCCCGCCAGGCGCCACCAGCGTAGTAGCGCGACCATGCGTTTTTTGTAACAACGCCGTAAAGCCGGTATCGATTACAGCCAGGCTTTT‑ATACTGCTCGCCTTCAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGGTGGTGa                                                                                                                                                                                                                           <  2:298347/244‑1 (MQ=255)                   aCGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAGGCTTTT‑ATACTGCTCGCCGTCAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTg                                                                                                                                                                                                                            >  1:342926/1‑250 (MQ=255)                          ccGCGTTGATCACTCGCCCCCGCGAAACGGGCCAGTTCGCCCACTACGCCGCCGCCAAGCCCCACCAGCGTGGTATCGCTACCATGCGTTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAGGCTTTT‑ATACTGCTCGCCGTCAGGTAGGATACCGCTCTCGACTTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACAt                                                                                                                                                                                                                    <  2:443846/251‑1 (MQ=255)                              gcTGATAACTCGCCGCCGCGAAGCCGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAGGCTTTT‑ATACTGCTCGCCGTCAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGAcc                                                                                                                                                                                                                >  1:247519/1‑251 (MQ=255)                               cTGATAACTCGCCGCCGCGAAGCCGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAGGCTTTT‑ATACTGCTCGCCGTCAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCt                                                                                                                                                                                                               >  1:259342/1‑251 (MQ=255)                                   tAACTCGCCGCCGCGAAGCCGGACAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAGGCTTTT‑ATACTGCTCGCCGTCAGGGAGGATAACGATATCGACGTTAACACCCGCCTGTTCAAGTACGCCGAGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACAAACATGACCTGCTc                                                                                                                                                                                                           >  2:435115/1‑251 (MQ=255)                                    aacccaccgcgacgtaGCCGGTAAGATCGCCCACTACTCCGCCGCCAAGCGCCTCCGGGGTGGTATCGCGCGCGTGCGGTTTGTGAAACAAGGCCGTAAAGCGGGTTTCGGGTACAGCCAGGCTTTT‑ATACTGCTCGCCGTGAGGGGGGATAACGGTAGCAACGTTAACGCCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTc                                                                                                                                                                                                           <  2:400037/235‑1 (MQ=255)                                    aaCTCGCCGCCGCGAAGCCGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGCCCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAGGCTTTT‑ATACTGCTCGCCGTCAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCg                                                                                                                                                                                                          <  1:203062/251‑1 (MQ=255)                                                                     cTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAGGCTTTT‑ATACTGCTCGCCGTCAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTc                                                                                                                                                                          >  1:32084/1‑250 (MQ=255)                                                                       aCGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAGGCTTTT‑ATACTGCTCGCCGTCAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTTCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCAtt                                                                                                                                                                       <  1:382899/251‑1 (MQ=255)                                                                        cgccgccgcCAAGCGCCCCCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAGGCTTTT‑ATACTGCTCGCCGTCAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTa                                                                                                                                                                      <  1:436633/251‑1 (MQ=255)                                                                           cgcnccaaaGCCCCACCAGCGTAATATCGCGACGTGGCTGTTTTTTTAACAACGCCGTAAAGACGGTATAGAGTCCAGCCAGTCTTTT‑ATACTGCTCGCCGCCAGGGAGGATAACGCTATCGACGTTAACCCCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCAGTaaac                                                                                                                                                                   <  2:405073/244‑1 (MQ=255)                                                                                    gcgcCACTAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCAAAAAGACGGCACCGCGTACAGCCAGGCTTTT‑ATACTGCTCGCCGTCAGGGAGGATAACGCTATCGACGTTCACCCACGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGATTTCGTTGGTGACAAAAATGACCTGCTCGCACGCTTTCAGCGGTAAGAATGAAGCTGGTTCATGAAACAAAcc                                                                                                                                                              >  1:289620/1‑247 (MQ=255)                                                                                         accaGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAGGCTTTT‑ATACTGCTCGCCGTCAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGCACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGAt                                                                                                                                                     >  1:36007/1‑251 (MQ=255)                                                                                                       gcgACCGTGCGGTTTTTGTAACAACGCCGTAAAGACGTTATCGAGTACAGCCAGGCTTTT‑ATACGGCTCGCCGTCAGGGAGGAAAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGTTAATTGGGTAAc                                                                                                                                       <  2:193021/251‑1 (MQ=255)                                                                                                             aTGCGTTTTTTGTAACACGGCCGTAAAGACGGTGTCGAGTACAGCCCGGCTTTT‑ATACTGCCCGCCGTCAGGGAGGAAAACCCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGtt                                                                                                                                 <  2:83171/251‑1 (MQ=255)                                                                                                                    ttttGTAACAACACCGTCAAGACCGTATCGAGTACAGCCAGGCTTTT‑CTACTGCTCGCCGTCAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTAGCCCGATCTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGTGATGGTAATTGGGTAACTACGTTCCCCg                                                                                                                            >  2:256163/1‑249 (MQ=255)                                                                                                                             accGCCGTAAAGACGGTATCGAGTACAGCCAGGCTTTT‑ATCCTGCTCGCCGTCAGGGAGGATACCGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACCAAc                                                                                                                                                                                                                      <  2:126435/148‑1 (MQ=255)                                                                                                                             aacGCCGTAAAGACGGTATCGAGTACAGCCAGGCTTTT‑ATACTGCTCGCCGTCAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACCAAc                                                                                                                                                                                                                      >  1:126435/1‑150 (MQ=255)                                                                                                                                         gcGGTACCGAGTACAGCCAGGCTTTT‑ATACTGCTCGCCGTCGGGGAGGATAACGCTATCGCCGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCAt                                                                                                     <  2:259342/250‑1 (MQ=255)                                                                                                                                         aCGGTATCGAGTACAGCCAGGCTTTT‑ATACTGCTCGCCGTCAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCAt                                                                                                     >  1:171644/1‑251 (MQ=255)                                                                                                                                         aCGGTATCGAGTACAGCCAGGCTTTT‑ATACTGCTCGCCGTAGGGGAGGATAACTCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCAt                                                                                                     <  2:95378/251‑1 (MQ=255)                                                                                                                                                     aCAGCCAGGCTTTT‑ATACTGCTCGCCGTCAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATcc                                                                                         >  1:55179/1‑251 (MQ=255)                                                                                                                                                           cagCTTTTCATACTGCGCCCCGTCTGGGTGGAAAACGCTATCGACGTAAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTACCTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACCtt                                                                                    <  2:359405/249‑1 (MQ=255)                                                                                                                                                                                  aaGGGAGGATAACGCTAGCGACGTTAACACACGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGCGAGTAACGACAATCCTCTCCATAACGCGACAGCCACCTTAATTACTGTACCCGCAGACGAg                                                              <  1:220216/249‑1 (MQ=255)                                                                                                                                                                                        ggacaccGCTATGGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACTACAATCCTCTCCATAACGCGACATCCACCTTAATTACTGTGCCCGCAGACGAGTGtatat                                                       <  2:223734/245‑1 (MQ=255)                                                                                                                                                                                                         ttAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACCTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTg                                      >  1:298391/1‑251 (MQ=255)                                                                                                                                                                                                         taaacccccGCCTGTTCAAGTACGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTGCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACCTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGtt                                       <  2:271531/244‑1 (MQ=255)                                                                                                                                                                                                                              aCGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACCTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGaataat                  <  1:350875/250‑1 (MQ=255)                                                                                                                                                                                                                              aCGCCGCGGACCTTATCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACCTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATc                 <  2:342926/251‑1 (MQ=255)                                                                                                                                                                                                                                            aTCGAAATACAGATGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACCTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGTATAATCTGGTTTGCAACCa    <  2:36007/250‑1 (MQ=255)                                                                                                                                                                                                                                             tCGAAATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACCTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACt  >  2:352674/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       CGTCGGGACTTGAATGAAACGGACACCGCGCTGATAACTCGCCGCCGCGAAGCCGGTCAGATCGCCCACTACGCCGCCGCCAAGCGCCACCAGCGTAGTATCGCGACCATGCGGTTTTTGTAACAACGCCGTAAAGACGGTATCGAGTACAGCCAGGCTTTT‑ATACTGCTCGCCGTCAGGGAGGATAACGCTATCGACGTTAACACCCGCCTGTTCAAGTACGCCGCGGACCTTATCGAGATACAGAGGAGCCAGGGTTTCGTTGGTGACCAACATGACCTGCTCGCCCGATTTCAGCGGTAAGAATGAAGCTGGTTCATTAAACAAACCAGATGCGATGGTAATTGGGTAACTACGTTCCCCGAGAGTAACGACAATCCTCTCCATAACGCGACATCCACCTTAATTACTGTACCCGCAGACGAGTGTATATAAAGCCAGAATTAGTTGCTTTCCAGCATGTGAATAATCTGGTTTGCAACCACT  >  NC_000913/3518100‑3518585

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