Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 710,687 T→C intergenic (+33/‑181) pgm → / → BW25113_RS03560 phosphoglucomutase (alpha‑D‑glucose‑1,6‑bisphosphate‑dependent)/lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273710,6870TC100.0% 66.2 / NA 20intergenic (+33/‑181)pgm/BW25113_RS03560phosphoglucomutase (alpha‑D‑glucose‑1,6‑bisphosphate‑dependent)/lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (7/13);  total (7/13)

GCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAA  >  NZ_CP009273/710600‑710761
                                                                                       |                                                                          
gCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGACCg                                                                          <  1:324867/90‑1 (MQ=255)
   gATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGACCGccc                                                                       <  1:201517/90‑1 (MQ=255)
                 cGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGACCGCCCTTTTTTACGTATGa                                                         >  2:439711/1‑90 (MQ=255)
                           tGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGACCGCCCTTTTTTACGTATGACAAACACAGa                                               >  1:695636/1‑90 (MQ=255)
                              tAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGACCGCCCTTTTTTACGTATGACAAACACAGAAtt                                            <  1:232781/90‑1 (MQ=255)
                                            aaaaCGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGACCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTAc                              <  2:725483/90‑1 (MQ=255)
                                            aaaaCGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGACCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTAc                              <  2:120759/90‑1 (MQ=255)
                                            aaaaCGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGACCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTAc                              <  2:147664/90‑1 (MQ=255)
                                               aCGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGACCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGa                                      >  1:287500/1‑79 (MQ=255)
                                               aCGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGACCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGa                                      <  2:287500/79‑1 (MQ=255)
                                               aCGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGACCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCt                           >  1:417058/1‑90 (MQ=255)
                                                   gTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGACCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATc                       <  1:297943/90‑1 (MQ=255)
                                                     aaaCACATTTAATAAAAAAAGGGCGGTCGCAAGACCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAg                     <  1:1066134/90‑1 (MQ=255)
                                                                 taaAAAAAGGGCGGTCGCAAGACCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGAt         <  1:96637/90‑1 (MQ=255)
                                                                  aaaaaaaGGGCGGTCGCAAGACCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCt                 >  2:783977/1‑81 (MQ=255)
                                                                  aaaaaaaGGGCGGTCGCAAGACCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCt                 <  1:783977/81‑1 (MQ=255)
                                                                  aaaaaaaGGGCGGTCGCAAGACCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATg        >  2:501954/1‑90 (MQ=255)
                                                                    aaaaaGGGCGGTCGCAAGACCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGt      >  1:831149/1‑90 (MQ=255)
                                                                        aGGGCGGTCGCAAGACCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCaa  <  2:775506/90‑1 (MQ=255)
                                                                        aGGGCGGTCGCAAGACCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCaa  <  1:439711/90‑1 (MQ=255)
                                                                                       |                                                                          
GCAGATTGAGAAAGAAGCGGTTGAGATTGTTAGCGAAGTTCTGAAAAACGCGTAAACACATTTAATAAAAAAAGGGCGGTCGCAAGATCGCCCTTTTTTACGTATGACAAACACAGAATTGCCTGATGCGCTACGCTTATCAGGCCTACGAGGATGGTGCAA  >  NZ_CP009273/710600‑710761

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