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

AAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAGCGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGGTGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGA  >  NC_000913/3937754‑3938239                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    aaaTCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGCACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGACGATGCAACGAAGGCGGCACAGGTACTGCATGAAAAAGGCATCCGTACTGTACTGATTACTTTAGGACGTCGTGGTGTATGGGCTAGCGGGAATGATGAAGGt                                                                                                                                                                                                                                              >  2:10664/1‑251 (MQ=255) aaaTCACCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTACGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGAGGACATTATTACGCCAACCGAACCGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGACGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTATGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAggg                                                                                                                                                                                                                                              >  2:462753/1‑250 (MQ=255)  aaTCTCCCATCAAGATAAGACGATCGGTGCGCTTCACCCGGTTCCGGCTCGCGAACTTCCTGACGGACTGCTGGCGCGGGTGGACATTATTACGCCAAACGAAACTGAAGCAGAAAAGCTCACCGTTATTCGTGTTGAAAATGATGAAGATGCAACGAAGTGGTCTCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATGACTTTATTAATTCGTTGTTGATGGGCTAGCGTGATTGATTAAGGTc                                                                                                                                                                                                                                             <  2:423393/251‑1 (MQ=255)   aTCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTAGTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCa                                                                                                                                                                                                                                            <  2:461957/251‑1 (MQ=255)                 tAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATTAAGGGCAGCGTGTTCCTGGAt                                                                                                                                                                                                                              >  2:344272/1‑251 (MQ=255)                        aTCGTTGCGCGTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGt                                                                                                                                                                                                                       <  1:109699/251‑1 (MQ=255)                          cGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGTTGTATGGGCTAGGGTGAATGATGAAGGTCAGCTCGTTCCTGGATTCCGGGTg                                                                                                                                                                                                                      <  2:163210/250‑1 (MQ=255)                          cGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTg                                                                                                                                                                                                                      <  2:128211/250‑1 (MQ=255)                               cgcTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAACCGGAAGCAGAAAAGCTCACCGGAATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGGTTCCGGGTGCAGGCt                                                                                                                                                                                                                >  2:170534/1‑251 (MQ=255)                                gcTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATGCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑g                                                                                                                                                                                                              <  1:54734/251‑1 (MQ=255)                                   tAACACGGCTCCTGCTCGCGAACTTCCGGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACTAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCCTGGAAAAAGTATCTGTACTGTACGGATTACTTTTTTAAGTCGAGGTGTACGGGTTAGCTTGGATGATGAAGGTCAGCGCGTGCCTGGATTCCGTGTGctggttgt                                                                                                                                                                                                              <  2:159241/249‑8 (MQ=255)                                    aaCCCGACTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCg                                                                                                                                                                                                           >  1:349207/1‑250 (MQ=255)                                                  gcgcgAACGTCCTAACGAACTGCTGGCGCTGGTTGACATTAGTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTTTTGAAAATGAGGAAGATGCAACGAAGGCGGCGCAGGTACTGCAGGAAAAAGGTATCCGTACTGTACTGATTGCTATAGGAAGTCGTGGGGTATGGGCTAGCGTGGATGATGAAGGTCATGGCGTTCCTGGTCTCCGGGTTCACGCT‑gatgggaccagtgctgccg                                                                                                                                                                                            <  2:220318/250‑19 (MQ=255)                                                    gcgAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCa                                                                                                                                                                                                                                            >  1:406202/1‑202 (MQ=255)                                                    gcgAACGTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAATGTATCCGTACTGTACTGAGTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGga                                                                                                                                                                                          <  1:458135/251‑1 (MQ=255)                                                     cgAACTTCCTTACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTTACCGTGATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTTCTGATTACTTTAGTAAGTCGTGGTGTATGGGCTAGCGGGAATGATGAAGTGCAGCGCGTTCCTGGATTCCGGGTGCAGGCTTGT‑GATACCATTTCTGCCGgag                                                                                                                                                                                         <  2:428839/251‑1 (MQ=255)                                                       aaCTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGAt                                                                                                                                                                                       <  2:458313/251‑1 (MQ=255)                                                        aCTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATa                                                                                                                                                                                      <  2:422023/251‑1 (MQ=255)                                                           tCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTACCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATAc                                                                                                                                                                                     <  2:355431/249‑1 (MQ=255)                                                              tGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATtgctgc                                                                                                                                                                                              >  1:394989/1‑237 (MQ=255)                                                              tGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCGGGATTCCGGGTGCAGGCT‑GTCGATACCATtgctgc                                                                                                                                                                                              <  2:394989/237‑1 (MQ=255)                                                               gACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACGGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTa                                                                                                                                                                                <  2:206827/250‑1 (MQ=255)                                                                  gaactgccagcgcaggaggccATTATTACGCCAAACGAACCGGAAGCAGAAAAGCTCCCCGGTCTTCGTGTTGCACATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTAACCGAACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCAAGCGGGATTGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTAGATACCATTGTTGCAGGAGATACCTTTAACgg                                                                                                                                                                            >  2:39650/21‑251 (MQ=255)                                                                       gctgGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATGCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGTTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGc                                                                                                                                                                         <  1:175777/249‑1 (MQ=255)                                                                                ggtggACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATTTGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACg                                                                                                                                                               <  1:264550/250‑1 (MQ=255)                                                                                     aCATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTg                                                                                                                                                         <  2:202489/251‑1 (MQ=255)                                                                                          attaCGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTgg                                                                                                                                                     >  2:276631/1‑250 (MQ=255)                                                                                               gCCAAACGAAACGGAAACAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGCAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGgaagaa                                                                                                                                                >  2:92168/1‑250 (MQ=255)                                                                                                 cAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGaaaaa                                                                                                                                             <  1:46397/251‑1 (MQ=255)                                                                                                       aaacgGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGCAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTg                                                                                                                                       >  2:228734/1‑251 (MQ=255)                                                                                                       aaacgCAAGCAGAAAAGCTCCCCGGAATTCGTGTTGAAAATGATGAAGATGCAACGAAGGAGGCGCAGGTACTGCATGAACAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGAGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGGTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATAACTTTAACGGTGCGTTAATCAAGGCATTGCTGGAACAAAAACCAtcg                                                                                                                                       >  2:235966/1‑249 (MQ=255)                                                                                                        aacgGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGAGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTAGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGCAAAACCAttcc                                                                                                                                      >  2:311215/1‑249 (MQ=255)                                                                                                                cAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCg                                                                                                                              >  1:217262/1‑251 (MQ=255)                                                                                                                   aaaaGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGAtt                                                                                                                           >  2:160699/1‑251 (MQ=255)                                                                                                                              ggTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCCTTGCCAGAGGCGATTCGTTTTGccc                                                                                                                 >  2:4334/1‑250 (MQ=255)                                                                                                                                              atgatgAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATTAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCg                                                                                                                              <  1:32595/221‑1 (MQ=255)                                                                                                                                              atgatgAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAg                                                                                                                                                                                                                                           >  2:392141/1‑113 (MQ=255)                                                                                                                                              atgatgAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAg                                                                                                                                                                                                                                           <  1:392141/113‑1 (MQ=255)                                                                                                                                                 atgaAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGcc                                                                                             >  1:70426/1‑251 (MQ=255)                                                                                                                                                   gaAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCg                                                                                            <  1:160699/250‑1 (MQ=255)                                                                                                                                                         gcaacGAAGGCGGCGCAGGTACGGCATGAAAAAGGTATCCGTACTGTACTGATTACGTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCGTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGAGGCGTTTTGCCCATGCTGCCGTTGCGATTGCCGTAacac                                                                                      <  2:428869/246‑1 (MQ=255)                                                                                                                                                                                    gAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAGCCTTCCGTACCGt                                                           <  1:138073/250‑1 (MQ=255)                                                                                                                                                                                               ccGTACTGTACTGATGTCTTTAGGAAGGCGTGGTGTATGGGCTAGCGTGAATGAGGAATGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATGCCTTTAACGGTGCGTTAATCACTGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTGGCCCATGCTGCCGCTGCGATTGCCGTAACACGGAAAGGCGCACAACCTTCCGTACCGTGGCGTGAagag                                                <  1:3207/250‑1 (MQ=255)                                                                                                                                                                                                           gATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCAtt                                     <  1:4334/249‑1 (MQ=255)                                                                                                                                                                                                            aTTACTGTAGTAAGTCGTGGTGTATGGGCTAGCGTGAATGATTAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGTT‑GTCGATACCATTGCTGCCGGAGAGCCCTTTAACGGTGCGTTAATCACGGCAATGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCTTACCGTGGCGTGAAGAGATCGACGCAttttt                                  <  1:39650/251‑1 (MQ=255)                                                                                                                                                                                                             ttACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCAttttt                                  >  1:353524/1‑250 (MQ=255)                                                                                                                                                                                                               aCTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGAGTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTAACGGGGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCTCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGGACCGTGGCGTGAAGAGATCGACGCATTTTTAg                                <  2:219824/250‑1 (MQ=255)                                                                                                                                                                                                                    aGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGc                          <  1:228734/251‑1 (MQ=255)                                                                                                                                                                                                                      gAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCag                        <  1:74733/251‑1 (MQ=255)                                                                                                                                                                                                                             tggtgAATGGGATAGCGAGACTGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTAGTTTTGCCCATGCTGCCGCTGCGATTGCAGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCAattttcacaaggcagaggggac                 >  2:359348/1‑234 (MQ=255)                                                                                                                                                                                                                                    tGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATc                                             >  1:295920/1‑216 (MQ=255)                                                                                                                                                                                                                                    tGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTGTTGCCCATGCTGCCGCTGCGAGTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCTTGGCTTGAAGAGATc                                             <  2:295920/216‑1 (MQ=255)                                                                                                                                                                                                                                          aGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACaa     <  1:36280/250‑1 (MQ=255)                                                                                                                                                                                                                                            cGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATg   >  1:313674/1‑250 (MQ=255)                                                                                                                                                                                                                                             gTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGa  >  1:28920/1‑250 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    AAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAGCGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGGTGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCT‑GTCGATACCATTGCTGCCGGAGATACCTTTAACGGTGCGTTAATCACGGCATTGCTGGAAGAAAAACCATTGCCAGAGGCGATTCGTTTTGCCCATGCTGCCGCTGCGATTGCCGTAACACGTAAAGGCGCACAACCTTCCGTACCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGA  >  NC_000913/3937754‑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: ‑