Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,201,626 A→G pseudogene (34/786 nt) BW25113_RS05980 → baseplate J/gp47 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,201,6260AG100.0% 33.4 / NA 10pseudogene (34/786 nt)BW25113_RS05980baseplate J/gp47 family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (5/5);  total (5/5)

CTGTCTGTGCCTGAGGGGAAACTTTTACGCCTACAAAGTGAAAGCATTTGGCGAAGTGGCTGAACTGCTGCCCGTCGATCCCGGCTGTGTGGTATATGCGCTGGGAAGGTGTCAGCGATGGCCTGAAGGTGACCGCCGGGAGTGTTATTCAGCGCGATGACCTGGTGCAGTACACGACAACTGACGATGCAACCAGCTCCGGTGGTGTCCTGCGCGTGCCGATCGCCTGCTCAAGTGCAGGTGCGGTCGGTA  >  NZ_CP009273/1201501‑1201752
                                                                                                                             |                                                                                                                              
ctgtctgtGCCTGAGGGGAAACTTTTACGCCTACAAAGTGAAAGCATTTGGCGAAGTGGCTGAACTGCTGCCCGTCGATCCCGGCTGTGTGGTATATGCGCTGGGAAGGTGTCAGCGATGGCCTGGAGGTGACCGCCgg                                                                                                                   <  1:265359/139‑1 (MQ=255)
   tctgtGCCTGAGGGGAAACTTTTACGCCTACAAAGTGAAAGCATTTGGCGAAGTGGCTGAACTGCTGCCCGTCGATCCCGGCTGTGTGGTATATGCGCTGGGAAGGTGTCAGCGATGGCCTGGAGGTGACCGCCGGGAg                                                                                                                >  2:256124/1‑139 (MQ=255)
       tgCCTGAGGGGAAACTTTTACGCCTACAAAGTGAAAGCATTTGGCGAAGTGGCTGAACTGCTGCCCGTCGATCCCGGCTGTGTGGTATATGCGCTGGGAAGGTGTCAGCGATGGCCTGGAGGTGACCGCCGGGAgtgt                                                                                                             >  2:87331/1‑138 (MQ=255)
                 gAAACTTTTACGCCTACAAAGTGAAAGCATTTGGCGAAGTGGCTGAACTGCTGCCCGTCGATCCCGGCTGTGTGGTATATGCGCTGGGAAGGTGTCAGCGATGGCCTGGAGGTGACCGCCGGGAGTGTTATTCAgcgc                                                                                                   >  2:425993/1‑138 (MQ=255)
                 gAAACTTTTACGCCTACAAAGTGAAAGCATTTGGCGAAGTGGCTGAACTGCTGCCCGTCGATCCCGGCTGTGTGGTATATGCGCTGGGAAGGTGTCAGCGATGGCCTGGAGGTGACCGCCGGGAGTGTTATTCAGcgcg                                                                                                  <  2:134650/139‑1 (MQ=255)
                        ttACGCCTACAAAGTGAAAGCATTTGGCGAAGTGGCTGAACTGCTGCCCGTCGATCCCGGCTGTGTGGTATATGCGCTGGGAAGGTGTCAGCGATGGCCTGGAGGTGACCGCCGGGAGTGTTATTCAGCGCGATGAcc                                                                                            <  1:247472/138‑1 (MQ=255)
                                                                             aTCCCGGCTGTGTGGTATATGCGCTGGGAAGGTGTCAGCGATGGCCTGGAGGTGACCGCCGGGAGTGTTATTCAGCGCGATGACCTGGTGCAGTACACGACAACTGACGATGCAACCAGCTCCGGTGGTGTCCTGcgcg                                      <  1:109533/139‑1 (MQ=255)
                                                                                                 gcgcTGGGAAGGTGTCAGCGATGGCCTGGAGGTGACCGCCGGGAGTGTTATTCAGCGCGATGACCTGGTGCAGTACACGACAACTGACGATGCAACCAGCTCCGGTGGTGTCCTGCGCGTGCCGATCGCCTGCTCAAGt                  >  1:362068/1‑139 (MQ=255)
                                                                                                      gggAAGGTGTCAGCGATGGCCTGGAGGTGACCGCCGGGAGTGTTATTCAGCGCGATGACCTGGTGCAGTACACGACAACTGACGATGCAACCAGCTCCGGTGGTGTCCTGCGCGTGCCGATCGCCTGCTCAAGTGCAgg             >  2:360941/1‑139 (MQ=255)
                                                                                                                 aGCGATGGCCTGGAGGTGACCGCCGGGAGTGTTATTCAGCGCGATGACCTGGTGCAGTACACGACAACTGACGATGCAACCAGCTCCGGTGGTGTCCTGCGCGTGCCGATCGCCTGCTCAAGTGCAGGTGCGGTCGGTa  <  1:87331/139‑1 (MQ=255)
                                                                                                                             |                                                                                                                              
CTGTCTGTGCCTGAGGGGAAACTTTTACGCCTACAAAGTGAAAGCATTTGGCGAAGTGGCTGAACTGCTGCCCGTCGATCCCGGCTGTGTGGTATATGCGCTGGGAAGGTGTCAGCGATGGCCTGAAGGTGACCGCCGGGAGTGTTATTCAGCGCGATGACCTGGTGCAGTACACGACAACTGACGATGCAACCAGCTCCGGTGGTGTCCTGCGCGTGCCGATCGCCTGCTCAAGTGCAGGTGCGGTCGGTA  >  NZ_CP009273/1201501‑1201752

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