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%	131.5 / NA	38	G235D (GGT→GAT)	rbsK	ribokinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/21);  total (17/21)

CGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAGCGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGGTGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAG  >  NC_000913/3937758‑3938242                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      cGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAgc                                                                                                                                                                                                                                            <  1:250706/251‑1 (MQ=255)               aGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTc                                                                                                                                                                                                                              <  1:45767/251‑1 (MQ=255)                   tATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCggg                                                                                                                                                                                                                          >  1:218678/1‑251 (MQ=255)                    aTCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCgg                                                                                                                                                                                                                           <  1:208666/249‑1 (MQ=255)                      cGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGc                                                                                                                                                                                                                       <  1:150734/251‑1 (MQ=255)                         tGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAgg                                                                                                                                                                                                                    >  1:26677/1‑251 (MQ=255)                         tGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAg                                                                                                                                                                                                                     <  2:43506/250‑1 (MQ=255)                                     ggctccggctcGCGAACTGCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCa                                                                                                                                                                                                        <  2:77377/251‑1 (MQ=255)                                               cgcgAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGcc                                                                                                                                                                                                <  2:130981/249‑1 (MQ=255)                                                 cgAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGgag                                                                                                                                                                                            <  2:466/251‑1 (MQ=255)                                                  gAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGgaga                                                                                                                                                                                           <  1:45587/251‑1 (MQ=255)                                                              gAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACgg                                                                                                                                                                               <  1:221538/251‑1 (MQ=255)                                                               aaCTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGCAGTCGTGGTGTATGGGCTCGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCAGGAGATACCTTTAACGGt                                                                                                                                                                              >  2:252133/1‑251 (MQ=255)                                                                       gcgcTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAAt                                                                                                                                                                      <  2:53262/251‑1 (MQ=255)                                                                          cTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCAc                                                                                                                                                                   >  2:234253/1‑251 (MQ=255)                                                                                             cAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGAGTAGAACCTGCTGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGACAAAGGTATCCGTACTGTACTGCTTACTTTAGGAAGTCGTGGTGTATGGGCTCGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATCCCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGaaaa                                                                                                                                                 >  1:151545/1‑250 (MQ=255)                                                                                                             aGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCg                                                                                                                                 >  2:102407/1‑250 (MQ=255)                                                                                                              gAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTgagcttttc                                                                                                                                                                     <  2:19397/213‑9 (MQ=255)                                                                                                              gAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTAGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTgagcttttc                                                                                                                                                                     >  1:19397/1‑205 (MQ=255)                                                                                                                aaaGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTc                                                                                                                             >  2:164913/1‑251 (MQ=255)                                                                                                                    cTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGttt                                                                                                                         >  2:78492/1‑251 (MQ=255)                                                                                                                     tCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGtttt                                                                                                                        >  2:238164/1‑251 (MQ=255)                                                                                                                                            gatgaAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAGACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGc                                                                                                 >  1:289528/1‑251 (MQ=255)                                                                                                                                            gatgaAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTg                                                                                                  <  1:247475/250‑1 (MQ=255)                                                                                                                                             atgaAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGcc                                                                                                <  1:81092/251‑1 (MQ=255)                                                                                                                                                   aTGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAac                                                                                           >  1:227192/1‑250 (MQ=255)                                                                                                                                                     gcaacGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGAGTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACg                                                                                        <  1:31189/247‑1 (MQ=255)                                                                                                                                                      caacGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGTTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTGTAACGGTGCGTTAAGCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTGTGAGGATGCTGCCGCTGCGATTGCCGTAACACGt                                                                                       <  2:202146/248‑1 (MQ=255)                                                                                                                                                                     aGGTATTGCATGAAAAAGGCATCCGTACTGTACTGATTACTTTAGGAAGTCGAGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGCAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGGAAAGGCGCACAACCt                                                                        >  1:62191/1‑251 (MQ=255)                                                                                                                                                                        tACTGAATGAAAACGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTcc                                                                     >  2:129339/1‑251 (MQ=255)                                                                                                                                                                                         aTCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGa                                                       >  1:265889/1‑248 (MQ=255)                                                                                                                                                                                           ccGTACTGTACTGATTACTTTAGGAAGTCGTGGAGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAag                                                     >  1:304270/1‑248 (MQ=255)                                                                                                                                                                                             gtactgtactGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATc                                                <  1:261416/251‑1 (MQ=255)                                                                                                                                                                                             gaactgtactgATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAag                                                     >  2:58970/3‑246 (MQ=255)                                                                                                                                                                                                           aCTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAg                                   <  2:246960/250‑1 (MQ=255)                                                                                                                                                                                                              ttAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGAc                                 <  1:274283/249‑1 (MQ=255)                                                                                                                                                                                                                aGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCGTCCGTACCGTGGCGTTAAGAGATCGACGCATTTGTAGACAgg                              <  1:267387/250‑1 (MQ=255)                                                                                                                                                                                                                                           gAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAg  <  2:102158/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      CGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAGCGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGGTGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAG  >  NC_000913/3937758‑3938242

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