Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,498,429 T→G R2S (AGA→AGC BW25113_RS22240 ← IS30‑like element IS30 family transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,498,4290TG100.0% 98.5 / NA 31R2S (AGA→AGCBW25113_RS22240IS30‑like element IS30 family transposase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (17/14);  total (17/14)

CACTGAAGCCTGTTCCGTTCTTCCATAGTTCAAAAACAGAGGCTTTTTCCTCTGCTGTAAATGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTT  >  NZ_CP009273/4498362‑4498563
                                                                   |                                                                                                                                      
cACTGAAGCCTGTTCCGTTCTTCCATAGTTCAAAAACAGAGGCTTTTTCCTCTGCTGTAAATGTTCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTg                                                                 >  1:496269/1‑139 (MQ=14)
  cTGAAGCCTGTTCCGTTCTTCCATAGTTCAAAAACAGAGGCTTTTTCCTCTGCTGTAAATGTTCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCt                                                               <  2:167803/139‑1 (MQ=18)
     aaGCCTGTTCCGTTCTTCCATAGTTCAAAAACAGAGGCTTTTTCCTCTGCTGTAAATGTTCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCtttt                                                            >  1:232351/1‑139 (MQ=21)
      aGCCTGTTCCGTTCTTCCATAGTTCAAAAACAGAGGCTTTTTCCTCTGCTGTAAATGTTCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCttttt                                                           <  2:735394/139‑1 (MQ=21)
            ttCCGTTCTTCCATAGTTCAAAAACAGAGGCTTTTTCCTCTGCTGTAAATGTTCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAata                                                     >  1:633663/1‑139 (MQ=21)
            ttCCGTTCTTCCATAGTTCAAAAACAGAGGCTTTTTCCTCTGCTGTAAATGTTCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAat                                                      <  1:444723/138‑1 (MQ=21)
               cGTTCTTCCATAGTTCAAAAACAGAGGCTTTTTCCTCTGCTGTAAATGTTCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAatatat                                                  <  2:56551/139‑1 (MQ=37)
                 ttcttcCATAGTTCAAAAACAGAGGCTTTTTCCTCTGCTGTAAATGTTCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATAtct                                                >  1:715940/1‑139 (MQ=38)
                 ttcttcCATAGTTCAAAAACAGAGGCTTTTTCCTCTGCTGTAAATGTTCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATAtct                                                <  1:110535/139‑1 (MQ=38)
                 ttcttcCATAGTTCAAAAACAGAGGCTTTTTCCTCTGCTGTAAATGTTCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATAtct                                                >  2:134182/1‑139 (MQ=38)
                 ttcttcCATAGTTCAAAAACAGAGGCTTTTTCCTCTGCTGTAAATGTTCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATAtct                                                >  1:479526/1‑139 (MQ=38)
                 ttcttcCATAGTTCAAAAACAGAGGCTTTTTCCTCTGCTGTAAATGTTCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATAtct                                                >  1:231651/1‑139 (MQ=38)
                 ttcttcCATAGTTCAAAAACAGAGGCTTTTTCCTCTGCTGTAAATGTTCGGCTCATTCAAAAAACCTCCGCAACCCCATGTCTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATAtct                                                <  2:543646/139‑1 (MQ=255)
                      ccATAGTTCAAAAACAGAGGCTTTTTCCTCTGCTGTAAATGTTCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAat                                                      >  1:331658/1‑128 (MQ=25)
                      ccATAGTTCAAAAACAGAGGCTTTTTCCTCTGCTGTAAATGTTCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAat                                                      <  2:331658/128‑1 (MQ=25)
                      ccATAGTTCAAAAACAGAGGCTTTTTCCTCTGCTGTAAATGTTCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCt                                            >  2:11840/1‑138 (MQ=255)
                               aaaaaCAGAGGCTTTTTCCTCTGCTGTAAATGTTCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTATATATCTCGCTCAAGGCGAGCt                                 <  1:659052/139‑1 (MQ=255)
                                aaaaCAGAGGCTTTTTCCTCTGCTGTAAATGTTCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCt                                 >  1:152238/1‑139 (MQ=255)
                                      gagGCTTTTTCCTCTGCTGTAAATGTTCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCAtt                            >  1:256864/1‑138 (MQ=255)
                                      gagGCTTTTTCCTCTGCTGTAAATGTTCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCAtt                            >  2:429220/1‑138 (MQ=255)
                                       agGCTTTTTCCTCTGCTGTAAATGTTCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTa                          <  1:496610/139‑1 (MQ=255)
                                             tttCCTCTGCTGTAAATGTTCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGcc                     <  1:134182/138‑1 (MQ=255)
                                                    tgctgTAAATGTTCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGc               <  2:13961/137‑1 (MQ=255)
                                                     gctgTAAATGTTCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTa                          >  1:612467/1‑125 (MQ=255)
                                                     gctgTAAATGTTCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTa                          <  2:612467/125‑1 (MQ=255)
                                                                tCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAGTTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGcc                     <  2:124548/119‑1 (MQ=255)
                                                                tCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGcc                     >  1:124548/1‑119 (MQ=255)
                                                                tCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTg          <  1:72011/130‑1 (MQ=255)
                                                                tCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTg          >  2:72011/1‑130 (MQ=255)
                                                                tCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAAttt  >  2:527983/1‑138 (MQ=255)
                                                                tCGGCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAAttt  >  1:114796/1‑138 (MQ=255)
                                                                   |                                                                                                                                      
CACTGAAGCCTGTTCCGTTCTTCCATAGTTCAAAAACAGAGGCTTTTTCCTCTGCTGTAAATGTTCGTCTCATTCAAAAAACCTCCGCAACCCCATGTTTTCACATAACTGTTGCGTTGACCAATTGAATCTACAGTTGCTTTTTTTAATATATCTCGCTCAAGGCGAGCTTCATTTAACGCCTTACGCAGTTGCAGAATTT  >  NZ_CP009273/4498362‑4498563

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: