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%	75.7 / NA	28	L59L (CTC→CTT)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/15);  total (13/15)

TCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATAT  >  NC_000913/4638113‑4638584                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                       tCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTa                                                                                                                                                                                                                                 <  1:247743/249‑1 (MQ=255)    caTAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCt                                                                                                                                                                                                                            <  1:176622/251‑1 (MQ=255)                                 tCGTATCCTGCTGCCATTGCAAATGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGAGTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGTTGCAGGaa                                                                                                                                                                                                 <  2:193601/249‑1 (MQ=255)                                      tCCTGCTGCCATTGCAAAGGAGAAGACTATGATGAAATCCCCCCTGTTCCGGAAAATGACCACCTTGTTTGGTGAAGTATTGCTGTTGTTGATTCCGGTAATGATGCTTCGGCAAAGGATTGTCGAAGGTGCAGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGATGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATAAGATGCAGGAAGAGGAt                                                                                                                                                                                           >  1:244464/1‑250 (MQ=255)                                         tgctgcAATTGCAAAGGAGAAGACTATGTTGACATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGCTAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTTGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCCAAAACTCGTTGTGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATaccc                                                                                                                                                                                       >  1:72179/1‑248 (MQ=255)                                             gcCATTGCAAAGTAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGGACGTGCCGATTACCGTAGCGATGTGGAAGATGCGGTTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATCCGGTGCAGGAGGAGGTTAACACCgg                                                                                                                                                                                   >  2:78373/1‑250 (MQ=255)                                                                            ccccccTGTTCTGGAAAATGACCTCCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTg                                                                                                                                                     <  1:185522/250‑1 (MQ=255)                                                                                 cTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTAcc                                                                                                                                               >  1:82029/1‑251 (MQ=255)                                                                                      cTGGAAAATGACCACCTTGTTTGGTGCAGTCTTGATGTTGTTGATTCCGATCATGCTGATTCGGCAAGTGATTGTCGCACGTGCCGATTCCCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACAAGCGGGCCGCAAAAACTCGTTGGGCCGCTTCTCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGt                                                                                                                                          >  2:246886/1‑251 (MQ=255)                                                                                        ggaaaatgacccccttgtttgttGCATTATCGCTGGTGTTGATTCCTATAATGCTGAGTCGGCAAGTGATTGTCGAACGTGCCGATTGCCTTAGCGATGTGGAAGATGCGATTCGCCAAATTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGGGGAGCGGAAACGAAGGTTGATCCATTTTTGTTTAGCTggg                                                                                                                                           <  2:185190/233‑3 (MQ=255)                                                                                                     cttGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCAttgatggttgatg                                                                                                                           >  2:61258/3‑251 (MQ=255)                                                                                                                   tATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGCAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGAGGAAACGAAGATTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTg                                                                                                             >  1:232786/1‑251 (MQ=255)                                                                                                                      tGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACAGTAGCGATGTGGAAGCTGCGATTCGCCAACGTACCAGAGGGCCGCAACAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATCCGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACATGAGTCATTGATGGTTAATGGCACTCAGAACGTGga                                                                                                           >  2:269266/1‑250 (MQ=255)                                                                                                                                            gCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTtat                                                                                     >  1:60806/1‑250 (MQ=255)                                                                                                                                                        agTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCt                                                                        <  2:198728/250‑1 (MQ=255)                                                                                                                                                                    aCGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCTGTCAGGTCTGGCACAGTGAtt                                                            <  1:213367/251‑1 (MQ=255)                                                                                                                                                                        gccgATTACCGGAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCTGTCAGGTCTGGCACAGTGATTTAAc                                                        <  2:13467/248‑1 (MQ=255)                                                                                                                                                                                   tAGCGATGTGGACGATGCGATTCGCCAAAGTACCAGCGGGCCGCCCAACCTCGTTGATCCGCTTATCGCTATTCCTGTGAACGAGCTTTATACGGTGCAGGAAGAGGATAACACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCGCAGTGCTTTAACGCTAAAAGCCg                                             >  2:135482/1‑251 (MQ=255)                                                                                                                                                                                     gggATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGGTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGTGATTGATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGAt                                           <  1:251785/249‑1 (MQ=255)                                                                                                                                                                                      cGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGAtt                                          <  2:127481/251‑1 (MQ=255)                                                                                                                                                                                        aTGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCAGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGATTTCTACGGTGCAGGACGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGCTAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTc                                        >  1:186689/1‑251 (MQ=255)                                                                                                                                                                                             gAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATTCCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGt                                   <  1:78373/251‑1 (MQ=255)                                                                                                                                                                                                             aaaGTACCAGCGTGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAATAACGCAAGATAGGGATTTATACCGGTCAGGTCAGGTACAGGGATGTAACGCTAAAACCCGATTGCGATGTTTCGCGTCTTAGCGag                   <  2:180330/251‑2 (MQ=255)                                                                                                                                                                                                                   ccAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAAc             <  1:269266/251‑1 (MQ=255)                                                                                                                                                                                                                    taGCGGGACGCAAAAACTCGTAGGGCCGCTTATCGATATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATCAAACCGTGGAGCGGAACCGACGATTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGACTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAAcg            >  2:17374/2‑251 (MQ=255)                                                                                                                                                                                                                          gCCGCAAAAACTCGTTGGGCCGCTTATCGCTAGTCCTGTGACCGAGATTTATACGGTGCAGGAAGAGGATAAAACCGGGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTATCGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCa       >  1:111556/1‑250 (MQ=255)                                                                                                                                                                                                                            cGCAAAAACGCGGTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAat    <  2:64695/251‑1 (MQ=255)                                                                                                                                                                                                                              cAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTGGCGAGCTCAACGCGCCAAatat  <  2:82029/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                       TCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATAT  >  NC_000913/4638113‑4638584

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: ‑