Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,634,151 Δ1 bp coding (164/534 nt) rrrQ ← lysozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,634,1500G.100.0% 105.3 / NA 25coding (165/534 nt)rrrQlysozyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (10/15);  total (10/15)

CCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAA  >  NZ_CP009273/1634062‑1634235
                                                                                        |                                                                                     
ccTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACtgtgc                                                                                      <  2:660833/90‑2 (MQ=255)
  ttATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACtgtgcac                                                                                    >  1:1135722/1‑87 (MQ=255)
  ttATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACtgtgcac                                                                                    >  2:998859/1‑87 (MQ=255)
  ttATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACtgtgcac                                                                                    >  2:1228811/1‑87 (MQ=255)
   tATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACtgtgcacc                                                                                   <  2:363072/90‑5 (MQ=255)
            cAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGAt                                                                         <  2:615465/90‑1 (MQ=255)
            cAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGAt                                                                         >  1:94964/1‑90 (MQ=255)
              aTGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATgg                                                                       >  2:345427/1‑90 (MQ=255)
              aTGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATgg                                                                       <  1:626613/90‑1 (MQ=255)
                   gTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCat                                                                  <  1:1228811/90‑1 (MQ=255)
                   gTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCat                                                                  >  2:696772/1‑90 (MQ=255)
                   gTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCat                                                                  >  2:111153/1‑90 (MQ=255)
                   gTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCat                                                                  <  1:696772/90‑1 (MQ=255)
                      gACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATg                                                               >  2:359763/1‑90 (MQ=255)
                           ggTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGa                                                          >  1:30641/1‑90 (MQ=255)
                                 aTTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATc                                                    <  2:335123/90‑1 (MQ=255)
                                   ttttCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATcgc                                                  <  2:94964/90‑1 (MQ=255)
                                   ttttCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATcgc                                                  <  2:858906/90‑1 (MQ=255)
                                   ttttCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATcgc                                                  <  1:591649/90‑1 (MQ=255)
                                   ttttCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATcgc                                                  <  2:720063/90‑1 (MQ=255)
                                     ttCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCgg                                                <  2:1100665/90‑1 (MQ=255)
                                     ttCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCgg                                                <  1:495612/90‑1 (MQ=255)
                                            gACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCa                                         >  2:81737/1‑90 (MQ=255)
                                                gTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATtgt                                     >  1:534995/1‑90 (MQ=255)
                                                       atTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTa                              <  1:386375/90‑1 (MQ=255)
                                                          gggAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACct                           >  2:45721/1‑90 (MQ=255)
                                                          gggAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACct                           >  2:277435/1‑90 (MQ=255)
                                                                            atccaCCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAAc         <  1:617261/90‑1 (MQ=255)
                                                                              ccaccaCTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTg       <  1:81737/90‑1 (MQ=255)
                                                                                   acTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCaa  <  2:30641/90‑1 (MQ=255)
                                                                                        |                                                                                     
CCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAA  >  NZ_CP009273/1634062‑1634235

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: