Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,071,788	C→T	A124A (GCG→GCA)	rutC ←	putative aminoacrylate peracid reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,071,788	0	C	T	100.0%	169.4 / NA	48	A124A (GCG→GCA)	rutC	putative aminoacrylate peracid reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/20);  total (28/20)

CTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTCGGCGTAGATTTCGTTAATCGCGGCGTAATTTTTCCAGTCGGTAATA  >  NC_000913/1071652‑1071923                                                                                                                                         |                                                                                                                                        cTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCa                                                                                                                                      >  2:837818/1‑140 (MQ=255) cTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCa                                                                                                                                      >  1:488431/1‑140 (MQ=255)    cTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAAAACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAAtt                                                                                                                                   >  2:371886/1‑140 (MQ=255)        aGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAACCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTgg                                                                                                                               <  1:315801/140‑1 (MQ=255)          cACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCg                                                                                                                             <  1:661488/140‑1 (MQ=255)           aCCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGa                                                                                                                            >  1:926295/1‑140 (MQ=255)               ccAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGAGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTc                                                                                                                        >  2:959550/1‑140 (MQ=255)                        gTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAgcgc                                                                                                               >  2:957617/1‑140 (MQ=255)                         tAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAg                                                                                                                  <  1:887707/136‑1 (MQ=255)                         tAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAg                                                                                                                  >  2:887707/1‑136 (MQ=255)                         tAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGcgcg                                                                                                              >  1:916536/1‑140 (MQ=255)                          aaCCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGt                                                                                                             <  2:519969/140‑1 (MQ=255)                          aaCCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGt                                                                                                             <  1:837818/140‑1 (MQ=255)                              aGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAgg                                                                                                         <  2:659473/140‑1 (MQ=255)                                 aaCTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGttt                                                                                                      >  2:572699/1‑140 (MQ=255)                                     gCCGCTACCCCCAAGACCCGAAATCAAAACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACAAGCGCGTCAGGTTTTAcc                                                                                                  <  1:371886/140‑1 (MQ=255)                                       cGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAg                                                                                                                  >  1:386907/1‑122 (MQ=255)                                       cGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAg                                                                                                                  <  2:386907/122‑1 (MQ=255)                                          tACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTAcc                                                                                                  >  2:593391/1‑135 (MQ=255)                                          tACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTAcc                                                                                                  <  1:593391/135‑1 (MQ=255)                                               ccAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGGGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACt                                                                                        <  2:916536/140‑1 (MQ=255)                                               ccAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACt                                                                                        <  1:957617/140‑1 (MQ=255)                                               ccAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACt                                                                                        <  2:204814/140‑1 (MQ=255)                                                 aaGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACt                                                                                        <  1:28066/138‑1 (MQ=255)                                                 aaGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACt                                                                                        >  2:28066/1‑138 (MQ=255)                                                              aaCACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATcgcg                                                                         <  1:654031/140‑1 (MQ=255)                                                                      tACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATcgc                                                                          >  2:369439/1‑131 (MQ=255)                                                                      tACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATcgc                                                                          <  1:369439/131‑1 (MQ=255)                                                                      tACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTa                                                                 >  1:895467/1‑140 (MQ=255)                                                                           gcgcATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCAcc                                                            >  2:347677/1‑140 (MQ=255)                                                                             gcATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCg                                                          >  1:585399/1‑140 (MQ=255)                                                                              cATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCgg                                                         >  1:926743/1‑140 (MQ=255)                                                                               aTCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGa                                                        >  2:85026/1‑140 (MQ=255)                                                                                      tAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACt                                                 <  2:926295/140‑1 (MQ=255)                                                                                         gggggAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTCgg                                              >  1:339721/1‑140 (MQ=255)                                                                                           gggAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTCGGCg                                            <  2:926743/140‑1 (MQ=255)                                                                                              aGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTCGGCGTAg                                         >  2:329806/1‑140 (MQ=255)                                                                                              aGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTCGGCGTAg                                         >  1:434956/1‑140 (MQ=255)                                                                                              aGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTCGGCGTAg                                         >  2:315015/1‑140 (MQ=255)                                                                                                   agagTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTCGGCGTAGATTTc                                    >  1:650397/1‑140 (MQ=255)                                                                                                   agagTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTCGGCGTAGATTTc                                    >  1:392091/1‑140 (MQ=255)                                                                                                          aaGTTTCATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTCGGCGTAGATTTCGTTAATc                             >  2:630809/1‑140 (MQ=255)                                                                                                                cATCGCGGCCTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTCGGCGTAGATTTCGTTAATCgcggcg                       >  2:709129/1‑140 (MQ=255)                                                                                                                         cTCACTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTCGGCGTAGATTTCGTTAATCGCGGCGTAATTTTTc              <  1:512901/140‑1 (MQ=255)                                                                                                                             cTTGGCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTCGGCGTAGATTTCGTTAATCGCGGCGTAATTTTTCCAGt          <  2:632270/140‑1 (MQ=255)                                                                                                                                 gCGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTCGGCGTAGATTTCGTTAATCGCGGCGTAATTTTTCCAGTCGGt      <  1:709129/140‑1 (MQ=255)                                                                                                                                  cGATATGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTCGGCGTAGATTTCGTTAATCGCGGCGTAATTTTTCCAGTCGGta     >  2:506181/1‑140 (MQ=255)                                                                                                                                     tatGTGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTCGGCGTAGATTTCGTTAATCGCGACGTAATTTTTCCAGTCGGtaata  >  1:676612/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                        CTGCTCCAGCACCGCCAGTTGCGGTAACCAGTAACTGCCGCTACCCCCAAGACCCGAAATCAACACCACTACGGGCGCATCAGCATAAGGGGGAGGTGAGAGTGAAAGTTTCATCGCGGCCTCACTTGGCGATATGCGCAATTGTGGCGATTTCCACCAGCGCGTCAGGTTTTACCAGTCCGCACTGAATGCAGAATCGCGCCGGTTTATCACCCGGAAAAAACTCGGCGTAGATTTCGTTAATCGCGGCGTAATTTTTCCAGTCGGTAATA  >  NC_000913/1071652‑1071923

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑