Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,632,437 (T)7→6 intergenic (‑71/+8) BW25113_RS08135 ← / ← BW25113_RS25620 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,632,4310T.100.0% 71.2 / NA 18intergenic (‑65/+14)BW25113_RS08135/BW25113_RS25620hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (5/13);  total (5/13)

GCCATTTGGATATTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAAAAAGCTCCCGAAGGAGCCTTCATTTTCAC  >  NZ_CP009273/1632346‑1632516
                                                                                     |                                                                                     
gCCATTTGGATATTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGAttttt                                                                                   <  1:331558/90‑1 (MQ=255)
  cATTTGGATATTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGAttttttc                                                                                 <  1:281193/90‑2 (MQ=255)
       ggATATTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCatgat                                                                           <  2:471146/90‑1 (MQ=255)
               aCGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCgt                                                                   <  1:258608/90‑1 (MQ=255)
               aCGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCgt                                                                   <  1:569568/90‑1 (MQ=255)
               aCGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCgt                                                                   <  2:285991/90‑1 (MQ=255)
               aCGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCgt                                                                   <  2:347293/90‑1 (MQ=255)
                    tgatgGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTAt                                                              <  1:126377/90‑1 (MQ=255)
                             tAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATc                                                     <  1:467790/90‑1 (MQ=255)
                                 tCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATc                                                 <  1:16071/90‑1 (MQ=255)
                                 tCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATc                                                 <  1:53707/90‑1 (MQ=255)
                                         tCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTAtgt                                                            <  1:356970/71‑1 (MQ=255)
                                         tCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTAtgt                                                            >  2:356970/1‑71 (MQ=255)
                                          cGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTa                                                                       <  1:219432/59‑1 (MQ=255)
                                          cGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTa                                                                       >  2:219432/1‑59 (MQ=255)
                                                          gctgGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAAAAAGCTc                        >  2:488545/1‑90 (MQ=255)
                                                          gctgGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAAAAAGCTc                        >  1:537130/1‑90 (MQ=255)
                                                            tgGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAAAAAGCTccc                      >  2:34143/1‑90 (MQ=255)
                                                            tgGTTCAGAGTTAAAAAAACTATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAAAAAGCTccc                      <  1:478636/90‑1 (MQ=255)
                                                                                tATGATTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAAAAAGCTCCCGAAGGAGCCTTCATTTTCAc  <  2:125848/90‑1 (MQ=255)
                                                                                     |                                                                                     
GCCATTTGGATATTCACGCATGATGGTGTTAACTCCAGTCATCGCTGGCACTACCACTGCTGGTTCAGAGTTAAAAAAACTATGATTTTTTTCATGATGTTACCGTAGTATGTGAGTATCCATCGAATAGACACCAAGCAAAAAAGCTCCCGAAGGAGCCTTCATTTTCAC  >  NZ_CP009273/1632346‑1632516

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: