Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,956,739	C→T	T470I (ACC→ATC)	ilvA →	threonine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,956,739	0	C	T	100.0%	167.6 / NA	47	T470I (ACC→ATC)	ilvA	threonine deaminase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/22);  total (25/22)

CGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCACCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCG  >  NC_000913/3956497‑3956970                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                        cGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTAc                                                                                                                                                                                                                                  <  2:71338/250‑1 (MQ=255)    agagCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCCTCGACTACTGgc                                                                                                                                                                                                                              >  2:324569/1‑251 (MQ=255)     gagCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGgcg                                                                                                                                                                                                                             >  1:348877/1‑251 (MQ=255)       gcgcAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGcgcg                                                                                                                                                                                                                           <  1:312129/251‑1 (MQ=255)         gcAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGt                                                                                                                                                                                                                          >  1:291279/1‑250 (MQ=255)           aaaGAAATTTTGCAGATGCACAACGACGGCGGCTACAGCGTGGTGGATCTCTCCGACGACGAAAGGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGGACTTGAACATGTCGTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACt                                                                                                                                                                                                                       <  1:62522/251‑1 (MQ=255)                       cAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCg                                                                                                                                                                                                            >  1:71936/1‑250 (MQ=255)                         gATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAAc                                                                                                                                                                                                         >  1:50157/1‑251 (MQ=255)                           tGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCCCCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTAAGGGCGCGTACTGGCGGCGTTCGCACtt                                                                                                                                                                                                       >  2:225670/1‑251 (MQ=255)                            gCTCAACTACGGCGGCTACAGCTTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTTCGCTATATGGTCGGCGGACGTCCAGCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCTCTCTGCTGCGCTTCCACAACACGCTGGGTACGTACTGGAACATTTCTTGGTTCCACTATCGGAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACtt                                                                                                                                                                                                       <  1:51522/250‑1 (MQ=255)                                   gacgGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACCTGCGCTATATGGTCGGCGGACGTCCCTCGCATCCGTTGCAGGAACCCCTCTACCGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTTTTTGTTCCACTATCGCAGCCATGGCCTCGACTACCGGCGCGTCCTTGCGGCGTTTTAACTTggagcccc                                                                                                                                                                                               >  2:321291/1‑245 (MQ=255)                                       gcggcTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGaa                                                                                                                                                                                           <  2:356065/251‑1 (MQ=255)                                             aCAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGAt                                                                                                                                                                                     >  1:277460/1‑251 (MQ=255)                                               aGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACt                                                                                                                                                                                                                                    <  1:222137/202‑1 (MQ=255)                                               aGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACt                                                                                                                                                                                                                                    >  2:222137/1‑202 (MQ=255)                                                    ggTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGaaa                                                                                                                                                                              <  2:163306/251‑1 (MQ=255)                                                    ggTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTACAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGGACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGCATTCGaa                                                                                                                                                                               >  2:107611/1‑250 (MQ=255)                                                                cgacgacgaAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGAGTTCGAAACCCGGCTGAATg                                                                                                                                                                  <  1:324569/251‑1 (MQ=255)                                                                 gacgacgaAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATg                                                                                                                                                                  >  1:176039/1‑250 (MQ=255)                                                                            ggCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACg                                                                                                                                                      <  1:345837/251‑1 (MQ=255)                                                                            ggCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACg                                                                                                                                                      >  2:395822/1‑251 (MQ=255)                                                                              cGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGACACCCGGCTGACTGAGCTGGGCTACGAt                                                                                                                                                    >  1:252788/1‑251 (MQ=255)                                                                                aaGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTTCAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTGCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCAAGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTg                                                                                                                                                  <  2:252788/251‑1 (MQ=255)                                                                                       acGTGCGCTATAGGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCATCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCacgac                                                                                                                                           <  2:297320/251‑1 (MQ=255)                                                                                            cgcTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAAc                                                                                                                                      >  2:212210/1‑251 (MQ=255)                                                                                                    ggTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAAcc                                                                                                                              <  2:348877/251‑1 (MQ=255)                                                                                                               gTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCTCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGTCTACGATTGCCTCGTCGAAACCAATAACCCGGCGTTCAgg                                                                                                                   <  1:158416/251‑1 (MQ=255)                                                                                                                      gCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCtttt                                                                                                            >  2:68790/1‑251 (MQ=255)                                                                                                                      gCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCtttt                                                                                                            >  2:56799/1‑251 (MQ=255)                                                                                                                                cAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTGCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTa                                                                                                  <  2:285718/251‑1 (MQ=255)                                                                                                                                           ctACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAAt                                                                                         <  1:183314/249‑1 (MQ=255)                                                                                                                                                gCTTCGAATTCCCGGAATCCCCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTTGCGGCGTTCGAACTTGGCGACCATGATCCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGAt                                                                                  <  2:81433/251‑1 (MQ=255)                                                                                                                                                gCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGCAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAAttac                                                                                      >  2:290163/1‑244 (MQ=255)                                                                                                                                                    cGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAgc                                                                               <  2:152240/250‑1 (MQ=255)                                                                                                                                                      aaTTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCCGGTTCTTTTTGGCGGGGTAGGGAAAAATGCATGAAAGcgct                                                                            >  1:251541/1‑251 (MQ=255)                                                                                                                                                           ccGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGGTAGGGAAAAAAGCCTGATAGcgcttcgct                                                                       >  1:257818/1‑251 (MQ=255)                                                                                                                                                           ccGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGAACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCCGATAGcgattcgct                                                                       >  1:274984/1‑251 (MQ=255)                                                                                                                                                               aaTCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAAACAATAACCCCGCGTTCAGGTTCTTTTTGGCGGGTTTGGGGAAAATGCCTgcacacgcttcccttctt                                                                   >  2:152635/1‑234 (MQ=255)                                                                                                                                                               aaTCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAAACGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGACTGATAGCGCCTCGCTTATc                                                                   >  1:55632/1‑251 (MQ=255)                                                                                                                                                                atccccGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGGTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCa                                                                  <  2:302535/247‑1 (MQ=255)                                                                                                                                                                       ggCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTTGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTAc                                                           <  2:352561/251‑1 (MQ=255)                                                                                                                                                                        gcgcgcTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTAcc                                                          <  2:81733/251‑1 (MQ=255)                                                                                                                                                                                       tCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCCCGACAACGTc                                            >  2:366/1‑250 (MQ=255)                                                                                                                                                                                             acacGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTtgtg                                     >  1:222704/1‑251 (MQ=255)                                                                                                                                                                                                                          ttGTTCCACTATCGCAGACATGGAATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACAGCATTTCGAGACCCGGCTGAATGAGCGGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTATGGAAAAATGCCTTATAGCGCTTCGCCTATCAGCCATACCCGTGCGACAACTTCATTTGTTGTTCGGCAAAATCTTCCAGAATTCCt           <  2:70916/248‑1 (MQ=255)                                                                                                                                                                                                                               ccACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCACTAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCCGAAAGCCCTTCGGTTATTAGGTCTACCCGCGGGGCAACCGCATTTTTGGGTAGGCAAAATCTTTCAGGATGGCACAAttc     >  2:313960/1‑248 (MQ=255)                                                                                                                                                                                                                                cACTATCGCAGCCATGGCATCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGGTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTGCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGGGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCg  <  1:170326/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                        CGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCACCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCG  >  NC_000913/3956497‑3956970

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