Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 904,547 G→A intergenic (‑80/‑46) ybjP ← / → ybjQ DUF3828 domain‑containing lipoprotein YbjP/putative heavy metal binding protein YbjQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913904,5470GA100.0% 64.2 / NA 19intergenic (‑80/‑46)ybjP/ybjQDUF3828 domain‑containing lipoprotein YbjP/putative heavy metal binding protein YbjQ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/5);  total (14/5)

GATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTG  >  NC_000913/904471‑904619
                                                                            |                                                                        
gATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCtaataa                                                                      >  2:168090/1‑81 (MQ=255)
             gCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTAATAATGCtttt                                                               >  2:168082/1‑75 (MQ=255)
                  cGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTAATAATGCTTTTGAGATCGAAgg                                                    >  1:168100/1‑81 (MQ=255)
                    ctctGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTAATAATGCTTTTGAGATCGAAGGCt                                                  >  1:168118/1‑81 (MQ=255)
                       tGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTAATAATGCTTTTGAGATCGAAGGCtt                                                 <  2:168064/79‑1 (MQ=255)
                           aGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTAATAATGCTTTTGAGATCGAAGGCTTAGCaaa                                           <  1:168065/81‑1 (MQ=255)
                             cGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTAATAATGCTTTTGAGATCGAAGGCTTAGCAAAc                                          >  2:168094/1‑80 (MQ=255)
                                tCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTAATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAg                                       >  2:168088/1‑80 (MQ=255)
                                 cTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTAATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAgg                                      <  1:168066/80‑1 (MQ=255)
                                   aCGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTAATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAg                                       >  2:168091/1‑77 (MQ=255)
                                       taataaTGATATCCTTTCAATAATAGCGTATCAGTCTAATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAgatcga                               >  2:168086/1‑81 (MQ=255)
                                             tGATATCCTTTCAATAATAGCGTATCAGTCTAATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGc                         <  1:168067/81‑1 (MQ=255)
                                              gATATCCTTTCAATAATAGCGTATCAGTCTAATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGc                         >  2:168083/1‑80 (MQ=255)
                                              gATATCCTTTCAATAATAGCGTATCAGTCTAATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGc                         >  1:168105/1‑80 (MQ=255)
                                                  tCCTTTCAATAATAGCGTATCAGTCTAATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAAttt                    <  2:168068/81‑1 (MQ=255)
                                                    cTTTCAATAATAGCGTATCAGTCTAATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTc                  >  2:168098/1‑81 (MQ=255)
                                                        cAATAATAGCGTATCAGTCTAATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTAc               >  2:168080/1‑80 (MQ=255)
                                                            aataGCGTATCAGTCTAATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAAccc          >  2:168108/1‑81 (MQ=255)
                                                                  gTATCAGTCTAATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAActc    >  2:168095/1‑81 (MQ=255)
                                                                    aTCAGTCTAATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTg  >  1:168069/1‑81 (MQ=255)
                                                                            |                                                                        
GATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTG  >  NC_000913/904471‑904619

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