Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,354	C→T	L59L (CTC→CTT)	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,354	0	C	T	100.0%	86.6 / NA	31	L59L (CTC→CTT)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/15);  total (16/15)

TAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTG‑CTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGC  >  NC_000913/4638139‑4638594                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 tagccTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTG‑CTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGc                                                                                                                                                                                                                >  1:164264/5‑250 (MQ=255)              ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTG‑CTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATaa                                                                                                                                                                                                   >  1:291218/1‑250 (MQ=255)              cagctgtCATTGCAAAGGAGAAGAATATGTTGACAACCCCCCTGTACTGTAATATGCCCACCTTGTTTGGTGTAG‑GTGGTCGTTTGTGGATTTCTATATTGCTTTTGCGGCAAGTGATTGTCTTACGTGCCGGTTACCGTAGCGATGGGGAAGTTGCGATTCACCAAAGTACCAGCGGGTCGCAAAAACGCGTTGGGCCTCGTAACTCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATaaa                                                                                                                                                                                                  <  2:202727/249‑1 (MQ=255)                gctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTG‑CTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAAc                                                                                                                                                                                                <  1:268299/251‑1 (MQ=255)                   gcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTG‑CTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGt                                                                                                                                                                                             <  1:381024/251‑1 (MQ=255)                      aTTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTG‑CTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTgg                                                                                                                                                                                           >  1:228434/1‑250 (MQ=255)                       ttGCAAAGGAGAAGACTATGTTGAAATCCGCCCTGTGCTTGAAAATGACCACCTTGTATGAGGCAGTATTG‑CTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTTGAACGTGCAGATTACCGGAGCGGTGTGGCAGATGCGATTCGCCAAAGTACCAGCGGGCCGAAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGATGTATACGGGGCAGGAAGAggtacaaaccggcagg                                                                                                                                                                                         >  2:351465/1‑237 (MQ=255)                          cAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTG‑CTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCgg                                                                                                                                                                                      >  1:284779/1‑251 (MQ=255)                                   agaCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTG‑CTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAgg                                                                                                                                                                             >  1:127580/1‑250 (MQ=255)                                          gTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTG‑CTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAGCCGTGGAGCGGAAACGAAGTTTTATcc                                                                                                                                                                      >  1:296917/1‑251 (MQ=255)                                             gAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTG‑CTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCTCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAgg                                                                                                                                                                                                             <  2:303665/209‑1 (MQ=255)                                             gAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTG‑CTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAgg                                                                                                                                                                                                             >  1:303665/1‑209 (MQ=255)                                                  ccccccTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTG‑CTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGAc                                                                                                                                                                                                                                 <  1:148538/184‑1 (MQ=255)                                                                     gtccaccttGTTTGGTGCAGTACTG‑CTGTGTTTGATTCCGTTAATGCTTATACGGCAAGTGATTGACGAACGTGCTGATGACCGTAGTGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGTGCCGCAAAAACTCGTTGGGTCGCTTATCGCGATTCCTGTGAACGAGCTTTATACGGTGCATGTAGAGGATATAACCGTGGATCGGAAACGTAGATTTATCCATTTTTGGTTACCGGAGTCAttgtta                                                                                                                                            <  2:151729/242‑4 (MQ=255)                                                                           cttGTTTGGTGCAGTATTG‑CTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCAttgatggttgatg                                                                                                                                     >  1:378130/3‑251 (MQ=255)                                                                                                       gattcctataagtttgtatcgtcaagTGATTGTCGAATGTGCCGATTCCCTTAGCGATGAGGCAGATGCGATCCGCCAAAGGACCAGCGGGCCGCAAAAACGCGTTGGTCCGCTTATCGCTATTCCTGAGACCGAGCTTTATACGGTGCAGGAAGAGGATAACACCGTGGAGCTGAAACGAAGTTTTATCCATTTTTGGTTACCTGATTCATTGATGGTTGATGGCAATCAGAACGTTGAAGATCGCAAg                                                                                                           <  2:143313/225‑1 (MQ=255)                                                                                                                  tGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTAACGTAGCGATGTGGTAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTTCAGGGAGAGGATAAAACCGTGGAGCGGAAACGAAATTTTCTCCCTTTTTGGTTACCTGGGTCATGGTTGGTTGATGGCAATCAGAAAGTGGAAGAACGCAAGATaagggtgac                                                                                                >  2:351234/1‑242 (MQ=255)                                                                                                                      gATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATAcc                                                                                            <  1:349761/250‑1 (MQ=255)                                                                                                                      gATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCg                                                                                           >  1:63570/1‑251 (MQ=255)                                                                                                                            gcaagTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAgg                                                                                     >  1:274922/5‑251 (MQ=255)                                                                                                                                      gTCGAACGTGCCGATTACCGTAGCGATGGGGAAGATTCGATTCGCTAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCTGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAAATAGGGATTTATACCGGTCAGGTCTGGCACAg                                                                           <  2:161739/251‑1 (MQ=255)                                                                                                                                                 cgATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACgc                                                                >  2:287739/2‑250 (MQ=255)                                                                                                                                                     gaccgtagctagttggatgggtcGATTCGCCAAAGTAACAGACGGACGTAAAAAATCGTGTGGCCGATTATCGCTATTCCTGTGACCGAGCTGTATACGGTGCAGGAAGAGGATAAAACCGTGGCGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTTATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGTGATTTATACCGGTCAGGTCTGGCACAGGGATTTAACGCTaa                                                             <  2:258411/228‑1 (MQ=255)                                                                                                                                                             cGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGAtt                                                    <  1:287739/251‑1 (MQ=255)                                                                                                                                                              gATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGAttt                                                   <  1:120951/251‑1 (MQ=255)                                                                                                                                                                                tgcctaaGTACCAGCGGTACGCAAAAACGCGTTGGGCTGCTTTTCGATATTCCTGTGTCCGAGCGTTATACGTTTCAGGAAGAGTATAAAACCGGTGAGCGGAAACTAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGACCGTGGAATAACTCAATATAGGGATTGATACCGTTCAGTTCTGGCACATTGATTTAACGCTAAAAGCCGAGTTCTATGTTTCTCGTCtga                                  <  2:249687/245‑3 (MQ=255)                                                                                                                                                                                    aaaGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGag                             <  1:44320/251‑2 (MQ=255)                                                                                                                                                                                            cggtgtccGCAAAAACTCGGTTGGCCGTTTTGCGCTATTCCTTTGTCTGAGCGTTATACGGTTCAGGAAGAGGATAAAACTGTGGTGCTGAAACGAAGTTTTATCCATTTTTGGTTACCTTAGTCAGTGATTGTTGATGGCATTCAGAACGTGGAAGAACGCAAGATAGCGACTTATGCCGGTCAGGTCTTTCTTAGTGCTTTAACGCTAAAAGTTGATTTAGATTTTTCGCGTCTTTGCGATCTCAAcgt                     <  2:127580/245‑2 (MQ=255)                                                                                                                                                                                                          aCTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCatct       >  1:193414/1‑248 (MQ=255)                                                                                                                                                                                                            tCGTTGGGCAGCTTATCGCTATTCCTGTGACCGCGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACGTGAGTCATTGATGGTTGAGGGATATCAGAACGTGGATGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCTATTACGATGTTTCGCGTCTTAGCGGGCTCAACGCGCAAAATATCATCtt      >  2:93614/1‑250 (MQ=255)                                                                                                                                                                                                                tGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGc  >  1:328199/1‑250 (MQ=255)                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 TAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTG‑CTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGC  >  NC_000913/4638139‑4638594

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