Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,103,871 (A)7→6 intergenic (+71/+78) BW25113_RS15410 → / ← yghD tRNA‑Phe/GspM family type II secretion system protein YghD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,103,8650A.100.0% 67.3 / NA 17intergenic (+65/+84)BW25113_RS15410/yghDtRNA‑Phe/GspM family type II secretion system protein YghD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (7/10);  total (7/10)

GGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTT  >  NZ_CP009273/3103776‑3103949
                                                                                         |                                                                                    
ggTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGa                                                                                      <  2:497602/90‑1 (MQ=255)
 gTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGaa                                                                                     <  1:165334/90‑1 (MQ=255)
      aTTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGaaaaaat                                                                                <  1:417155/90‑2 (MQ=255)
                cGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATaa                                                                     <  2:355629/90‑1 (MQ=255)
                                ttAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGc                                                                           <  1:235844/68‑1 (MQ=255)
                                ttAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGc                                                                           >  2:235844/1‑68 (MQ=255)
                                ttAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTAc                                                     >  1:329207/1‑90 (MQ=255)
                                ttAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTAc                                                     >  2:30365/1‑90 (MQ=255)
                                      aCCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGc                                               >  1:250922/1‑90 (MQ=255)
                                      aCCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGc                                               <  1:400055/90‑1 (MQ=255)
                                       cccGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGcc                                              >  2:490801/1‑90 (MQ=255)
                                         cGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCAc                                            >  1:134660/1‑90 (MQ=255)
                                                     gggTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGt                                <  2:289212/90‑1 (MQ=255)
                                                     gggTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGt                                <  1:235498/90‑1 (MQ=255)
                                                          tttGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAAttt                           <  1:22747/90‑1 (MQ=255)
                                                                 ttGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGat                    <  2:250922/90‑1 (MQ=255)
                                                                     aTCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTa                <  2:176325/90‑1 (MQ=255)
                                                                               aaCCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGAt      <  2:368368/90‑1 (MQ=255)
                                                                               aaCCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGAt      <  2:408439/90‑1 (MQ=255)
                                                                                   ttCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACtt  >  2:13510/1‑90 (MQ=255)
                                                                                         |                                                                                    
GGTTAGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTGATACTT  >  NZ_CP009273/3103776‑3103949

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