Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,777,187 T→C pseudogene (1184/1374 nt) BW25113_RS25705 ← DUF5507 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,777,1870TC100.0% 90.5 / NA 25pseudogene (1184/1374 nt)BW25113_RS25705DUF5507 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/13);  total (12/13)

TCCTCAGAGGAACTATCATCAAAGTCATCATCCATAAATAAAATATCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTTCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACGTTCCCCAGTTCACATATTGAATATCATAGTTACTGAAGTTCTGT  >  NZ_CP009273/2777067‑2777318
                                                                                                                        |                                                                                                                                   
tcCTCAGAGGAACTATCATCAAAGTCATCATCCATAAATAAAATATCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAAt                                                                                                                    >  1:101705/1‑138 (MQ=255)
     agagGAACTATCATCAAAGTCATCATCCATAAATAAAATATCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGAtttt                                                                                                              >  1:13870/1‑139 (MQ=255)
     agagGAACTATCATCAAAGTCATCATCCATAAATAAAATATCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGAtt                                                                                                                <  2:101705/137‑1 (MQ=255)
           aCTATCATCAAAGTCATCATCCATAAATAAAATATCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTa                                                                                                        >  2:468970/1‑139 (MQ=255)
           aCTATCATCAAAGTCATCATCCATAAATAAAATATCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTa                                                                                                        >  1:304868/1‑139 (MQ=255)
               tcatcaAAGTCATCATCCATAAATAAAATATCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTAct                                                                                                      <  1:426996/137‑1 (MQ=255)
                    aaaGTCATCATCCATAAATAAAATATCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCttgtt                                                                                                <  2:304868/138‑1 (MQ=255)
                                         aataTCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATa                                                                          >  2:3194/1‑139 (MQ=255)
                                           tatCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATaca                                                                        <  1:408728/139‑1 (MQ=255)
                                               atcGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACgg                                                                     <  1:443477/138‑1 (MQ=255)
                                                        tGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACa                                                           <  1:3194/139‑1 (MQ=255)
                                                          ccACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACAt                                                          <  1:468970/138‑1 (MQ=255)
                                                             cGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACgg                                                                     <  2:544656/124‑1 (MQ=255)
                                                             cGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACgg                                                                     >  1:544656/1‑124 (MQ=255)
                                                                  gTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAAt                                                  >  2:9046/1‑138 (MQ=255)
                                                                  gTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAAt                                                  >  1:241840/1‑138 (MQ=255)
                                                                                   attattaCGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACGTTCCCCAGTTcac                                 <  2:241840/138‑1 (MQ=255)
                                                                                      attaCGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACGTTCCCCAGTTCACata                              <  2:13870/138‑1 (MQ=255)
                                                                                      attaCGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACGTTCCCCAGTTCACata                              <  2:385438/138‑1 (MQ=255)
                                                                                                           aCTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTActc                                                                                                     <  2:531775/46‑1 (MQ=255)
                                                                                                           aCTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTActc                                                                                                     >  1:531775/1‑46 (MQ=255)
                                                                                                               aaaCCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACg                                              <  2:127337/97‑1 (MQ=255)
                                                                                                               aaaCCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACg                                              >  1:127337/1‑97 (MQ=255)
                                                                                                               aaaCCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACGTTCCCCAGTTCACATATTGAATATCATAGTTACTGAAGTTCt    >  1:241986/1‑139 (MQ=255)
                                                                                                                 aCCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACGTTCCCCAGTTCACATATTGAATATCATAGTTACTGAAGTTCTGt  >  1:155786/1‑139 (MQ=255)
                                                                                                                        |                                                                                                                                   
TCCTCAGAGGAACTATCATCAAAGTCATCATCCATAAATAAAATATCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTTCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACGTTCCCCAGTTCACATATTGAATATCATAGTTACTGAAGTTCTGT  >  NZ_CP009273/2777067‑2777318

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