Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,897,890	C→T	L129L (CTG→TTG)	cbrB →	putative inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,897,890	0	C	T	100.0%	92.5 / NA	26	L129L (CTG→TTG)	cbrB	putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/14);  total (12/14)

CGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCTGGTGGCATAGGTGGCGTCGTTATCACCGAGATCTATTGTGCAGTTATTGTACATATTAAGGGCATGGCTTCCTCGGAGTTGTTTGAAAACATTCTTTCTGGTGACAGTCTCGTTGTCCGCATCATTCCTGCATTGCTGGCAGGTGTGGTGATGAGCAGAATCATTACCCGTCTACCCGGATTGGATATTTCATGTCCTGAAACAGACTCTTTAAGTTAAGCGGGATACTTTATCTTTGGGCTACTCAAAAGCAGACAGGATGTTTCTATGACTCAAAATATCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCG  >  NC_000913/3897657‑3898115                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                  cGATACCTACGGGAGTTATTCCTGCGTTGTTACCCGGTGTAACGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAAAACCTTTATCGTTGCCTGGCTGGTTGCATAGGTGGCTTCGTT‑TCACCGAGATCTATTGTGCAGTTATTGTACATATTAAGGGCATGGCTTCCTCGGAGTTGTTTGAAAACATTCTTTCTGGTGACAGTCTCGTTGTCCGCATCATTCCTGCAGTGTTGGCAGGTGTGGtgatg                                                                                                                                                                                                                  <  1:139315/250‑1 (MQ=255) cGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCTGGTGGCATAGGTGGCGTCGTTATCACCGAGATCTATTGTGCAGTTATTGTACATATTAAGGGCATGGCTTCCTCGGAGTTGTTTGAAAACATTCTTTCTGGTGACAGTCTCGTTGTCCGCATCATTCCTGCATTGTTGGCAGGTGTGGtgatg                                                                                                                                                                                                                  <  2:131581/251‑1 (MQ=255)   aTAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCTGGTGGCATAGGTGGCGTCGTTATCACCGAGATCTATTGTGCAGTTATTGTACATATTAAGGGCATGGCTTCCTCGGAGTTGTTTGAAAACATTCTTTCTGGTGACAGTCTCGTTGTCCGCATCATTCCTGCATTGTTGGCAGGTGTGGtgatg                                                                                                                                                                                                                  >  1:162210/1‑249 (MQ=255)    tAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCTGGTGGCATAGGTGGCGTCGTTATCACCGAGATCTATTGTGCAGTTATTGTACATATTAAGGGCATGGCTTCCTCGGAGTTGTTTGAAAACATTCTTTCTGGTGACAGTCTCGTTGTCCGCATCATTCCTGCATTGTTGGCAGGTGTGGGGATGAGc                                                                                                                                                                                                               >  2:104046/1‑251 (MQ=255)                        gCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCTGGTGGCATAGGTGGCGTCGTTATCACCGAGATCTATTGTGCAGTTATTGTACATATTAAGGGCATGGCTTCCTCGGAGTTGTTTGAAAACATTCTTTCTGGTGACAGTCTCGTTGTCCGCATCATTCCTGCATTGTTGGCAGGTGTGGTGATGAGCAGAATCATTACCCGTCTAcc                                                                                                                                                                                           <  1:361332/251‑1 (MQ=255)                                  ccGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCTGGTGGCATAGGTGGCGTCGTTATCACCGAGATCTATTGTGCAGTTATTGTACATATTAAGGGCATGGCTTCCTCGGAGTTGTTTGAAAACATTCTTTCTGGTGACAGTCTCGTTGTCCGCATCATTCCTGCATTGTTGGCAGGTGTGGTGATGAGCAGAATCATTACCCGTCTACCCGGATTGGat                                                                                                                                                                                 <  1:63743/251‑1 (MQ=255)                                             gTTGCCTGTCTGCCGGAAAAGATCGGGTCACACACACGTTATCGTTGTCTGGCTGGTGGCATAGGTGGCGTCGTTATCACCGAGATCTATTGTGCAGTTATTGTACATATTAAGGGCATGGCTTCCTCGGAGTTGTTTGAAAACATTCTTTCTGGTGACAGTCTCGTTGTCCGCATCATTCCTGCATTGTTGGTAGGTGTGGTGATGAGCAGAATCATTACCCGTCTACCCGGATTGGATATTTCATGTcc                                                                                                                                                                      >  1:333066/1‑251 (MQ=255)                                                            gAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCTGGTGGCATAGGTGGCGTCGTTATCACCGAGATCTATTGTGCAGTTATTGTACATATTAAGGGCATGGCTTCCTCGGAGTTGTTTGAAAACATTCTTTCTGGTGACAGTCTCGTTGTCCGCATCATTCCTGCATTGTTGGCAGGTGTGGTGATGAGCAGAATCATTACCCGTCTACCCGGATTGGATATTTCATGTCCTGAAACAGGCTCttt                                                                                                                                                       <  1:230322/251‑1 (MQ=255)                                                            gAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCTGGTGGCATAGGTGGCGTCGTTATCACCGAGATCTATTGTGCAGTTATTGTACATATTAAGGGCATGGCTTCCTCGGAGTTGTTTGAAAACATTCTTTCTGGTGACAGTCTCGTTGTCCGCATCATTCCTGCATTGTTGGCAGGTGTGGTGATGAGCAGAATCATTACCCGTCTACCCGGATTGGATATTTCATGTCCTGAAACAGACTCttt                                                                                                                                                       <  1:194005/251‑1 (MQ=255)                                                                         cacaGAAACGTTATCGTTGTCTGGCTGGTGGCATAGGTGGCGTCGTTATCACCGAGATCTATTGTGCAGTTATTGTACATATTAAGGGCATGGCTTCCTCGGAGTTGTTTGAAAACATTCTTTCTGGTGACAGTCTCGTTGTCCGCATCATTCCTGCATTGTTGGCAGGTGTGGTGATGAGCAGAATCATTACCCGTCTACCCGGATTGGATATTTCATGTCCTGAAACAGACTCTTTAAGTTAAGCGGGa                                                                                                                                          <  1:359059/251‑1 (MQ=255)                                                                            aGAAACGTTATCGTTGTCTGGCTGGTGGCATAGGTGGCGTCGTTATCACCGAGATCTATTGTGCAGTTATTGTACATATTAAGGGCATGGCTTCCTCGGAGTTGTTTGAAAACATTCTTTCTGGTGACAGTCTCGTTGTCCGCATCATTCCTGCATTGTTGGCAGGTGTGGTGATGAGCAGAATCATTACCCGTCTACCCGGATTGTATCTTTCATGTCCTGAAACAGACTCTTTAATTTAAGCGGGCTAc                                                                                                                                       >  2:207055/1‑251 (MQ=255)                                                                                    tATCGTTGTCTGGCTGGTGGCATAGGTGGCGTCGTTATCACCGAGATCTATTGTGCAGTTATTGTACATATTAAGGGCATGGCTTCCTCGGAGTTGTTTGAAAACATTCTTTCTGGTGACAGTCTCGTTGTCCGCATCATTCCTGCATTGTTGGCAGGTGTGGTGATGAGCAGAATCATTACCCGTCTACCCGGATTGGATATTTCATGTCCTGAAACAGACTCTTTAAGTTAAGCGGGATACTTTATCtt                                                                                                                               >  1:106837/1‑251 (MQ=255)                                                                                        gTTGTCTGGCTGGTGGCATAGGTGGCGTCGTTATCACCGAGATCTATTGTGCAGTTATTGTACATATTAAGGGCATGGCTTCCTCGGAGTTGTTTGAAAACATTCTTTCTGGCGACAGTCTCGTTGTCCGCATCATTCCTGCATTGTTGGCAGGTGTGGTGATGAGCAGAATCATTACCCGTCTACCCGGATTGGATATTTCATGTCCTGAAACAGACTCTTTAAGTTAAGTGGGATAATTTATCTTTggg                                                                                                                           >  1:36568/1‑251 (MQ=255)                                                                                                    gtggCATAGGTGGCGTCGTTATCACCGAGATCTATTGTGCAGTTATTGTACATATTAAGGGCATGGCTTCCTCGGAGTTGTTTGAAAACATTCTTTCTGGTGACAGTCTCGTTGTCCGCATCATTCCTGCATTGTTGGCAGGTGTGGTGATGAGCAGAATCATTACCCGTCTACCCGGATTGGATATTTCATGTCCTGAAACAGACTCTTTAAGTTAAGCGGGCTACTTTATCTTTGGGCTACTCAAAAGc                                                                                                               >  1:181791/1‑251 (MQ=255)                                                                                                                                gaTCTATTGTGCAGTTATTGTACATATTAAGGGCATGGCTTCCTCGGAGTTGTTTGAAAACATTCTTTCTGGTGACAGTCTCGTTGTCCGCATCATTCCTGCATTGTTGGCAGGTGTGGTGATGAGCAGAATCATTACCCGTCTACCCGGATTGGATATTTCATGTCCTGAAACAGACTCTTTAAGTTAAGCGGGATACTTTATCTTTGGGCTACTCAAAAGCAGACAGGATGTTTCTATGACTCAAAat                                                                                    >  1:251813/1‑250 (MQ=255)                                                                                                                                     aTTGTGCAGTTATTGTACATATTAAGGGCATGGCTTCCTCGGAGTTGTTTGAAAACATTCTTTCTGGTGACAGTCTCGTTGTCCGCATCATTCCTGCATTGTTGGCAGGTGTGGTGATGAGCAGAATCATTACCCGTCTACCCGGATTGGATATTTCATGTCCTGAAACAGACTCTTTAAGTTAAGCGGGATACTTTATCTTTGGGCTACTCAAAAGCAGACAGGATGTTTCTATGACTCAAAata                                                                                   <  1:63892/246‑1 (MQ=255)                                                                                                                                           cAGTTATTGTACATATTAAGGGCATGGCTTCCTCGGAGTTGTTTGAAAACATTCTTTCTGGTGACAGTCTCGTTGTCCGCATCATTCCTGCATTGTTGGCAGGTGTGGTGATGAGCAGAATCATTACCCGTCTACCCGGATTGGATATTTCATGTCCTGAAACAGACTCTTTAAGTTAAGCGGGATACTTTATCTTTGGGCTACTCAAAAGCAGACAGGATGTTTCTATGACTCAAAATATCAGGCCGTTa                                                                        <  2:70775/251‑1 (MQ=255)                                                                                                                                               tATTGTACATATTAAGGGCATGGCTTCCTCGGAGTTGTTTGAAAACATTCTTTCTGGTGACAGTCTCGTTGTCCGCATCATTCCTGCATTGTTGGCAGGTGTGGTGATGAGCAGAATCATTACCCGTCTACCCGGATTGGATATTTCATGTCCTGAAACAGACTCTTTAAGTTAAGCGGGATACTTTATCTTTGGGCTACTCAAAAGCAGACAGGATGTTTCTATGACTCAAAATATCAGGCCGTTAccc                                                                     <  1:230846/250‑1 (MQ=255)                                                                                                                                                   gtaaataatAAGGGCATGGCTCCCTCGGAGTTGTTTGAAAACATTCTTTCTGGTGACAGTCTCGATGTCCGCATCATTCCTGCATTGTTGGCAGGTGGGGTGATGAGCAGAATCATTACCCGTCTACCCGGATTGGAAATTTCATGTCCTGAAACAGACTCTTTAGGTTAAGCGGGATACTTTATCTTGGGGATACTCAAAAGCAGCCAGGATGTTTCTATGACTCAAAATATAAAGCCGTTCACCCAAtt                                                                >  2:204261/9‑251 (MQ=255)                                                                                                                                                         atTAAGGGCATGGCTTCCTCGGAGTTGTTTGAAAACATTCTTTCTGGTGACAGTCTCGTTGTCCGCATCATTCCTGCATTGTTGGCAGGTGTGGTGATGAGCAGAATCATTACCCGTCTACCCGGATTTGATATTTCATGTCCTGAAACAGACTCTTTAAGTTAAGCGGGATACTTTATCTTTGGGCTACTCAAAAGCAGACAGGATGTTTCTATGACTCAAAATATCAGGCCGTTTCCCCAATTCAAata                                                          <  2:33168/251‑1 (MQ=255)                                                                                                                                                                  aTGGCTTCCTCGGAGTTGTTTGAAAACATTCTTTCTGGTGACAGTCTCGTTGTCCGCATCATTCCTGCATTGTTGGCATGTGTGGTGATGAGCAGAATCATTACCCGTCTACCCGGATTGGATATTTCATGTCCTGAAACAGACTCTTTAAGTTAAGCGGGATACTTTATCTTTGGGCTACTCAAAAGCAGACAGGATGTTTCTATGACTCAAAATATCAGGCCGTTACCCCAATTCAAATATCATCCCaa                                                 <  1:293119/251‑1 (MQ=255)                                                                                                                                                                             ggAGTTGTTTGAAAACATTCTTTCTGGTGACAGTCTCGTTGTCCGCATCATTCCTGCATTGTTGGCAGGTGTGGTGCTGAGCAGAATCATTACCCGTCTACCCGGATTGGATATTTCATGTCCTGAAACAGACTCTTTAAGTTAAGCGGGATACTTTATCTTTGGGCTACTCAAAAGCAGACAGGATGTTTCTATGACTCAAAATATCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGaa                                       >  1:382594/1‑250 (MQ=255)                                                                                                                                                                              gAGTTGTTTGAAAACATTCTTTCTGGTGACAGTCTCGTTGTCCGCATCATTCCTGCATTGTTGGCAGGTGTGGTGATGAGCAGAATCATTACCCGTCTACCCGGATTGGATATTTCATGTCCTGAAACAGACTCTTTAAGTTAAGCGGGATACTTTATCTTTGGGCTACTCAAAAGCAGACAGGATGTTTCTATGACTCAAAATATCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAAc                                     >  1:224018/1‑251 (MQ=255)                                                                                                                                                                                          aaCATTCTTTCTGGTGACAGTCTCGTTGTCCGCATCATTCCTGCATTGTTGGCAGGTGTGGTGATGAGCAGAATCATTACCCGTCTACCCGGATTGGATATTTCATGTCCTGAAACAGACTCTTTAAGTTAAGCGGGATACTTTATCTTTGGGCTACTCAAAAGCAGACAGGATGTTTCTATGACTCAAAATATCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTg                          >  1:284270/1‑250 (MQ=255)                                                                                                                                                                                            cATTCTTTCTGGTGACAGTCTCGTTGTCCGCATCATTCCTGCATTGTTGGCAGGTGTGGTGATGAGCAGAATCATTACCCGTCTACCCGGATGGGATATTTCATGTCCTGAAACAGACTCTTTAAGTTAAGCGGGATACTTTATCTTTGGGCTACTCAAAAGCAGACAGGATGTTTCTATGACTCAAAATATCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAAc                       <  1:393384/251‑1 (MQ=255)                                                                                                                                                                                                                 cGTTGGCCGCATCATTCCTGCATGTTTGGCAGGTGTGGAGATGAGCAGAATCATTACCCGTCTACCCGGATTGGATATTTCATGTCCTGAAACAGACTCTTTAAGTTAAGCGGGATACTTTATCTTTGGGCTACTCAAAAGCAGACAGGATGTTTCTATGACTCAAAATATCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCg  <  2:66874/251‑1 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                  CGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCTGGTGGCATAGGTGGCGTCGTTATCACCGAGATCTATTGTGCAGTTATTGTACATATTAAGGGCATGGCTTCCTCGGAGTTGTTTGAAAACATTCTTTCTGGTGACAGTCTCGTTGTCCGCATCATTCCTGCATTGCTGGCAGGTGTGGTGATGAGCAGAATCATTACCCGTCTACCCGGATTGGATATTTCATGTCCTGAAACAGACTCTTTAAGTTAAGCGGGATACTTTATCTTTGGGCTACTCAAAAGCAGACAGGATGTTTCTATGACTCAAAATATCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCG  >  NC_000913/3897657‑3898115

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