Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,335,018 G→A L543L (CTG→TTG)  yfaA ← DUF2138 family protein, putative host defense protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,335,0180GA100.0% 81.5 / NA 25L543L (CTG→TTG) yfaADUF2138 family protein, putative host defense protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/10);  total (15/10)

TTCGACATTCAGCATTTCGCTATGGGCAACAACACAACACAGCCAGCAAATCAGCGCCAGCAGCCCGTGCCTCATAGCGGCTGCCAGGTTATGGGTAGCCACTGCCAGGCGGCACCGGGTTCCATCTGGGCCAGCTTCATGACATAACGCGGTTGTTGAGATAAAGCGTCCAGCTTCGGCATTAATAAAGTTTGTGCGGCGTTATAAAAAACCGGTTCGAGATTCTTCGGCAGACTGGTCAACGTTTCGTTACGCAGCAGTTTCGCTA  >  NC_000913/2334885‑2335152
                                                                                                                                     |                                                                                                                                      
ttCGACATTCAGCATTTCGCTATGGGCAACAACACAACACAGCCAGCAAATCAGCGCCAGCAGCCCGTGCCTCATAGCGGCTGCCAGGTTATGGGTAGCCACTGCCAGGCGGCACCGGGTTCCATCTGGGCCAACTTCa                                                                                                                                   >  2:518729/1‑139 (MQ=255)
        tCAGCATTTCGCTATGGGCAACAACACAACACAGCCAGCAAATCAGCGCCAGCAGCCCGTGCCTCATAGCGGCTGCCAGGTTATGGGTAGCCACTGCCAGGCGGCACCGGGTTCCATCTGGGCCAACTTCATGACATaa                                                                                                                           >  2:218788/1‑139 (MQ=255)
          aGCATTTCGCTATGGGCAACAACACAACACAGCCAGCAAATCAGCGCCAGCAGCCCGTGCCTCATAGCGGCTGCCAGGTTATGGGTAGCCACTGCCAGGCGGCACCGGGTTCCATCTGGGCCAACTTCATGACATAAcg                                                                                                                         <  1:525169/139‑1 (MQ=255)
                 cGCTATGGGCAACAACACAACACAGCCAGCAAATCAGCGCCAGCAGCCCGTGCCTCATAGCGGCTGCCAGGTTATGGGTAGCCACTGCCAGGCGGCACCGGGTTCCATCTGGGCCAACTTCATGACATAACGCGgttgt                                                                                                                  >  2:466085/1‑139 (MQ=255)
                       gggCAACAACACAACACAGCCAGCAAATCAGCGCCAGCAGCCCGTGCCTCATAGCGGCTGCCAGGTTATGGGTAGCCACTGCCAGGCGGCACCGGGTTCCATCTGGGCCAACTTCATGACATAACGCGGTTGTTGAGAt                                                                                                            >  1:302334/1‑139 (MQ=255)
                                cacaacacaGCCAGCAAATCAGCGCCAGCAGCCCGTGCCTCATAGCGGCTGCCAGGTTATGGGTAGCCACTGTCAGGCGGCACCGGGTTCCATCTGGGCCAACTTCATGACATAAcg                                                                                                                         <  1:449791/117‑1 (MQ=255)
                                cacaacacaGCCAGCAAATCAGCGCCAGCAGCCCGTGCCTCATAGCGGCTGCCAGGTTATGGGTAGCCACTGTCAGGCGGCACCGGGTTCCATCTGGGCCAACTTCATGACATAAcg                                                                                                                         >  2:449791/1‑117 (MQ=255)
                                  caacacaGCCAGCAAATCAGCGCCAGCAGCCCGTGCCTCATAGCGGCTGCCAGGTTATGGGTAGCCACTGCCAGGCGGCACCGGGTTCCATCTGGGCCAACTTCATGACATAACGCGGTTGTTGAGATAAAGCGTCCAg                                                                                                 >  2:539300/1‑139 (MQ=255)
                                   aacacaGCCAGCAAATCAGCGCCAGCAGCCCGTGCCTCATAGCGGCTGCCAGGTTATGGGTAGCCACTGCCAGGCGGCACCGGGTTCCATCTGGGCCAACTTCATGACATAAcg                                                                                                                         <  1:18029/114‑1 (MQ=255)
                                   aacacaGCCAGCAAATCAGCGCCAGCAGCCCGTGCCTCATAGCGGCTGCCAGGTTATGGGTAGCCACTGCCAGGCGGCACCGGGTTCCATCTGGGCCAACTTCATGACATAAcg                                                                                                                         >  2:18029/1‑114 (MQ=255)
                                        agccagcAAATCAGCGCCAGCAGCCCGTGCCTCATAGCGGCTGCCAGGTTATGGGTAGCCACTGCCAGGCGGCACCGGGTTCCATCTGGGCCAACTTCATGACATAACGCGGTTGTTGAGATAAAGCGTCCAGCTTCgg                                                                                           >  2:520667/1‑139 (MQ=255)
                                                        ccAGCAGCCCGTGCCTCATAGCGGCTGCCAGGTTATGGGTAGCCACTGCCAGGCGGCACCGGGTTCCATCTGGGCCAACTTCATGACATAACGCGgttgtt                                                                                                                 >  1:256817/1‑101 (MQ=255)
                                                        ccAGCAGCCCGTGCCTCATAGCGGCTGCCAGGTTATGGGTAGCCACTGCCAGGCGGCACCGGGTTCCATCTGGGCCAACTTCATGACATAACGCGgttgtt                                                                                                                 <  2:256817/101‑1 (MQ=255)
                                                         cagcagCCCGTGCCTCATAGCGGCTGCCAGGTTATGGGTAGCCACTGCCAGGCGGCACCGGGTTCCATCTGGGCCAACTTCATGACATAACGCGGTTGTTGAGATAAAGCGTCCAGCTTCGGCATTAATAAAGTTtgtg                                                                          <  2:302334/139‑1 (MQ=255)
                                                          agcagcCCGTGCCTCATAGCGGCTGCCAGGTTATGGGTAGCCACTGCCAGGCGGCACCGGGTTCCATCTGGGCCAACTTCATGACATAACGCGGTTGTTGAGATAAAGCGTCCAGCTTCGGCATTAATAAAGTTTGTgc                                                                         <  1:466085/139‑1 (MQ=255)
                                                                                cTGCCAGGTTATGGGTAGCCACTGCCAGGCGGCACCGGGTTCCATCTGGGCCAACTTCATGACATAACGCGGTTGTTGAGATAAAGCGTCCAGCTTCGGCATTAATAAAGTTTGTGCGGCGTTATAAAAAACCGGTTCg                                                   >  2:108704/1‑139 (MQ=255)
                                                                                               tAGCCACTGCCAGGCGGCACCGGGTTCCATCTGGGCCAACTTCATGACATAACGCGGTTGTTGAGATAAAGCGTCCAGCTTCGGCATTAATAAAGTTTGTGCGGCGTTATAAAAAACCGGTTCGAGATTCTTCGGCAGa                                    <  1:539300/139‑1 (MQ=255)
                                                                                                aGCCATTGCCAGGCGGCACCGGGTTCCATCTGGGCCAACTTCATGACATAACGCGGTTGTTGAGATAAAGCGTCCAGCTTCGGCATTAATAAAGTTTGTGCGGCGTTATAAAAAACCGGTTCGAGATTCTTCGGCAGAc                                   >  1:563327/1‑139 (MQ=255)
                                                                                                 gCCACTGCCAGGCGGCACCGGGTTCCATCTGGGCCAACTTCATGACATAACGCGGTTGTTGAGATAAAGCGTCCAGCTTCGGCATTAATAAAGTTTGTGCGGCGTTATAAAAAACCGGTTCGAGATTCTTCGGCAGACt                                  >  1:323132/1‑139 (MQ=255)
                                                                                                   cACTGCCAGGCGGCACCGGGTTCCATCTGGGCCAACTTCATGACATAACGCGGTTGTTGAGATAAAGCGTCCAGCTTCGGCATTAATAAAGTTTGTGCGGCGTTATAAAAAACCGGTTCGAGATTCTTCGGCAGACTgg                                <  2:323132/139‑1 (MQ=255)
                                                                                                                      gTTCCATCTGGGCCAACTTCATGACATAACGCGGTGGTTGAGATAAAGCGTCCAGCTTCGGCATTAATAAAGTTTGTGCGGCGTTATAAAAAACCGGTTCGAGATTCTTCGGCAGACTGGTCAACGTTTCGTTACgcag             <  1:520667/139‑1 (MQ=255)
                                                                                                                                ggCCAACTTCATGACATAACGCGGTTGTTGAGATAAAGCGTCCAGCTTCGGCATTAATAAAGTTTGTGCGGCGTTATAAAAAACCGGTTCGAGATTCTTCGGCAGACTgg                                <  2:40763/110‑1 (MQ=255)
                                                                                                                                ggCCAACTTCATGACATAACGCGGTTGTTGAGATAAAGCGTCCAGCTTCGGCATTAATAAAGTTTGTGCGGCGTTATAAAAAACCGGTTCGAGATTCTTCGGCAGACTgg                                >  1:40763/1‑110 (MQ=255)
                                                                                                                                 gCCAACTTCATGACATAACGCGGTTGTTGAGATAAAGCGTCCAGCTTCGGCATTAATAAAGTTTGTGCGGCGTTATAAAAAACCGGTTCGAGATTCTTCGGCAGACTGGTCAACGTTTCGTTACGCAGCAGTTTCGCTa  >  2:204643/1‑139 (MQ=255)
                                                                                                                                 gCCAACTTCATGACATAACGCGGTTGTTGAGATAAAGCGTCCAGCTTCGGCATTAATAAAGTTTGTGCGGCGTTATAAAAAACCGGTTCGAGATTCTTCGGCAGACTGGTCAACGTTTCGTTACGCAGCAGTTTCGCTa  >  1:539239/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                      
TTCGACATTCAGCATTTCGCTATGGGCAACAACACAACACAGCCAGCAAATCAGCGCCAGCAGCCCGTGCCTCATAGCGGCTGCCAGGTTATGGGTAGCCACTGCCAGGCGGCACCGGGTTCCATCTGGGCCAGCTTCATGACATAACGCGGTTGTTGAGATAAAGCGTCCAGCTTCGGCATTAATAAAGTTTGTGCGGCGTTATAAAAAACCGGTTCGAGATTCTTCGGCAGACTGGTCAACGTTTCGTTACGCAGCAGTTTCGCTA  >  NC_000913/2334885‑2335152

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