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_0009131,073,7320CT100.0% 62.3 / NA 21A94T (GCG→ACG) rutApyrimidine oxygenase, FMN‑dependent
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (9/12);  total (9/12)

GCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGG  >  NC_000913/1073598‑1073860
                                                                                                                                      |                                                                                                                                
gccagTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGc                                                                                                                              >  1:330511/1‑139 (MQ=255)
                 cGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGaa                                                                                                             <  2:513566/139‑1 (MQ=255)
                 cGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGaa                                                                                                             <  1:431108/139‑1 (MQ=255)
                 cGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGGGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCTCCATCAAGGTGaa                                                                                                             <  1:389855/139‑1 (MQ=255)
                     aaaaCGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAAGGTTAAGGTGGCAGCAGTGGCGTAAATCTGAAAGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGAc                                                                                                         >  1:523051/1‑139 (MQ=255)
                      aaaCGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACt                                                                                                        <  1:200868/139‑1 (MQ=255)
                        aCGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCa                                                                                                      <  1:741273/139‑1 (MQ=255)
                                                   tGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGCCACGGCCGTCAGTCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTg                                                                           >  1:627670/1‑139 (MQ=255)
                                                            aCGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCaaa                                                                   <  2:78567/138‑1 (MQ=255)
                                                            aCGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCaaa                                                                   >  1:78567/1‑138 (MQ=255)
                                                                      aTTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGtt                                                        <  2:20708/139‑1 (MQ=255)
                                                                      aTTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGtt                                                        >  1:162149/1‑139 (MQ=255)
                                                                      aTTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGCCACGGCCGTCAGTCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGtt                                                        <  2:627670/139‑1 (MQ=255)
                                                                               gggcaCGGTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCg                                               >  1:13935/5‑139 (MQ=255)
                                                                                        aaGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGcga                                      >  2:249020/1‑139 (MQ=255)
                                                                                                agcagTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGaa                                                                                                             >  1:489970/1‑60 (MQ=255)
                                                                                                agcagTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCAGCCATCAAGGTGaa                                                                                                             <  2:489970/60‑1 (MQ=255)
                                                                                                        cGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACt                                                                                                        <  1:284533/57‑1 (MQ=255)
                                                                                                        cGTAAATCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACt                                                                                                        >  2:284533/1‑57 (MQ=255)
                                                                                                           aaaTCTGAATGCGCGAGGTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTc                   <  2:162149/139‑1 (MQ=255)
                                                                                                                            gTCACGGCCGTCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATgg  <  2:13935/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                
GCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGG  >  NC_000913/1073598‑1073860

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