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% 96.9 / NA 27pseudogene (1184/1374 nt)BW25113_RS25705DUF5507 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/14);  total (13/14)

TTCTCAACGGGATCATCCTCAGAGGAACTATCATCAAAGTCATCATCCATAAATAAAATATCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTTCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACGTTCCCCAGTTCACATATTGAATATCATAGTTACTGAAGTTCT  >  NZ_CP009273/2777052‑2777316
                                                                                                                                       |                                                                                                                                 
ttCTCAACGGGATCATCCTCAGAGGAACTATCATCAAAGTCATCATCCATAAATAAAATATCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGc                                                                                                                                >  1:119976/1‑139 (MQ=255)
                              tcatcaAAGTCATCATCCATAAATAAAATATCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTActc                                                                                                   <  2:185387/138‑1 (MQ=255)
                                atcaAAGTCATCATCCATAAATAAAATATCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCttg                                                                                                >  2:338407/1‑139 (MQ=255)
                                       tcatcaTCCATAAATAAAATATCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAg                                                                                         <  2:174670/139‑1 (MQ=255)
                                                        aataTCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGAt                                                                                                               <  1:117559/100‑1 (MQ=255)
                                                        aataTCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGAt                                                                                                               >  2:117559/1‑100 (MQ=255)
                                                        aataTCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATa                                                                        >  1:168478/1‑139 (MQ=255)
                                                          tatCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATaca                                                                      <  1:366027/139‑1 (MQ=255)
                                                               tcGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTcac                                                                                                                            <  2:407641/80‑1 (MQ=255)
                                                               tcGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTcac                                                                                                                            >  1:407641/1‑80 (MQ=255)
                                                               tcGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAg                                                                 >  2:113024/1‑139 (MQ=255)
                                                                       tGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACa                                                         <  2:192479/139‑1 (MQ=255)
                                                                         ccACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATc                                                       <  1:272605/139‑1 (MQ=255)
                                                                         ccACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTTTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATc                                                       >  1:332214/1‑138 (MQ=255)
                                                                            cGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGGATCAACATCaaa                                                    >  2:40802/1‑139 (MQ=255)
                                                                                 gTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCAt                                                                                                        <  1:357028/82‑1 (MQ=255)
                                                                                 gTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCAt                                                                                                        >  2:357028/1‑82 (MQ=255)
                                                                                 gTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAAt                                                >  1:317715/1‑138 (MQ=255)
                                                                                        gTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACGTTcc                                        <  2:317715/139‑1 (MQ=255)
                                                                                                   ttattaCGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGGATCAACATCAAACAATGACGTTCCCCAGTTcaca                              >  1:269474/1‑138 (MQ=255)
                                                                                                      ttaCGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACGTTCCCCAGTTCACATAtt                          <  1:113024/139‑1 (MQ=255)
                                                                                                      ttaCGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACGTTCCCCAGTTCACATAtt                          <  2:467515/139‑1 (MQ=255)
                                                                                                                  aaTACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACGTTCCCCAGTTCACATATTGAATATCATAGt              <  2:82026/139‑1 (MQ=255)
                                                                                                                  aaTACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGGATCAACATCAAACAATGACGTTCCCCAGTTCACATATTGAATATCATAGt              <  2:102593/139‑1 (MQ=255)
                                                                                                                          aCTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGGATCAACATCAAACAATGACGTTCCCCAGTTCACata                            <  1:398566/117‑1 (MQ=255)
                                                                                                                          aCTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGGATCAACATCAAACAATGACGTTCCCCAGTTCACata                            >  2:398566/1‑117 (MQ=255)
                                                                                                                              aaaCCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACGTTCCCCAGTTCACATATTGAATATCATAGTTACTGAATTTCt  >  2:406122/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                 
TTCTCAACGGGATCATCCTCAGAGGAACTATCATCAAAGTCATCATCCATAAATAAAATATCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTTCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACGTTCCCCAGTTCACATATTGAATATCATAGTTACTGAAGTTCT  >  NZ_CP009273/2777052‑2777316

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: