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% 149.2 / NA 42pseudogene (773/986 nt)BW25113_RS21020DUF2713 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/16);  total (26/16)

CCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAA  >  NZ_CP009273/4250888‑4251154
                                                                                                                                     |                                                                                                                                     
ccTACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAAt                                                                                                                                  >  2:334377/1‑139 (MQ=255)
            tATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATg                                                                                                                       <  1:362035/138‑1 (MQ=255)
             aTCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATgg                                                                                                                      <  2:565700/138‑1 (MQ=255)
                   aaCAGTGTATATTAATGAAAATATTGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGTTATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTg                                                                                                               >  2:557938/1‑139 (MQ=255)
                        tgtATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGAt                                                                                                            >  1:311978/1‑137 (MQ=255)
                                  tGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTAt                                                                                                                  >  1:156371/1‑121 (MQ=255)
                                  tGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTAt                                                                                                                  <  2:156371/121‑1 (MQ=255)
                                  tGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGt                                                                                                     <  2:582193/134‑1 (MQ=255)
                                  tGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGt                                                                                                     >  1:582193/1‑134 (MQ=255)
                                  tGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACAt                                                                                                >  2:197780/1‑139 (MQ=255)
                                  tGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACAt                                                                                                >  1:366599/1‑139 (MQ=255)
                                      aataatGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTgg                                                                                            >  1:359298/1‑139 (MQ=255)
                                                 cTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGa                                                                                 <  1:422418/139‑1 (MQ=255)
                                                       gAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTc                                                                           >  2:161048/1‑139 (MQ=255)
                                                              cAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCt                                                                     >  1:23860/1‑138 (MQ=255)
                                                               aGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTg                                                                    <  1:197780/138‑1 (MQ=255)
                                                               aGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGa                                                                   >  2:308068/1‑139 (MQ=255)
                                                                 tttGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGATGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGata                                                                 >  2:496877/1‑139 (MQ=255)
                                                                                             agaagTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGcac                                      >  1:166536/1‑138 (MQ=255)
                                                                                               aagTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGcaccac                                   <  2:366599/139‑1 (MQ=255)
                                                                                                    ccTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATc                                                                       <  1:380721/98‑1 (MQ=255)
                                                                                                    ccTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATc                                                                       >  2:380721/1‑98 (MQ=255)
                                                                                                      tGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGcac                                      <  2:452665/129‑1 (MQ=255)
                                                                                                      tGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGcac                                      >  1:452665/1‑129 (MQ=255)
                                                                                                       gggAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGcacca                                    >  1:387996/1‑130 (MQ=255)
                                                                                                       gggAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGcacca                                    <  2:387996/130‑1 (MQ=255)
                                                                                                        ggAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCgaga                                              <  2:437204/119‑1 (MQ=255)
                                                                                                        ggAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCgaga                                              >  1:437204/1‑119 (MQ=255)
                                                                                                            agaTGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGt                                                                                                     <  2:122071/60‑1 (MQ=255)
                                                                                                            agaTGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGt                                                                                                     >  1:122071/1‑60 (MQ=255)
                                                                                                            agaTGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCtt                      >  1:124899/1‑139 (MQ=255)
                                                                                                              aTGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGa                                                                         <  2:198621/86‑1 (MQ=255)
                                                                                                              aTGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGa                                                                         >  1:198621/1‑86 (MQ=255)
                                                                                                              aTGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTcg                    >  2:353936/1‑139 (MQ=255)
                                                                                                              aTGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTcg                    >  2:241769/1‑139 (MQ=255)
                                                                                                              aTGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTcg                    >  2:146299/1‑139 (MQ=255)
                                                                                                                   gCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAg               >  2:445822/1‑139 (MQ=255)
                                                                                                                      tCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCt                                                                            <  2:415713/75‑1 (MQ=255)
                                                                                                                      tCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCt                                                                            >  1:415713/1‑75 (MQ=255)
                                                                                                                      tCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCg            >  1:316305/1‑139 (MQ=255)
                                                                                                                              aTGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGAcc     <  2:11531/138‑1 (MQ=255)
                                                                                                                              aTGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGAcc     <  1:353936/138‑1 (MQ=255)
                                                                                                                                ggggggGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCGCCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCaa  >  1:153436/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                     
CCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAA  >  NZ_CP009273/4250888‑4251154

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