Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,150,316	G→A	E219K (GAG→AAG)	yghA →	NADP(+)‑dependent aldehyde reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,150,316	0	G	A	100.0%	100.6 / NA	36	E219K (GAG→AAG)	yghA	NADP(+)‑dependent aldehyde reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/12);  total (24/12)

CCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGGAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATG  >  NC_000913/3150088‑3150553                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              ccGGATATTGCAGCCATCCCCAGCGAACAGATTCAAAAGACCTTGGCCATTAACTTGTTCGCGCTGTGCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCATTCCAGGCATACCAGCCAAGTCCGTATTTACTGGACTATGCGGCGACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAata                                                                                                                                                                                                                          <  2:143088/250‑1 (MQ=255) ccGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAatat                                                                                                                                                                                                                         >  1:227244/1‑251 (MQ=255)     atatTACAGACCTCACCAGCGAAAAGTTTCACAAGACCTTTGCCATTCACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCACTGCTACCGACAGGTGCAAGTATCCTCAACACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGACGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTc                                                                                                                                                                                                                     >  2:230366/1‑251 (MQ=255)           cAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCa                                                                                                                                                                                                               <  2:109714/251‑1 (MQ=255)                   aCCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGAt                                                                                                                                                                                                       >  1:319442/1‑251 (MQ=255)                           aCAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAg                                                                                                                                                                                               <  1:30241/251‑1 (MQ=255)                                  taaaaGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGAGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGc                                                                                                                                                                                         >  2:16322/2‑250 (MQ=255)                                           tttGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCcgg                                                                                                                                                                               >  1:25447/1‑251 (MQ=255)                                                cATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCTGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGAc                                                                                                                                                                          <  1:199143/251‑1 (MQ=255)                                                  ttAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACa                                                                                                                                                                         >  1:130188/1‑250 (MQ=255)                                                         tttCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAgg                                                                                                                                                                 >  1:355293/1‑251 (MQ=255)                                                                   ttCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGCCAGCACTGCAAATTTCCGGCGGACAAACGAAGGATAAGATccc                                                                                                                                                       >  1:207717/1‑251 (MQ=255)                                                                              ccAGGAAGCGTTCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCACGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCATTTTGGTc                                                                                                                                            <  2:25447/251‑1 (MQ=255)                                                                                   aaGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGACGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGACTATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACACACGCAGGATAAGATCCCGCAGTTTGGTCAGCa                                                                                                                                        >  2:18127/1‑250 (MQ=255)                                                                                    aGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGGCAATCCCGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCTGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGAGCTGGACATCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCaaa                                                                                                                                      <  2:326477/251‑1 (MQ=255)                                                                                        aacccccTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGCAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAACCAGGTAGCGAAGAAAGGAATTCGGGTGAATATTGTCGCGCCAGGCCCGCACTGGACAGCCATGCAAATTTCCGGCGGAAAAACGCAGGATAAAATCCCGCAGTTTGGTCTGCAAACGcc                                                                                                                                  >  2:223161/3‑251 (MQ=255)                                                                                           ccccTGCTACCGAAAGGAGAACGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGCTTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGAt                                                                                                                               >  2:29595/1‑251 (MQ=255)                                                                                            cccTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATg                                                                                                                              >  1:290311/1‑251 (MQ=255)                                                                                                   aCCGAAAGGTGCCAGTATCATCAACACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGCACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGAAAACGCCGATGAAACGTg                                                                                                                       >  2:276282/1‑251 (MQ=255)                                                                                                     cGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCg                                                                                                                     >  1:305296/1‑251 (MQ=255)                                                                                                                      accaccaatgcGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGGCTATGCGGCTACGAAGGCGGCGATTCTGAACTGCAGCCGTGGCTTGGCAAATCAGGTCGCGAAGAAAGGTATTCGGGTGAATTTTGTCGCGCCAGGCCCGATCGGGACAGCACTGCAAATTTCCGGCGGGCAAACGCGGGATAAGTTCCCGCAGGTTGGTCAGCAAACGCCGATGAAACGTTCGGGGCAGCCGGCGGAACt                                                                                                    <  2:164841/241‑1 (MQ=255)                                                                                                                                 gTCAATCCAGGCATACAAGCCAAGTACGCATTTACTGGACAATGAGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAACCGCAGGATAAGCTCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGtat                                                                                         >  2:232818/1‑251 (MQ=255)                                                                                                                                   cAATCCAGTCCTGCCATCCAAGTCCGCATGTACTGGACTATGCGGATACGAAGGCGGCGATTCTGAACGACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGGGAATATTGTCGCGCCAGGCCCGATCGTGACAGCACTGCAAAGTTCCGGCTTACAAACGCAGGATAATATCCCGCAGTGTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCTGACCTGGTCCCTGtata                                                                                        <  2:27041/250‑1 (MQ=255)                                                                                                                                    aaTCCAGGAATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAGTATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAAAGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATg                                                                                      >  1:94381/1‑251 (MQ=255)                                                                                                                                                     cAAGTCCGCATTTACTGGACTCTGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCCAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAACCGCCGATGAAACGTGAGGGGCCACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAggc                                                                       >  2:59599/1‑247 (MQ=255)                                                                                                                                                        gTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAgg                                                                    >  1:144358/1‑249 (MQ=255)                                                                                                                                                                       aCTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCTATGAAACGTGCGGGGCAACCGGCGGAACTTGCCCCTGTATATGTGTATCTGGCAAGCCAGGAGTCGAGCTACGTCAcc                                                   <  1:383823/251‑1 (MQ=255)                                                                                                                                                                       aCTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCAcc                                                   >  1:134400/1‑251 (MQ=255)                                                                                                                                                                                     aaGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACgg                                     >  1:237820/1‑251 (MQ=255)                                                                                                                                                                                          ggcgATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACTTCACCGCAGAAGTGCACGGCgtgg                                >  2:313236/1‑250 (MQ=255)                                                                                                                                                                                            cgATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGGGTgc                              >  2:289980/1‑251 (MQ=255)                                                                                                                                                                                                        ccAGCCGTGTCTTGGCAAAGCCGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCCGTCCCGTGCTGGACAGCACGGCAAATTTCCGTCGGACAAACGCAGTTTAATATCCCGCAGTTTGTTCAGCTAATGCCTAAGATACGTGCGGGGCAGCCGTCGGAACTGGCCCCTGTATTTGTTTATCTGTGAAGTCAGTAGTCAAGCTACGTCTCCGCTGAAGTGCGCGTCGTGTGCTTCGTCGAGCAt                  <  2:366029/250‑1 (MQ=255)                                                                                                                                                                                                          agcatgggCTTGGCACAACAGGTCGCGAAGAAAGGGATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGAGCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCAttt                <  2:219886/245‑1 (MQ=255)                                                                                                                                                                                                               tGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGGGCATTTAGGtt           >  1:114360/1‑251 (MQ=255)                                                                                                                                                                                                                     ggCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTaaaaa      <  2:114072/250‑1 (MQ=255)                                                                                                                                                                                                                        aaaaCAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATg  >  2:327039/1‑251 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              CCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGGAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATG  >  NC_000913/3150088‑3150553

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