Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,481,187	T→C	N884S (AAC→AGC)	valS ←	valine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,481,187	0	T	C	100.0%	168.9 / NA	48	N884S (AAC→AGC)	valS	valine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (35/13);  total (35/13)

GGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCTGCATCCGCGCTGCAACCACGCAGCAGCAGCTCCAGCGGTTTGCCCGGCGCGATGTTCATTTCTGCACGGATGTTACGT  >  NC_000913/4480954‑4481417                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                       ggCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCACTCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATc                                                                                                                                                                                                                       >  2:185469/1‑251 (MQ=255)      tCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATATCAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCggg                                                                                                                                                                                                                    <  2:345122/249‑1 (MQ=255)       cGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATc                                                                                                                                                                                                                 >  2:288039/1‑251 (MQ=255)                    gtgtTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCAc                                                                                                                                                                                                    >  2:218384/1‑251 (MQ=255)                           gATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCgat                                                                                                                                                                                             <  2:292101/251‑1 (MQ=255)                               aCAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCTAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATc                                                                                                                                                                                         >  2:252714/1‑251 (MQ=255)                               aCAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATc                                                                                                                                                                                         >  2:77063/1‑251 (MQ=255)                                 aGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCtt                                                                                                                                                                                       >  2:174034/1‑251 (MQ=255)                                    gcggcgCTAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCCGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGt                                                                                                                                                                                    >  2:238452/1‑251 (MQ=255)                                    gcggcgATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGATCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGt                                                                                                                                                                                    <  2:258440/251‑1 (MQ=255)                                         gATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACgg                                                                                                                                                                               >  1:164700/1‑251 (MQ=255)                                                       gTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCttt                                                                                                                                                                  <  1:288039/250‑1 (MQ=255)                                                           aaTCAGGTTCGCTTTCGCTTCCGCATAGCCGTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCGTTGTCAt                                                                                                                                                             <  2:308395/251‑1 (MQ=255)                                                               aGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCggc                                                                                                                                                         >  1:17081/1‑251 (MQ=255)                                                                    cgctttcgcttCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCggcaggc                                                                                                                                                     >  1:135286/1‑250 (MQ=255)                                                                         tcgcttCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCAc                                                                                                                                                <  1:371860/250‑1 (MQ=255)                                                                                 gCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACt                                                                                                                                       >  1:86502/1‑251 (MQ=255)                                                                                          tCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGa                                                                                                                                <  2:164700/249‑1 (MQ=255)                                                                                                   tcACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGAGAACGACTTCGTTCGGCAGTTTGTTATCGATACGGGTGATTTCACCTTAAATCTTCGCCACTTATTTAGACAGACGAGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGGGATACTTTCCAGACGCGACAgggg                                                                                                                     >  1:365504/1‑250 (MQ=255)                                                                                                     aCGCTCTTTCGCGATGACCGATTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCAGGGttt                                                                                                                   >  2:136067/1‑251 (MQ=255)                                                                                                       gCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGCCAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGCTACTTTCCAGACTCGCCAGGGCTTGc                                                                                                                 >  2:800/1‑251 (MQ=255)                                                                                                       gCTCTTTCGCGATGACCGCTTCCGGAGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTCCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGGCCTTTGTCATCGGCAGGCAGCACGGGGATACTTTCCAGACACGCCAGGGTTTGc                                                                                                                 >  2:815/1‑251 (MQ=255)                                                                                                            ttCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAgg                                                                                                              <  1:105921/249‑1 (MQ=255)                                                                                                                        gCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGGTTTc                                                                                                >  2:71045/1‑251 (MQ=255)                                                                                                                          ttCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCAt                                                                                              >  2:118957/1‑251 (MQ=255)                                                                                                                            ccGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCAtt                                                                                             <  1:185469/250‑1 (MQ=255)                                                                                                                                tGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTaagc                                                                                        >  2:175595/1‑248 (MQ=255)                                                                                                                                                  ttcgttcgCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCTGCATCCgcgc                                                                      >  1:150924/1‑251 (MQ=255)                                                                                                                                                     gttcgCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGCTCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATAATTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGCCGTTCTGCATCCGCCCTGc                                                                   >  2:282939/1‑251 (MQ=255)                                                                                                                                                         gCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCTGCATCCGCGCTGCAAcc                                                               >  1:344633/1‑251 (MQ=255)                                                                                                                                                           cAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCTGCATCCGCGCTGCAACCAc                                                             >  1:14757/1‑251 (MQ=255)                                                                                                                                                                   tctcGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCTGCATCCGCGCTGCAACCACgcagcagc                                                     >  2:48163/1‑251 (MQ=255)                                                                                                                                                                    ctcGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCa                                                                                                                         >  1:339136/1‑182 (MQ=255)                                                                                                                                                                    ctcGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCa                                                                                                                         <  2:339136/182‑1 (MQ=255)                                                                                                                                                                      cGATCCGGCTGAGTTCACCTTCAATCTTCGCCGATTCTTTCACCAGGCGCGCCAGCTCATCTTCTTGGCTGATGAGGACAGCCAGCGGGATCAGCTGCTCTGCACCGTCGATGATCGTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCGGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGGTCTGCATCCGCTCGTCCACCACGcaacagcag                                                   <  2:334170/250‑1 (MQ=255)                                                                                                                                                                          aCGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCTGCATCCGCGCTGCAACCACGCAGCAGCAGCTccc                                              >  2:69834/1‑250 (MQ=255)                                                                                                                                                                             gCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCTGCATCCGCGCTGCAACCACGCAGCAGCAGCTCCAGCg                                           <  1:218384/251‑1 (MQ=255)                                                                                                                                                                                  tttCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCTGCATCCGCGCTGCAACCACGCAGCAGCAGCTCCAGCGGTTTg                                      >  1:267883/1‑251 (MQ=255)                                                                                                                                                                                      aCCTTCCATCTTCGCCACTTCTTTCGACAGACGCCCCAGCTCATCTTCTTTGCTGCTGAGGACAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGAACTTTGTCATCGGCAGGCAGAACGGTGATACTTTCCATAAGCGACAGGGTTTGCAGGACGCCACGGTTTTCCTTGACGCGACGTTCAGCATCCGCGCTGCAACCACGCAGCAGCAGCTCCAGCGGTTTGcccc                                  >  2:160137/1‑250 (MQ=255)                                                                                                                                                                                      aCCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCTGCATCCGCGCTGCAACCACGCAGCAGCAGCTCCAGCGGTTTGCCCg                                  >  1:154745/1‑251 (MQ=255)                                                                                                                                                                                         ttCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCTGCATCCGCGCTGCAACCACGCCGCAGCAGCTCCAGCGGTTTGCCCGgcg                               >  1:130831/1‑251 (MQ=255)                                                                                                                                                                                               cTTCGCCACTTGTTTCGCCAGACGCGCCATCTCATCTTCTTTGCTGATCAGGCCAGCCATTGGGATCAGCAGCTGTGCACCGTTGATGATGTTCGTACCGGAAACCTGACCTTTGTCATCGGCAGGCAGCACGGTGATAGTTTCCAGACGCGGCAGGGTTTGCAGGGATCCTTTGTTTTTATTGTCGCGACGGTCGGCCTCCGCTCTCCAGCCTCTCAGCAGCCTCTCCAGCGGGTTGCCGGGCGCGAt                           <  2:161432/249‑1 (MQ=255)                                                                                                                                                                                               cTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCTGCACCCGCGCTGCAACCACGCAGCAGCAGCTCCAGCGGCTTGACCGCCGCGacgt                         >  2:82837/1‑248 (MQ=255)                                                                                                                                                                                                 tCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCTGCATCCGCGCTGCAACCACGCAGCAGCAGCTCCAGCGGTTTGCCCGGCGCGATGTTc                       >  1:234798/1‑251 (MQ=255)                                                                                                                                                                                                       cttcttTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCTGCATCCGCGCTGCAACCACGCAGCAGCAGCTCCAGCGGTTTGCCCGGCGCGATGTTCATTTCt                 >  1:208277/1‑251 (MQ=255)                                                                                                                                                                                                          cttTCGCCAGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGAACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCTGCATCCGCGCTGCAACCACGCAGCAGCAGCTCCAGCGGTTTGCCCGACGCGATGTTCATTTCTGc               >  1:331839/1‑250 (MQ=255)                                                                                                                                                                                                                  aGACGCGCCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCTGCATCCGCGCTGCAACCACGCAGCAGCAGCTCCAGCGGTTTGCCCGGCGCGCTGTTCATTTCTGCACGGATGtt      >  1:254454/1‑251 (MQ=255)                                                                                                                                                                                                                      gcgcCAGCTCATCTTCTTTGCTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCTGCATCCGCGCTGCAACCACGCAGCAGCAGCTCCAGCGGTTTGCCCGGCGCGATGTTCATTTCTGCACGGATGTTACGt  >  2:434/1‑251 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                       GGCCTTCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGCTTCCGCATAGCCTTCCAGCTTCTCACGCTCTTTCGCGATGACCGCTTCCGGTGCGCGGGCGACAAAGCCTTCGTTCGCCAGTTTGTTCTCGATACGGCTGATTTCACCTTCAATCTTCGCCACTTCTTTCGCCAGACGCGCCAGCTCATCTTCTTTGTTGATGAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCTGCATCCGCGCTGCAACCACGCAGCAGCAGCTCCAGCGGTTTGCCCGGCGCGATGTTCATTTCTGCACGGATGTTACGT  >  NC_000913/4480954‑4481417

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