Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,156,839	G→A	intergenic (+50/‑11)	argB → / → argH	acetylglutamate kinase/argininosuccinate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,156,839	0	G	A	100.0%	119.6 / NA	36	intergenic (+50/‑11)	argB/argH	acetylglutamate kinase/argininosuccinate lyase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/16);  total (20/16)

GCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGGCGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGAAGATGTTCGCGCCAGGCCACAACA  >  NC_000913/4156613‑4157073                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           gCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGGCGCTGGATGCGGCCCGCGCGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTgggg                                                                                                                                                                                                                    <  1:166532/251‑1 (MQ=255)          cTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGGCGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGttt                                                                                                                                                                                                           >  1:74560/1‑251 (MQ=255)             gttGAGCAGGGCATTTTTACTTACTGCATGTTAGTGAAAGTGAACGCGGCGCTGGATGCGGCCCGCAGGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGGCATGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTAc                                                                                                                                                                                                        <  2:74546/250‑1 (MQ=255)               tGAGCATGGCATTATTGCTGACGGCATGATAGTGAAATTGAACGCGGCGATGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCGGAGGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGTGTTTTAcc                                                                                                                                                                                                       <  1:344032/250‑1 (MQ=255)                 cgCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGGCGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTAAAAATAAGAAAACAGAGTTATAGCAATTTGGGGCGGGCGTTTTAccc                                                                                                                                                                                                      >  2:235671/2‑249 (MQ=255)                 aGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGGCGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTAccc                                                                                                                                                                                                      >  1:277770/1‑249 (MQ=255)                               aCTGACGGCATGATAGTGAAAGTGAACGCGGCGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATTGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACg                                                                                                                                                                                      >  1:348240/1‑251 (MQ=255)                               aCTGACGGCATGATAGTGAAAGTGAACGCGGCGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAAc                                                                                                                                                                                       <  1:234579/250‑1 (MQ=255)                                 tGACGGCATGATAGTGAAAGTGACCGCGGCGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGGAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAGGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGGGTTTTACCCCGGCAGCAGATCAACggg                                                                                                                                                                                    >  2:81508/1‑250 (MQ=255)                                  gACGGCATGATAGTGAAAGTGAACGCGGCGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACgg                                                                                                                                                                                     <  1:121141/249‑1 (MQ=255)                                              gtgaaagtaaaCGCGGCGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTAATGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGAACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTACAACAATTCa                                                                                                                                                                        >  1:227904/1‑250 (MQ=255)                                                         cgcgGCGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTgc                                                                                                                                                            >  1:72301/1‑251 (MQ=255)                                                          gcggcgCAGGATGTGGCCCGCCCGGTGGGCCGCCCGGTAGATATCGCCTCCTGGCGTCGTGCGGAGCAGCTTCCGGCACTGTTTAACTGGATGCCGATGGGGACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTgcg                                                                                                                                                           <  2:186622/251‑1 (MQ=255)                                                          ccggggcggcatgtggCCCGCACGCTGTGACGTCAGGTAGATATCGCCTCCTGGCGTCAGGCGGAGCAGCTTCCGTCACTGTTTAAGTGTATGCTGAGTGGTACGCGGATGTTAGCGTAAGTATGGTTGGCCGGAGGCGCAGCTGTCAGGCATTGAAGTTCAAAATAAGAAAACAGAGTGATTGCACTTTAGGGCGGGCGTTTTACCCAGGCAGCAGATCAACTGTTGAAACAATTCCACGACACACTgcg                                                                                                                                                           <  2:68264/237‑1 (MQ=255)                                                               gcTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAAAGACTCACTGCGCTTTg                                                                                                                                                      >  1:161679/1‑251 (MQ=255)                                                                                     ggCCGTCCGGTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAgg                                                                                                                                 >  1:269315/1‑250 (MQ=255)                                                                                             ggTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGGGTTTTACCCAGGCAGCAGATCAACGGTTCAAACCATTCCACGACTCACTGCGCTTTGATTACCg                                                                                                                                               >  2:250909/1‑228 (MQ=255)                                                                                             ggTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCg                                                                                                                                               <  1:250909/228‑1 (MQ=255)                                                                                                              ggCGTCATGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTcc                                                                                                       <  2:269315/251‑1 (MQ=255)                                                                                                                 gTCATGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCttt                                                                                                                                                       <  2:55052/200‑1 (MQ=255)                                                                                                                 gTCATGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCttt                                                                                                                                                       >  1:55052/1‑200 (MQ=255)                                                                                                                         gggcagcTTCCGGCACTGGTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGTGCAGCTTTCCTGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCGGTTCCAAAGCCCTGGt                                                                                            <  2:348240/249‑1 (MQ=255)                                                                                                                               cTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATAAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTTTTGGCTCTGTGGCCCGGTCCAAAGCCCTGGTCCCggg                                                                                      >  2:145778/1‑250 (MQ=255)                                                                                                                                       aCTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCgtgt                                                                              >  1:179415/1‑251 (MQ=255)                                                                                                                                            ttAACGGGATGCAGATGGGTACGCGGATTTTAGCTTAATTTTTGTTGGCCGGAGTCGTAGTTTTCCGGCATTGAATTACAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGTGGTTTACCCAGGCAGCAGATCAATGGATCAAACAAGTCAACGACTCACTGCGCTTTGATTACCTTCTGGCGGAGAAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTAACTGTATGCGTGTTAAc                                                                          <  2:99809/250‑1 (MQ=255)                                                                                                                                                             ggTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCG‑CATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCAATGCGATTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCCGGTCCAAAGCCCTGGGCCCGGTTGGCGTGTTAACCGGAGAAGGGCCGgcggc                                                        >  1:4792/1‑248 (MQ=255)                                                                                                                                                                  gcgGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTgg                                                   <  1:145778/251‑1 (MQ=255)                                                                                                                                                                   cgGATTTTAGCTTATGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGa                                                  <  1:318754/251‑1 (MQ=255)                                                                                                                                                                                 aGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCCGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACg                                    >  1:83308/1‑251 (MQ=255)                                                                                                                                                                                     ttgttgGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACgtgt                                 <  1:84668/250‑1 (MQ=255)                                                                                                                                                                                       gttgGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAACCAGAGTTATGGCACTTTGGGGCGGGCGGTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGATTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCACGGGAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGAGCTGAACGTGTTgg                              >  1:302362/1‑250 (MQ=255)                                                                                                                                                                                             cGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTAACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGa                          >  1:241985/1‑249 (MQ=255)                                                                                                                                                                                              ggAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGa                          >  1:11663/1‑248 (MQ=255)                                                                                                                                                                                              ggAGGCGCAGCTTTACGGCATTGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTAACGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCACGGTAAGCTTGTTAACCGCAGCAGAGCAGGCGCACCTGGAAGAGGCGCTGAACGTGTTGCTGGa                          >  1:262697/1‑248 (MQ=255)                                                                                                                                                                                                                   tGAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGAAGATGTTCGCGCCAGGCCacaac   <  1:37371/250‑1 (MQ=255)                                                                                                                                                                                                                    gAATTTCAAAATAAGAAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGAAGATGTTCGCGCCAGGCCacaaca  >  1:339474/1‑250 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           GCAGAACAACTGATTGAGCAGGGCATTATTACTGACGGCATGATAGTGAAAGTGAACGCGGCGCTGGATGCGGCCCGCACGCTGGGCCGTCCGGTAGATATCGCCTCCTGGCGTCATGCGGAGCAGCTTCCGGCACTGTTTAACGGTATGCCGATGGGTACGCGGATTTTAGCTTAAGTTTTGTTGGCCGGAGGCGCAGCTTTCCGGCATTGAATTTCAAAATAAGGAAACAGAGTTATGGCACTTTGGGGCGGGCGTTTTACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGAAGATGTTCGCGCCAGGCCACAACA  >  NC_000913/4156613‑4157073

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