Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,180,586 (A)7→6 intergenic (+333/+77) yjaZ → / ← thiH heat‑shock protein C/2‑iminoacetate synthase ThiH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,180,5800A.100.0% 38.2 / NA 10intergenic (+327/+83)yjaZ/thiHheat‑shock protein C/2‑iminoacetate synthase ThiH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (3/7);  total (3/7)

TTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCG  >  NZ_CP009273/4180492‑4180660
                                                                                        |                                                                                
ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaa                                                                                 >  2:23019/1‑90 (MQ=255)
            aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGc                                                                    <  1:64952/90‑1 (MQ=255)
                                                   tCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTc                             <  2:269240/90‑1 (MQ=255)
                                                                caTTCCTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTAc                <  1:213657/90‑1 (MQ=255)
                                                                    ccTGCTGACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACCTTc            >  1:346843/1‑90 (MQ=255)
                                                                         tgACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACa       >  2:149566/1‑90 (MQ=255)
                                                                          gACGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACAt      <  2:303903/90‑1 (MQ=255)
                                                                            cGGCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTc    <  2:307277/90‑1 (MQ=255)
                                                                              gCACTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCg  <  2:92844/90‑1 (MQ=255)
                                                                                 cTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCg  >  1:45137/1‑87 (MQ=255)
                                                                                 cTACCATAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCg  <  2:45137/87‑1 (MQ=255)
                                                                                        |                                                                                
TTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCG  >  NZ_CP009273/4180492‑4180660

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