Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,073,996 Δ1 bp coding (16/1149 nt) rutA ← pyrimidine oxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,073,9930G.100.0% 67.0 / NA 17coding (19/1149 nt)rutApyrimidine oxygenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (4/13);  total (4/13)

CAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGT  >  NC_000913/1073914‑1074068
                                                                               |                                                                           
cagccGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACgg                                                                            <  2:197363/81‑1 (MQ=255)
  gccGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACgggc                                                                          <  1:197364/81‑2 (MQ=255)
        gttgCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGCGCTGc                                                                    <  2:197365/80‑1 (MQ=255)
         ttgCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGCGCTGCAt                                                                  <  2:197366/81‑1 (MQ=255)
                tAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGCGCTGCATCCTGCa                                                            >  2:197396/1‑80 (MQ=255)
                             gCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGCGCTGCATCCTGCATATCCTTTTCAGcc                                              <  1:197367/81‑1 (MQ=255)
                             gCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGCGCTGCATCCTGCATATCCATTTCAGcc                                              <  2:197368/81‑1 (MQ=255)
                                      catcatcaACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTg                                      <  1:197369/80‑1 (MQ=255)
                                       atcatcaACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGc                                    <  2:197370/81‑1 (MQ=255)
                                             atccTCTCTTCGTCTCGTAAAGTGAAAGTCATACGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGtt                               <  2:197371/78‑1 (MQ=255)
                                                   tctTCGTCTCGTAAAGTGAAAGTCAGACGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAg                        >  1:197394/1‑81 (MQ=255)
                                                     ttCGTCTCGTAAAGTGAAAGTCAGACGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCaaa                         <  1:197372/78‑1 (MQ=255)
                                                       agTCTCGTAAAGTGAAAGTCAGACGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGc                       <  2:197373/77‑1 (MQ=255)
                                                        gTCTCGTAAAGTGAAAGTCAGACGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGAt                   >  1:197395/1‑81 (MQ=255)
                                                        gTCTCGTAAAGTGAAAGTCAGACGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGAt                   <  2:197374/81‑1 (MQ=255)
                                                             gTAAAGTGAAAGTCAGACGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAg              <  2:197375/81‑1 (MQ=255)
                                                                   tGAAAGTCAGACGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGttttt         >  2:197393/1‑80 (MQ=255)
                                                                     aaaGTCAGACGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTaa       <  1:197376/80‑1 (MQ=255)
                                                                         tCAGACGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGt  <  1:197377/81‑1 (MQ=255)
                                                                               |                                                                           
CAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCTCTCTTCGTCTCGTAAAGTGAAAGTCAGACGGGGCGCTGCATCCTGCATATCCTTTTCAGCCGCGTATTGGCTTGTTTGCAAAGCGGATGCCAGTTTTTAAAAAGT  >  NC_000913/1073914‑1074068

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