Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,251,021 A→G pseudogene (773/986 nt) BW25113_RS21020 → DUF2713 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,251,0210AG100.0% 88.1 / NA 25pseudogene (773/986 nt)BW25113_RS21020DUF2713 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/11);  total (14/11)

ACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACC  >  NZ_CP009273/4250891‑4251151
                                                                                                                                  |                                                                                                                                  
aCGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGaa                                                                                                                            >  2:51259/1‑139 (MQ=255)
             aGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTAt                                                                                                               <  2:36434/139‑1 (MQ=255)
                     tgtATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGAt                                                                                                         >  2:25640/1‑137 (MQ=255)
                               tGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGAt                                                                                                         <  1:177269/127‑1 (MQ=255)
                               tGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGAt                                                                                                         >  2:177269/1‑127 (MQ=255)
                                   aataatGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATgg                                                                                                                   <  2:329491/113‑1 (MQ=255)
                                   aataatGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATgg                                                                                                                   >  1:329491/1‑113 (MQ=255)
                                   aataatGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTgg                                                                                         >  2:105805/1‑139 (MQ=255)
                                            agCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAg                                                                                >  1:375737/1‑139 (MQ=255)
                                                 ccTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGt                                                                                                  >  2:7894/1‑116 (MQ=255)
                                                 ccTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGt                                                                                                  <  1:7894/116‑1 (MQ=255)
                                                    gAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTc                                                                        >  2:126304/1‑139 (MQ=255)
                                                           cAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTg                                                                                          >  2:302696/1‑114 (MQ=255)
                                                           cAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTg                                                                                          <  1:302696/114‑1 (MQ=255)
                                                           cAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCt                                                                  >  1:420941/1‑138 (MQ=255)
                                                                       ttATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCa                                                     <  2:317522/139‑1 (MQ=255)
                                                                              aaCGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCg                                              <  1:105805/139‑1 (MQ=255)
                                                                               aCGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCga                                             <  1:316511/139‑1 (MQ=255)
                                                                                 gTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCgaga                                           <  2:421461/139‑1 (MQ=255)
                                                                                              gTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaa                              >  1:243841/1‑139 (MQ=255)
                                                                                                           aTGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTcg                 >  1:62282/1‑139 (MQ=255)
                                                                                                           aTGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTcg                 >  2:123917/1‑139 (MQ=255)
                                                                                                                     tGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATTTGGt                                                                                        >  1:407387/1‑58 (MQ=255)
                                                                                                                     tGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGt                                                                                        <  2:407387/58‑1 (MQ=255)
                                                                                                                           aTGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGAcc  <  2:62282/138‑1 (MQ=255)
                                                                                                                              ggggggAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGtt       <  2:231273/130‑1 (MQ=255)
                                                                                                                              ggggggAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGtt       >  1:231273/1‑130 (MQ=255)
                                                                                                                                  |                                                                                                                                  
ACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACC  >  NZ_CP009273/4250891‑4251151

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