Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 442,756 G→A G403D (GGC→GAC)  thiI → tRNA s(4)U8 sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913442,7560GA100.0% 145.5 / NA 46G403D (GGC→GAC) thiItRNA s(4)U8 sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/26);  total (20/26)

GAAGCGGAAGAAGAGAAGTTCGACTTCAGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGGCTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGATTGATGTGGTTTCTCTGCCGTTCTATAAACTGAGCACCAAATTTGGCGATCTC  >  NC_000913/442629‑442886
                                                                                                                               |                                                                                                                                  
gAAGCGGAAGAAGAGAAGTTCGACTTCAGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGa                                                                                                                         <  1:752691/139‑1 (MQ=255)
      gaagaaGAGAAGTTCGACTTCAGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAacgacg                                                                                                                   <  2:900381/139‑1 (MQ=255)
         gaagagaagTTCGACTTCAGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGa                                                                                                                <  1:550510/139‑1 (MQ=255)
           agagaagTTCGACTTCAGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATc                                                                                                              <  2:295149/139‑1 (MQ=255)
              gaagTTCGACTTCAGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTc                                                                                                           <  1:993571/139‑1 (MQ=255)
              gaagTTCGACTTCAGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTc                                                                                                           <  1:764831/139‑1 (MQ=255)
                 gTTCGACTTCAGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGat                                                                                                        >  1:393798/1‑139 (MQ=255)
                      aCTTCAGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGATATCCg                                                                                                   <  1:557716/139‑1 (MQ=255)
                                 ctcGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATg                                                                                        >  2:35316/1‑139 (MQ=255)
                                  tcGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGatat                                                                                                      <  2:973748/124‑1 (MQ=255)
                                  tcGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGatat                                                                                                      >  1:973748/1‑124 (MQ=255)
                                    gATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAAc                                                                                     <  1:805743/139‑1 (MQ=255)
                                            ggTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGGGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGAt                                                                             >  1:79832/1‑139 (MQ=255)
                                              ttGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATaa                                                                           >  2:210132/1‑139 (MQ=255)
                                               tGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAg                                                                          <  1:35316/139‑1 (MQ=255)
                                                       cGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGaa                                                                  >  1:639590/1‑139 (MQ=255)
                                                         aaTAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATc                                                                                                              >  1:586170/1‑93 (MQ=255)
                                                         aaTAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATc                                                                                                              <  2:586170/93‑1 (MQ=255)
                                                                   atatCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGatat                                                                                                      >  2:33915/1‑91 (MQ=255)
                                                                   atatCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGatat                                                                                                      <  1:33915/91‑1 (MQ=255)
                                                                   atatCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGat                                                                                                        >  1:406508/1‑89 (MQ=255)
                                                                   atatCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGat                                                                                                        <  2:406508/89‑1 (MQ=255)
                                                                   atatCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGAt                                                      >  2:1059576/1‑139 (MQ=255)
                                                                         gcgAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGa                                                                                                                >  2:536916/1‑75 (MQ=255)
                                                                         gcgAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGa                                                                                                                <  1:536916/75‑1 (MQ=255)
                                                                                     agcagACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGATTGATGTGGTTTCTCTGcc                                    <  1:210132/139‑1 (MQ=255)
                                                                                       cagACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGATTGATGTGGTTTCTCTGCCGt                                  <  2:639590/139‑1 (MQ=255)
                                                                                       cagACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGATTGATGTGGTTTCTCTGCCGt                                  <  2:79832/139‑1 (MQ=255)
                                                                                            cGAGCAGGAAGTGGTGGAAGTGGAAACCGTCACTGACTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGATTGATGTGGTTTCTCTg                                      <  2:643210/130‑1 (MQ=255)
                                                                                            cGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGATTGATGTGGTTTCTCTg                                      >  1:643210/1‑130 (MQ=255)
                                                                                             gAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGAt                                                                             <  1:1008964/90‑1 (MQ=255)
                                                                                             gAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGAt                                                                             >  2:1008964/1‑90 (MQ=255)
                                                                                             gAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGATTGATGTGGTTTCTCTGcc                                    >  2:559774/1‑131 (MQ=255)
                                                                                             gAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGATTGATGTGGTTTCTCTGcc                                    <  1:559774/131‑1 (MQ=255)
                                                                                               gCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGATTGATGTGGTTTCTCTGCCGTTCTATaaa                          >  2:347081/1‑139 (MQ=255)
                                                                                               gCAGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGATTGATGTGGTTTCTCTGCCGTTCTATaaa                          <  1:1059576/139‑1 (MQ=255)
                                                                                                 aGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGATTGATGTGGTTTCTCTg                                      <  2:359443/125‑1 (MQ=255)
                                                                                                 aGGAAGTGGTGGAAGTGGAAACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGATTGATGTGGTTTCTCTg                                      >  1:359443/1‑125 (MQ=255)
                                                                                                         gtggaagtggaaACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGATATCCGTTc                                                                                                <  2:21025/59‑1 (MQ=255)
                                                                                                         gtggaagtggaaACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGATATCCGTTc                                                                                                >  1:21025/1‑59 (MQ=255)
                                                                                                         gtggaagtggaaACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAg                                                                                   >  1:953438/1‑72 (MQ=255)
                                                                                                         gtggaagtggaaACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAg                                                                                   <  2:953438/72‑1 (MQ=255)
                                                                                                               gtggaaACCGTCAATGACTTCGGCCCGAACGACGTGATCCTc                                                                                                           >  2:1019702/1‑42 (MQ=255)
                                                                                                               gtggaaACCGTCAATGACTTCGGCCCGAACGACGTGATCCTc                                                                                                           <  1:1019702/42‑1 (MQ=255)
                                                                                                                  gaaACCGTCAATGACTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGATTGATGTGGTTTCTCTGCCGTTCTATAAACTGAGCACCAAATTTGGCg       <  2:393798/139‑1 (MQ=255)
                                                                                                                       cGTCAATGACTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGCAAGTCGAAGGGATTGATGTGGTTTCTCTGACGTTCTATAAACTTAGCACCAAATTTGGCGAtctc  >  2:511659/1‑139 (MQ=255)
                                                                                                                               |                                                                                                                                  
GAAGCGGAAGAAGAGAAGTTCGACTTCAGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGGCTTCGGCCCGAACGACGTGATCCTCGATATCCGTTCTATCGATGAACAGGAAGATAAGCCACTGAAAGTCGAAGGGATTGATGTGGTTTCTCTGCCGTTCTATAAACTGAGCACCAAATTTGGCGATCTC  >  NC_000913/442629‑442886

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: