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

TCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTT  >  NC_000913/1073603‑1073864                                                                                                                                  |                                                                                                                                     tCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCg                                                                                                                              >  1:300177/1‑138 (MQ=255) tCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCg                                                                                                                              <  2:300177/138‑1 (MQ=255) tCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGc                                                                                                                             >  2:38658/1‑139 (MQ=255)   aCGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCa                                                                                                                           <  2:916457/139‑1 (MQ=255)      aGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCa                                                                                                                        <  2:429931/139‑1 (MQ=255)             cGCCAAAACTCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGaa                                                                                                                 <  2:261918/139‑1 (MQ=255)                  aaaCGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACt                                                                                                            <  2:396087/139‑1 (MQ=255)                           ccAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTtgt                                                                                                   >  2:230685/1‑139 (MQ=255)                             agagATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGa                                                                                                 <  1:1116893/139‑1 (MQ=255)                                    gAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGa                                                                                          >  1:24597/1‑139 (MQ=255)                                           tGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTtgt                                                                                                   >  2:637734/1‑123 (MQ=255)                                           tGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTtgt                                                                                                   <  1:637734/123‑1 (MQ=255)                                                 cGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTgcc                                                                             >  2:995077/1‑139 (MQ=255)                                                       tACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACt                                                                                                            >  1:87044/1‑102 (MQ=255)                                                       tACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACt                                                                                                            <  2:87044/102‑1 (MQ=255)                                                          gggCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGcc                                                                    <  1:597971/139‑1 (MQ=255)                                                            gCGACGATGGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCAc                                                                  <  2:328602/139‑1 (MQ=255)                                                                  aTTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGtt                                                            <  1:230685/139‑1 (MQ=255)                                                                  aTTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGtt                                                            <  1:443039/139‑1 (MQ=255)                                                                  aTTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGtt                                                            <  1:476333/139‑1 (MQ=255)                                                                                   tAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGcg                                           >  1:768495/1‑139 (MQ=255)                                                                                    aaGGTGGGAGCAGTGGCGTAAATCTGACGGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCCGTTTTGCCGCCAAAGCCTCGCAGTTTGATCATCGACAGGGcga                                          <  2:751214/139‑1 (MQ=255)                                                                                    aaGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGcga                                          >  1:1132544/1‑139 (MQ=255)                                                                                      ggTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATcc                                                                                              >  2:1292769/1‑85 (MQ=255)                                                                                      ggTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATcc                                                                                              <  1:1292769/85‑1 (MQ=255)                                                                                       gTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGtgaaat                                       <  2:24597/139‑1 (MQ=255)                                                                                       gTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGcgaaat                                       <  1:38658/139‑1 (MQ=255)                                                                                       gTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGcgaaat                                       >  2:629131/1‑139 (MQ=255)                                                                                           cagcagTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGc                                                                <  2:735417/110‑1 (MQ=255)                                                                                           cagcagTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGc                                                                >  1:735417/1‑110 (MQ=255)                                                                                                           cTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCa                                                                                                                           <  1:900686/35‑1 (MQ=255)                                                                                                           cTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCa                                                                                                                           >  2:900686/1‑35 (MQ=255)                                                                                                                   gcgAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGAc                                                                                                             >  1:901686/1‑41 (MQ=255)                                                                                                                   gcgAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGAc                                                                                                             <  2:901686/41‑1 (MQ=255)                                                                                                                     gAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACaa         <  1:1192447/139‑1 (MQ=255)                                                                                                                     gAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACaa         <  1:629131/139‑1 (MQ=255)                                                                                                                     gAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACaa         <  1:111137/139‑1 (MQ=255)                                                                                                                            cGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCttt  >  2:117051/1‑139 (MQ=255)                                                                                                                                  |                                                                                                                                     TCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTT  >  NC_000913/1073603‑1073864

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