Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,073,732	C→T	A94T (GCG→ACG)	rutA ←	pyrimidine oxygenase, FMN‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,073,732	0	C	T	100.0%	185.5 / NA	58	A94T (GCG→ACG)	rutA	pyrimidine oxygenase, FMN‑dependent
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/34);  total (24/34)

GCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAA  >  NC_000913/1073598‑1073857                                                                                                                                       |                                                                                                                              gccagTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGc                                                                                                                           >  1:1377527/1‑139 (MQ=255)      tCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGc                                                                                                                      >  1:713908/1‑139 (MQ=255)                  cgacAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGaa                                                                                                          <  1:774118/136‑1 (MQ=255)                  cGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTATCGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGaa                                                                                                          <  2:1258196/139‑1 (MQ=255)                  cGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGaa                                                                                                          <  2:713908/139‑1 (MQ=255)                   gCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAAc                                                                                                         <  2:778818/139‑1 (MQ=255)                    ccAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACg                                                                                                        >  2:313835/1‑139 (MQ=255)                      aaaaCGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGAc                                                                                                      >  1:8192/1‑139 (MQ=255)                      aaaaCGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGAc                                                                                                      <  2:1097560/139‑1 (MQ=255)                        aaCGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCa                                                                                                                 <  1:1204132/126‑1 (MQ=255)                        aaCGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCa                                                                                                                 >  2:1204132/1‑126 (MQ=255)                                     gaTGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATcc                                                                                       <  1:1563008/139‑1 (MQ=255)                                     gaTGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATcc                                                                                       <  1:1321591/139‑1 (MQ=255)                                       tGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCa                                                                                     >  1:431972/1‑139 (MQ=255)                                          aGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGaa                                                                                  >  1:905533/1‑139 (MQ=255)                                           gTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAAc                                                                                 <  1:1617913/139‑1 (MQ=255)                                              gATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCTATCAAGGTGaa                                                                                                          >  2:1499143/1‑111 (MQ=255)                                              gATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCTATCAAGGTGaa                                                                                                          <  1:1499143/111‑1 (MQ=255)                                              gATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCa                                                                                                                    <  2:495801/101‑1 (MQ=255)                                              gATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCa                                                                                                                    >  1:495801/1‑101 (MQ=255)                                                tGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGCTAAGGTGGCAGCAGTGGCGTAAATCTGCATGCGCGAGGTAACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGt                                                                            >  2:967563/1‑139 (MQ=255)                                                      cGGCCATACGGGCGACGATTGCCGGAGGTAACGTTGAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGaa                                                                                                          <  2:967889/103‑1 (MQ=255)                                                      cGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGaa                                                                                                          >  1:967889/1‑103 (MQ=255)                                                            tACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCaaa                                                                >  2:589666/1‑139 (MQ=255)                                                               gggCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGcc                                                             <  1:589666/139‑1 (MQ=255)                                                                ggCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCa                                                            <  2:1307731/139‑1 (MQ=255)                                                                      gaTTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGt                                                      >  1:427438/1‑139 (MQ=255)                                                                       aTTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGTACTCAGTTTTGCCGCCAAAGCCACGCAGtt                                                     <  1:1371473/139‑1 (MQ=255)                                                                       aTTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGtt                                                     <  2:375553/139‑1 (MQ=255)                                                                       aTTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGtt                                                     <  2:1628441/139‑1 (MQ=255)                                                                       aTTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGtt                                                     <  2:1387863/139‑1 (MQ=255)                                                                       aTTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGtt                                                     <  1:313835/139‑1 (MQ=255)                                                                            cGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGatc                                                <  2:1377527/139‑1 (MQ=255)                                                                                    aCGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGcc                                                                                                                     <  1:399393/62‑1 (MQ=255)                                                                                    aCGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGcc                                                                                                                     >  2:399393/1‑62 (MQ=255)                                                                                        tAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGcg                                    <  1:1065486/139‑1 (MQ=255)                                                                                         aaGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAG‑CCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGcgaa                                  <  2:1416238/139‑1 (MQ=255)                                                                                            gTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGcgaaat                                <  1:195520/139‑1 (MQ=255)                                                                                            gTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGcgaaat                                <  2:90632/139‑1 (MQ=255)                                                                                               gcagcaGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGcgaaatcga                             >  2:1458963/1‑139 (MQ=255)                                                                                                      tGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACt                                                                                                     >  1:849740/1‑60 (MQ=255)                                                                                                      tGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACt                                                                                                     <  2:849740/60‑1 (MQ=255)                                                                                                      tGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAg                                                                                                 >  2:486625/1‑64 (MQ=255)                                                                                                      tGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAg                                                                                                 <  1:486625/64‑1 (MQ=255)                                                                                                      tGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAggg                                      >  2:780732/1‑123 (MQ=255)                                                                                                      tGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTa                      >  1:1362231/1‑139 (MQ=255)                                                                                                      tGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTATGATCATCGACAggg                                      <  1:780732/123‑1 (MQ=255)                                                                                                        gCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGc                                                                                                                           <  2:1085885/36‑1 (MQ=255)                                                                                                        gCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGc                                                                                                                           >  1:1085885/1‑36 (MQ=37)                                                                                                           tAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCt                 <  1:1458963/139‑1 (MQ=255)                                                                                                                cTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCt            <  2:512593/139‑1 (MQ=255)                                                                                                                   aaTGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCtttt         >  2:1084937/1‑139 (MQ=255)                                                                                                                         cgAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGcaca   <  1:1438646/139‑1 (MQ=255)                                                                                                                         cgAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGcaca   >  2:1075375/1‑139 (MQ=255)                                                                                                                         cgAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGcaca   >  1:108146/1‑139 (MQ=255)                                                                                                                          gAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACaa  <  1:43620/139‑1 (MQ=255)                                                                                                                             gTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACt                                                                                                     >  2:318686/1‑37 (MQ=255)                                                                                                                             gTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACt                                                                                                     <  1:318686/37‑1 (MQ=255)                                                                                                                                       |                                                                                                                              GCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAA  >  NC_000913/1073598‑1073857

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