Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,240	G→A	L21L (TTG→TTA)	creD →	putative inner membrane protein CreD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,240	0	G	A	100.0%	166.1 / NA	45	L21L (TTG→TTA)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/25);  total (20/25)

GGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGA  >  NC_000913/4638008‑4638461                                                                                                                                                                                                                                         |                                                                                                                                                                                                                              ggTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGAtt                                                                                                                                                                                                             >  1:187460/1‑251 (MQ=255)      gCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCa                                                                                                                                                                                                        <  1:444528/251‑1 (MQ=255)       cGTTCGTCAGTGAGTTCGCCCGTTTGTTTAACGGCTAAGTCACCCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCTACTTCACCGTCACTTCACATAGCCTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCATTATTGCTGTTGTTAATTCCGATTATGCTGAGTCGGCAg                                                                                                                                                                                                       <  2:335357/251‑1 (MQ=255)        gTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAgg                                                                                                                                                                                                      >  1:219911/1‑251 (MQ=255)        gTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAgg                                                                                                                                                                                                      <  1:277790/251‑1 (MQ=255)           cGTCAGTGAGGTCGCCCGTTTGTTTCACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTACTTCCGATAATGCTGATTCGGCAGGTg                                                                                                                                                                                                    >  2:441517/1‑250 (MQ=255)                        gCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGAGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATCGCTTCAAATTCTTCCCACATAGTCTTCGTATACTGCTGCCATTGCAAAGGAGAAGACTATGATGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATCATGCTGATTCGGCAGGTGATTGTCGAACGTGc                                                                                                                                                                                      >  2:519902/1‑251 (MQ=255)                                  ttAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCg                                                                                                                                                                             <  1:111070/250‑1 (MQ=255)                                    aaCGGCGAAGTCAAGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAACTTCAACTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAACGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCga                                                                                                                                                                            >  2:538747/1‑248 (MQ=255)                                                     gcgcAACGTGCAGGAAGGGTGCGTGCTGTCCTCGCTTCGACTTCACCGTCACTGCACATAGCTTCAAATTCTTCCCACATAGTCTGCGTGTCCTGCTGCCATTGCAAAGGAGAAGACTAGGTTGAAATCCCCCCTGTGCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATGCCGAGAAGGCTGATGCGGCAGTTGATTGTCGAACGTGCTTATTACCGTAGCGCTGTGCAGGATGCGa                                                                                                                                                         <  2:231704/251‑1 (MQ=255)                                                                aGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGt                                                                                                                                              <  1:279753/251‑1 (MQ=255)                                                                               tGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGCAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCa                                                                                                                                >  1:462622/1‑250 (MQ=255)                                                                                ggCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCaa                                                                                                                               >  1:31719/1‑250 (MQ=255)                                                                                 gCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCATGTGATTGGCGAACGTGCTGATTGCCGTAGCGATGTGGAAGATGCGATTCGCCAAATTACCAGCGGGCCGCaaaa                                                                                                                             <  2:219911/251‑1 (MQ=255)                                                                                  ccTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCaaaaa                                                                                                                            <  2:405822/251‑1 (MQ=255)                                                                                            aCTTCCCCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGc                                                                                                                  <  1:370108/251‑1 (MQ=255)                                                                                                cACCGTCACTTCACATAGCTTAAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCt                                                                                                              <  2:83988/251‑1 (MQ=255)                                                                                                  ccGTCACTTCACATAGCCTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTc                                                                                                             >  1:456017/1‑250 (MQ=255)                                                                                                    gTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATc                                                                                                          <  1:519902/251‑1 (MQ=255)                                                                                                                     cAAATTTTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCg                                                                                          >  1:146434/1‑250 (MQ=255)                                                                                                                     cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCg                                                                                          <  1:387309/250‑1 (MQ=255)                                                                                                                     cAAATTCTACCCCCATAGTCTTCGTATCCTGCTGCCAGTGCAAAGGAGTAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCGGTTGGGTGCAGTAGTGCTGTGGTTAATTCCGATAATGCTTATTCGGCAGGTGATTGTCGATCGTGCTGTGTACCGTAGCGATGTGGAAGATGCGAGTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCGCATAGCTATGCCTGTGACCg                                                                                          <  2:252605/250‑1 (MQ=255)                                                                                                                        aTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCt                                                                                      >  1:329781/1‑251 (MQ=255)                                                                                                                             tcCCACATAGTCTTCGTAACCTGCTGTCATTGCAAAGGAGAAGACTATGTTGAAAGCCCCCCTGTTCTGGAAAATGTCTATCCTGTTGGGTGCAGTATTTCTGTTGTTAATTCCGATACTGCTGCTTCTGCAGGTGATTGTCGAACGTGCTGATTACCGTTGCGATGTGGAAGTTGCGACTCTCCAAAGTACCATCAGGCCGCACAAACTCGTGGGACCGCTCATCGCTCGTCCTGTGTCCTGGCTTtat                                                                                  <  2:80890/250‑1 (MQ=255)                                                                                                                                cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTTCCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTGTATACgg                                                                              <  2:45178/251‑1 (MQ=255)                                                                                                                                   aTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGc                                                                           >  1:450064/1‑251 (MQ=255)                                                                                                                                      gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAgg                                                                        <  1:357906/251‑1 (MQ=255)                                                                                                                                        cTTCGCATCCTGCTGCCATTGCAAAGGAGAAGACTCTGTTGAAATCCCCCGTGTTATGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTACGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTAGTCCTGTGACAGAGCCTTATCCGGTGCAGGa                                                                       >  1:232958/1‑250 (MQ=255)                                                                                                                                             tATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTCTACGGTGCAGGAAGAgg                                                                  >  2:187031/1‑250 (MQ=255)                                                                                                                                                    ctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAAcc                                                          >  1:361267/1‑251 (MQ=255)                                                                                                                                                    ctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAAATCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAAcc                                                          >  2:377818/1‑251 (MQ=255)                                                                                                                                                     tgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCg                                                         >  1:438957/1‑251 (MQ=255)                                                                                                                                                         aTTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGACTTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGCTGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGAACGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTgg                                                      >  1:529001/1‑250 (MQ=255)                                                                                                                                                                 ggAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAAc                                             <  2:187460/251‑1 (MQ=255)                                                                                                                                                                  gagaaaaatATGTTGCAATCCACCCTGTTCTGGACACTGACTAGCCTGTTTGGGGCAGTATTTCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAAGTGCTGATTACCGTAGCGATGTGGAAGAAGCAATTCGACCAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCGGCTCATCGCTATTCCTGTAACCGAGCTTTATACGGTGCAGGTAGAGGATAAAACCGTGGAGCGGAAACg                                            >  2:461496/9‑251 (MQ=255)                                                                                                                                                                    gaagaCTATGTTGAAAGCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGaa                                          <  2:361267/251‑1 (MQ=255)                                                                                                                                                                     aagaCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGAATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGGGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAACACCGTGGAGCGGAAACGacg                                         >  1:61666/1‑249 (MQ=255)                                                                                                                                                                      agaCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGGGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGt                                        <  2:456017/251‑1 (MQ=255)                                                                                                                                                                           aTGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCATGAAGAGGATAAAACCGTGGAGCGGAAACGAAGtttt                                     <  1:441517/249‑1 (MQ=255)                                                                                                                                                                                gAAATTCCCCCGGTTCTGCAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTGGTTAATTCCGATAATGCTGAGTCGGCAGGTGATTGACGAACGTGCTGATTCCCGTAGCGATGTGGAAGATGCGTTTCGCCATAGTCCCAGCTCCCCGCAAAAGCTCGTTGGGTCTCTCATCGCTATTCCTTTGACCGCGCTTTATACGGTGCAGGAAGCGGCTAATTCCGTGCAGCGGAAACGAACGTTTATCCAt                               <  2:140184/250‑1 (MQ=255)                                                                                                                                                                                            gTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACGCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTg                  <  2:31719/251‑1 (MQ=255)                                                                                                                                                                                                 ggAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAg                <  1:187031/248‑1 (MQ=255)                                                                                                                                                                                                 ggAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGAGGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGGTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAg                <  2:146434/248‑1 (MQ=255)                                                                                                                                                                                                       tGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGATGGGCCGCTCATCGCTATTCCTGTGACCGCGCTTTATACGGTGCCGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGGCAttgatg       >  2:459862/1‑251 (MQ=255)                                                                                                                                                                                                            aGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCAttgatggtgga  <  2:130524/251‑4 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                              GGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGA  >  NC_000913/4638008‑4638461

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