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% 65.9 / NA 17coding (165/534 nt)rrrQlysozyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (4/13);  total (4/13)

TTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTC  >  NZ_CP009273/1634064‑1634233
                                                                                      |                                                                                   
ttATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACtgtgcac                                                                                  <  2:257775/90‑4 (MQ=255)
    cACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCg                                                                             <  1:205254/90‑1 (MQ=255)
     aCGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGa                                                                            >  2:8214/1‑90 (MQ=255)
       gCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACa                                                                          <  1:8214/90‑1 (MQ=255)
          cAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGAt                                                                       <  1:282238/90‑1 (MQ=255)
              ggCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTc                                                                   <  1:263749/90‑1 (MQ=255)
               gCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTcc                                                                  >  1:329896/1‑90 (MQ=255)
                               aTTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATc                                                  <  1:173728/90‑1 (MQ=255)
                                   ttCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCgg                                              <  1:221340/90‑1 (MQ=255)
                                              gTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATtgt                                   >  2:213281/1‑90 (MQ=255)
                                                ttCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATtgtgt                                 <  2:193851/90‑1 (MQ=255)
                                                        gggAAAAACGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACct                         >  2:24620/1‑90 (MQ=255)
                                                                cGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTc                 <  1:139223/90‑1 (MQ=255)
                                                                cGGTTTTTCCATCCACCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTc                 <  1:93646/90‑1 (MQ=255)
                                                                        ccatccaCCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGaa         <  1:117983/90‑1 (MQ=255)
                                                                          atccaCCACTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAAc       <  1:361581/90‑1 (MQ=255)
                                                                               ccacTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTc  <  1:353475/90‑1 (MQ=255)
                                                                               ccacTGTGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTc  <  1:213281/90‑1 (MQ=255)
                                                                                      |                                                                                   
TTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTC  >  NZ_CP009273/1634064‑1634233

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