Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,503,125 G→A *223* (TAG→TAA ydcL → lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,503,1250GA100.0% 89.6 / NA 29*223* (TAG→TAAydcLlipoprotein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/16);  total (13/16)

GTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAGTTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCAATGTATTGGCG  >  NC_000913/1503000‑1503253
                                                                                                                             |                                                                                                                                
gTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGa                                                                                                                           >  2:856025/1‑133 (MQ=255)
gTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGa                                                                                                                           <  1:856025/133‑1 (MQ=255)
                              aaTAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGa                                                                                                                           >  2:576869/1‑103 (MQ=255)
                              aaTAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGa                                                                                                                           <  1:576869/103‑1 (MQ=255)
                                   cGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTa                                                                                  >  1:793802/1‑139 (MQ=255)
                                        gTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAAcc                                                                             >  1:1188085/1‑139 (MQ=255)
                                           cGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCa                                                                          <  2:859827/139‑1 (MQ=255)
                                           cGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCa                                                                          <  1:278688/139‑1 (MQ=255)
                                           cGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCa                                                                          >  2:1181267/1‑139 (MQ=255)
                                                 tGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                      >  1:241998/1‑89 (MQ=255)
                                                 tGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                      <  2:241998/89‑1 (MQ=255)
                                                                     cAAGTTATTGATGACATGGCGACCGATGCCACCATTTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGgcgc                                                <  1:781880/139‑1 (MQ=255)
                                                                               atgaCATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCAc                                      <  1:1181267/139‑1 (MQ=255)
                                                                                   cATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                      >  2:699565/1‑55 (MQ=255)
                                                                                   cATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                      >  2:739498/1‑55 (MQ=255)
                                                                                   cATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                      <  1:739498/55‑1 (MQ=255)
                                                                                   cATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                      <  1:699565/55‑1 (MQ=255)
                                                                                   cATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATgg                                                                                                         <  2:1060165/68‑1 (MQ=255)
                                                                                   cATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATgg                                                                                                         >  1:1060165/1‑68 (MQ=255)
                                                                                      ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCa                                                                                                                    >  1:490750/1‑54 (MQ=255)
                                                                                      ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCa                                                                                                                    <  2:490750/54‑1 (MQ=255)
                                                                                      ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCa                               <  2:1188085/139‑1 (MQ=255)
                                                                                      ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCa                               >  1:689586/1‑139 (MQ=255)
                                                                                                     caTGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTACAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCtgttgt                >  2:879663/1‑139 (MQ=255)
                                                                                                       tGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTc              <  2:1296700/139‑1 (MQ=255)
                                                                                                        gTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCa             <  2:689586/139‑1 (MQ=255)
                                                                                                                  aaCAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAAcccc                                                                           <  2:1272956/67‑1 (MQ=255)
                                                                                                                  aaCAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAAcccc                                                                           >  1:1272956/1‑67 (MQ=255)
                                                                                                                   aCAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCAATGTATTGGCg  <  1:1098213/139‑1 (MQ=255)
                                                                                                                             |                                                                                                                                
GTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAGTTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCAATGTATTGGCG  >  NC_000913/1503000‑1503253

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