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%	198.1 / NA	61	A94T (GCG→ACG)	rutA	pyrimidine oxygenase, FMN‑dependent
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/29);  total (32/29)

AGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCAC  >  NC_000913/1073597‑1073855                                                                                                                                        |                                                                                                                            agccagTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAg                                                                                                                          <  1:720703/139‑1 (MQ=255)       tCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGc                                                                                                                    >  2:1496469/1‑139 (MQ=255)                     ccAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGa                                                                                                         >  2:326204/1‑136 (MQ=255)                     ccAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGa                                                                                                         <  1:326204/136‑1 (MQ=255)                       aaaaCGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGAc                                                                                                    >  2:602587/1‑139 (MQ=255)                                                 tGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGt                                                                          >  1:312603/1‑139 (MQ=255)                                                  ggTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGtt                                                                         >  2:1490820/1‑139 (MQ=255)                                                             tACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCaaa                                                              >  2:752798/1‑139 (MQ=255)                                                              aCGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAg                                                             >  2:1216909/1‑139 (MQ=255)                                                                   cgacgaTTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACg                                                        <  2:676843/139‑1 (MQ=255)                                                                    gacgaTTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGc                                                       >  2:1061236/1‑139 (MQ=255)                                                                        aTTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGaa                                                                                                        <  1:278765/86‑1 (MQ=255)                                                                        aTTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGaa                                                                                                        >  2:278765/1‑86 (MQ=255)                                                                        aTTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGtt                                                   <  1:1496469/139‑1 (MQ=255)                                                                        aTTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGtt                                                   <  2:603452/139‑1 (MQ=255)                                                                        aTTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGtt                                                   <  2:776542/139‑1 (MQ=255)                                                                        aTTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGtt                                                   <  1:1386467/139‑1 (MQ=255)                                                                        aTTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGtt                                                   <  1:1216909/139‑1 (MQ=255)                                                                                   tAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGAc                                        <  1:752798/139‑1 (MQ=255)                                                                                      cGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACt                                                                                                   <  1:770900/77‑1 (MQ=255)                                                                                      cGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACt                                                                                                   >  2:770900/1‑77 (MQ=255)                                                                                         tAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCa                                                                >  2:704962/1‑109 (MQ=255)                                                                                         tAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCa                                                                <  1:704962/109‑1 (MQ=255)                                                                                            ggTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAg                                                                                                             <  1:587071/61‑1 (MQ=255)                                                                                            ggTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAg                                                                                                             >  2:587071/1‑61 (MQ=255)                                                                                                gcagcaGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGcgaaatcga                           <  1:1490820/139‑1 (MQ=255)                                                                                                 cagcagTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGTTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGcgaaatcgaa                          >  1:1404434/1‑139 (MQ=255)                                                                                                 cagcagTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCa                                                                                                 <  1:1591705/68‑1 (MQ=255)                                                                                                 cagcagTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCa                                                                                                 >  2:1591705/1‑68 (MQ=255)                                                                                                 cagcagTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGcgaaatcgaa                          >  1:667244/1‑139 (MQ=255)                                                                                                 cagcagTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGcgaaatcgaa                          >  2:355082/1‑139 (MQ=255)                                                                                                 cagcagTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGcgaaatcgaa                          >  1:447228/1‑139 (MQ=255)                                                                                                  agcagTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGccc                                                                                                                       >  2:903251/1‑45 (MQ=255)                                                                                                  agcagTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGccc                                                                                                                       <  1:903251/45‑1 (MQ=255)                                                                                                       tGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCa                                                          <  1:335029/101‑1 (MQ=255)                                                                                                       tGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCa                                                          >  2:335029/1‑101 (MQ=255)                                                                                                          cGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTg                 >  1:514615/1‑139 (MQ=255)                                                                                                                   gAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGtt                                                   >  2:1029592/1‑96 (MQ=255)                                                                                                                   gAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGtt                                                   <  1:1029592/96‑1 (MQ=255)                                                                                                                   gAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCttt        <  1:1061236/139‑1 (MQ=255)                                                                                                                   gAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCttt        >  1:172160/1‑139 (MQ=255)                                                                                                                    aaTGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGa                                                                                        >  1:729836/1‑58 (MQ=255)                                                                                                                    aaTGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGa                                                                                        <  2:729836/58‑1 (MQ=255)                                                                                                                    aaTGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCAc                                                         <  2:1659668/89‑1 (MQ=255)                                                                                                                    aaTGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCAc                                                         >  1:1659668/1‑89 (MQ=255)                                                                                                                    aaTGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCtttt       >  1:190769/1‑139 (MQ=255)                                                                                                                     aTGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGTTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCC‑CTTTTTg     <  2:1404434/139‑1 (MQ=255)                                                                                                                     aTGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCttt        <  1:531078/137‑1 (MQ=255)                                                                                                                     aTGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCttt        >  2:531078/1‑137 (MQ=255)                                                                                                                      tGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGcgaaatcga                           <  2:1603000/117‑1 (MQ=255)                                                                                                                      tGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGcgaaatcga                           >  1:1603000/1‑117 (MQ=255)                                                                                                                         gcgAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGcac  <  2:447228/139‑1 (MQ=255)                                                                                                                         gcgAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGcac  >  2:1369948/1‑139 (MQ=255)                                                                                                                          cgAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGaa                                                                                                        >  1:820726/1‑36 (MQ=255)                                                                                                                          cgAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGaa                                                                                                        <  2:820726/36‑1 (MQ=255)                                                                                                                          cgAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTgccgc                                                                  >  1:713624/1‑74 (MQ=255)                                                                                                                          cgAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTgccgc                                                                  <  2:713624/74‑1 (MQ=255)                                                                                                                                   ggCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTtgt                                                                                          <  1:1364218/41‑1 (MQ=255)                                                                                                                                   ggCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTtgt                                                                                          >  2:1364218/1‑41 (MQ=255)                                                                                                                                   ggCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGcgaaat                              >  2:1090179/1‑101 (MQ=255)                                                                                                                                   ggCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGcgaaat                              <  1:1090179/101‑1 (MQ=255)                                                                                                                                        |                                                                                                                            AGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCAC  >  NC_000913/1073597‑1073855

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