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% 60.1 / NA 15intergenic (+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 . (4/11);  total (4/11)

TCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTG  >  NZ_CP009273/3103784‑3103943
                                                                                 |                                                                              
tCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGaaaaaatta                                                                        <  1:130631/90‑4 (MQ=255)
        cGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATaa                                                               <  2:236239/90‑1 (MQ=255)
        cGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATaa                                                               <  1:82773/90‑1 (MQ=255)
                        ttAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGaaaaaat                                                                          <  2:338674/64‑2 (MQ=255)
                        ttAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGaaaaaat                                                                          >  1:338674/1‑63 (MQ=255)
                        ttAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTAc                                               >  1:234894/1‑90 (MQ=255)
                              aCCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGc                                         <  1:74595/90‑1 (MQ=255)
                                 cGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCAc                                      >  2:103546/1‑90 (MQ=255)
                                        aaGGCGGGTTTTTGCTTTTTGAACTGACAATAACCTTTCCGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACg                               <  2:227891/90‑1 (MQ=255)
                                         aGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACga                              >  2:187629/1‑90 (MQ=255)
                                             gggTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGt                          <  1:59481/90‑1 (MQ=255)
                                             gggTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGt                          <  2:199039/90‑1 (MQ=255)
                                                               cTGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAAt        <  1:343248/90‑1 (MQ=255)
                                                                tGACAATAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAAtc       <  1:41565/90‑1 (MQ=255)
                                                                    aaTAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCt   <  2:281722/90‑1 (MQ=255)
                                                                     aTAACCTTCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTg  <  2:248901/90‑1 (MQ=255)
                                                                           ttCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACg                               <  2:33122/55‑1 (MQ=255)
                                                                           ttCACGAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACg                               >  1:33122/1‑55 (MQ=255)
                                                                                 |                                                                              
TCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACGTTAACGAGAGTAATTTTATTGATATTAATCTCCTG  >  NZ_CP009273/3103784‑3103943

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