Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,325,369 T→G pseudogene (2423/4515 nt) BW25113_RS11635 ← alpha‑2‑macroglobulin family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,325,3690TG100.0% 87.6 / NA 26pseudogene (2423/4515 nt)BW25113_RS11635alpha‑2‑macroglobulin family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (13/13);  total (13/13)

GTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGATTTCTGGTTGCAGCGCGTAGATCATTTCATCGACCACGCCTACCGTTAGC  >  NZ_CP009273/2325235‑2325503
                                                                                                                                      |                                                                                                                                      
gTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGGCg                                                                                                                                      <  1:320061/137‑1 (MQ=255)
      cGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGGCGTAGCTGa                                                                                                                               >  2:75927/1‑138 (MQ=255)
       ggTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGGCGTAGCTGAt                                                                                                                              <  2:397411/138‑1 (MQ=255)
            tgagtgaCGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGGCGTAGCTGATAAACGa                                                                                                                        <  2:120245/139‑1 (MQ=255)
            tgagtgaCGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGGCGTAGCTGATAAACGa                                                                                                                        <  1:75927/139‑1 (MQ=255)
                   cGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGGCGTAGCTGATAAACGACAGACTg                                                                                                                 >  1:345503/1‑139 (MQ=255)
                   cGGCATCCATGCCGCGGGATCCACCTCTTCACGCCGTGGACGGTCAAGCATTTTTACTCGCCGCTCGCTGCGGGTAGTTGCGCCATGCGCACCCGGCTCGCTTAAGAGCGCCTGGGCGTAGCTGATAAACGACAGACTg                                                                                                                 >  2:21658/1‑139 (MQ=255)
                    ggCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGGCGTAGCTGATAAACGACAGACTgg                                                                                                                <  1:539756/139‑1 (MQ=255)
                              ccGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGGCGTAGCTGATAAACGACAGACTGGAGCTGGTACg                                                                                                      >  1:169052/1‑139 (MQ=255)
                                  ggTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGGCGTAGCTGATAAACGACAGACTGGAGCTGGTACGcaca                                                                                                  <  2:441241/139‑1 (MQ=255)
                                  ggTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGGCGTAGCTGATAAACGACAGACTGGAGCTGGTACGcaca                                                                                                  <  2:345503/139‑1 (MQ=255)
                                           ctctTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGAGCCTGGGCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGc                                                                                         >  2:95451/1‑139 (MQ=255)
                                              ttCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGGCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGcccc                                                                                      >  2:560842/1‑139 (MQ=255)
                                                       tGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGGCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTa                                                                                            <  1:465760/124‑1 (MQ=255)
                                                       tGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGGCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTa                                                                                            >  2:465760/1‑124 (MQ=255)
                                                        ggACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGGCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAg                                                                             >  1:298104/1‑138 (MQ=255)
                                                                              gccgcTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGGCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGgcgc                                                      >  2:180899/1‑139 (MQ=255)
                                                                                cgctcgctGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGGCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGa                                                    >  1:582444/1‑139 (MQ=255)
                                                                                                          cgcAACCGGCTCGCTGGAGAGCGCCTGGGCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGATTTCTGGTTGCAGCGCGTAGATCAtt                          >  1:429869/1‑139 (MQ=255)
                                                                                                             aaCCGGCTCGCTGGAGAGCGCCTGGGCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGATTTCTGGTTGCAGCGCGTAGATCATTTCa                       <  1:351417/139‑1 (MQ=255)
                                                                                                                  gttCGCTGGAGAGCGCCTGGGCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGATTTCTGGTTGCAGCGCGTAGATCATTTCATCGAc                  <  2:293074/137‑1 (MQ=255)
                                                                                                                  gCTCGCTGGAGAGCGCCTGGGCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGATTTCTGGTTGCAGCGCGTAGATCATTTCATCGAc                  >  2:101124/1‑139 (MQ=255)
                                                                                                                       cTGGAGAGCGCCTGGGCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGATTTCTGGTTGCAGCGCGTAGATCATTTCATCGACCACGc             <  1:101124/139‑1 (MQ=255)
                                                                                                                               cgcCTGGGCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGATTTCTGGTTGCa                                        >  2:32341/1‑104 (MQ=255)
                                                                                                                               cgcCTGGGCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGATTTCTGGTTGCa                                        <  1:32341/104‑1 (MQ=255)
                                                                                                                                  cTGGGCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGATTTCTGGTTGCAGCGCGTAGATCATTTCATCGACCACGCCTACCGTTAGc  <  2:348986/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                      
GTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCAGGCGCAACCGGCTCGCTGGAGAGCGCCTGGTCGTAGCTGATAAACGACAGACTGGAGCTGGTACGCACATTGTTACGCCCCAGCGGATAGAAAAATTTGCCGATATTCGGCGCGATTTCTGGTTGCAGCGCGTAGATCATTTCATCGACCACGCCTACCGTTAGC  >  NZ_CP009273/2325235‑2325503

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