Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,199,864 G→A pseudogene (239/246 nt)
pseudogene (8/1137 nt)
BW25113_RS26105 →
BW25113_RS05965 →
helix‑turn‑helix domain‑containing protein
terminase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,199,8640GA100.0% 133.5 / NA 34pseudogene (239/246 nt)
pseudogene (8/1137 nt)
BW25113_RS26105
BW25113_RS05965
helix‑turn‑helix domain‑containing protein
terminase large subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/22);  total (12/22)

GAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCC  >  NZ_CP009273/1199731‑1199988
                                                                                                                                     |                                                                                                                            
gAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCTTCGCCAGGGTGATGATtcat                                                                                                                         <  2:222361/139‑1 (MQ=255)
gAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATtcat                                                                                                                         <  1:115385/139‑1 (MQ=255)
              gCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACtgtg                                                                                                             <  2:510607/137‑1 (MQ=255)
              gCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGcc                                                                                                           >  1:42491/1‑139 (MQ=255)
                  tCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGTCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCgtgg                                                                                                       <  2:736471/139‑1 (MQ=255)
                                   ccGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAg                                                                                      <  1:389290/139‑1 (MQ=255)
                                   ccGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAg                                                                                      <  1:623724/139‑1 (MQ=255)
                                     gTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCa                                                                                    <  1:353459/139‑1 (MQ=255)
                                               cgcgGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGAt                                                                          <  1:426946/139‑1 (MQ=255)
                                               cgcgGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGAt                                                                          <  2:42491/139‑1 (MQ=255)
                                               cgcgGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGAt                                                                          <  2:319889/139‑1 (MQ=255)
                                                             ccGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAg                                                                                      >  2:668069/1‑113 (MQ=255)
                                                             ccGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAg                                                                                      <  1:668069/113‑1 (MQ=255)
                                                              cGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTa                                                                  >  2:104465/1‑132 (MQ=255)
                                                              cGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTa                                                                  <  1:104465/132‑1 (MQ=255)
                                                                tCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGAt                                                         <  2:529618/139‑1 (MQ=255)
                                                                      cGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTAc                                                   <  1:790807/139‑1 (MQ=255)
                                                                                  ggCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGAtt                                                                                                                            >  1:413853/1‑54 (MQ=255)
                                                                                  ggCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGAtt                                                                                                                            <  2:413853/54‑1 (MQ=255)
                                                                                              aaGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACtgtg                                                                                                             <  2:43839/57‑1 (MQ=255)
                                                                                              aaGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACtgtg                                                                                                             >  1:43839/1‑57 (MQ=255)
                                                                                              aaGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACgccagc                            >  1:382370/1‑138 (MQ=255)
                                                                                              aaGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACgccagc                            >  2:21163/1‑138 (MQ=255)
                                                                                                     gTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACgccagcca                          <  1:674086/133‑1 (MQ=255)
                                                                                                     gTTAACGCGTAAGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACgccagcca                          >  2:674086/1‑133 (MQ=255)
                                                                                                                aGGCGCGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTac         >  1:637860/1‑139 (MQ=255)
                                                                                                                  gcgcGTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTacca       <  1:351640/139‑1 (MQ=255)
                                                                                                                     cgTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGa                                                                                                   >  2:274847/1‑44 (MQ=255)
                                                                                                                     cgTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGa                                                                                                   <  1:274847/44‑1 (MQ=255)
                                                                                                                     cgTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATg                                                        >  2:504454/1‑87 (MQ=255)
                                                                                                                     cgTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATg                                                        <  1:504454/87‑1 (MQ=255)
                                                                                                                     cgTCGCCAGGGTGATGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCAccg    <  1:471025/139‑1 (MQ=255)
                                                                                                                                  aTGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCAccgcc  >  2:756118/1‑128 (MQ=255)
                                                                                                                                  aTGATTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCAccgcc  <  1:756118/128‑1 (MQ=255)
                                                                                                                                     |                                                                                                                            
GAAGGTGTGTGCTGGCCATCAATTGAAACCATTGCCCGTCAGATTGGCGCGGGGATGAGTACCGTCAGAACGGCTATCGCACGGCTGGAAGCAGAAGGCTGGTTAACGCGTAAGGCGCGTCGCCAGGGTGATGGTTCATCACCCCACTGTGCCGTGGTGGATGAATATCACGAGCACGCCACAGATGCGCTTTACACCACGATGCTTACCGGGATGGGGGCGCGACGCCAGCCACTGATGTGGGCCATTACCACCGCC  >  NZ_CP009273/1199731‑1199988

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