Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,003,303	A→G	R292R (CGT→CGC)	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,303	0	A	G	100.0%	119.4 / NA	33	R292R (CGT→CGC)	rarD	putative chloramphenical resistance permease
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/13);  total (20/13)

TACCGTTAGATGTTGGGATAAATATACGTAACATAAATTTTACATCCTTGTATGAGTCTCCGGTCAGCATGGCAATATGCCCACTCTCATGCAAGAGTCGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAA  >  NC_000913/4003160‑4003447                                                                                                                                                |                                                                                                                                                 tACCGTTAGATGTTGGGATAAATATACGTAACATAAATTTTACATCCTTGTATGAGTCTCCGGTCAGCATGGCAATATGCCCACTCTCATGCAAGAGTCGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTGCGCt                                                                                                                                              <  1:80396/148‑1 (MQ=255)            gTTGGGATAAATATACGTAACATAAATTTTACATCCTTGTATGAGTCTCCGGTCAGCATGGCAATATGCCCACTCTCATGCAAGAGTCGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTGCGCTGCGTATAAATCGca                                                                                                                                <  2:466941/151‑1 (MQ=255)                    aaaTATACGTAACATAAATTTTACATCCTTGTATGAGTCTCCGGTCAGCATGGCAATATGCCCACTCTCATGCAAGAGTCGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTGCGCTGCGTATAAATCGCATCCATcac                                                                                                                        >  2:140978/1‑151 (MQ=255)                        ataCGTAACATAAATTTTACATCCTTGTATGAGTCTCCGGTCAGCATGGCAATATGCCCACTCTCATGCAAGAGTCGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTGCGCTGCGTATAAATCGCATCCATCACaaa                                                                                                                     >  1:356262/1‑150 (MQ=255)                           cGTAACATAAATTTTACATCCTTGTATGAGTCTCCGGTCAGCATGGCAATATGCCCACTCTCATGCAAGAGTCGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTGCGCTGCGTATAAATCGCATCCATCACAAAAAtt                                                                                                                 >  1:582374/1‑151 (MQ=255)                              aaCATAAATTTTACATCCTTGTATGAGTCTCCGGTCAGCATGGCAATATGCCCCCTCTCATGCAAGAGTCGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTGCGCTGCGTATAAATCGCATCCATCTCAAAAAtt                                                                                                                 <  2:82918/148‑1 (MQ=255)                              aaCATAAATTTTACATCCTTGTATGAGTCTCCGGTCAGCATGGCAATATGCCCACTCTCATGCAAGAGTCGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTGCGCTGCGTATAAATCGCATCCATCACAAAAATTGcc                                                                                                              >  1:445410/1‑151 (MQ=255)                                     aTTTTACATCCTTGTATGAGTCTCCGGTCAGCATGGCAATATGCCCACTCTCATGCAAGAGTCGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTGCGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGcc                                                                                                        >  2:50710/1‑150 (MQ=255)                                       tttACATCCTTGTATGAGTCTCCGGTCAGCATGGCAATATGCCCACTCTCATGCAAGAGTCGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTGCGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCAcc                                                                                                     >  1:288521/1‑151 (MQ=255)                                          aCATCCTTGTATGAGTCTCCGGTCAGCATGGCAATATGCCCACTCTCATGCAAGAGTCGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTGCGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCaa                                                                                                  <  1:126931/151‑1 (MQ=255)                                                                          aaTATGCCCACTCTCATGCAAGAGTCGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTGCGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCAcc                                                                  >  1:454049/1‑151 (MQ=255)                                                                          aaTATGCCCACTCTCATGCAAGAGTCGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTGCGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCAcc                                                                  <  1:140981/151‑1 (MQ=255)                                                                                cccACTCTCATGCAAGAGTCGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTGCGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGttt                                                            >  1:530353/1‑151 (MQ=255)                                                                                   aCTCTCATGCAAGAGTCGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTGCGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTc                                                         <  1:225824/151‑1 (MQ=255)                                                                                        cATGCAAGAGTCGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTGCGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCAt                                                    >  1:464955/1‑151 (MQ=255)                                                                                           gCAAGAGTCGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTGCGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATaaa                                                 <  2:399162/151‑1 (MQ=255)                                                                                                agTCGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTGCGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTc                                            <  2:464947/151‑1 (MQ=255)                                                                                                         atTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTGCGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCAcagcc                                      >  1:114836/1‑148 (MQ=255)                                                                                                           tttttCAGAATATATTTATTTTTTATTTGGACGTTCTGCGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCAcagccagc                                   >  1:63980/1‑149 (MQ=255)                                                                                                            ttttCAGAATATATTTATTTTTTATTTGGACGTTCTGCGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAgg                                >  2:296925/1‑151 (MQ=255)                                                                                                                cAGAATATATTTATTTTTTATTTGGACGTTCTGCGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACa                            <  2:253832/151‑1 (MQ=255)                                                                                                                      atatTTATTTTTTATTTGGACGTTCTGCGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCg                      >  2:406960/1‑151 (MQ=255)                                                                                                                      atatTTATTTTTTATTTGGACGTTCTGCGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCg                      >  2:597662/1‑151 (MQ=255)                                                                                                                      atatTTATTTTTTATTTGGACGTTCTGCGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCg                      >  1:440713/1‑151 (MQ=255)                                                                                                                            attttttatttGGACGTTCTGCGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCggg                 <  2:107659/150‑1 (MQ=255)                                                                                                                            attttttatttGGACGTTCTGCGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCggg                 <  1:583332/150‑1 (MQ=255)                                                                                                                              tttttatttGGACGTTCTGCGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGcc               >  1:176989/1‑150 (MQ=255)                                                                                                                                  tatttGGACGTTCTGCGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGt          >  1:611633/1‑151 (MQ=255)                                                                                                                                   atttGGACGTTCTGCGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGGACATCAGCGTCGGGCCAATGTa         <  2:36628/151‑1 (MQ=255)                                                                                                                                   atttGGACGTTCTGCGCTCCGGATAAATCGCATCCATCTTAAAAATCGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGTCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGt          >  2:355134/1‑150 (MQ=255)                                                                                                                                     ttGGACGTTCTGCGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACt       >  2:593652/1‑151 (MQ=255)                                                                                                                                      tGGACGTTCTGCGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACt       <  1:314876/150‑1 (MQ=255)                                                                                                                                           gTTCTGCGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGaaa  >  2:126788/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                 TACCGTTAGATGTTGGGATAAATATACGTAACATAAATTTTACATCCTTGTATGAGTCTCCGGTCAGCATGGCAATATGCCCACTCTCATGCAAGAGTCGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAA  >  NC_000913/4003160‑4003447

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