Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,207,461	G→A	A24A (GCC→GCT)	purH ←	bifunctional AICAR transformylase/IMP cyclohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,207,461	0	G	A	100.0%	142.4 / NA	42	A24A (GCC→GCT)	purH	bifunctional AICAR transformylase/IMP cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/25);  total (17/25)

GTTAACAACCACCATATCGATAGGCTGGATCTGATGTTCTTCCATAATGGCATCGTCCTGGCCGCGACGGCCCAGAATGCCACCATGTACTTTCGGATGCAGGGTCTTCACGCGTCCATCCATCATCTCCGGGAAACCGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCAC  >  NC_000913/4207221‑4207691                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                       gTTAACAACCACCATATCGATAGGCTGGATCTGATGTTCTTCCATAATGGCATCGTCCTGGCCGCGCCGGCCCAGAATGCCACCATGTACTTTCTGATGCAGGGTCTTCACGCGTCCATCCATCATCTCCGGGAAACCGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTcga                                                                                                                                                                                                                              <  1:31793/251‑1 (MQ=255)      caacCACCATATCGATAGGCTGGATCTGATGTTCTTCCATAATGGCATCGTCCTGGCCGCGACGGCCCAGAATGCCACCATGTACTTTCGGATGCAGGGTCTTCACGCGTCCATCCATCATCTCCGGGAAACCGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGAt                                                                                                                                                                                                                          >  1:207557/1‑250 (MQ=255)            ccaTATCGATAGGCTGGCTCTGATGTTCTTCCATAATGGCATCGTCCTGGACGCGACGGCCCAGAATGCCACCATGTACTTTCGGATGCAGGGTCTTCACGCGTCCATCCATCATCTCCGGGAAACCGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGCAAGTGCCTGAGCGAATTCGACGATACCGGCt                                                                                                                                                                                                                   >  2:77482/1‑251 (MQ=255)                atCGATAGGCTGGATCTGATGTTCTTCCATAATGGCATCGTCCTGGCCGCGACGGCCCAGAATGCCACCATGTACTTTCGGATGCAGGGTCTTCACGCGTCCATCCATCATCTCCGGGAAACCGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTg                                                                                                                                                                                                                <  1:47098/250‑1 (MQ=255)                    aTAGGCTGGATCTGATGTTCTTCCATAATGGCATCGTCCTGGCCGCGACGGGACAGAATGCCACCATGTACTTTCGGATGCAGGGTCTTCACGCGTCCATCCATCATCTACGG‑AAACCGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGAAGTACCCCCTGTAGACCGCAGCTCCACACCTGGTGCGGAAAGCGACTGAGCGAATTCGACGATACGGGCTTTGTCAGa                                                                                                                                                                                                           >  1:222642/1‑250 (MQ=255)                       ggCTGGATCTGATGTTCTTCCATAATGGCATCGTCCTGGCCGCGACGGCCCAGAATGCCACCATGTACTTTCGGATGCAGGGTCTTCACGCGTCCATCCATCATCTCCGGGAAACCGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACAGGCTTTGTCAGAAac                                                                                                                                                                                                        >  1:194056/1‑251 (MQ=255)                            gATCTGATGTTCTTCCATAATGGCATCGTCCTGGCCGCGACGGCCCAGAATGCCACCATGTACTTTCGGATGCAGGGTCTTCACGCGTCCATCCATCATCTCCGGGAAACCGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTg                                                                                                                                                                                                    >  2:356824/1‑250 (MQ=255)                                 gATGTTCTTCCATAATGGCATCGTCCTGGCCGCGACGGCCCAGAATGCCACCATGTACTTTCGGATGCAGGGTCTTCACGCGTCCATCCATCATCTCCGGGAAACCGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGa                                                                                                                                                                                                   <  1:61511/246‑1 (MQ=255)                                 gATGTTCTTCCATAATGGCATCGTCCTGGCCGCGACGGCCCAGAATGCCACCATGTACTTTCGGATGCAGGGTCTTCACGCGTCCATCCATCATCTCCGGGAAACCGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCaga                                                                                                                                                                                              <  2:272896/251‑1 (MQ=255)                                            aTAATGGCATCGTCCTGGCCGCGACGGCCCAGAATGCCACCATGTACTTTCGGATGCAGGGTCTTCACGCGTCCATCCATCATCTCCGGGAAACCGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGATCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGAc                                                                                                                                                                                   >  1:50735/1‑251 (MQ=255)                                             tAATGGCATCGTCCTGGCCGCGACGGCACAGAATGCCACCATGTACTTTCGGATGCAGGGTCTTCACGCGTCCATCCATCATCTCCGGGACACCGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAGTTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGNACt                                                                                                                                                                                  >  2:378136/1‑251 (MQ=255)                                                    aTCGTCCTGGCCGCGACGGCCCAGAATGCCACCATGTACTTTCGGATGCAGGGTCTTCGCGCGTCCATCCATCATCTCCGGGAAACCGGTGTAATCGGAAACTTCGGTTACCGTCAGACCCTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGgacgac                                                                                                                                                                           <  1:404474/251‑1 (MQ=255)                                                            ggCCGCGACGGCCCAGAATGCCACCATGTACTTTCGGATGCAGGGTCTTCACGCGTCCATCCATCATCTCCGGGAAACCGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACgttgtt                                                                                                                                                                     <  1:148684/249‑1 (MQ=255)                                                                 cgACGGCCCAGAATGCCACCATGTACTTTCGGATGCAGGGTCTTCACGCGTCCATCCATCATCTCCGGGAAACCGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGt                                                                                                                                                              <  1:240233/251‑1 (MQ=255)                                                                    cGGCCCAGAATGCCACCATGTACTTTCGGATGCAGGGTCTTCCCGCGTCCATCCATCATCTCCGGGAAACCGGTGTAAGCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTGTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTaaa                                                                                                                                                           <  1:43725/251‑1 (MQ=255)                                                                        ccAGAATGCCACCATGTACTTTCGGATGCAGGGTCTTCACGCGTCCATCCATCATCTCCGGGAAACCGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATccc                                                                                                                                                       >  1:357126/1‑251 (MQ=255)                                                                           gAATGCCACCATGTACTTTCGGAGGCAGGGTCTTCACGCGTCCATCCATCATCTCCGGGAAACCGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTg                                                                                                                                                    <  2:207557/251‑1 (MQ=255)                                                                                   ccaTGTACTTTCGGATGCAGGGTCTTCACGCGTCCATCCATCATCTCCGGGAAACCGGTGTAATCGGAAACTTCTGTTACCGTCAGACCTTTTTCTGCTAACAGGCGGGCAGTGCCCCCTGTAGACAGCAGCTCCTCACCGCGTGCGGAAAGTACCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACGCTGAGCAGAGCGCGGCGGACGGGACTACGTTGTTGCATGGTAAATCCCCTGGATTTGAc                                                                                                                                            <  2:270059/251‑1 (MQ=255)                                                                                        tACTTTCGGATGCAGGGTCTTCACGCGTCCATCCATCATCTCCGGGAAACCGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTAtt                                                                                                                                        >  1:209098/1‑250 (MQ=255)                                                                                            ttCGGATGCAGGGTCTTCACGCGTCCATCCATCATCTCCGGGAAACCGGTGTAATCGGAAACGTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACag                                                                                                                                    <  2:283944/250‑1 (MQ=255)                                                                                            tnCGGATGCAGGGTCTTCACGCGTCCATCCATCATCTCCGGGAAACCGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACa                                                                                                                                     <  2:371066/249‑1 (MQ=255)                                                                                                 aTGCAGGGTCTTCACGCGTCCATCCATCATCTCCGGGAAACCGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCg                                                                                                                              <  1:406971/251‑1 (MQ=255)                                                                                                        gTCTTCACGCGTCCATCCATCATCTCCGGGAAACCGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTg                                                                                                                       <  1:239603/251‑1 (MQ=255)                                                                                                           ttCACGCGTCCATCCATCCTCGACGGGAAACCGCAGTACTCGAAACCTTCGGACACCGACAGACCTCTTTCTGCTAAACGACCGGCAGTGCCCCCTCAAAACCGCCGCTCAACACCGCCTGCGCAAAGTGCCTGAGCGAGTTCGACGATACCGGCTTTCTCAGAAACACTCAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTACATCCCCTGCATTTGACTATTACACAGAGCGTTacctncat                                                                                                                    >  2:377041/1‑244 (MQ=255)                                                                                                                gcgTCCATCCATCATCTCCGGGAAACCGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTc                                                                                                               <  1:198178/251‑1 (MQ=255)                                                                                                                     catccatcATCTCCGGGAAACCGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAGTTTTTCGCGa                                                                                                           >  1:128738/1‑250 (MQ=255)                                                                                                                                 ccGGGAAACCGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCt                                                                                                <  2:231750/249‑1 (MQ=255)                                                                                                                                     gAAACCGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTTCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGcc                                                                                          <  1:325123/251‑1 (MQ=255)                                                                                                                                         ccGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGCCTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCt                                                                                       >  2:347056/1‑250 (MQ=255)                                                                                                                                         ccGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTa                                                                                      <  1:150249/251‑1 (MQ=255)                                                                                                                                          cGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTa                                                                                      >  1:129160/1‑250 (MQ=255)                                                                                                                                                aaTCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGc                                                                               <  2:357126/251‑1 (MQ=255)                                                                                                                                                    ggAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCGTcc                                                                           <  1:390374/251‑1 (MQ=255)                                                                                                                                                           tCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTAttt                                                                     <  2:312573/250‑1 (MQ=255)                                                                                                                                                                          ccTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTaa                                                     <  2:128738/251‑1 (MQ=255)                                                                                                                                                                           cTTTTTCTGATAACAGACGGGCCGTGCCCCCTGTAGACAGCAGCTCAACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAAc                                                    >  2:160190/1‑251 (MQ=255)                                                                                                                                                                                 cTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACCCCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATCCCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGANAAc                                              >  2:376306/1‑251 (MQ=255)                                                                                                                                                                                            cGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAAc                                   <  1:109322/251‑1 (MQ=255)                                                                                                                                                                                                      ccccTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAAt                         <  2:209098/251‑1 (MQ=255)                                                                                                                                                                                                              gACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCCGAGCGCGGCAGACTGGACGACGTTGCTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCACGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAACACGTTTGCGCAACGCTCGCGAATTTTTctct                 >  1:260753/1‑251 (MQ=255)                                                                                                                                                                                                                        ccACACCGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATCTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATgg         <  1:87785/249‑1 (MQ=255)                                                                                                                                                                                                                             ccGCGTGCGGAAAGTGCCTGAGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATcac  >  2:108339/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                       GTTAACAACCACCATATCGATAGGCTGGATCTGATGTTCTTCCATAATGGCATCGTCCTGGCCGCGACGGCCCAGAATGCCACCATGTACTTTCGGATGCAGGGTCTTCACGCGTCCATCCATCATCTCCGGGAAACCGGTGTAATCGGAAACTTCGGTTACCGGCAGACCTTTTTCTGCTAACAGACGGGCAGTGCCCCCTGTAGACAGCAGCTCCACACCGCGTGCGGAAAGTGCCTGGGCGAATTCGACGATACCGGCTTTGTCAGAAACACTGAGCAGAGCGCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCTAACGGGGCATCCTTATTTTTCGCCCGCATTGTAACGAAAACGTTTGCGCAACGCTCGCGAATTTTTCTCTTTCAATGGTGATCAC  >  NC_000913/4207221‑4207691

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