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% 104.6 / NA 33*223* (TAG→TAAydcLlipoprotein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/16);  total (17/16)

GTTATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAGTTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCAATGTATTGG  >  NC_000913/1502991‑1503251
                                                                                                                                      |                                                                                                                              
gTTATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTcc                                                                                                                            >  2:935785/1‑139 (MQ=255)
     cAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGAcg                                                                                                                       <  1:935785/139‑1 (MQ=255)
       aaGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGAcgcg                                                                                                                     <  2:1328857/139‑1 (MQ=255)
                                aGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATc                                                                                            >  2:351445/1‑139 (MQ=255)
                                      gAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTgc                                                                                      >  2:608630/1‑139 (MQ=255)
                                      gAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTgc                                                                                      >  1:45173/1‑139 (MQ=255)
                                            cGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTa                                                                                >  1:1215969/1‑139 (MQ=255)
                                                 gTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATgg                                                                                                       >  2:847935/1‑111 (MQ=255)
                                                 gTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATgg                                                                                                       <  1:847935/111‑1 (MQ=255)
                                                            gCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAg                                                                >  1:410000/1‑139 (MQ=255)
                                                              tttCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAggg                                                              <  1:1420005/139‑1 (MQ=255)
                                                              tttCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAggg                                                              <  1:351445/139‑1 (MQ=255)
                                                              tttCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAggg                                                              <  1:608630/139‑1 (MQ=255)
                                                                tCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCg                                                            <  2:45173/139‑1 (MQ=255)
                                                                 cGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGgcg                                                                                                     <  2:694987/97‑1 (MQ=255)
                                                                 cGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGgcg                                                                                                     >  1:694987/1‑97 (MQ=255)
                                                                      atatTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGata                                                      <  2:1215969/139‑1 (MQ=255)
                                                                                 gTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTgc                                                                                      <  2:1049597/96‑1 (MQ=255)
                                                                                 gTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTgc                                                                                      >  1:1049597/1‑96 (MQ=255)
                                                                                        atgaCATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCg                                                            <  2:1435573/115‑1 (MQ=255)
                                                                                        atgaCATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCg                                                            >  1:1435573/1‑115 (MQ=255)
                                                                                                       aTGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAg                            <  1:761894/132‑1 (MQ=255)
                                                                                                       aTGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAg                            >  2:761894/1‑132 (MQ=255)
                                                                                                       aTGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTg                     >  2:1164532/1‑139 (MQ=255)
                                                                                                                tGTTTGACGTTAACGAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCg                                                                                          >  1:915853/1‑61 (MQ=255)
                                                                                                                tGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCg                                                                                          <  2:915853/61‑1 (MQ=255)
                                                                                                                          tAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCAATGTATTgg  >  2:984735/1‑139 (MQ=255)
                                                                                                                           aaCAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCttt                                                                                                 >  2:1112580/1‑43 (MQ=255)
                                                                                                                           aaCAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCttt                                                                                                 <  1:1112580/43‑1 (MQ=255)
                                                                                                                           aaCAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTg                                                                                       >  2:251209/1‑53 (MQ=255)
                                                                                                                           aaCAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTg                                                                                       <  1:251209/53‑1 (MQ=255)
                                                                                                                            aCAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCa                             >  2:619667/1‑110 (MQ=255)
                                                                                                                            aCAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCa                             <  1:619667/110‑1 (MQ=255)
                                                                                                                                      |                                                                                                                              
GTTATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAGTTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCAATGTATTGG  >  NC_000913/1502991‑1503251

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