Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	603,208	G→A	C417Y (TGC→TAC)	pheP →	phenylalanine:H(+) symporter PheP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	603,208	0	G	A	100.0%	163.5 / NA	44	C417Y (TGC→TAC)	pheP	phenylalanine:H(+) symporter PheP
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/23);  total (21/23)

GTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCG  >  NC_000913/602968‑603435                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                    gTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTAGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTc                                                                                                                                                                                                                           <  1:191706/251‑1 (MQ=255)         gtgCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTAGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCAt                                                                                                                                                                                                                   >  2:120157/1‑251 (MQ=255)               aTTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTAGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGAtttt                                                                                                                                                                                                             >  1:13127/1‑251 (MQ=255)                ttAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTAGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTtg                                                                                                                                                                                                            >  1:164556/1‑251 (MQ=255)                           aTGCTTTCCGGAGCGATCACTTCGCTGGTAGTGTTAATCAACTATCTGCTGACGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGACGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGAtgtg                                                                                                                                                                                                 >  1:48701/1‑251 (MQ=255)                               tttCCGGAGCGATCACTTCGCTGGTAGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACg                                                                                                                                                                                             >  2:334559/1‑251 (MQ=255)                                      aGCGATCACTTCGCTGGTAGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACgatggatg                                                                                                                                                                                      >  2:257272/1‑251 (MQ=255)                                        cGATCACTTCGCTGGTAGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGgatgat                                                                                                                                                                                    <  2:164556/251‑1 (MQ=255)                                        cGATCACTTCGCTGGTAGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGGGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGgatgat                                                                                                                                                                                    >  2:180753/1‑251 (MQ=255)                                             aCTTCGCTGGTAGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACgatggat                                                                                                                                                                                       >  1:335503/1‑243 (MQ=255)                                             aCTTCGCTGGTAGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACgatggat                                                                                                                                                                                       <  2:335503/243‑1 (MQ=255)                                               ttCGCTGGGAGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGCTGTGCACGATGGATGATATGCGCt                                                                                                                                                                             >  2:355234/1‑251 (MQ=255)                                                 cgcgtgtagtgtTAATCAACTTTCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCtt                                                                                                                                                                            <  2:23404/242‑1 (MQ=255)                                                              aaTCAACTATCTGCTGCCGCCAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTTGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGAACGATGGTTGATATGCGCTTGTCAGCGATCctgc                                                                                                                                                              >  2:139283/1‑251 (MQ=255)                                                               aTCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCctgct                                                                                                                                                             >  1:33255/1‑251 (MQ=255)                                                               aTCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCctgct                                                                                                                                                             <  1:58343/251‑1 (MQ=255)                                                                            ggcCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGAtt                                                                                                                                                <  2:33255/250‑1 (MQ=255)                                                                             gcCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGAt                                                                                                                                                 <  1:120157/249‑1 (MQ=255)                                                                                  aaaaaGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTc                                                                                                                                          >  2:204787/1‑251 (MQ=255)                                                                                        ggTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGttt                                                                                                                                    <  2:13127/250‑1 (MQ=255)                                                                                                         gcgcTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCt                                                                                                                   <  1:254556/251‑1 (MQ=255)                                                                                                                                     aCTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACa                                                                                       <  1:72606/251‑1 (MQ=255)                                                                                                                                     aCTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCATTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTAc                                                                                        >  1:57588/1‑250 (MQ=255)                                                                                                                                         gATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGAttt                                                                                                                                                                                                              <  1:188776/128‑1 (MQ=255)                                                                                                                                         gATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGAttt                                                                                                                                                                                                              >  2:188776/1‑128 (MQ=255)                                                                                                                                             aTGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCCGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACa                                                                                >  2:3536/1‑250 (MQ=255)                                                                                                                                               gATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCAttt                                                                                                                           >  1:55348/1‑205 (MQ=255)                                                                                                                                               gATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCAttt                                                                                                                           <  2:55348/205‑1 (MQ=255)                                                                                                                                                          cgcATCTGCGTTGTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGGCAAcc                                                                  <  2:385561/251‑1 (MQ=255)                                                                                                                                                           gcATCTGCGTTGTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACATTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACAGTGCCTTCCTCGGCATGAGTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCGGCCGGTGTGGATTGTATTCCTGTGTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCa                                                                 <  2:402501/251‑1 (MQ=255)                                                                                                                                                                 gCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATc                                                               <  1:334559/247‑1 (MQ=255)                                                                                                                                                                            cAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCAt                                                <  1:174495/251‑1 (MQ=255)                                                                                                                                                                                 aTGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCAtt                                           >  2:322218/1‑251 (MQ=255)                                                                                                                                                                                             gggCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCgg                               <  1:14577/251‑1 (MQ=255)                                                                                                                                                                                                 gTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTAGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGtttt                           <  1:3536/251‑1 (MQ=255)                                                                                                                                                                                                   gAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATTGCATTTAAAACGCTGCGTCGGAAATAATTCATTCACTCTACATCCGAGAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATGCGATGAGAGCGGtttttt                         <  2:41458/251‑1 (MQ=255)                                                                                                                                                                                                     aaCACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGttttttt                        >  1:242249/1‑250 (MQ=255)                                                                                                                                                                                                      acacAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTaa                      <  2:158368/251‑1 (MQ=255)                                                                                                                                                                                                                   cgcTGCTCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAAt         >  1:77762/1‑251 (MQ=255)                                                                                                                                                                                                                       gctCTATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGc       <  2:77762/249‑1 (MQ=255)                                                                                                                                                                                                                          ctATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCAcccc  >  1:340194/1‑250 (MQ=255)                                                                                                                                                                                                                          ctATCCGTTCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCAcc    <  1:257272/249‑1 (MQ=255)                                                                                                                                                                                                                                  tCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTAt                                                                                                                                  >  2:330316/1‑115 (MQ=255)                                                                                                                                                                                                                                  tCGGCAACTATCTCTACATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTGGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTAt                                                                                                                                  <  1:330316/115‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                    GTCGCGGTGTGCCGATTAACTCGCTGATGCTTTCCGGAGCGATCACTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCG  >  NC_000913/602968‑603435

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