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

TTCTCAACGGGATCATCCTCAGAGGAACTATCATCAAAGTCATCATCCATAAATAAAATATCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTTCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACGTTCCCCAGTTCACATATTGAATATCATAGTTACTGAAGTTCTGTC  >  NZ_CP009273/2777052‑2777319
                                                                                                                                       |                                                                                                                                    
ttCTCAACGGGATCATCCTCAGAGGAACTATCATCAAAGTCATCATCCATAAATAAAATATCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGc                                                                                                                                   >  2:162594/1‑139 (MQ=255)
          gATCATCCTCAGAGGAACTATCATCAAAGTCATCATCCATAAATAAAATATCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGttt                                                                                                                         >  2:219366/1‑139 (MQ=255)
                   cAGAGGAACTATCATCAAAGTCATCATCCATAGATAAAATATCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGAttt                                                                                                                >  1:662147/1‑139 (MQ=255)
                        gAACTATCATCAAAGTCATCATCCATAAATAAAATATCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCAt                                                                                                           <  1:308075/139‑1 (MQ=255)
                                                          tatCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATaca                                                                         <  2:122729/139‑1 (MQ=255)
                                                            tcatcGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCTCTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATacac                                                                        <  1:707984/138‑1 (MQ=255)
                                                              atcGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACgg                                                                      <  1:162594/138‑1 (MQ=255)
                                                              atcGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACgg                                                                      <  2:662147/138‑1 (MQ=255)
                                                                                 gTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGgc                                                                                       >  2:613701/1‑102 (MQ=255)
                                                                                 gTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGgc                                                                                       <  1:613701/102‑1 (MQ=255)
                                                                                 gTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAAt                                                   >  1:446668/1‑138 (MQ=255)
                                                                                 gTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATg                                                  <  2:446668/139‑1 (MQ=255)
                                                                                       aGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACGtt                                             <  1:219366/138‑1 (MQ=255)
                                                                                                           aTCAGTCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACGTTCCCCAGTTCACATATTGAata                        >  1:302099/1‑139 (MQ=255)
                                                                                                                tCAATACTCCACTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACGTTCCCCAGTTCACATATTGAATATCa                     >  1:106737/1‑137 (MQ=255)
                                                                                                                          aCTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCAt                                                                                                           >  2:264229/1‑41 (MQ=255)
                                                                                                                          aCTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCAt                                                                                                           <  1:264229/41‑1 (MQ=255)
                                                                                                                          aCTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACg                                                                       >  2:7335/1‑77 (MQ=255)
                                                                                                                          aCTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACg                                                                       <  1:7335/77‑1 (MQ=255)
                                                                                                                          aCTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTAt                                                                 >  2:345171/1‑83 (MQ=255)
                                                                                                                          aCTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTAt                                                                 <  1:345171/83‑1 (MQ=255)
                                                                                                                          aCTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACaa                                                    <  1:272152/96‑1 (MQ=255)
                                                                                                                          aCTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACaa                                                    >  2:272152/1‑96 (MQ=255)
                                                                                                                          aCTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACGTTCCCCAGTTCACata                               <  2:493945/117‑1 (MQ=255)
                                                                                                                          aCTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACGTTCCCCAGTTCACata                               >  1:493945/1‑117 (MQ=255)
                                                                                                                          aCTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACGTTCCCCAGTTCACATATTGAATATCATAGTTACTg            <  1:224118/136‑1 (MQ=255)
                                                                                                                          aCTTAAACCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACGTTCCCCAGTTCACATATTGAATATCATAGTTACTg            >  2:224118/1‑136 (MQ=255)
                                                                                                                              aaaCCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACGTTCCCCAGTTCACATATTGAATATCATAGt                 <  2:370533/127‑1 (MQ=255)
                                                                                                                              aaaCCGTTTCCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACGTTCCCCAGTTCACATATTGAATATCATAGt                 >  1:370533/1‑127 (MQ=255)
                                                                                                                                    tttcCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACGTTCCCCAGTTCACATATTGAATATCATAGTTACTGAAGTTCTGTc  >  1:278940/1‑136 (MQ=255)
                                                                                                                                       |                                                                                                                                    
TTCTCAACGGGATCATCCTCAGAGGAACTATCATCAAAGTCATCATCCATAAATAAAATATCATCGAATGGTGCCACGCCCGTGATGAGTTTTATTTTATTATTACGATCAGTCAATACTCCACTTAAACCGTTTTCGCTCACAGGTTTTAATGATTTTTCATTACTCTTGTTGTAAGCAGGCGCATTAAAAATACACGGAGTATCAACATCAAACAATGACGTTCCCCAGTTCACATATTGAATATCATAGTTACTGAAGTTCTGTC  >  NZ_CP009273/2777052‑2777319

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: