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,0210AG98.2% 189.7 / ‑3.5 55pseudogene (773/986 nt)BW25113_RS21020DUF2713 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (32/22);  total (33/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

ACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTT  >  NZ_CP009273/4250891‑4251146
                                                                                                                                  |                                                                                                                             
aCGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGaa                                                                                                                       >  1:207529/1‑139 (MQ=255)
     gATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAAtttt                                                                                                                   >  2:481057/1‑138 (MQ=255)
         tATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATg                                                                                                               <  2:186523/138‑1 (MQ=255)
         tATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATg                                                                                                               <  1:235341/138‑1 (MQ=255)
         tATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATg                                                                                                               <  2:293564/138‑1 (MQ=255)
           tCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGc                                                                                                             >  1:679594/1‑138 (MQ=255)
             aGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTAt                                                                                                          <  1:32078/139‑1 (MQ=255)
                  cAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGa                                                                                                     <  2:562508/139‑1 (MQ=255)
                  cAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGa                                                                                                     <  2:547579/139‑1 (MQ=255)
                   aGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTg                                                                                                                  <  2:610025/125‑1 (MQ=255)
                   aGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTg                                                                                                                  >  1:610025/1‑125 (MQ=255)
                             aaTGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTca                                                                                           >  1:534446/1‑138 (MQ=255)
                               tGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACAt                                                                                        >  2:839522/1‑139 (MQ=255)
                               tGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACAt                                                                                        >  1:836703/1‑139 (MQ=255)
                               tGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACAt                                                                                        >  2:466948/1‑139 (MQ=255)
                               tGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACAt                                                                                        >  1:160100/1‑139 (MQ=255)
                                 aaaaTAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCt                                                                                      >  2:389994/1‑139 (MQ=255)
                                   aataatGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTgg                                                                                    >  2:502754/1‑139 (MQ=255)
                                   aataatGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTgg                                                                                    >  2:672948/1‑139 (MQ=255)
                                                aCCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGt                                                                                   >  2:176027/1‑127 (MQ=255)
                                                aCCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGt                                                                                   <  1:176027/127‑1 (MQ=255)
                                                 ccTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACg                                                                      <  1:502754/139‑1 (MQ=255)
                                                         cTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATcc                                                              >  2:552211/1‑139 (MQ=255)
                                                           cAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAAt                                                                                                                          <  2:434783/77‑1 (MQ=255)
                                                           cAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTca                                                                                           >  1:152228/1‑108 (MQ=255)
                                                           cAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTca                                                                                           <  2:152228/108‑1 (MQ=255)
                                                           cAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCt                                                             >  1:7110/1‑138 (MQ=255)
                                                           cAGTTTGAAAAATTATCTCAACGCGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAAt                                                                                                                          >  1:434783/1‑77 (MQ=255)
                                                                            tcAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATg                                           >  1:659945/1‑139 (MQ=255)
                                                                              aaCGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCg                                         <  2:7110/139‑1 (MQ=255)
                                                                              aaCGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCg                                         <  2:342575/139‑1 (MQ=255)
                                                                              aaCGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCg                                         <  1:552211/139‑1 (MQ=255)
                                                                                 gTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCgaga                                      <  2:521263/139‑1 (MQ=255)
                                                                                     tatGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGAt                                                                                               >  2:430334/1‑78 (MQ=255)
                                                                                     tatGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGAt                                                                                               <  1:430334/78‑1 (MQ=255)
                                                                                         aagaagTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGcac                              <  2:679594/139‑1 (MQ=255)
                                                                                            aagTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGcacca                            >  1:114808/1‑138 (MQ=255)
                                                                                              gTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaa                         >  1:748665/1‑139 (MQ=255)
                                                                                              gTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaa                         >  2:807863/1‑139 (MQ=255)
                                                                                               tACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaat                        >  2:537173/1‑139 (MQ=255)
                                                                                                 ccTGGGAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaatt                       >  1:279000/1‑138 (MQ=255)
                                                                                                    gggAGAGATGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaattaat                    >  2:818160/1‑138 (MQ=255)
                                                                                                           aTGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCt                                                                    <  2:791093/83‑1 (MQ=255)
                                                                                                           aTGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCt                                                                    >  1:791093/1‑83 (MQ=255)
                                                                                                           aTGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTCcg            >  2:155948/1‑139 (MQ=255)
                                                                                                           aTGAAGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTCcg            >  1:776544/1‑139 (MQ=255)
                                                                                                              aaGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTCCGCGa         <  2:651064/139‑1 (MQ=255)
                                                                                                              aaGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTCCGCGa         <  2:776544/139‑1 (MQ=255)
                                                                                                               aGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGcaccac                           <  1:594216/120‑1 (MQ=255)
                                                                                                               aGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGcaccac                           >  2:594216/1‑120 (MQ=255)
                                                                                                               aGCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTCCGCGaa        <  1:818160/139‑1 (MQ=255)
                                                                                                                gCGTCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTCCGCGAAg       >  1:300759/1‑139 (MQ=255)
                                                                                                                   tCTGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTCCGCGAAGTCg    >  1:770717/1‑139 (MQ=255)
                                                                                                                     tGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTCCGCGAAGTCGtt  >  2:743696/1‑139 (MQ=255)
                                                                                                                     tGGCTGATGGGGGGGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTCCGCGAAGTCGtt  >  2:618649/1‑139 (MQ=255)
                                                                                                                                  |                                                                                                                             
ACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGTGATATGAAGAAGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTT  >  NZ_CP009273/4250891‑4251146

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