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

CAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTA  >  NC_000913/1073596‑1073865                                                                                                                                         |                                                                                                                                      cagccagTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCa                                                                                                                                     <  2:397874/139‑1 (MQ=255)               gTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAg                                                                                                                       >  1:372468/1‑139 (MQ=255)                         aaaCGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACt                                                                                                             <  1:545978/139‑1 (MQ=255)                          aaCGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTc                                                                                                            >  1:1015906/1‑139 (MQ=255)                                cgccAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGtt                                                                                                      <  2:179478/139‑1 (MQ=255)                                 gccAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTtg                                                                                                     >  1:865374/1‑139 (MQ=255)                                           gAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGa                                                                                           >  2:603612/1‑139 (MQ=255)                                                gATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCa                                                                                      >  2:295030/1‑139 (MQ=255)                                                gATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCa                                                                                      >  2:1006457/1‑139 (MQ=255)                                                 aTGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAg                                                                                     >  1:1007192/1‑139 (MQ=255)                                                        cGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTgcc                                                                              >  1:560871/1‑139 (MQ=255)                                                              tACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCaaa                                                                        >  1:220952/1‑139 (MQ=255)                                                                   gCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGCATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCAc                                                                   <  2:865374/139‑1 (MQ=255)                                                                   gCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCAc                                                                   >  2:691919/1‑139 (MQ=255)                                                                    cgacgaTTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACg                                                                  <  1:242812/139‑1 (MQ=255)                                                                     gacgaTTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGc                                                                 <  1:691919/139‑1 (MQ=255)                                                                     gacgaTTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGc                                                                 <  2:560871/139‑1 (MQ=255)                                                                       cgaTTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAg                                                               >  1:62091/1‑139 (MQ=255)                                                                               ggAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGatca                                                       <  1:155596/139‑1 (MQ=255)                                                                                     aaCGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACa                                                 >  2:272976/1‑139 (MQ=255)                                                                                                 gcagcaGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGa                                                   >  2:180682/1‑125 (MQ=255)                                                                                                 gcagcaGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGa                                                   <  1:180682/125‑1 (MQ=255)                                                                                                             tAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCt                         >  2:1049738/1‑139 (MQ=255)                                                                                                                    gAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCttt                  >  1:93919/1‑139 (MQ=255)                                                                                                                         cgcgAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGca             >  2:554140/1‑139 (MQ=255)                                                                                                                           cgAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGaa                                                                                                                  >  1:126777/1‑36 (MQ=255)                                                                                                                           cgAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGaa                                                                                                                  <  2:126777/36‑1 (MQ=255)                                                                                                                           cgAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGcaca           >  2:465286/1‑139 (MQ=255)                                                                                                                            gAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACaa          >  1:738739/1‑139 (MQ=255)                                                                                                                            gAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACaa          <  1:603612/139‑1 (MQ=255)                                                                                                                            gAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACaa          <  2:1002227/139‑1 (MQ=255)                                                                                                                                  aCGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCtt    >  2:193620/1‑139 (MQ=255)                                                                                                                                    ggCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTa  <  1:733126/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                      CAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTA  >  NC_000913/1073596‑1073865

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