Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,079,433	T→C	V46A (GTG→GCG)	yiiE →	CopG family putative transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,079,433	0	T	C	100.0%	181.5 / NA	51	V46A (GTG→GCG)	yiiE	CopG family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/27);  total (24/27)

GGCGATGAATACGGTTTTTCTTCATTTATCAGAAGAGGCAATTAAACGGCTGAACAAGCTTCGGGGCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGTGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAATGACAAAAATTTTTCATTCTTGAATATAAAAAACAGATGCCCTTATTCTGGTATTAATACAAGGCTGTT  >  NC_000913/4079293‑4079575                                                                                                                                             |                                                                                                                                               ggCGATGAATACGGTTTTTCTTCATTTATCAGAAGAGGCAATTAAACGGCTGAACAAGCTTCGGGGCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCa                                                                                                                                      <  1:542895/151‑1 (MQ=255)   cGATGAATACGGTTTTTCTTCATTTATCAGAAGAGGCAATTAAACGGCTGAACAAGCTTCGGGGCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCaaa                                                                                                                                    <  1:995876/151‑1 (MQ=255)    gATGAATACGGTTTTTCTTCATTTATCAGAAGAGGCAATTAAACGGCTGAACAAGCTTCGGGGCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCaa                                                                                                                                     >  1:734783/1‑149 (MQ=255)       gAATACGGTTTTTCTTCATTTATCAGAAGAGGCAATTAAACGGCTGAACAAGCTTCGGGGCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAggg                                                                                                                                 >  2:487639/1‑150 (MQ=255)       gAATACGGATTTTCTTCATTTATCAGAAGAGGCAATTAAACGGCTGAACAAGCTTCGGGGCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAg                                                                                                                                   >  2:152033/1‑148 (MQ=255)           aCGGTTTTTCTTCATTTATCAGAAGAGGCAATTAAACGGCTGAACAAGCTTCGGGGCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGaggca                                                                                                                            <  1:827954/151‑1 (MQ=255)           aCGGTTTTTCTTCATTTATCAGAAGAGGCAATTAAACGGCTGAACAAGCTTCGGGGCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGaggc                                                                                                                             <  1:763009/150‑1 (MQ=255)                  ttcttcATTTATCAGAAGAGGCAATTAAACGGCTGAACAAGCTTCGGGGCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCaaa                                                                                                                      <  2:426383/150‑1 (MQ=255)                   tcttcATTTATCAGAAGAGGCAATTAAACGGCTGAACAAGCTTCGGGGCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCaaaa                                                                                                                     <  1:418120/150‑1 (MQ=255)                    cttcATTTATCAGAAGAGGCAATTAAACGGCTGAACAAGCTTCGGGGCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAgg                                                                                                                   >  1:400874/1‑151 (MQ=255)                       cATTTATCAGAAGAGGCAATTAAACGGCTGAACAAGCTTCGGGGCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAggggg                                                                                                                <  1:753660/151‑1 (MQ=255)                               agaagaGGCAATTAAACGGCTGAACAAGCTTCGGGGCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGtt                                                                                                        <  2:423157/151‑1 (MQ=255)                                 aagaGGCAATTAAACGGCTGAACAAGCTTCGGGGCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTgg                                                                                                      >  2:771349/1‑151 (MQ=255)                                  agagGCAATTAAACGGCTGAACAAGCTTCGGGGCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGCTAGTTGgt                                                                                                     >  1:209756/1‑151 (MQ=255)                                      gCAATTAAACGGCTGAACAAGCTTCGGGGCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTtg                                                                                                 >  2:1064825/1‑151 (MQ=255)                                            aaaCGGCTGAACAAGCTTCGGGGCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAgggg                                                                                                                 <  2:5674/129‑1 (MQ=255)                                            aaaCGGCTGAACAAGCTTCGGGGCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAgggg                                                                                                                 >  1:5674/1‑129 (MQ=255)                                                  cTGAACAAGCTTCGGGGCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAgtgt                                                                                      >  2:110785/1‑150 (MQ=255)                                                     aaCAAGCTTCGGGGCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAgg                                                                                  >  1:818510/1‑151 (MQ=255)                                                      aCAAGCTTCGGGGCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGa                                                                                 <  2:743315/151‑1 (MQ=255)                                                       cAAGCTTCGGGGCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAgg                                                                                  >  2:694344/1‑149 (MQ=255)                                                        aaGCTTCGGGGCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGc                                                                               >  1:437212/1‑151 (MQ=255)                                                             tCGGGGCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATaa                                                                           <  1:818408/150‑1 (MQ=255)                                                             tCGGGGCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATa                                                                            <  2:437186/149‑1 (MQ=255)                                                               ggggCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAgg                                                                         >  1:574409/1‑150 (MQ=255)                                                                 ggCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAAt                                                                      <  1:951700/151‑1 (MQ=255)                                                                 ggCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAAt                                                                      <  1:681049/151‑1 (MQ=255)                                                                  gCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAATg                                                                     <  1:694379/151‑1 (MQ=255)                                                                     ggCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAATGAc                                                                   <  2:553416/150‑1 (MQ=255)                                                                      gCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAATGACaa                                                                 <  2:762459/151‑1 (MQ=255)                                                                                tCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAATGACAAAAATTTTTc                                                        >  2:568054/1‑150 (MQ=255)                                                                                       cTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAATGACAAAAATTTTTCATTCTTg                                                 >  2:106104/1‑150 (MQ=255)                                                                                         gCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAATGACAAAAATTTTTCATTCtt                                                  <  2:186197/147‑1 (MQ=255)                                                                                         gCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAATGACAAAAATTTTTCATTCTTGa                                                <  1:109557/149‑1 (MQ=255)                                                                                            aTTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAATGACAAAAATTTTTCATTCTTGAatat                                            >  1:179049/1‑150 (MQ=255)                                                                                            aTTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAATGACAAAAATTTTTCATTCTTGAAtata                                           <  1:949995/151‑1 (MQ=255)                                                                                               ttACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAATGACAAAAATTTTTCATTCTTGAATATaaa                                         >  2:255014/1‑150 (MQ=255)                                                                                                 aCGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAATGACAAAAATTTTTCATTCTTGAATATaaaaa                                       <  2:962887/150‑1 (MQ=255)                                                                                                 aCGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAATGACAAAAATTTTTCATTCTTGAATATaaaaa                                       <  1:294964/150‑1 (MQ=255)                                                                                                   gcgAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAATGACAAAAATTTTTCATTCTTGAATATAAAAAAc                                     >  2:498383/1‑150 (MQ=255)                                                                                                     gAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGCGCGTCAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGAAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAATGACAAAAATTTTTCATTCTTGAATATAAAAAAc                                     >  2:671689/1‑148 (MQ=255)                                                                                                                  aaTATCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAATGACAAAAATTTTTCATTCTTGAATATAAAAAACAGATGCCCTTATTCTg                     >  1:318935/1‑151 (MQ=255)                                                                                                                   atatCTGGAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAATGACAAAAATTTTTCATTCTTGAATATAAAAAACAGATGCCCTTATTCTgg                    <  1:1003615/151‑1 (MQ=255)                                                                                                                         ggAACGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAATGACAAAAATTTTTCATTCTTGAATATAAAAAACAGATGCCCTTATTCTGGTAtt                <  2:1013887/149‑1 (MQ=255)                                                                                                                            aCGTCAGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAATGACAAAAATTTTTCATTCTTGAATATAAAAAACAGATGCCCTTATTCTGGTATTAATAc           >  1:85386/1‑151 (MQ=255)                                                                                                                             cGTCAGAAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAATGACAAAAATTTTTCATTCTTGAATATAAAAAACAGATGCCCTTATTCTGTTATTAATAc           >  2:331222/1‑150 (MQ=255)                                                                                                                                 aGCAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAATGACAAAAATTTTTCATTCTTGAATATAAAAAACAGATGCCCTTATTCTGGTATTAATACAAGGc      <  2:897401/151‑1 (MQ=255)                                                                                                                                   cAATTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAATGACAAAAATTTTTCATTCTTGAATATAAAAAACAGATGCCCTTATTCTGGTATTAATACAAGGCTg    <  2:286365/151‑1 (MQ=255)                                                                                                                                    aaTTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAATGACAAAAATTTTTCATTCTTGAATATAAAAAACAGATGCCCTTATTCTGGTATTAATACAAGGCTGt   >  1:668858/1‑151 (MQ=255)                                                                                                                                     aTTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAATGACAAAAATTTTTCATTCTTGAATATAAAAAACAGATGCCCTTATTCTGGTATTAATACAAGGCTGtt  >  1:547088/1‑151 (MQ=255)                                                                                                                                     aTTCCCGGCGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAATGACAAAAATTTTTCATTCTTGAATATAAAAAACAGATGCCCTTATTCTGGTATTAATACAAGGCTGtt  <  1:885393/151‑1 (MQ=255)                                                                                                                                             |                                                                                                                                               GGCGATGAATACGGTTTTTCTTCATTTATCAGAAGAGGCAATTAAACGGCTGAACAAGCTTCGGGGCTGGCGTAAGGTTTCACGATCTGCAATTTTACGCGAAGCGGTAGAGCAATATCTGGAACGTCAGCAATTCCCGGTGCGTAAAGCAAAGGGAGGCAGGCAAAAGGGGGAGGTAGTTGGTGTTGATGATCAGTGTAAGGAGCATAAGGAATGACAAAAATTTTTCATTCTTGAATATAAAAAACAGATGCCCTTATTCTGGTATTAATACAAGGCTGTT  >  NC_000913/4079293‑4079575

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