Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,003,844	A→G	V112A (GTG→GCG)	rarD ←	putative chloramphenical resistance permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,003,844	0	A	G	100.0%	151.9 / NA	44	V112A (GTG→GCG)	rarD	putative chloramphenical resistance permease
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/25);  total (19/25)

AAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCG  >  NC_000913/4003696‑4003985                                                                                                                                                     |                                                                                                                                              aaaaCTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGcc                                                                                                                                             <  2:627917/151‑1 (MQ=255)  aaaCTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCa                                                                                                                                            <  1:262202/151‑1 (MQ=255)    aCTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGCTCGCCAGc                                                                                                                                          >  1:83232/1‑151 (MQ=255)       aaaTGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCggg                                                                                                                                       <  1:326374/151‑1 (MQ=255)       aaaTGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTAGCCAg                                                                                                                                           >  2:607131/1‑147 (MQ=255)        aaTGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGt                                                                                                                                      >  2:506170/1‑151 (MQ=255)           gCCAGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTta                                                                                                                                    >  2:645432/1‑150 (MQ=255)           gCCAGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTta                                                                                                                                    >  2:536948/1‑150 (MQ=255)            ccAGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTtaat                                                                                                                                  >  2:607403/1‑151 (MQ=255)             cAGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTtaata                                                                                                                                 <  2:459400/151‑1 (MQ=255)              aGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTtaataa                                                                                                                                <  2:456474/151‑1 (MQ=255)                  ccAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGt                                                                                                                             >  2:283644/1‑150 (MQ=255)                          ataataGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCaa                                                                                                                       >  2:433078/1‑148 (MQ=255)                                  tAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTcc                                                                                                            >  1:95725/1‑151 (MQ=255)                                   aGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCa                                                                                                           <  2:239601/151‑1 (MQ=255)                                       aaCCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCat                                                                                                       <  1:105175/151‑1 (MQ=255)                                          cAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCatat                                                                                                     >  1:202511/1‑150 (MQ=255)                                          cAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATg                                                                                                    >  2:428725/1‑151 (MQ=255)                                             aaGCCCACAGCGGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTg                                                                                                 <  1:480291/151‑1 (MQ=255)                                                    caGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTc                                                                                          >  1:350209/1‑151 (MQ=255)                                                     aGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCa                                                                                         <  2:732299/151‑1 (MQ=255)                                                         ggACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGc                                                                                     >  2:455348/1‑151 (MQ=255)                                                           aCTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGccc                                                                                   <  2:721772/151‑1 (MQ=255)                                                             tAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGccc                                                                                   >  1:663671/1‑149 (MQ=255)                                                              aaCACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCa                                                                                  <  2:258825/149‑1 (MQ=255)                                                              aaCACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAg                                                                                 <  1:428674/150‑1 (MQ=255)                                                                 acGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATaa                                                                             <  2:12368/151‑1 (MQ=255)                                                                        atatCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAg                                                                      <  2:11106/151‑1 (MQ=255)                                                                              ccAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGtt                                                                >  2:34260/1‑151 (MQ=255)                                                                                          ccagccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAg                                                    <  2:143176/151‑1 (MQ=255)                                                                                           cagccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGc                                                   >  2:117890/1‑151 (MQ=255)                                                                                                aTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCAc                                                 <  1:450997/148‑1 (MQ=255)                                                                                                  tGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCa                                             <  1:243071/150‑1 (MQ=255)                                                                                                       ccGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACt                                       <  2:350243/151‑1 (MQ=255)                                                                                                                     cgcCGAGGAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAAt                         >  1:487370/1‑151 (MQ=255)                                                                                                                           ggAAAATCATCCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCtg                   <  1:653932/151‑1 (MQ=255)                                                                                                                             aaaaTCATCCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCtgtg                 <  2:619466/151‑1 (MQ=255)                                                                                                                               aaTCATCCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGc               <  2:692999/151‑1 (MQ=255)                                                                                                                                atcatcCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCg              <  2:706162/151‑1 (MQ=255)                                                                                                                                  catcCCCAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGt             <  2:188994/150‑1 (MQ=255)                                                                                                                                     ccccAGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTg          >  1:749603/1‑150 (MQ=255)                                                                                                                                         aGCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCa     <  2:303495/151‑1 (MQ=255)                                                                                                                                          gCACAATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGTCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAg    <  2:563550/151‑1 (MQ=255)                                                                                                                                            acaATGTTCGCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCg  >  1:3527/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                              AAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCG  >  NC_000913/4003696‑4003985

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