Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,937,997	G→A	G235D (GGT→GAT)	rbsK →	ribokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,937,997	0	G	A	100.0%	190.4 / NA	57	G235D (GGT→GAT)	rbsK	ribokinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/33);  total (24/33)

CGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAGCGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGGTGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGA  >  NC_000913/3937751‑3938239                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   cgAAAATCGCCCATCAAAATAAGACTACCGTTGCGCTTAACCCGGCTCCGGCTCGCTAACTTCCTGACGAACTTCTGGCTCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTTTTGAAAATGATGAAGATGCTACGAAGGCGTCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTTATTACTTTAGTAAGTCTTGGTGTATTGGCTAGCGTGAATGATGaa                                                                                                                                                                                                                                                <  2:316518/251‑1 (MQ=255)   aaaaTCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAgg                                                                                                                                                                                                                                              <  1:393958/251‑1 (MQ=255)                  aaTAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTgg                                                                                                                                                                                                                               >  1:255391/1‑251 (MQ=255)                   aTAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGaa                                                                                                                                                                                                                                                <  1:288431/233‑1 (MQ=255)                   aTAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGAATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGaaagatcggaagggcgtcgg                                                                                                                                                                                                                              >  2:288431/1‑233 (MQ=255)                             cGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGGCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGc                                                                                                                                                                                                                    <  2:153817/251‑1 (MQ=255)                                  cgcTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCt                                                                                                                                                                                                               <  2:77493/251‑1 (MQ=255)                                       aaCCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGCAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGAGCAGGCTGTCg                                                                                                                                                                                                           >  1:357084/1‑250 (MQ=255)                                          ccggctccggctCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATGCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGAGTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATAc                                                                                                                                                                                                       <  1:38060/251‑1 (MQ=255)                                          ccggcgccggctCTCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTTAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGAttc                                                                                                                                                                                                       <  2:175981/251‑3 (MQ=255)                                                      cgcgAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGAGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCgg                                                                                                                                                                                           >  2:336018/1‑251 (MQ=255)                                                        cgAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGgag                                                                                                                                                                                         <  1:208903/251‑1 (MQ=255)                                                            cTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATAc                                                                                                                                                                                     >  2:274999/1‑251 (MQ=255)                                                                   aCGAACTGCGAGCGCTGGTGGTCATTATGACGCCAAACTAAAAGGAAGCAGAAAAGCGCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGATGGCGGCGCTGGGACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGACTCCGGGTGCAGGCTGTCGAGACCATTGCTGCCGgagttcccattaa                                                                                                                                                                               <  2:105867/250‑11 (MQ=255)                                                                     gAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGAAGCAAAGCTCACCGGTATTCGTGTTGACAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACgg                                                                                                                                                                            >  1:267261/1‑251 (MQ=255)                                                                             ggCGCTGGTGGACATTATTACGCCAAACGACACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTa                                                                                                                                                                     >  2:84627/1‑250 (MQ=255)                                                                             ggCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTa                                                                                                                                                                     >  1:234661/1‑250 (MQ=255)                                                                               cgcTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATc                                                                                                                                                                  <  1:392258/251‑1 (MQ=255)                                                                                    gtggACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCTATACCATTGCTGCCGGAGATTCCTTTAACGGTGCGTTAATCACGGc                                                                                                                                                             <  2:73181/251‑1 (MQ=255)                                                                                     tggACATTATTACGCCACACGAAACGGAAGCAGAAAAGCTCACCGGCATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGAACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGcc                                                                                                                                                            >  2:308588/1‑250 (MQ=255)                                                                                      ggACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCa                                                                                                                                                            <  2:124088/250‑1 (MQ=255)                                                                                       gACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCAtt                                                                                                                                                          <  2:131779/251‑1 (MQ=255)                                                                                                         gaaacgGAAGCCGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTCTCCGTACTGTACTGATTACTTTAGGCAGTCGTGGAGTATGGGCTAGCGTGAATGATGACGGTAAGCGCGTTCCTGGCTTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAACCCAtt                                                                                                                                        >  2:270682/1‑251 (MQ=255)                                                                                                              ggAAGAAGACAAGCTCACCGGTATTCGTGTTGAATATGATGAAGATGCAACGACGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTAATGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGACTGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGcc                                                                                                                                     >  2:260440/1‑249 (MQ=255)                                                                                                               gAAGCAGAAAAGCTCCCCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCAt                                                                                                                                         >  1:303473/1‑244 (MQ=255)                                                                                                               gAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCAt                                                                                                                                         <  2:303473/244‑1 (MQ=255)                                                                                                                aaGCAGGAAAGCTCTCCGTGATTCGTGTTGAAAATGATGAAGAGGCAACGAATGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTGCTGTACTGATTACTTTAGGAAGTCTGGGTGGATGGGCGAGCTTTAATGATGAAGGTCAGCGCGTTCCTGGATTCCAGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTATAACGGTGCGTTAATCACGGCAGTGCTGGAAGAAAAACCATTTCCagag                                                                                                                                 <  2:128066/251‑1 (MQ=255)                                                                                                                      aaaaACTCAACGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGAGGTGTATGGACTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGCAGAAAAACCATGGCCAGAGACGAtt                                                                                                                           >  2:302092/1‑251 (MQ=255)                                                                                                                           cTCACCGGTATTCGTGTTGAAAATGATGCAGATGCAACGAAGACGGCGCAGGTACTGAATGAAAAAGGAATCCGTACTGTACTGATTACTTTAGGAAGTAGTGGTGTATGGGCTAGCGGGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGACGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGACAGAGGCGATTCGttt                                                                                                                      >  1:105187/1‑251 (MQ=255)                                                                                                                              cccGGTATTCGTGTTGAAAATGATGAAGAGGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCGGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGc                                                                                                                   <  2:274678/250‑1 (MQ=255)                                                                                                                                            tGAAAATGATGAAGATGCAACGAAGGCGGCGCAGTTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTTAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCAGTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCt                                                                                                      <  2:195269/250‑1 (MQ=255)                                                                                                                                                      gaAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGCGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGt                                                                                           >  2:184306/1‑251 (MQ=255)                                                                                                                                                        aGATGCACCGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTa                                                                                          >  1:312591/1‑250 (MQ=255)                                                                                                                                                         gATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAac                                                                                        <  1:394081/251‑1 (MQ=255)                                                                                                                                                                cGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTaaa                                                                                 <  1:94471/251‑1 (MQ=255)                                                                                                                                                                  aaGGCGGCTCAGGGACTGCATGAAAAAGGTATCCGTACGGTACTGATTACTTTAGGATGTCTTGGTGTAGGGGCTAGCGGGAATGATGAAGTTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTGAACGGTGCGTTAATCCCGGCATTGCTGGAAGAAAAACCAGTGCTGGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACTTAAAgg                                                                               <  2:365316/251‑1 (MQ=255)                                                                                                                                                                      cggcgCATGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGca                                                                           <  1:394698/251‑1 (MQ=255)                                                                                                                                                                          gcAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTCACACGTAAAGGCGCACAAc                                                                       >  2:23656/1‑251 (MQ=255)                                                                                                                                                                                       gAAAAAGGGATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTAGGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGGTAATCACGGCATTGCTTGCAGAAAAACAATTGCCAGAGGCGATTGGTTTTGCCCATGCTGCCGCTGAGATTGCCGTAAGACGTAAAGGCGCACAACCTACCGAACCGTg                                                          >  2:90515/1‑251 (MQ=255)                                                                                                                                                                                        aaaaaGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTGAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTga                                                         <  2:24876/251‑2 (MQ=255)                                                                                                                                                                                           aaGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGt                                                      >  2:34709/1‑251 (MQ=255)                                                                                                                                                                                                  ccGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAag                                                  >  2:57028/1‑248 (MQ=255)                                                                                                                                                                                                   cggactgtactgATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGAt                                              <  1:194684/248‑1 (MQ=255)                                                                                                                                                                                                      actgtactgATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTGTCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGGATTGCGGGAAGAAAATCCATTGCCAGAGGCGATTCTTTTTGCCCATGCTGCCGCTGCGATTGACGTAACACGTAAAGGCGCACAACCGGTCGTACCGTGGCGTGAAgagagcg                                            <  2:232884/250‑4 (MQ=255)                                                                                                                                                                                                        tgtactgATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGTAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGc                                                                                                                   <  2:266989/177‑1 (MQ=255)                                                                                                                                                                                                        tgtactgATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGTAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGc                                                                                                                   >  1:266989/1‑177 (MQ=255)                                                                                                                                                                                                               aTTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTaaa                                                                                 <  1:59448/204‑1 (MQ=255)                                                                                                                                                                                                               aTTACTTTAGGAAGGCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTaaa                                                                                 >  2:59448/1‑204 (MQ=255)                                                                                                                                                                                                                 tACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTa                                 <  1:57028/250‑1 (MQ=255)                                                                                                                                                                                                                         gAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCag                        <  2:312591/251‑1 (MQ=255)                                                                                                                                                                                                                          aaGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCaga                       <  1:402240/251‑1 (MQ=255)                                                                                                                                                                                                                                 ggtgTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACg                >  2:147060/1‑251 (MQ=255)                                                                                                                                                                                                                                         ggCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCt         >  1:221144/1‑250 (MQ=255)                                                                                                                                                                                                                                          gcgaGCGGGAATGATGAAGGGCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTAc       <  2:255391/248‑1 (MQ=255)                                                                                                                                                                                                                                          gCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTTCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGGGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCCTGGCTAc       <  2:228951/251‑1 (MQ=255)                                                                                                                                                                                                                                              gCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGACGTAACACGTAAAGGCGCACAACCTTCCGTGCCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATg   <  2:234661/251‑1 (MQ=255)                                                                                                                                                                                                                                               cGTGAATGATGAAGGTCAGCGCGTTCGTGGATTCCGGGTGCAGGCTGTCGAGACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGGTTTGCCCAGGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAAGCTTCCGTACCGTGGCGGGAAGAGAACGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGa  <  1:180334/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   CGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAGCGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGGTGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGA  >  NC_000913/3937751‑3938239

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