Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 606,147 G→A A18V (GCT→GTT)  ybdJ ← DUF1158 domain‑containing protein YbdJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913606,1470GA100.0% 39.9 / NA 16A18V (GCT→GTT) ybdJDUF1158 domain‑containing protein YbdJ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/7);  total (9/7)

GCCGCGAAGTACGGCGATAATCGACGTTTAATCCGCCAGCGAGAACCAGCGTCGCCAGATAAAGCGCAGAACAAAATACTCAATAGCGCCCAGCACTAAAAACCACAGACAAAACAATAAAGTGTAAAGCTGACTAAGATCCATCAGATGGAACATGGTCACCAGTTTTTGTGCCAGCGCCAGCCCCAGTGCGGGGGCGGGCAGCAGCAGGCAAGAGAGGAAAGCCATCAGCAAAATGCCTGCTGCGGTGGTCAAGGTTTCTAAAGGGTGTTTCATAACATTGTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACG  >  NC_000913/605924‑606370
                                                                                                                                                                                                                               |                                                                                                                                                                                                                               
gCCGCGAAGTACGGCGATAATCGACGTTTAATCCGCCAGCGAGAACCAGCGTCGCCAGATAAAGCGCAGAACAAAATACTCAATAGCGCCCAGAACTAACAACCACAGACAACACAATAAAGTGTAAAGCTGACTAAGATCCATCAGATGGAACATGGTCACCAGTTTTTGTGCCAGCGCCAGCCCCAGTGCGGGGGCGGGCAGCAGCAGGCAAGAGAGGAAAACCATCa                                                                                                                                                                                                                           >  1:218214/1‑230 (MQ=255)
     gAAGTACGGCGATAATCGACGTTTAAACCGCCAGCGAGAACCAGCGTCGCCAGATAAAGCGCAGAACAAAATACTCAATAGCGCCCAGCACTAAAAACCACAGACAAAACAATAAAGTGTAAAGCTGACTAAGATCCATCAGATGGAACATGGTCACCAGTGTTTGTGCCAGGGCGAGCCCCAATGCGGGGGCGGGCAGCAGCAGGCAATAGAGGAAAACCATCAGCAAAATGCCTGCTGCGGTGGTCAAg                                                                                                                                                                                                 <  1:140830/251‑1 (MQ=255)
                                                    cGCCAGATAAAGCGCAGAACAAACTACTCAATACCGCCCAGCACTAAAACCACAGCCAAAAACATAAAGTGTAAAGCTGACTAAGCTCCATCAGATGGACCATGGTCACCAGTTTTTGTGCCAGCGCCAGCACCAGTGAGGGGGCGGGCAGCAGCAGGCAAGAGAGGAAAACCATCAGCAAAATGCCTGCTGCGGTGGTCAAGGTTTCTAAAGGGAGTTTCATAACATTGTTAAATGTAAGTTAAAACAcc                                                                                                                                                 >  1:139002/1‑251 (MQ=255)
                                                                          aaTACTCAATAGCGCCCAGCACTAAAAACCACAGACAAAACAATAAAGTGTAAAGCTGACTAAGATCCATCAGATGGAACATGGTCACCAGTTTTTGTGCCAGCGCCAGCCCCAGTGCGGGGGCGGGCAGCAGCAGGCAAGAGAGGAAAACCATCAGCAAAATGCCTGCTGCGGTGGTCAAGGTTTCTAAAGGGTGTTTCATAACATTGTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTg                                                                                                                            <  1:54850/251‑1 (MQ=255)
                                                                          aaTACTCAATAGCGCCCAGCACTAAAAACCACAGACAAAACAATAAAGTGTAAAGCTGACTAAGATCCATCAGATGGAACATGGTCACCAGTTTTTGTGCCAGCGCCAGCCCCAGTGCGGGGGCGGGCAGCAGCAGGCAAGAGAGGAAAACCATCAGCAAAATGCCTGCTGCGGTGGTCAAGGTTTCTAAAGGGTGTTTCATAACATTGTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTg                                                                                                                            <  1:54859/251‑1 (MQ=255)
                                                                            tACTCAATAGCGCCCAGCACTAAAAACCACAGACAAAACAATAAAGTGTAAAGCTGACTAAGATCCATCAGATGGAACATGGTCACCAGTTTTTGTGCCAGCGCCAGCCCCCGTGCGGGGGCGGGCAGCAGCAGGCAAGAGAGGAAAACCATCAGCAAAATGCCTGCTGCGGTGGTCAAGGTTTCTAAAGGGTGTTTCATAACATTGTTAAATGTAAGTTAAAACACCATTGTCAg                                                                                                                                         >  1:183759/1‑236 (MQ=255)
                                                                                  aTAGCGCCCAGCACTAAAAACCACAGACAAAACAATACACTGTAAAGCTGACTAAAATCCATCAGATGGAACATGGTCACCAGTTTTTGTGCCAGCGCCAGCCCCAGTGCGGGGGCGGGCAGCCGCAGGCAAGAGAGGAAAACCATCAGCAAAATGCCTGCTGCGGTGGTCAAGGTTTCTAAAGGGTGTTTCATAACATTGtt                                                                                                                                                                    >  1:329312/1‑203 (MQ=255)
                                                                                       gcCCAGCACTACAAAACACAGACAAAACAATCAAGTGTAAAGCTGACTAAGATCCATCAGATGGAACATGGTCACCAGTTTTTGTGCCAGCGCCAGCCCCAGTGCGGGGGCGGGCAGCAGCAGGCAAGAGAGGAAAACCATCAGCAAACTGCCTGCTGCGGTGGTCAAGGTTTCTAAAGGGTGTTTCAt                                                                                                                                                                             >  1:253210/1‑189 (MQ=255)
                                                                                        nccAGCACTAAAAACCACAGACAAAACAATAAAGTGTAAAGCTGACTAAGATCCATCAGATGGAACATGGTCACCAGTTTTTGTGCCAGCGCCAGCCCCAGTGCGGGGGCGGGCAGCAGCAGGCAAGAGAGGAAAACCATCAGCAAAATGCCTGCTGCGGTGGTCAAGGTTTCTAAAGGGTGTTTCATAACATTGTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCg                                                                                                              >  1:78313/2‑251 (MQ=255)
                                                                                                 aaaaaCCACAGACAAAACAATAAAGTGTAAAGCTGACTAAGAGCCATCAGATGGAACATGGTCACCAGTTTTTGTGCCAGCGCCAGCCCCAGTGCGGGGGCGGGCAGCAGCAGGCAAGAGAGGAAAACCATCAGCAAAATGCCTGCTGCGGTGGTCAAGGTTTCTAAAGGGTGTTTCATAACATTGTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGTCTTAGTc                                                                                                      <  1:260299/250‑1 (MQ=255)
                                                                                                                                            ccATCAGATGGAACATGGTCACCATTTTTTGTGCCAGCGCCAGCCCCAGTGCGGTGGCGGTCAGCAGCATGCAAGAGAGGAAAACCATCAGCAAAATGCCTGCTTGGGTGGTCAAGGTGTCGTAAGGGTGTTTGATAACATTGTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGc                                                                                                                      <  1:149769/191‑1 (MQ=255)
                                                                                                                                                 aGATGGAACATGGTCACCAGTTTTTGTGCCAGCGCCAGCCCCAGTGCGGGGGCGGGCAGCAGCAGGCAAGAGAGGAAAACCATCAGCAAAATGCCTGCTGCGGTGGTCAAGGTTTCTAAAGGGTGTTTCATAACATTGTTAAATGTAAGTTAAAACACCATTGTc                                                                                                                                           <  1:23639/165‑1 (MQ=255)
                                                                                                                                                            ngTCACCAGTTTTTGTGCCAGCGCCAGCCCCAGTGCGGGGGCGGGCAGCAGCAGGCAAGAGAGGAAAACCATCAGCAAAATGCCTGCTGCGGTGGTCAAGGTTTCTAAAGGGTGTTTCATAACATTGTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGcgac                                          >  1:57061/2‑251 (MQ=255)
                                                                                                                                                                 ncAGTTTTTGTGCCAGCGCCAGCCCCAGTGCGGGGGCGGGCAGCAGCAGGCAAGAGAGGAAAACCATCAGCAAAATGCCTGCTGCGGTGGTCAAGGTTTCTAAAGGGTGTTTCATAACATTGTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAAt                                     >  1:93499/2‑251 (MQ=255)
                                                                                                                                                                                                    gcgggcAGCAGCAGGCAAGAGAGGAAAACCATCAGCAAAATGCCTGCTGCGGAGGTCAAGGTTTCTAAAGGGTGTTTCATAACATTGTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTgacg  >  1:230057/1‑251 (MQ=255)
                                                                                                                                                                                                                           gAAAACCATCAGCAAAATGCCTTCTGCTGTGGTCAAGGTTTCGAAAGGGTGGTTCATAACATTGTTAAATGTAAGTTAAAACGCCATTGTCGGGGATAGGCTTCTGGTAGTCAAGTCCCTGCTTTGTCACCGTCCCAGATCTCACAGTGCGTTTTCACCAGCCCAATCAGCGAGCCGCCATCgg                                              <  1:210618/184‑1 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                               
GCCGCGAAGTACGGCGATAATCGACGTTTAATCCGCCAGCGAGAACCAGCGTCGCCAGATAAAGCGCAGAACAAAATACTCAATAGCGCCCAGCACTAAAAACCACAGACAAAACAATAAAGTGTAAAGCTGACTAAGATCCATCAGATGGAACATGGTCACCAGTTTTTGTGCCAGCGCCAGCCCCAGTGCGGGGGCGGGCAGCAGCAGGCAAGAGAGGAAAGCCATCAGCAAAATGCCTGCTGCGGTGGTCAAGGTTTCTAAAGGGTGTTTCATAACATTGTTAAATGTAAGTTAAAACACCATTGTCAGGGATATTCTTCTGTAAGGCAATTCCCGGCTTAGTCACCGGCCCAGATCTCACAATGCTTTTTCACCAGCCCAATCAGCGAGCCGCCATCGGCGACGAAATCGCGCATTAGCTGCGCTTCATTCAGACCGCTGACG  >  NC_000913/605924‑606370

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